飞机带孔蒙皮局部优化设计.docx

-- -- -- -- 飞机带孔蒙皮局部应力优化设计 探究有效的方法，以减缓应力集中，提高构建的承载力量和使用寿命,对于工程实际问题格外重要。一、带中心孔平板应力集中优化设计模型参数正方形薄板 2m 探究有效的方法，以减缓应力集中，提高构建的承载力量和使用寿命, 对于工程实际问题格外重要。 一、带中心孔平板应力集中优化设计 模型参数 正方形薄板 2m×２m ,厚度 t=0.5m； 材料为硬铝，弹性模量 E＝7．2×1０１0Ｎ/m2,泊松比 μ=0.３2； 边界条件：薄板底边沿载荷方向的位移为０，薄板顶边受到均布载荷 q=1000 N/m2。 薄板圆孔应力集中探究 ２．1 无限大板圆孔应力集中的弹性力学理论 在平板宽度和圆孔直径的比值无限大时,在极坐标(?，?)中，应用弹性力学理论可 以 知 道 ： 沿圆孔边,在孔顶部和底部，即?=R 沿圆孔边,在孔顶部和底部，即?=R 和 ?=±90°，有最大应力 ,即应力 Ａnsｙs 模型计算 图错误!不能识别的开关参数。建模和网格划分 以ｒ＝0.１ｍ为例说明计算结果: 图错误!不能识别的开关参数。 r=0.1m 的应力云图 经过 anｓys 计算得到的结果如下： 半径 r／ｍ 0.５ 0.4 0.３ 0.2 0．1 0.０5 最大应力 5 3193 ２982 2852 应力集中系数 5.０1２ 4．074 3．52１ 3.193 2.９82 2.852 作出应力集中系数随半径变化的图线如下: 可见,应力集中系数随着半径的减小而渐渐减小,当半径 r０．1m 时，应力集中系数3。当孔径 r=０.1m,根本可以模拟无限大板的抱负状况了。 2.2 弹性力学理论验证 依据弹性力学理论可知,无限大板中心孔最大应力集中系数为 3,与孔径无关。下面通过连续减小孔径来验证这个结论： 半径ｒ/m 0.04 0．03 ０.02 ０.０1 ０.008 0.００6 最大应力 2９41 29３7 293３ 29３2 2９３1 2９30 应力集中系数 2.94１ ２.９37 2.93３ 2.93２ 2.93１ ２.930 作出应力集中系数随半径变化的图线如下： 可见,当半径足够小时,应力集中系数随着半径的变化根本不变,大约是 2．93，与理论值有肯定的误差,这可能是计算软件不够准确造成。在无视软件计算误差的状况下，可以认为验证 了弹性力学的结论。 平板圆孔应力集中优化探究 模型参数 正方形薄板 2m×２ｍ ，厚度 t＝2m;中心孔孔径r=0.1m; 材料为硬铝，弹性模量 E=7.２×1010Ｎ/m2,泊松比 μ=0．3２； 边界条件：薄板底边沿载荷方向的位移为 0，薄板顶边受到均布载荷ｑ=1００ ０ N／m2。 模型未优化前的应力云图如上文图 1 所示。 3．２ 单孔优化问题 工程实际中，常常通过在薄板上再打一个关心孔的方式来减小中心孔的应力集中。简洁知道， 小孔的位置和大小是影响优化效果的关键因素，并且小孔的孔径和与大孔孔心的距离之间存 在耦合。下面对角度、孔径、孔心距三个因素分开争论。 角度影响 〔此时小孔孔径r=0．05m，与中心孔孔心距为 0.5m) 以角度 120°时的计算结果说明过程: 图错误!未定义书签。角度 120°时的优化应力云图 Ａｎsys 计算结果列表如下： 半径 r/m 90° 105° １２0° 13 ５° 150° 1６5° １80° 最大应力 2945 2969 3 3０75 300９ 应力集中系数 2.945 2.969 3．045 3.１72 3.213 3.０7５ 3．００9 作出应力集中系数随角度变化的图线： 可见，当角度为 9０°时,应力集中最小，即优化效果最好。 得出结论：假设用打一个小孔的方法来优化时，应将小孔打在和大孔轴心连线与载荷平行的角度上。 孔间距的影响 (此时小孔在９０°位置,孔径为０.05m) 以孔间距０．25 为例说明计算过程: 图错误!不能识别的开关参数。 孔间距 0.2５优化时的应力云图 Anｓys 计算结果列表如下: 孔间距/m 0.５ 0.45 0．4 0.３5 ０.３ 0.25 0.2 最大应力 2945 2917 ２9１4 2９０6 290１ ２88８ 29０5 应力集中系数 2.9４5 ２ . ９ １7 ２．914 2.9０6 2.901 2.888 2.905 作出应力集中系数随孔间距的变化图线如下: 可见,当孔间距为 0．２5ｍ时,应力集中最小，即优化效果最好。 但由于此处数据间隔取得大,结果不准确，需要进一步细化步长,找出相对准确度的结果。 下面在孔间距为０.25 四周细化步长，计算结果列表如下: 孔 0.24 0. ２ ０.250 ０ .25 0.255 ０．257 0 ． 26 0.26 ０.2６ 0．2 0.2 间 ５ 47 3 ０