无碳小车设计说明书(一等奖作品)1..doc

(完整word版)无碳小车设计说明书(一等奖作品)1.第二届全国大学生工程训练综合能力竞赛无碳小车设计说明书 参赛者:王金卫 指导老师：刘吉兆 陈丰峰2011-1-16摘要第二届全国大学生工程训练综合能力竞赛命题主题为“无碳小车”。在设计小车过程中特别注重设计的方法，力求通过对命题的分析得到清晰开阔的设计思路；作品的设计做到有系统性规范性和创新性；设计过程中综合考虑材料 、加工 、制造成本等给方面因素。我们借鉴了参数化设计 、优化设计 、系统设计等现代设计发发明理论方法；采用了MATLAB、PROE等软件辅助设计。我们把小车的设计分为三个阶段：方案设计、技术设计、制作调试。通过每一阶段的深入分析、层层把关，是我们的设计尽可能向最优设计靠拢。方案设计阶段根据小车功能要求我们根据机器的构成（原动机构、传动机构、执行机构、控制部分、辅助部分）把小车分为车架 、原动机构 、传动机构 、转向机构 、行走机构 、微调机构六个模块，进行模块化设计.使其每个零件或结构件具有平衡性已达到减小摩擦.提高校车整体平衡的目的。分别针对每一个模块进行多方案设计，通过综合对比选择出最优的方案组合。我们的方案为：车架采用三角底板式、原动机构采用了锥形轴、传动机构采用齿轮或没有该机构、转向机构采用四连杆机构、行走机构采用单轮驱动实现差速、微调机构采用微调螺母螺杆。其中转向机构利用了调心轴承、关节轴承.圆锥滚子轴承。技术设计阶段我们先对方案建立数学模型进行理论分析，借助MATLAB分别进行了能耗规律分析、运动学分析、动力学分析、灵敏度分析。进而得出了小车的具体参数，和运动规律。接着应用PROE软件进行了小车的实体建模和部分运动仿真。在实体建模的基础上对每一个零件进行了详细的设计，综合考虑零件材料性能、加工工艺、成本等。小车大多的零件是标准件、可以购买，同时除部分要求加工精度高的部分需要特殊加工外，大多数都可以通过手工加工出来。对于塑料会采用自制的电锯切割。因为小车受力都不大，因此大量采用胶接，简化零件及零件装配。调试过程会通过微调等方式改变小车的参数进行试验，在试验的基础上验证小车的运动规律同时确定小车最优的参数。 关键字：无碳小车 参数化设计 软件辅助设计 微调机构 灵敏度分析目录摘要2一 绪论51.1本届竞赛命题主题51.2小车功能设计要求51.3小车整体设计要求61.4小车的设计方法6二 方案设计72.1车架102.2原动机构102.3传动机构112.4转向机构112.5行走机构132.6微调机构14三 技术设计153.1建立数学模型及参数确定163.1.1能耗规律模型163.1.2运动学分析模型183.1.3动力学分析模型233.1.4灵敏度分析模型253.1.5参数确定263.2零部件设计273.3整体设计303.3.1整体装配图303.3.2小车运动仿真分析31四 小车制作调试及改进324.1小车制作流程32详见工艺分析方案报告324.2小车调试方法324.3小车改进方法32五 评价分析335.1小车优缺点335.2自动行走比赛时的前行距离估计335.3改进方向33六 参考文献33七 附录347.1装配图347.2耗能分析程序387.3运动学分析程序397.4动力学分析程序417.5灵敏度分析程序43一 绪论1.1本届竞赛命题主题本届竞赛命题主题为“无碳小车”。命题与高校工程训练教学内容相衔接，体现综合性工程能力。命题内容体现“创新设计能力、制造工艺能力、实际操作能力和工程管理能力”四个方面的要求。1.2小车功能设计要求给定一重力势能，根据能量转换原理，设计一种可将该重力势能转换为机械能并可用来驱动小车行走的装置。该自行小车在前行时能够自动避开赛道上设置的障碍物（每间隔1米，放置一个直径20mm、高200mm的弹性障碍圆棒）。以小车前行距离的远近、以及避开障碍的多少来综合评定成绩。给定重力势能为5焦耳（取g=10m/s2），竞赛时统一用质量为1Kg的重块（5065 mm，普通碳钢）铅垂下降来获得，落差5002mm，重块落下后，须被小车承载并同小车一起运动，不允许掉落。要求小车前行过程中完成的所有动作所需的能量均由此能量转换获得，不可使用任何其他的能量形式。小车要求采用三轮结构（1个转向轮，2个驱动轮），具体结构造型以及材料选用均由参赛者自主设计完成。要求满足：小车上面要装载一件外形尺寸为6020 mm的实心圆柱型钢制质量块作为载荷，其质量应不小于750克；在小车行走过程中，载荷不允许掉落。转向轮最大外径应不小于30mm。1.3小车整体设计要求 小车设计过程中需要完成：机械设计、工艺方案设计、经济成本分析和工程管理方案设计。命题中的工程管理能力项要求综合考虑材料、加工、制造成本等各方面因素，提出合理的工程规划。设计能力项要求对参赛作品的设计具有创新性和规范性。命题中的制造工艺能力项以要求综合运用加工制造工艺知识的能力为主。1.4小车的设计方法小车的设计一定要做到目标明确，通过对命题的分析我们得到了比较清晰开阔的设计思路。作品的设计需要有系统性规范性和创新性。设计过程中需要综合考虑材料 、加工 、制造成本等给方面因素。小车的设计是提高小车性能的关键。在设计方法上我们借鉴了参数化设计 、优化设计 、系统设计等现代设计发发明理论方法。采用了MATLAB、PROE等软件辅助设计。下面是我们设计小车的流程（如图一）图一二 方案设计 通过对小车的功能分析小车需要完成重力势能的转换、驱动自身行走、自动避开障碍物。为了方便设计这里根据小车所要完成的功能将小车划分为五个部分进行模块化设计（车架 、原动机构 、传动机构 、转向机构 、行走机构 、微调机构）。为了得到令人满意方案，采用扩展性思维设计每一个模块，寻求多种可行的方案和构思。下面为我们设计图框（图二）图二在选择方案时应综合考虑功能、材料、加工、制造成本等各方面因素，同时尽量避免直接决策，减少决策时的主观因素，使得选择的方案能够综合最优。图三2.1车架车架不用承受很大的力，精度要求低。考虑到重量加工成本等，车架采用木材加工制作成三角底板式。可以通过回收废木材获得，已加工。2.2原动机构原动机构的作用是将重块的重力势能转化为小车的驱动力。能实现这一功能的方案有多种，就效率和简洁性来看绳轮最优。小车对原动机构还有其它的具体要求。1.驱动力适中，不至于小车拐弯时速度过大倾翻，或重块晃动厉害影响行走。2.到达终点前重块竖直方向的速度要尽可能小，避免对小车过大的冲击。同时使重块的动能尽可能的转化到驱动小车前进上，如果重块竖直方向的速度较大，重块本身还有较多动能未释放，能量利用率不高。3.由于不同的场地对轮子的摩擦摩擦可能不一样，在不同的场地小车是需要的动力也不一样。在调试时也不知道多大的驱动力恰到好处。因此原动机构还需要能根据不同的需要调整其驱动力。4.机构简单，效率高。基于以上分析我们提出了输出驱动力可调的绳轮式原动机构.和两级减速的齿轮结构。如下图四图四如上图我们可以通过改变绳子绕在绳轮上不同位置来改变其输出的动力。2.3传动机构传动机构的功能是把动力和运动传递到转向机构和驱动轮上。要使小车行驶的更远及按设计的轨道精确地行驶，传动机构必需传递效率高、传动稳定、结构简单重量轻等。1.不用其它额外的传动装置，直接由动力轴驱动轮子和转向机构，此种方式效率最高、结构最简单。在不考虑其它条件时这是最优的方式。2.带轮具有结构简单、传动平稳、价格低廉、缓冲吸震等特点但其效率及传动精度并不高。不适合本小车设计。3.齿轮具有效率高、结构紧凑、工作可靠、传动比稳定但价格较高。因此在第一种方式不能够满足要求的情况下优先考虑使用齿轮传动。2.4转向机构转向机构是本小车设计的关键部分，直接决定着小车的功能。转向机构也同样需要尽可能的减少摩擦耗能，结构简单，零部件已获得等基本条件，同时还需要有特殊的运动特性。能够将旋转运动转化为满足要求的来回摆动，带动转向轮左右转动从而实现拐弯避障的功能。能实现该功能的机构有：凸轮机构+摇杆、曲柄连杆+摇杆、曲柄摇杆、差速转弯等等。凸轮：凸轮是具有一定曲线轮廓或凹槽的构件，它运动时，通过高副接触可以使从动件获得连续或不连续的任意预期往复运动。优点：只需设计适当的凸轮轮廓，便可使从动件得到任意的预期运动，而且结构简单、紧凑、设计方便；缺点：凸轮轮廓加工比较困难。在本小车设计中由于：凸轮轮廓加工比较困难、尺寸不能够可逆的改变、精度也很难保证、重量较大、效率低能量损失大（滑动摩擦）因此不采用 曲柄连杆+摇杆 优点：运动副单位面积所受压力较小，且面接触便于润滑，故磨损减小，制造方便，已获得较高精度；两构件之间的接触是靠本身的几何封闭来维系的，它不像凸轮机构有时需利用弹簧等力封闭来保持接触。缺点：一般情况下只能近似实现给定的运动规律或运动轨迹，且设计较为复杂；当给定的运动要求较多或较复杂时，需要的构件数和运动副数往往比较多，这样就使机构结构复杂，工作效率降低，不仅发生自锁的可能性增加，而且机构运动规律对制造、安装误差的敏感性增加；机构中做平面复杂运动和作往复运动的构件所长生的惯性力难以平衡，在高速时将引起较大的振动和动载荷，故连杆机构常用于速度较低的场合。在本小车设计中由于小车转向频率和传递的力不大故机构可以做的比较轻，可以忽略惯性力，机构并不复杂，利用MATLAB进行参数化设计并不困难，加上个链接可以利用轴承大大减小摩擦损耗提高效率。对于安装误差的敏感性问题我们可以增加微调机构来解决。曲柄摇杆结构较为简单，但和凸轮一样有一个滑动的摩擦副，其效率低。其急回特性导致难以设计出较好的机构。差速转弯差速拐是利用两个偏心轮作为驱动轮，由于两轮子的角速度一样而转动半径不一样，从而使两个轮子的速度不一样，产生了差速。小车通过差速实现拐弯避障。差速转弯，是理论上小车能走的最远的设计方案。和凸轮同样，对轮子的加工精度要求很高，加工出来后也无法根据需要来调整轮子的尺寸。（由于加工和装配的误差是不可避免的）综合上面分析我们选择曲柄连杆+摇杆作为小车转向机构的方案。2.5行走机构 行走机构即为三个轮子，轮子又厚薄之分，大小之别，材料之不同需要综合考虑。有摩擦理论知道摩擦力矩与正压力的关系为 对于相同的材料为一定值。而滚动摩擦阻力，所以轮子越大小车受到的阻力越小，因此能够走的更远。但由于加工问题材料问题安装问题等等具体尺寸需要进一步分析确定。由于小车是沿着曲线前进的，后轮必定会产生差速。对于后轮可以采用双轮同步驱动，双轮差速驱动，单轮驱动。双轮同步驱动必定有轮子会与地面打滑，由于滑动摩擦远比滚动摩擦大会损失大量能量，同时小车前进受到过多的约束，无法确定其轨迹，不能够有效避免碰到障碍。双轮差速驱动可以避免双轮同步驱动出现的问题，可以通过差速器或单向轴承来实现差速。差速器涉及到最小能耗原理，能较好的减少摩擦损耗，同时能够实现满足要运动。单向轴承实现差速的原理是但其中一个轮子速度较大时便成为从动轮，速度较慢的轮子成为主动轮，这样交替变换着。但由于单向轴承存在侧隙，在主动轮从动轮切换过程中出现误差导致运动不准确，但影响有多大会不会影响小车的功能还需进一步分析。单轮驱动即只利用一个轮子作为驱动轮，一个为导向轮，另一个为从动轮。就如一辆自行车外加一个车轮一样。从动轮与驱动轮间的差速依靠与地面的运动约束确定的。其效率比利用差速器高，但前进速度不如差速器稳定，传动精度比利用单向轴承高。综上所述行走机构的轮子应有恰当的尺寸，可以如果有条件可以通过实验来确定实现差速的机构方案，如果规则允许可以采用单轮驱动。2.6微调机构一台完整的机器包括：原动机、传动机、执行机构、控制部分、辅助设备。微调机构就属于小车的控制部分。由于前面确定了转向采用曲柄连杆+摇杆方案，由于曲柄连杆机构对于加工误差和装配误差很敏感，因此就必须加上微调机构，对误差进行修正。这是采用微调机构的原因之一，其二是为了调整小车的轨迹（幅值，周期，方向等），使小车走一条最优的轨迹。微调机构可以采用下面两种方式微调螺母式杆、滑块式如图五图五 由于理论分析与实际情况有差距，只能通过理论分析得出较优的方案而不能得到最优的方案。因此我们设计了一种机构简单的小车，通过小部分的改动便可以改装成其它方案，再通过试验比较得到最优的小车。三 技术设计技术设计阶段的目标是完成详细设计确定个零部件的的尺寸。设计的同时综合考虑材料加工成本等各因素。3.1建立数学模型及参数确定通过对小车建立数学模型，可以实现小车的参数化设计和优化设计，提高设计的效率和得到较优的设计方案。充分发挥计算机在辅助设计中的作用。 3.1.1能耗规律模型为了简化分析，先不考虑小车内部的能耗机理。设小车内部的能耗系数为，即小车能量的传递效率为。小车轮与地面的摩阻系数为，理想情况下认为重块的重力势能都用在小车克服阻力前进上。则有为第i个轮子对地面的压力。为第i个轮子的半径。为第i个轮子行走的距离为小车总质量为了更全面的理解小车的各个参数变化对小车前进距离的变化下面分别从1.轮子与地面的滚动摩阻系数、2.轮子的半径、3.小车的重量、4.小车能量转换效率。四方面考虑。 通过查阅资料知道一般材料的滚动摩阻系数为0.1-0.8间。下图为当车轮半径分别为（222mm，70mm）摩阻系数分别为0.3，0.4，0.5.mm时小车行走的距离与小车内部转换效率的坐标图（图六）有上图六可知滚动摩阻系数对小车的运动影响非常显著，因此在设计小车时也特别注意考虑轮子的材料，轮子的刚度尽可能大，与地面的摩阻系数尽可能小。同时可看到小车为轮子提供能量的效率提高一倍小车前进的距离也提高一倍。因此应尽可能减少小车内部的摩擦损耗，简化机构，充分润滑。图七为当摩阻系数为0.5mm，车轮半径依次增加10mm时的小车行走的距离与小车内部转换效率的坐标图图六图七由图可知当小车的半径每增加1cm小车便可多前进1m到2m。因此在设计时应考虑尽可能增大轮子的半径。3.1.2运动学分析模型符号说明：驱动轮半径齿轮传动比驱动轮A与转向轮横向偏距驱动轮B与转向轮横向偏距驱动轴（轴2）与转向轮中心距离曲柄轴（轴1）与转向轮中心距离曲柄的旋转半径摇杆长连杆长轴的绳轮半径a、驱动：当重物下降时，驱动轴（轴2）转过的角度为，则有则曲柄轴（轴1）转过的角度小车移动的距离为（以A轮为参考）b、转向：当转向杆与驱动轴间的夹角为时，曲柄转过的角度为则与满足以下关： 解上述方程可得与的函数关系式c、小车行走轨迹只有A轮为驱动轮，当转向轮转过角度时，如图：则小车转弯的曲率半径为小车行走过程中，小车整体转过的角度当小车转过的角度为时，有d、小车其他轮的轨迹以轮A为参考，则在小车的运动坐标系中，B的坐标C的坐标在地面坐标系中，有整理上述表达式有： 为求解方程，把上述微分方程改成差分方程求解，通过设定合理的参数的到了小车运动轨迹如（图六）图六3.1.3动力学分析模型a、驱动如图：重物以加速度向下加速运动，绳子拉力为，有产生的扭矩，（其中是考虑到摩擦产生的影响而设置的系数。）驱动轮受到的力矩，曲柄轮受到的扭矩，为驱动轮A受到的压力，为驱动轮A提供的动力，有（其中是考虑到摩擦产生的影响而设置的系数）b、转向假设小车在转向过程中转向轮受到的阻力矩恒为，其大小可由赫兹公式求得，由于b比较小，故对于连杆的拉力，有c、小车行走受力分析设小车惯量为，质心在则此时对于旋转中心的惯量为（平行轴定理）小车的加速度为：整理上述表达式得：3.1.4灵敏度分析模型小车一旦设计出来在不改变其参数的条件下小车的轨迹就已经确定，但由于加工误差和装配误差的存在，装配好小车后可能会出现其轨迹与预先设计的轨迹有偏离，需要纠正。其次开始设计的轨迹也许并不是最优的，需要通过调试试验来确定最优路径，着同样需要改变小车的某些参数。为了得到改变不同参数对小车运行轨迹的影响，和指导如何调试这里对小车各个参数进行灵敏度分析。通过MATLAB编程得到幅值周期方向i-0.0117-0.09158528.135b176.5727-35.3795578.82R-0.316316.39132528.1437a11.465469-0.27592528.5547曲柄半径r123.71445-18.9437535.3565d0.040819-117.738528.1465转向杆的长-1.637693.525236527.5711连杆长度-176.955-196.268477.35613.1.5参数确定单位：m转向轮与曲柄轴轴心距 b=0.15;摇杆长c=0.06;驱动轮直径D=0.355;驱动轮A与转向轮横向偏距a1=0.08驱动轮B与转向轮横向偏距a2=0.08;驱动轴与转向轮的距离d=0.18;曲柄长r1=0.01347;绳轮半径r2=0.0063.2零部件设计需加工的零件：a驱动轴6061空心铝合金管。外径6mm 内径3mmb车轮聚甲醛板（POM板材）。厚度：8mm，规格尺寸：600*1200mm2.2可购买的标准件：a单向离合器轴承2个型號Bearing Number外型尺寸(mm)FC系列dDFC-6K(2)610bRBL关节轴承1个：SQ 5-RSc调心球轴承1个轴承代号dDB1355196d深沟球轴承1个型号内径(d)外径(D)R85zz58d圆柱直齿轮1对小齿轮：模数=1，齿数=15，外径=17mm，内孔=3mm，厚度：6.5mm大齿轮：模数=1，齿数=45，外径=47mm，内径=10mm，厚度=10mm材质：夹布塑料3.3整体设计3.3.1整体装配图3.3.2小车运动仿真分析为了进一步分析本方案的可行性，我们利用了proe和MATLAB进行了动态仿真，详见视频。四 小车制作调试及改进4.1小车制作流程详见工艺分析方案报告4.2小车调试方法小车的调试是个很重要的过程，有了大量的理论依据支撑，还必须用大量的实践去验证。小车的调试涉及到很多的内容，如车速的快慢，绕过障碍物，小车整体的协调性，小车前进的距离等。（1）小车的速度的调试：通过小车在指定的赛道上行走，测量通过指定点的时间，得到多组数据，从而得出小车行驶的速度，通过试验，发现小车后半程速度较快，整体协调性能不是太好，于是车小了绕绳驱动轴，减小过大的驱动力同时也增大了小车前进的距离。（2）小车避障的调试：虽然本组小车各个机构相对来说较简单，损耗能量较少，但是避障不是很好，但与此同时，小车由于设计时采用了多组微调机构，通过观察小车在指定赛道上行走时避障的特点，微调螺母，慢慢小车避障性能改善，并做好标记。4.3小车改进方法由于本组小车采用胶水黏贴各处，虽然少了许多的加工成本费用，也避免了能量的过多损耗，但小车会有时出现脱胶的现象，导致无法前进，于是想法改进，使小车能量损失减少，同时故障出现的次数减少，稳定性能较好，避障多，前进远。另外，本组采用微调机构，但通过计算编程发现要求精度非常高，改变0.001mm都可能使小车偏离原轨道，于是想法改进使小车精度降低，加工成本也减低。五 评价分析5.1小车优缺点优点：（1）小车机构简单，单级齿轮传动，损耗能量少，（2）多处采用微调机构，便于纠正轨迹，避开障碍物，（3）采用大的驱动轮，滚阻系数小，行走距离远，（4）采用磁阻尼，小车稳定性提高，不致使车速过快，缺点： 小车精度要求高，使得加工零件成本高，以及微调各个机构都很费时，避障稳定行差，时而偏左，时而偏右。5.2自动行走比赛时的前行距离估计通过理论与实践结合，小车行走距离（包括绕开障碍物）约20-25米。5.3改进方向小车最大的缺点是精度要求非常高，改进小车的精度要求，使能调整简单，小车便能达到很好的行走效果。六 参考文献七 附录7.1装配图7.2耗能分析程序clearclctic%符号定义%重物下降的高度h%小车行驶的路程s%内部能耗系数ypxln=10000;h=0.5;nn=1000;ypxl=linspace(0.5,1,n);R2=111/nn;R1=35/nn;m=1;g=9.8;mz=2; sgm=0.5/nn;for i=1:10 % sgm=(0.1*i+0.2)/nn; %mz=1.75-0.2+0.2*i;R1=R1+20/nn;R2=R2+20/nn;s=ypxl*m*h/(mz*(1/R1+2/R2)*sgm);s=s/1.045615886000699;plot(ypxl,s);hold ongrid onendplot(0.5,0);toc7.3运动学分析程序clearclctic%符号定义%重物下降的高度h%驱动轴转过角度sd2%驱动轴传动比ii%转向轮轴心距b%转向杆的长c%转向轮转过的角度af%驱动轮半径R%驱动轮A与转向轮横向偏距a1%驱动轮B与转向轮横向偏距a2%驱动轴与转向轮的距离d%小车行驶的路程s%小车x方向的位移x%小车y方向的位移y%轨迹曲率半径rou%曲柄半径r1%绳轮半径r2%参数输入n=1000;h=linspace(0,0.5,n);ii=3;b=0.15;R=0.111;%驱动轮A与转向轮横向偏距a1a1=0.08;%驱动轮B与转向轮横向偏距a2a2=0.08;%曲柄半径r1r1=0.01347;%绳轮半径r2r2=0.006;%驱动轴与转向轮的距离dd=0.18; %转向杆的长cc=0.06;l=sqrt(b2+r12)+(0.351)/1000; %算法g=-10;sd2=h/r2;sd1=sd2/ii+pi/2;C=l2-2*c2-r12.*(cos(sd1).2-(b-r1.*sin(sd1).2;A=2.*c.*(b-r1.*sin(sd1);B=-2*c2;af=asin(C./sqrt(A.2+B.2)-atan(B./A);format longrou=a1+(d)./(tan(af);s=sd2*R;ds=s(2)-s(1);dbd=ds./(rou);bd=cumsum(dbd);dy=ds*cos(bd);dx=-ds*sin(bd);x=cumsum(dx);y=cumsum(dy);xb=x-(a1+a2).*cos(bd);yb=y-(a1+a2).*sin(bd);xc=x-a1*cos(bd)-d*sin(bd);yc=y-a1*sin(bd)+d*cos(bd); plot(x,y,b,xb,yb,b,xc,yc,m); hold on grid on for i=1:9 t=0:0.01:2*pi; xy=0.01.*cos(t)-0.23; yy=0.01.*sin(t)+i; plot(xy,yy); hold on endtoc7.4动力学分析程序clearclcticn=1000;h=linspace(0,0.5,n);ii=3;b=0.15;R=0.111;%驱动轮A与转向轮横向偏距a1a1=0.08;%驱动轮B与转向轮横向偏距a2a2=0.08;%曲柄半径r1r1=0.01347;%绳轮半径r2r2=0.006;%驱动轴与转向轮的距离dd=0.18;%转向杆的长cc=0.06;l=sqrt(b2+r12)+(0.351)/1000;%算法g=-10;sd2=h/r2;sd1=sd2/ii+pi/2;C=l2-2*c2-r12.*(cos(sd1).2-(b-r1.*sin(sd1).2;A=2.*c.*(b-r1.*sin(sd1);B=-2*c2;af=asin(C./sqrt(A.2+B.2)-atan(B./A);format longrou=a1+(d)./(tan(af);s=sd2*R;ds=s(2)-s(1);dbd=ds./(rou);bd=cumsum(dbd);dy=ds*cos(bd);dx=-ds*sin(bd);x=cumsum(dx);y=cumsum(dy);xb=x-(a1+a2).*cos(bd);yb=y-(a1+a2).*sin(bd);xc=x-a1*cos(bd)-d*sin(bd);yc=y-a1*sin(bd)+d*cos(bd); toc %动力学分析%参数输入%重物质量m=1;%小车总质量mc=1.6+1;Nc=9.8*mc/3;%小车惯量rc=0.07;I=mc*rc2;a3=0.05;II=I+m*(rou-a1).2+a32);%传动效率lmd=0.5;%前轮半径RC=0.05;%前轮宽度B=2/1000;%弹性模量E1=100*1000000000;E2=150*1000000000;mu=0.2;%接触应力sgmc=sqrt(Nc/B/RC)/(2*pi*(1-mu2)/E1);bc=Nc/sgmc/2/B;%摩擦因素mucmuc=0.1;%摩擦力矩McMc=sgmc*muc*bc*B2/4;%摩阻系数sgm=0.5/1000;mMN=rou.*(m*9.8*r2*lmd-Nc*sgm)/R;K=rou.*m*r22*lmd/R2;NCNB=Nc*sgm.*sqrt(rou-a1).2+d2)/RC+Nc*sgm*(rou-a1-a2);RIA=II./rou;NRA=NCNB*R./rou;aa=(mMN-NCNB)./(K+RIA); plot(y,aa)hold on 7.5灵敏度分析程序cleartic%符号定义%重物下降的高度h%驱动轴转过角度sd2%驱动轴与圆柱凸轮轴传动比ii%转向轮与圆柱凸轮轴心距b%转向杆的长c%转向轮转过的角度af%驱动轮半径R%驱动轮A与转向轮横向偏距a1%驱动轮B与转向轮横向偏距a2%驱动轴与转向轮的距离d%小车行驶的路程s%小车x方向的位移x%小车y方向的位移y%轨迹曲率半径rou%曲柄半径r1%绳轮半径r2%参数输入n=10000;h=linspace(0,0.5,n);ii=3;b=0.15;R=0.111;%驱动轮A与转向轮横向偏距a1a1=0.08;%驱动轮B与转向轮横向偏距a2a2=0.08;%曲柄半径r1r1=0.01347;%绳轮半径r2r2=0.006;%驱动轴与转向轮的距离dd=0.18;%转向杆的长cc=0.06;l=sqrt(b2+r12)+(0.351)/1000;aa=zeros(1,8);kk=zeros(3,8);A1=zeros(9,4);ddc=0.000001;aa(1,1)=ii;aa(1,2)=b;aa(1,3)=R;aa(1,4)=a1;aa(1,5)=r1;aa(1,6)=r2;aa(1,7)=c;aa(1,8)=l;%算法for i=1:9 if i1 aa(1,i-1)=aa(1,i-1)+ddc; endii=aa(1,1);b=aa(1,2);R=aa(1,3);a1=aa(1,4);r1=aa(1,5);r2=aa(1,6);c=aa(1,7);l=aa(1,8); g=-10;sd2=h/r2;sd1=sd2/ii+pi/2;C=l2-2*c2-r12.*(cos(sd1).2-(b-r1.*sin(sd1).2;A=2.*c.*(b-r1.*sin(sd1);B=-2*c2;af=asin(C./sqrt(A.2+B.2)-atan(B./A);format longrou=a1+(d)./(tan(af);s=sd2*R;ds=s(2)-s(1);dbd=ds./(rou);bd=cumsum(dbd);dy=ds*cos(bd);dx=-ds*sin(bd);x=cumsum(dx);y=cumsum(dy);plot(x,y)grid onhold onfor j=fix(6.5*n/9):fix(8.5*n/9) if x(j)=min(x(fix(6.5*n/9):fix(8.5*n/9) A1(i,1)=x(j); A1(i,2)=y(j); endif x(j)=max(x) A1(i,3)=x(j); A1(i,4)=y(j);endendif i1aa(1,i-1)=aa(1,i-1)-ddc;endendfor i=2:9 kk(1,i-1)=(A1(i,1)-A1(1,1)/ddc/14;%幅值 kk(2,i-1)=(A1(i,4)-A1(1,4)/ddc/4;%波长 kk(3,i-1)=A1(i,3)/ddc/8/8;%角度endtockk52