1、汽车空气悬架建模与仿真分析1汽车空气悬架建模与仿真分析作者 :刘宜东 指导教师:宋宇(安徽农业大学经济技术学院 合肥 07530087)摘要：悬架是汽车的重要总成之一，它对汽车的平顺性和操纵稳定性有很大的影响。在汽车悬架系统的设计和开发过程中，其运动学、动力学和采用控制策略的计算分析占有十分重要的地位。空气弹簧具有优良的弹性特性，用在车辆悬挂装置中不仅能大大改善车辆的动力性能，显著提高车辆的运行舒适度；还能降低汽车振动频率和车轮动载荷，使其获得良好的行驶平顺性、操纵稳定性和行驶安全性，减小高速行驶车辆对路面的破坏。本文就汽车空气弹簧悬架及空气悬架系统的特点，对整车的影响通过简单计算，再把计算的
2、结果输入计算机，由计算机模拟出图形，我们将得到研究结果。关键词：汽车空气悬架、空气弹簧、建模、计算机仿真1 引言30 年代初，美国法尔斯通轮胎和橡胶公司第一次真正把空气弹簧用于汽车工业。哈维法尔斯通在其好友亨利福特一世和托马斯阿瓦爱迪生的技术支持下，研制出了空气柱形式的空气弹簧悬架系统。于是在 1934 年就诞生了 AIREDE 空气弹簧。51 年前，美国纽威安柯洛克国际公司（Neway Anchorlok lnternational）成立时即作为一家架车悬架系统的生产厂家，为公路和非公路行驶的重型机车设计和制造钢板弹簧悬架系统。由于纽威在重型车辆市场上取得了成功，后来就向高速公路车辆悬架系统
3、方向发展。35 年前，纽威向市场上投放了世界上第一种实际应用的空气悬架系统。从此以后，纽威开发出一系列空气悬架产品，应用于世界各地的客车、载货车和架车。纽威提供的空气悬架产品约占北美和欧洲用于客车、载货车和架车市场的 70。2 发展趋势汽车空气悬架建模与仿真分析2随着高档客车制造技术的引进以及人们对舒适性要求的提高，加上国家对客车等级划分的标准要求，空气悬架才开始在我国逐步应用起来。目前空气悬架主要集中应用在高等级客车上，但是受多方面因素的制约，空气悬架的配置率仍然很低，基本上还属于“导入”阶段。中国是最新的前沿阵地，正在把钢板弹簧更换为空气悬架弹簧。空气悬架发展的历史经验告诉我们，引入空气悬
4、架的国家一般是首先将其用于客车，随后就向载货车和架车方向发展，中国也会有这样的发展过程。3 空气悬架系统简述悬架是连接车身和车轮之间一切传力装置的总称，主要由弹簧（如钢板弹簧、螺旋弹簧、空气弹簧、扭杆等） 、减振器和导向机构三部分组成。当汽车在不同路面上行驶时，由于悬架系统实现了车身和车轮之间的弹性支承，有效地降低了车身与车轮的振动，从而改善了汽车行驶的平顺性和操纵稳定性。3.1 空气悬架的特点空气悬架系统是以空气弹簧为弹性元件，利用气体的可压缩性实现其弹性作用的压缩气体的气压能够随载荷和道路条件变化而进行自动调节，不论满载还是空载，整车高度不会变化，可以大大提高乘坐的舒适性。空气弹簧的运动性
5、能特点是：负载能力可调；弹性系数随负载变化；负载变化肘,固有频率几乎不变；固有频率较低。这些特点决定了空气悬架具有以下优点：(1)乘坐更舒适安全；(2)改善车辆的行驶平顺性；(3)延长轮胎和制动片的使用寿命；(4)负载变化时车身高度不变；(5)减少电气、空调、排气系统、车桥、车身和底盘的维修成本；(6)减少对道路的冲击，保护路面，降低高速公路的维修费用；(7)延长车辆的使用寿命并增加折旧值。3.2 空气悬架的结构形式空气悬架的主要组成部分除了空气弹簧组件、减振阻尼装置之外，还有导向机构、高度控制阀、横向稳定器和缓冲限位块等。汽车空气悬架建模与仿真分析3空气弹簧是一种利用橡胶气囊内部压缩空气的反
6、力作为弹性恢复力的新型弹性元件，它利用橡胶的弹性和空气压力获得综合具有吸振、减振、隔振、防噪和缓冲等功能。作为隔振支承装置使用时，空气弹簧具有优良的弹性特性，与普通的钢制弹簧和橡胶制弹簧相比，它具有以下特点：(1)空气弹簧具有非线性特性，可将其特性曲线设计成理想形状；(2)空气弹簧质量轻，内摩擦小，对高频振动有很好的隔振、消声能力：(3)空气弹簧的刚度和承载能力可以通过调节橡胶气囊内的压力来调整；(4)负载变化时，固有频率几乎不变，且固有频率较低；(5)空气弹簧具有较高的疲劳寿命，其疲劳寿命可达 300 万次以上，实际使用寿命可达 5 年以上。如今，空气弹簧因具有良好的减振和隔振性能而被广泛应
7、用于汽车、铁道车辆以及舰船、航空、工业等领域的动力机械、电子设备、仪器仪表、化工机械等各种设备当中。空气弹簧的结构可以设计成多种类型。根据橡胶气囊工作时的变形方式，空气弹簧可分为囊式空气弹簧、膜式空气弹簧两种。（如图所示）囊式空气弹簧 膜式空气弹簧囊式空气弹簧由夹有帘线的橡胶制成的气囊和密闭在其中的压缩空气构成。气囊上下盖板将空气封于囊内。根据橡胶气囊曲数的不同，囊式空气弹簧可以分为单曲、双曲和多曲囊式空气弹簧。气囊各段之间镶有金属轮缘，目的是承受内压张力。囊式空气弹簧有效面积变化率较大，弹簧刚度较大，振动频率也较高。为了获得低的振动频率，可设置辅助气室。但过大的辅助气室对降低振动频率的效果不
8、是很明显，因此在设计时辅助气室的容积最大不超过原气囊的三倍。增加气囊曲数时，由于气囊的变形可由各个曲面平均分担，因此曲数越多，有效直径变化率就越小。可见增加气囊曲数会降低囊式空气弹簧的刚度，降低弹簧的振动频率。膜式空气弹簧的结构是在盖板和底座之问放置一圆柱形橡胶气巍，通过气囊挠曲变形实现整体伸缩。工作时，膜式空气弹簧橡胶气囊沿活塞面发生变形，通汽车空气悬架建模与仿真分析4过橡胶气囊的卷曲变形实现整体伸缩。膜式空气弹簧有效面积的变化率小，因此，膜式空气弹簧在辅助气室较小的情况下，也可得到较低的自振频率。在萁正常工作范围内，膜式空气弹簧刚度变化要比囊式小，弹性曲线更为理想，固有频牢更低些。同时膜式
9、窄气弹簧可以通过改变活塞底座的形状和利用活塞底庳的空心内腔作为辅助气室来优化其刚度特性，从而获得那想的非线性特性，膜式空气弹簧尺寸小，便于布置，刚而多用于小轿车上。3.3 空气弹簧的工作原理空气悬架工作原理就是用空气压缩机形成压缩空气，并将压缩空气送到弹簧和减振器的空气室中，以此来改变车辆的高度。在前轮和后轮的附近设有车高传感器，按车高传感器的输出信号，微机判断出车身高度的变化，再控制压缩机和排气阀，使弹簧压缩或伸长，从而起到减振的效果。 空气悬架给予了汽车更多的灵性。当你在高速行驶时悬架可以变硬来提高车身的稳定性；而长时间在低速不平的路面行驶时，控制单元会使悬架变软来提高车子的舒适性。3.4
10、 空气悬架对整车性能的影响1. 空气悬架为刚度可变的非线性悬架。当簧载质量变化时，刚度随之变化，以保持空载和满载时车身高度相同，悬架固有频率基本不变。根据需要，可以选择不同的气囊工作高度，获得理想的固有频率，从而得到良好的行驶平顺性。2. 空气悬架质量轻，弹簧刚度低，高速行驶时，轮胎与地面的附着能力强，制动距离短；转向时，过多转向和不足转向倾向减小，转向稳定性强，提高了整车的操纵稳定性。3. 空气弹簧内的空气压力直接反映了簧载质量，可取空气压力作为信号，控制制动缸内的气压，来控制制动时的制动力，更好地保证了行驶安全性。4. 可通过给空气弹簧气囊充气或放气来调节车身高度。在平坦的路面上，降低车身
11、高度，保持空气阻力系数为最佳值，可以减小油耗或在功率不变的情况下获得最大车速。4 空气弹簧悬架的建模及分析由于汽车是一个复杂的振动系统，为了便于分析解决问题，常对其进行简化。在预测汽车的振动响应情况时，首先应该建立一个尽可能真实反映系统的数学模汽车空气悬架建模与仿真分析5型。选择模型时考虑因素的全面与否同预测结果的精确度有很大的关系。建立理论模型的一般法则是把系统刚度较大的某一部分看作是具有集中质量的刚体，而把刚度较小的那一部分当作弹性连接元件，然后将这些刚体和连接元件组合起来，就形成一个机械振动体系。通常认为理论模型的自由度数目越多，预测的结果就越精确。但由于多自由度系统的计算要求测定的相关
12、参数就多，而在车辆设计阶段许多参数并不能准确测定，会给计算结果带来很大的误差。另就研究汽车平顺性而言，由路面不平激起的各种振动成分的作用也有大小和主次之分。在试验的基础上得知车身地板水平方向变化数值就很小，对平顺性影响不到 3。通常在研究平顺性时，暂时忽略汽车的侧向运动，着重考虑车辆垂直平面内的运动。4.1 建模在研究车辆振动性能的模型中，根据研究方便，我们简化成1/4悬架模型。汽车行使时受到路面阶跃为X 0，轮胎的缓冲作用令其弹性刚度为K 1， 空气弹簧刚度为K 2，减振器阻尼系数为Cd，根据研究的方便及需要在本文中简化成单质量系统模型。简化的结果如下图1 空气弹簧悬架的简化模型它由车身质量
13、M 2， 车质量轮M 1，车轮刚度系数K 1，减振阻尼器阻尼系数Cd的悬架组成。其中X 0为地面阶跃，X 2为车身位移，X 1为轮胎位移。4.2 振动微分方程的获得对模型的分析，利用牛顿第二定理，得到系统运动的微分方程：对车身而言： (1)201122012 )()() MPAXMPCkXd 汽车空气悬架建模与仿真分析6对轮胎而言： (2)()()() 2110121021 XKMPAXMPCdKX 对空气弹簧内部，设其截面积不变；则(3)2102Vpnek( Pv+Pa)Vm=(Pro+Pa)Vom (4)令 y 0=x2, y1=x2, y2=x1, y3=x1则 y 0=y1y1= 20
14、12201 )()()(MPAXMPCkd y2=y3y3=- )()()()( 21101210 XKK其中上式中：M1轮胎的质量，KgM2车身的质量，KgX0地面对车的阶跃，mX1轮胎的位移，mX2车身的位移，mCd阻尼系数，NsmP空气弹簧的气压，PPr0静平衡位置时的相对气压，PaV0静平衡位置时的相对体积，m 3V任意位置的相对容积，m 3,我们近似认为v 1+v2M空气弹簧内部缓冲系数。一般为M=1.31.38，缓慢时为M=1，坏路为M=1.4。A0空气弹簧的截面积 ,m 2V1-空气弹簧的主室体积，m 3V2-空气弹簧的副室体积，m 3n-空气热力指数，对静态时n=1，对动态时n
15、=1.314。5微分方程的求解5.1 式 my+c(y-q)+k(y-q)=0（1）的解析解直接求解上式的解析解由高等数学微分方程的求解可知此方程的解由自由振动奇次方程的解之和组成。汽车空气悬架建模与仿真分析7若令 ， mcn2k20则奇次方程为：y+zny+ 02 y=0 0为系统圆频率.而阻尼对运动的影响取决于 n 和 0的比值 。 为阻尼比： kcn20汽车悬架的阻尼比 的比值通常在 0.25 左右属于小阻尼,此时方程的解为:)sin(120atAeym分析此解可知,有阻尼自由振动时,质量 m 以有阻固有频率 120nr振动,其振幅按 me衰减。以上仅为方程的解析解。在通常的科学计算中，
16、大量的微分方程无法用解析的方法求解，且没有真正的应用价值。特别是随着计算机的迅速发展和广泛的应用，人们越来越认识到科学计算是科学研究的第三种方法，特别是工科类的大学生更应该具备这方面的知识和能力。为适应时代工科大学生的要求，故可对方程（1）进行科学计算。5.2 方程（1）的数值求解常微分方程是讨论一些典型方程解析解的基本方法。然而，在生产实际和科学研究中遇到的微分方程比较复杂，在很多情况下，都不可能给出方程的解析表达式，又因计算量太大而不实用；有时，即使是一些已经有了求解的基本方法的典型方程，但在实际使用时也是有困难的。以上情况都说明：用求解析式的基本方法来计算微分方程的数值解往往是不适宜的，
17、基本上是很难办到的。在实际问题中，对于求解微分方程，一般只要求解在若干上的近似值或者解的便于计算的近似表达式（只要满足规定的精度）即可，就像本课题为了求得最后的振动方程，只要求解若干个点最后连成光滑的曲线即可。5.2.1 选用解法把式（1）进行整理可得：y”+ y+ q=0.此方程是一个二阶微分方程，mck可由多种方法对其求数值解，改进的欧拉方法和四阶龙格库塔法比较合适，又汽车空气悬架建模与仿真分析8因改进的欧拉方法比较适用于本身不太光滑的曲线，故在本次设计中可选用经典的解法：四阶龙格库塔法。5.2.2 四阶龙格库塔法的基本思想从研究商 hxyii)(1(开始，由微分方程中值定理 )()(xy
18、iii （1,其中：1. 01时，此时为过阻尼，汽车振动不明显。汽车空气悬架建模与仿真分析10在本文中，我们分别对以上01两种情况下汽车悬架的振动情况进行仿真输出，对比各种情况下在受到相同阶跃输入时候，那种情况更理想，从而达到汽车舒适性和平顺性的目的。6.1 数据的处理利用设计出来的C语言程序计算出来数据处理后分别将车轮位移、车轮速度、车身位移、车身速度、车身加速度、车轮加速度转化为test. txt、test2. txt、test3. txt、test4. txt、test5. txt、test6.txt的数据文件存放在计算机磁盘中，以供Matlab绘图时候调用。在Matlab调用数据之前，
19、应先将以上的数据文件保存到Matlab程序文件夹的work文件里面，这样才能直接调用。在调用时候，将每个数据转化为矩阵形式，与时间参数一一对应（即t=0:0.01:1.00）进行绘图，从而的出车轮位移图、车轮速度图、车身位移图、车身加速度图、车身速度图。6.2 输出图形的结果本次仿真设计的结果是对悬架系统在汽车车轮受到阶跃输入时候，对车身的位移，速度，加速度和对车轮的位移和速度进行形象的图形输出。通过改变汽车阻尼比进行对比，具体结果如下：图1 阻尼比较小的情况，车身位移图汽车空气悬架建模与仿真分析11图2 阻尼比较小的情况下，车身速度图图3 阻尼比较小的情况下，车轮位移图图4 阻尼比较小的情况
20、下，车身加速度图汽车空气悬架建模与仿真分析12图5 阻尼比较小的情况下，车轮加速度图图6 阻尼比较大时，车身加速度图。在阻尼比大时，车身的加速度很快的降下来图7 在阻尼比大，车轮受到较大地面冲击时，也能在较短的时间内变小，使车轮在短时间平稳下来汽车空气悬架建模与仿真分析13图8 在阻尼比大，车身受到较大地面冲击时，也能在较短的时间内变小，使车身在短时间平稳下来。6.3 仿真结果分析通过对以上各图形的相互对比分析，汽车空气悬架的阻尼比在处于 01情况下的时候，当车轮受到地面阶跃输入的时候，车身的位移速度都变化的太快，给人身的刺激也就比较大，从而不能达到乘坐舒适性和驾驶平顺性的目的。所以汽车在大部
21、分情况下，阻尼比一般处于01情况下。总结及建议本文以客车空气悬架系统为研究对象，先对空气悬架系统的结构进行分析，再把空气悬架系统特点对整车的影响通过简单计算，得出计算结果后输入计算机，由计算机模拟出图形，我们得到研究结果。由于时间仓促和客观条件限制，再加上作者水平有限，所以在本文研究工作基础上还有诸多可以完善之处。现提出以下几点建议：（1）对悬架的仿真目的是为了寻找一个最合适的阻尼系数以达到汽车的最佳舒适性和平顺性，由于路面和外界干扰因素对车身的干扰都是未知的，人们无法用汽车空气悬架建模与仿真分析14什么模型来表示它，所以为了达到最佳舒适性和平顺性，悬架的阻尼系数也不应该是一个定值，只有通过不
22、断的改变悬架的阻尼系数，同时通过对不同的路面冲击的模拟，通过计算机的输出的结果，我们找出接近理想的仿真曲线，才能达到最理想的状态，这是悬架的控制，才能了解空气弹簧的优越性。（2）在空气悬架控制系统结构及控制算法的设计和研究上，国内大多停留在计算机仿真上。因此在条件允许的情况下，综合考虑整车的平顺性和操纵稳定性，对各种控制策略下的主动空气悬架系统进行实车试验，验证所研究控制策略的正确性和可行性，为进一步实用化打下基础。致 谢行文至此，我的这篇论文已接近尾声；岁月如梭，我四年的大学时光也即将敲响结束的钟声。离别在即，站在人生的又一个转折点上，心中难免思绪万千，一种感恩之情油然而生。感谢学院关心我们
23、毕业论文的老师和领导们，学院的计算机实验室给我的课题带来了很多方便；感谢在班里同学和朋友，感谢你们在我遇到困难的时候帮助我，给我支持和鼓励，感谢你们。 特别感谢我的指导老师宋宇，在本课题中给予我悉心指导，从一开始对空气弹簧的完全陌生到最后的恍然大悟，我学到了很多知识。宋老师给我的鼓励与指引，使我能够克服重重困难，将论文完成，在此谨向宋宇老师致以诚挚的谢意和崇高的敬意。谢谢！ 参 考 文 献汽车空气悬架建模与仿真分析151陈家瑞汽车构造(下)M北京：机械工业出版社，2005：199-2432徐士良编著.常用算法程序集清华大学出版社,1996 年3刘波客车电控空气悬架系统及其发展趋势J客车技术，2
24、007(6)：5-84何渝生编著.汽车控制理论基础及应用.重庆大学出版社,1995 年5陈中亮编著.Turbo科学计算与绘图.华东理工大学出版社,1996 年6方瑞华，解跃青，雷雨成空气悬架理论及其关键技术J同济大学学报:自然科学版，2003，31(9)：1072-10767严彦从 大型客车空气悬架结构特点及设计中的若干问题1997年8张建文,庄德军,林逸,王望予,刘宏伟汽车用空气弹簧悬架系统综述2002年 9ROD M V C，PRZLAD V A unable fuzzy logic controller for vehicle active suspension systemsJ Fuz
25、zy Sets and Systems,1997,85:11-2010Gieck JackRide on air：A history of air suspensionMNew York：11John Wood roffe Heavy Turuck Suspension Dynamics:Methods for Evaluacting SUSOENSION Road Friendliness and RIDE onality C.SAE paper.NO.962152Automotive Air Suspension Modeling and Simulation Author: liuyid
26、ong teacher: Sung Yu(The agriculture university hefei 07530087)Abstract Suspension is one of the important assembly car, it for car getting smooth-going and stability has very big effect. In the automobile suspension system design and development process, the kinematics and dynamics and control stra
27、tegy of calculation and analysis by occupies very important position. Air spring has excellent elastic properties, used in vehicle suspension devices can not only can greatly improve 汽车空气悬架建模与仿真分析16vehicle dynamic performance significantly improve the operation vehicle comfort; Still can reduce vehi
28、cle vibration frequency and wheel dynamic load, make its win good ride comfort and sex, handling stability and high-speed driving safety, reducing vehicle of automobi damage. This paper automobile air spring suspension and air suspension system characteristic, the influence of vehicle through a simp
29、le calculation, and then the calculation results into the computer and the graphic by computer simulation, we will get results. Keywords: the air suspension, air spring, modeling, computer simulation 附录编程模块 grkt1.c#include “stdlib.h“void grkt1(t,y,n,h,k,z)int n,k;double t,h,y,z; extern void grkt1f()
30、;int i,j,l;double a4,tt,*b,*d;b=malloc(n*sizeof(double);d=malloc(n*sizeof(double);a0=h/2.0; a1=a0;a2=h; a3=h;for (i=0; i=n-1; i+) zi*k=yi;for (l=1; l=k-1; l+) grkt1f(t,y,n,d);for (i=0; i=n-1; i+) bi=yi;for (j=0; j=2; j+) for (i=0; i=n-1; i+) yi=zi*k+l-1+aj*di;bi=bi+aj+1*di/3.0;汽车空气悬架建模与仿真分析17tt=t+aj
31、;grkt1f(tt,y,n,d);for (i=0; i=n-1; i+)yi=bi+h*di/6.0;for (i=0; i=n-1; i+)zi*k+l=yi;t=t+h;free(b); free(d);return;汽车车身速度的编程程序：#include “stdio.h“#include “grkt1.c“#include “math.h“main()int i,j;FILE *fp;double t,h,y4,z4101;y0=0.0;y1=0.0;y2=0.0;y3=0.0;t=0.0;h=0.01;grkt1(t,y,4,h,101,z);printf(“n“);if (f
32、p=fopen(“test3.txt“,“w“)=NULL)printf (“cannot open filen“);exit(0);for (i=0;i=100;i+)t=i*h;printf(“t=%5.2fn“,t);fprintf(fp,“%6.2fn“,z3i);for (j=0;j=3;j+)printf(“y(%d)=%e“,j,zji);printf(“n“);void grkt1f(t,y,n,d)int n;double t,y,d;t=t;n=n;d0=y1;汽车空气悬架建模与仿真分析18d1=-200*(y1-y3)-150*(y0-y2);d2=y3;d3=-100*
33、(y3-y1)+250*(y2-y0)-4000*(y2-sin(t);return;汽车的车身加速度的编程：#include “stdio.h“#include “grkt1.c“#include “math.h“main()int i,j;FILE *fp;double t,h,y4,z4101;y0=0.0;y1=0.0;y2=0.0;y3=0.0;t=0.0;h=0.01;grkt1(t,y,4,h,101,z);printf(“n“);if (fp=fopen(“test5.txt“,“w“)=NULL)printf (“cannot open filen“);exit(0);for
34、 (i=0;i=100;i+)t=i*h;printf(“t=%5.2fn“,t);fprintf(fp,“%6.2fn“,(z3i-z3I-1)/0.01);for (j=0;j=3;j+)printf(“y(%d)=%e“,j,zji);printf(“n“);void grkt1f(t,y,n,d)int n;double t,y,d;t=t;n=n;d0=y1;d1=-200*(y1-y3)-150*(y0-y2);d2=y3;d3=-100*(y3-y1)+250*(y2-y0)-4000*(y2-sin(t);return;汽车车轮速度的编程#include “stdio.h“#i
35、nclude “grkt1.c“汽车空气悬架建模与仿真分析19#include “math.h“main()int i,j;FILE *fp;double t,h,y4,z4101;y0=0.0;y1=0.0;y2=0.0;y3=0.0;t=0.0;h=0.01;grkt1(t,y,4,h,101,z);printf(“n“);if (fp=fopen(“test2.txt“,“w“)=NULL)printf (“cannot open filen“);exit(0);for (i=0;i=100;i+)t=i*h;printf(“t=%5.2fn“,t);fprintf(fp,“%6.2fn
36、“,z1i);for (j=0;j=3;j+)printf(“y(%d)=%e“,j,zji);printf(“n“);void grkt1f(t,y,n,d)int n;double t,y,d;t=t;n=n;d0=y1;d1=-200*(y1-y3)-150*(y0-y2);d2=y3;d3=-100*(y3-y1)+250*(y2-y0)-4000*(y2-sin(t);return;汽车车轮的加速度编程#include “stdio.h“#include “grkt1.c“#include “math.h“main()int i,j;FILE *fp;double t,h,y4,z4
37、101;y0=0.0;y1=0.0;y2=0.0;y3=0.0;t=0.0;h=0.01;grkt1(t,y,4,h,101,z);汽车空气悬架建模与仿真分析20printf(“n“);if (fp=fopen(“test6.txt“,“w“)=NULL)printf (“cannot open filen“);exit(0);for (i=0;i=100;i+)t=i*h;printf(“t=%5.2fn“,t);fprintf(fp,“%6.2fn“,(z1i-z1i-1)/0.01);for (j=0;j=3;j+)printf(“y(%d)=%e“,j,zji);printf(“n“)
38、;void grkt1f(t,y,n,d)int n;double t,y,d;t=t;n=n;d0=y1;d1=-200*(y1-y3)-150*(y0-y2);d2=y3;d3=-100*(y3-y1)+250*(y2-y0)-4000*(y2-sin(t);return;汽车车身的位移的编程#include “stdio.h“#include “grkt1.c“#include “math.h“main()int i,j;FILE *fp;double t,h,y4,z4101;y0=0.0;y1=0.0;y2=0.0;y3=0.0;t=0.0;h=0.01;grkt1(t,y,4,h,
39、101,z);printf(“n“);if (fp=fopen(“test3.txt“,“w“)=NULL)printf (“cannot open filen“);exit(0);for (i=0;i=100;i+)t=i*h;printf(“t=%5.2fn“,t);fprintf(fp,“%6.2fn“,z2i);汽车空气悬架建模与仿真分析21for (j=0;j=3;j+)printf(“y(%d)=%e“,j,zji);printf(“n“);void grkt1f(t,y,n,d)int n;double t,y,d;t=t;n=n;d0=y1;d1=-200*(y1-y3)-15
40、0*(y0-y2);d2=y3;d3=-100*(y3-y1)+250*(y2-y0)-4000*(y2-sin(t);return;汽车车轮的位移的编程#include “stdio.h“#include “grkt1.c“#include “math.h“main()int i,j;FILE *fp;double t,h,y4,z4101;y0=0.0;y1=0.0;y2=0.0;y3=0.0;t=0.0;h=0.01;grkt1(t,y,4,h,101,z);printf(“n“);if (fp=fopen(“test1.txt“,“w“)=NULL)printf (“cannot op
41、en filen“);exit(0);for (i=0;i=100;i+)t=i*h;printf(“t=%5.2fn“,t);fprintf(fp,“%6.2fn“,z0i);for (j=0;j=3;j+)printf(“y(%d)=%e“,j,zji);printf(“n“);void grkt1f(t,y,n,d)汽车空气悬架建模与仿真分析22int n;double t,y,d;t=t;n=n;d0=y1;d1=-200*(y1-y3)-150*(y0-y2);d2=y3;d3=-100*(y3-y1)+250*(y2-y0)-4000*(y2-sin(t);return;Mabtl
42、ab 车轮位移绘图程序：fip=fopen(test1.txt,r);y=fscanf(fip,%g);fclose(fip);t=0:0.01:1.00;plot(t,y);grid on;车轮速度绘图程序fip=fopen(test2.txt,r);y=fscanf(fip,%g);fclose(fip);t=0:0.01:1.00;plot(t,y);grid on; 汽车车身位移绘图程序fip=fopen(test3.txt,r);y=fscanf(fip,%g);fclose(fip);t=0:0.01:1.00;plot(t,y);grid on;汽车车身速度绘图程序fip=fopen(test4.txt,r);y=fscanf(fip,%g);fclose(fip);t=0:0.01:1.00;plot(t,y);grid on;汽车车身加速度绘图程序fip=fopen(test5.txt,r);y=fscanf(fip,%g);fclose(fip);t=0:0.01:1.00;plot(t,y);grid on;汽车空气悬架建模与仿真分析23汽车车轮加速度绘图程序fip=fopen(test6.txt,r);y=fscanf(fip,%g);fclose(fip);t=0:0.01:1.00;plot(t,y);grid on;
