发明名称 一种外部承受高压的U型波纹管优化设计方法
摘要 本发明涉及一种外部承受高压的U型波纹管优化设计方法,采用基于遗传算法的自适应的优化设计方法,首先通过建立波纹管优化设计的数学模型计算波纹管几何参数、对波纹管几何参数进行优化,最后根据有限元软件对优化结果进行稳定性分析,实现了对外部承压的波纹管的优化设计,为外部承压的波纹管设计提供了一种快速可靠的方法;试验表明:经过本发明优化设计的承外压的波纹管的计算结果与试验结果吻合性很好,误差在10%以内,是一种有效可靠的承受外压的波纹管优化设计方法。
申请公布号 CN103473410B 申请公布日期 2016.06.01
申请号 CN201310404365.6 申请日期 2013.09.06
申请人 北京宇航系统工程研究所;中国运载火箭技术研究院 发明人 孙海鹏;曾斌;李德权;孙法国;史刚;余锋;石朝锋
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 中国航天科技专利中心 11009 代理人 范晓毅
主权项 一种外部承受高压的U型波纹管优化设计方法,其特征在于:包括如下步骤:步骤(一)、计算波纹管几何参数(1)、建立波纹管优化设计的数学模型以波纹管的内外径尺寸、波距、波宽、层数、波数以及原材料厚度为优化变量,以强度评价系数、临界平面失稳压力以及允许外径最大尺寸为约束条件,建立波纹管优化设计的数学模型,即在下列ST所述的条件下,变化Φ中的各个参数,使得Kx最小,所述数学模型如下:<maths num="0001"><math><![CDATA[<mrow><mi>Min</mi><mo>:</mo><msub><mi>K</mi><mi>x</mi></msub><mo>,</mo><msub><mi>K</mi><mi>x</mi></msub><mo>=</mo><mn>1.7</mn><mfrac><mrow><mi>DE</mi><msup><mi>&delta;</mi><mn>3</mn></msup><mi>Z</mi></mrow><mrow><msup><mi>h</mi><mn>3</mn></msup><msub><mi>C</mi><mi>f</mi></msub></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000378561230000011.GIF" wi="661" he="133" /></maths>Φ=(D,d,t,a,n)<maths num="0002"><math><![CDATA[<mrow><mi>ST</mi><mo>:</mo><mi>erro</mi><mn>1</mn><mo>=</mo><mfrac><msub><mi>&sigma;</mi><mn>1</mn></msub><mrow><msub><mi>C</mi><mi>w</mi></msub><mo>[</mo><mi>&sigma;</mi><mo>]</mo></mrow></mfrac><mo>&le;</mo><mn>1</mn><mo>,</mo><mi>erro</mi><mn>2</mn><mo>=</mo><mfrac><msub><mi>&sigma;</mi><mn>2</mn></msub><mrow><msub><mi>C</mi><mi>w</mi></msub><mo>[</mo><mi>&sigma;</mi><mo>]</mo></mrow></mfrac><mo>&le;</mo><mn>1</mn><mo>,</mo><mi>erro</mi><mn>3</mn><mo>=</mo><mfrac><msub><mi>&sigma;</mi><mi>p</mi></msub><mrow><msub><mi>C</mi><mi>m</mi></msub><mo>[</mo><mi>&sigma;</mi><mo>]</mo></mrow></mfrac><mo>&le;</mo><mn>1</mn></mrow>]]></math><img file="FDA0000378561230000012.GIF" wi="1332" he="133" /></maths>P<sub>si</sub>≥pD≤D<sub>max</sub>其中:K<sub>x</sub>为波纹管轴向刚度,E为弹性模量,波高<img file="FDA0000378561230000013.GIF" wi="167" he="111" />C<sub>f</sub>为位移引起的波纹管子午向薄膜应力系数,δ为原材料厚度,Φ为优化设计变量,erro1、erro2、erro3为强度评价系数,P<sub>si</sub>为临界平面失稳压力,p为波纹管外部承受压力,D<sub>max</sub>为允许外径最大尺寸,D为外径尺寸,d为内径尺寸,t为波距,a为波宽,n为波数,Z为层数,σ为内应力,σ<sub>p</sub>为内压力引起的最大等效应力,C<sub>w</sub>为焊接系数;C<sub>m</sub>为材料强度系数;(2)、优化计算波纹管几何参数(a)、根据初始定义的外径尺寸D,内径尺寸d,波距t,波宽a,波数n和层数Z计算初始内应力σ、初始的内压力引起的最大等效应力σ<sub>p</sub>和初始的波纹管轴向刚度K<sub>x</sub>,再将焊接系数C<sub>w</sub>和计算得到的所述初始内应力σ、初始的内压力引起的最大等效应力σ<sub>p</sub>、初始的波纹管轴向刚度K<sub>x</sub>代入下式计算强度评价系数erro1、erro2、erro3和平面失稳临界压力P<sub>si</sub>;<maths num="0003"><math><![CDATA[<mrow><mi>erro</mi><mn>1</mn><mo>=</mo><mfrac><msub><mi>&sigma;</mi><mn>1</mn></msub><mrow><msub><mi>C</mi><mi>w</mi></msub><mo>[</mo><mi>&sigma;</mi><mo>]</mo></mrow></mfrac><mo>,</mo><mi>erro</mi><mn>2</mn><mo>=</mo><mfrac><msub><mi>&sigma;</mi><mn>2</mn></msub><mrow><msub><mi>C</mi><mi>w</mi></msub><mo>[</mo><mi>&sigma;</mi><mo>]</mo></mrow></mfrac><mo>,</mo><mi>erro</mi><mn>3</mn><mo>=</mo><mfrac><msub><mi>&sigma;</mi><mi>p</mi></msub><mrow><msub><mi>C</mi><mi>m</mi></msub><mo>[</mo><mi>&sigma;</mi><mo>]</mo></mrow></mfrac><mo>,</mo><msub><mi>p</mi><mi>si</mi></msub><mo>=</mo><mfrac><mrow><mn>0.34</mn><mi>&pi;</mi><msub><mi>K</mi><mi>x</mi></msub></mrow><mrow><msup><mi>n</mi><mn>2</mn></msup><mi>t</mi></mrow></mfrac></mrow>]]></math><img file="FDA0000378561230000021.GIF" wi="1278" he="135" /></maths>(b)、将计算得到的强度评价系数erro1、erro2、erro3和平面失稳临界压力P<sub>si</sub>根据步骤(1)中所述的ST条件进行判断,若不满足条件,则在设定的范围内对外径尺寸D,内径尺寸d,波距t,波宽a,波数n和层数Z进行调整,返回步骤(a),若满足条件,进入步骤(c);(c)、根据如下公式计算波纹管轴向刚度K<sub>x</sub>,并作为满足约束条件的优化目标值:<maths num="0004"><math><![CDATA[<mrow><msub><mi>K</mi><mi>x</mi></msub><mo>=</mo><mn>1.7</mn><mfrac><mrow><mi>DE</mi><msup><mi>&delta;</mi><mn>3</mn></msup><mi>Z</mi></mrow><mrow><msup><mi>h</mi><mn>3</mn></msup><msub><mi>C</mi><mi>f</mi></msub></mrow></mfrac></mrow>]]></math><img file="FDA0000378561230000022.GIF" wi="325" he="127" /></maths>(d)、在设定范围内变化外径尺寸D,内径尺寸d,波距t,波宽a,波数n和层数Z,根据步骤(a)中的公式重新计算强度评价系数erro1、erro2、erro3和平面失稳临界压力P<sub>si</sub>,若满足步骤(1)中所述的ST条件,则返回步骤(c),计算波纹管轴向刚度<img file="FDA0000378561230000023.GIF" wi="83" he="62" />并与步骤(c)中的K<sub>x</sub>进行比较,二者中取最小值作为新的优化目标值;(e)、所述新的优化目标值对应的外径尺寸D,内径尺寸d,波距t,波宽a,波数n作为最终优化的波纹管设计尺寸参数;步骤(二)、对波纹管进行平面稳定性有限元分析,具体方法如下:(3)、利用有限元软件根据步骤(e)得到的波纹管设计尺寸参数建立波纹管模型,并为波纹管模型定义材料属性,所述材料属性包括弹性模量、屈服强度、破坏强度以及延伸率;(4)、设置边界条件,波纹管直边段一边为固支约束,另一边为位移约束;(5)、定义接触,定义波纹管各波之间以及各层之间为自接触约束;(6)、添加波纹管承受的外压力以及波纹管受到的位移载荷;(7)、有限元离散,将增量本构方程与应力应变演化本构方程一起组成离散增量列式,进行有限元求解;(8)、利用有限元软件进行仿真计算,得到波纹管的应力值与波距的位移值;(9)、将稳定性判定准则对波纹管的应力值与波距的位移值进行判读分析,若波纹管波谷最大变形量超过15%,则不满足要求,返回步骤(一),重新计算,若满足要求,则所述步骤(e)得到的波纹管设计尺寸参数作为最终设计依据。
地址 100076 北京市丰台区南大红门路1号内35栋