学期飞机结构损伤与维修课程设计报告(2).docxVIP

飞机结构疲劳强度与维修课程设计报告专业：飞行器制造工程学号：　　100146　　学生姓名：　　　所属学院：　航空工程学院　指导教师：　　　于洪　　　　二〇一三年六月一、题目——程序块载荷谱下的疲劳寿命估算图1为一中心孔LY12-CZ铝合金板，板宽W=50mm，孔直径D=8mm。名义应力谱见表1，试求其疲劳寿命。图1 中心孔板表1 载荷谱表2aLY12-CZ铝合金KT=2时的S-N曲线（单位：kg/mm2）表2bLY12-CZ铝合金KT=3时的S-N曲线（单位：kg/mm2）表2cLY12-CZ铝合金KT=4时的S-N曲线（单位：kg/mm2）表2dLY12-CZ铝合金KT=5时的S-N曲线（单位：kg/mm2）图2 中心孔板理论应力集中系数二、解答流程（1）对表1数据进行单位转换。（2）对于、，取为横坐标，为纵坐标，进行三次样条插值。得一个对应插值。（3）对于重复步骤（2），得包含步骤（2）中在内共9个对应插值。（4）以为横坐标，以步骤（3）中9个插值Sa为纵坐标进行插值，得一条、下的曲线。（5）对于N=、，重复步骤（2）、（3）、（4）得包含步骤（4）中曲线在内共10条下的与10个疲劳寿命N相对应的插值曲线。（6）对于表1中第1级中的=3.62在步骤（5）的10条的插值曲线中依次找出与10个疲劳寿命N相对应的10个值。（7）以10个疲劳寿命值N为横坐标，以步骤（6）中相对应的10个值为纵坐标进行插值，得出一条=3.62下的曲线。（8）根据表1，与=3.62相对应的值为5.87。在步骤六得出的曲线中查找出相对应的N值，至此，得出，=3.62，=5.87下对应的N值。（9）重复步骤（2）至（8），得出2、3、4、5、6、7、8级下对应的N值。（10）根据线性累积损伤理论进行疲劳寿命估算。三、题目解答（1）数据处理（a）首先注意到，题目中表1和表2中数据的单位不同。表1单位为，表2单位为。需要将表一中的数据进行转换。单位换算公式如下：所以，表1里面均应除以9.8，才与表格2数据匹配。处理结果如下：表3a未经单位转换的载荷谱级数Smax/MpaSmin/MpaSm/MpaSa/Mpa循环次数n193-2235.557.5845021411980616723220-6478142584318-12198.5219.52534227184.5157.576193096.596.5130717641108.567.5982885-5913725180表3b经单位转换的载荷谱级数Smax(Kg/mm2)Smin(Kg/mm2)Sm(Kg/mm2)Sa(Kg/mm2)循环次数n19.49-2.243.625.878450214.391.948.16645-6.537.9614.4958432.45-12.3510.0522.402534.902.7618.8316.077619.690.009.859964.1811.076.8998288.67-6.021.337.355180（b）由题目知，中心孔LY12-CZ铝合金板的板宽W=50mm，孔直径D=8mm。，在图2中找出对应的。（2）表格绘制（a）根据表2中a，b，c，d各个表格的数据，对给定的疲劳寿命，可以以为横坐标，以为纵坐标，画出对应四个值的离散的点。用matlab作图如下所示：图3N=时不同KT值下的离散点（程序见附录(1)）（b）在matlab中对进行三次样条插值，得时离散的点，然后对进行三次样条插值，求得平滑的曲线图如下所示：图4aN=时据离散点进行插值得出的曲线（程序见附录(2)）然后依据相同的步骤，对N取、做出同样的图。如下所示： 从左至右S形排列N依次为、横坐标为，纵坐标为，图例请参考图4a 图4b对应各N值时据离散点进行插值得出的曲线（程序见附录(3)）（3）求取S-N曲线（a）根据经单位转换后的载荷谱，对级数1中的值在图4b中在所有对应的曲线中找出N分别为、、时的值。（值是在matlab中精确对应的，不是在图中目测得到，图片只是便于理解，所有参数均在程序中得到）结果如下：表4N离散点N29.560323.377815.664313.538711.7274N9.91248.26107.19655.99625.4821（b）对当时所得到的值进行三次样条插值，可得曲线。曲线如下图所示：图5时插值得出的曲线（程序见附录(4)）（c）由应力谱可知，当时，。因此，可以在上图曲线中查得时对应的疲劳寿命N值为3162300次。（d）重复上述步骤a、b、c，依次获得应力谱中2、3、4、5、6、7、8级对应的疲劳寿命N值。结果如下：表5各级对应的疲劳寿命N（表中的数据可根据附录（4）中程序得出）级数1234疲劳寿命N31623