北航最新-飞行器设计-课程设计报告

飞机带孔蒙皮局部应力优化报告专业：飞行器设计学号：39051623姓名：黄星指导老师：张铮2012年9月25日一、设计课程题目飞机带孔蒙皮局部应力优化设计二、研究对象飞机带孔蒙皮三、设计目的综合运用有关基础理论、专业知识和实际经验，独立地解决专业范围内比较简单的具有典型性的设计任务，为毕业设计以及毕业后在专业工作解决更全面而复杂的技术问题打好基础。四、研究内容１、矩形板和孔的位置与形状：设计说明：在一定载荷P下，构件宽度、孔径和空边应力集中系数的关系：在载荷、板宽和孔径都不变的条件下，沿板构件的纵轴线再打一个孔，孔的位置和孔径大小对原孔孔边应力集中系数的影响；进一步，可以再打第二个孔、第三个孔…再进一步，孔可以不打在纵轴线上，如何设计孔的位置和孔径大小？２、梯形板形状：设计说明：当载荷不变，板构件形状改变时(如错误!未找到引用源。所示)，一个孔及多个孔在考虑上述应力集中条件下的设计，其中，板构件的宽端尺寸不变时，窄端尺寸与应力集中系数的关系？宽孔径孔径３、双向载荷长圆孔：设计说明：如板构件受到双向拉力，纵向载荷是横向载荷的2倍（这是机舱段机壳常规的受载情况），原圆孔改为长圆孔（即原圆孔沿横向直径隔开，加入一等宽矩形段，如错误!未找到引用源。所示，这是机窗的基本形式），如何设计孔径和矩形边长，实现长圆孔周边等周向（切向）应力（或基本等切向应力）？五、实验环境ANSYS13有限元分析软件，模拟真实条件的应力状态。软件所设的各种参数：单元类型：QUAD8NODE183单元设置：PLANESTRSＷ／ＴＨＫ设定杨氏模量：E=2*105μ=0.3板及孔的长度单位为mm应力单位为MPa长六、实验过程与结果（一）矩形板构件：１、模拟无限大平板模型为１0０ｘ２００孔位于中心（0，0），初始孔径大小20加载：底边约束Ｙ方向的约束，自由端加载-１的均布载荷孔径大小为自变量，从20开始往下逐渐减小，仔细观察构件的应力分布图及读取孔边最大应力值中心孔半径2015108642应力集中系数3.750483.38843.165353.108753.06663.036643.02039因为半径小于6时，应力集中系数的变化率小于1%，故近似认为r小于等于6时，孔径对圆孔应力的影响忽然不计，此时可把100*200的平板看作是无限大的。而且r=6时，应力集中系数为3.0666，非常接近理论值，说明实验时的网格划的足够密，之后的模拟实验都是在此基础上进行的。要想减小孔边应力集中系数，我们可以采取改变中心孔的应力场分布的方法。从单孔的应力云图可以看出，在孔边沿载荷方向应力会出现很大的降低，我们从中可以得到启发：在中心孔的附近打孔，而且我们能够预测，在竖直方向打孔的效果将会是最好的。下面是各种打孔情况的实验。划分网格时，每一个孔边平均分为160份，矩形板长边平均分为40份，短边平均分为20份，既保证孔边1附加孔在斜方向2附加孔在水平方向3附加孔在竖直方向上面的实验很好的符合了我们预测，只有在竖直方向打孔才能很好的降低中心孔边的应力集中系数，显然，由对称的思想可知，在中心孔竖直方向对称打孔可以更好的减小中心孔的应力集中系数，所以下面我们具体讨论附加孔对称分布在竖直方向时，附加孔孔径及与中心孔孔心距对中心孔的应力集中系数的影响。对称打孔应力云图实验采用ANSYS的优化设计算法，设置两个参数附加孔半径R和孔心距Y，设置R的范围为（1-5.9），因为孔径太小时对中心孔应力集中系数的影响可以忽略，但是附加孔径大于6的话，由一个孔的模拟实验我们可以得出孔径越大，应力集中系数也越大的规律，此时附加孔的应力集中将比中心孔严重，不符合题目要求。Y的范围为（7.5-36）Y太大时由圣维南原理可知附加孔对中心孔的影响也可不计。Z为中心孔边与附加孔边的最小距离。设置将每次实验的最大应力赋给变量PAR，PAR的收敛精度设置为0.005采用一阶算法，循环次数设置为30。第一次优化结果如下：第一次优化应力集中系数与孔心距及孔半径的关系图第一次优化应力集中系数、孔心距及孔半径与循环系列号的关系图第一次优化各参数数据N123456789Z619.044316.299115.789716.580314.56587.653.64363.6412R34.66875.43595.56335.14874.75125.7045.02644.9418Y1529.71327.73527.35327.72925.31719.35414.6714.583PAR2.88972.80692.76042.76442.75752.7372.6352.51992.5428N101112131415161718Z3.4833.43383.52563.68363.30153.09572.16094.78894.534R4.9364.91324.88144.83544.81054.77834.73915.01314.971Y14.41914.34714.40714.51914.11213.87412.915.80215.505PAR2.54112.54872.5522.5562.55272.58022.59272.54082.5426N1920212223242526Z4.38394.14213.48292.89062.3392.59752.6292-1.5963R4.94414.90594.87715.09645.2435.08455.03785.1528Y15.32815.04814.3613.98713.58213.68213.6679.5565PAR2.55162.56042.55882.50352.51122.50592.50582.4067从上表可以看出，第26次循环的结果为最佳，但是发现此时中心孔已经被附加孔破坏，也不符合要求。增加一状态变量：中心孔边与附加孔边的最小距离Z。设置Z的范围为（0.2-24），这样就能保证附加孔不破坏中心孔了。进行第二次优化，优化数据如下：第二次优化计算结果第二次优化应力集中系数与孔心距及孔半径的关系图第二次优化应力集中系数、孔心距及孔半径与循环系列号的关系图第二次优化各参数数据N12345678Z612.362912.847212.784612.515920.085612.7676.7209R33.62713.62583.57043.49914.57245.84.9931Y1521.9922.47322.35522.01530.65824.56717.714PAR2.88972.82152.82422.83452.85322.832.752.5695N910111213141516Z6.28455.51565.47074.60074.40784.2934.59774.5051R4.88854.83044.80734.74434.72824.6834.96235.1329Y17.17316.34616.27815.34515.13614.97615.5615.638PAR2.57592.57882.58852.58162.5872.60142.55192.5538N1718192021222324Z2.01312.07912.11732.06991.55140.74950.33238.059R5.29295.36695.46775.54615.67265.38755.27175.058Y13.30613.44613.58513.61613.22412.13711.60419.117PAR2.51882.52432.52812.54092.54722.49932.47982.6258N2526272829303132Z8.08986.16684.97015.1715.13125.10765.0235R5.16624.88824.76394.7214.65184.63644.6105Y19.25617.05515.73415.89215.78315.74415.634PAR2.61392.57412.59092.60062.60722.61832.619从图和表中我们可以看出22、23组循环系列得到的结果很接近，且为所有最优系列的最小值。因为我们不能穷尽所有点而得到最佳解，故我们认为最在区域R(5.27175.3875);Y(11.60412.137)，此时能获得最优的应力集中系数。通过选取在些区域里的点进行实验，发现应力集中系数也在上述范围内，与我们的结论很好的相符。部分命令流如下CYL4,0,0,6/设置中心孔位置与大小CYL4,0,Y,R/设置附加孔位置与大小PLNSOL,S,EQV,0,1.0/显示等效应力云图*GET,PAR,PLNSOL,0,MAX/取出最大应力值并赋给目标变量PAROPVAR,R,DV,1,5.9,,/设置R的范围为1-5.9OPVAR,Y,DV,7.5,36,,/设置Y的范围为7.5-36OPVAR,Z,SV,0.2,24,,/设置Z的范围为0.2-24OPVAR,PAR,OBJ,,,0.005,/设置PAR收敛精度为0.005（二）梯形板问题在减少梯形窄边宽度的时候会引起梯形板上圆孔边应力增大，但是却能有效减小板的重量，当这个应力增量还在我们应许范围内的时候，我们却能减少大量的材料，从而达到使板重量减轻，我们用重量相对于初始情况下的减小量与应力相对初始情况下的增加量之比来衡量结果的优化程度，该值越大则表明增加相同的应力的情况下减重越明显。首先对方形板的情况进行实验，获得初始值。平面板的大小为200*200，孔半径为10.网格划分：孔边平均分为200份，与孔相连的边按0.08的比例各分为80份，其它各边则平均分为80份。网格如下图方形板网格初始应力云图梯形窄边长为128时应力云图当梯形窄边长一直减小时，我们可以预测最大应力将出现在窄边处而不是在孔边，这种情况也不是我们想要的。当窄边为46时，最大应力在窄边下，如下图：梯形窄边长为46时应力云图故使用ANSYS优化设计时，设计变量X(窄边的半宽度)应该限制在24-100之间，为了保证所加载荷不变，故在窄边加的压力值设为变量-100/X。为了不使窄边宽度很小的情况下其网格过密，设置窄边网格数为变量40*X/100。为了得到重量相对于初始情况下的减小量与应力相对初始情况下的增加量之比，我们设置变量：PAR/为最大应力值；S/为相对初始情况下的面积减少值；T/面积减少值与初始面积之比；DPAR/应力增大值；W/重量相对于初始情况下的减小量与应力相对初始情况下的增加量之比。采用零阶算法，步长为1，进行优化。得到71组数据，如下表。梯形板优化数据表X2526272829303132PAR4.07894.05154.0233.99424.0294.00233.9683.9936W1.705271.747691.796311.850191.731871.77681.847571.74909X3334353637383940PAR3.94793.95813.94723.95883.91193.89283.87883.8964W1.853691.795791.800521.739291.853371.887341.905791.8148X4142434445464748PAR3.89843.88143.87353.9073.82923.83653.83493.8239W1.777971.802991.798291.661871.893521.831811.803641.81002X4950515253545556PAR3.80543.84113.8383.77323.76423.7543.76113.7379W1.846651.680251.656951.868321.868311.873811.801981.86576X5758596061626364PAR3.72363.70653.70123.69573.6893.67333.66163.6603W1.891841.934851.917091.899431.8876