主权项 |
一种带孔薄壁旋转曲面结构的孔洞优化设计方法,其特征在于包括以下步骤:(a)建立旋转母线方程,确定极坐标系下旋转面上任一点径向坐标与轴向坐标之间的关系;设z轴为旋转轴,旋转母线位于y-z平面上,R表示极坐标系下旋转面上任一点的径向坐标;采用解析表达建立旋转母线方程表达式如下: <mrow> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mi>F</mi> <mrow> <mo>(</mo> <mi>y</mi> <mo>,</mo> <mi>z</mi> <mo>)</mo> </mrow> <mo>=</mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mi>x</mi> <mo>=</mo> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>则R与z之间的关系由(2)式确定F(sign(y)R(z),z)=0 (2)采用参数拟合方式建立的旋转母线方程表达式如下: <mrow> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mi>y</mi> <mo>=</mo> <munderover> <mi>Σ</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>N</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>v</mi> <mo>)</mo> </mrow> <mo>·</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mi>z</mi> <mo>=</mo> <munderover> <mi>Σ</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>N</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>v</mi> <mo>)</mo> </mrow> <mo>·</mo> <msub> <mi>z</mi> <mi>i</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mn>0</mn> <mo>≤</mo> <mi>v</mi> <mo>≤</mo> <mn>1</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>则R与z之间的关系由(4)式确定R(z)=|y(v(z))| (4)(b)建立给定旋转曲面的参数方程: <mrow> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mi>x</mi> <mo>=</mo> <mi>R</mi> <mrow> <mo>(</mo> <mi>z</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mi>cos</mi> <mo>[</mo> <msub> <mi>θ</mi> <mn>0</mn> </msub> <mo>+</mo> <mrow> <mo>(</mo> <msub> <mi>θ</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>θ</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mi>s</mi> <mo>/</mo> <msub> <mi>s</mi> <mn>0</mn> </msub> <mo>]</mo> </mtd> </mtr> <mtr> <mtd> <mi>y</mi> <mo>=</mo> <mi>R</mi> <mrow> <mo>(</mo> <mi>z</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mi>sin</mi> <mo>[</mo> <msub> <mi>θ</mi> <mn>0</mn> </msub> <mo>+</mo> <mrow> <mo>(</mo> <msub> <mi>θ</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>θ</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mi>s</mi> <mo>/</mo> <msub> <mi>s</mi> <mn>0</mn> </msub> <mo>]</mo> </mtd> </mtr> <mtr> <mtd> <mi>z</mi> <mo>=</mo> <msub> <mi>H</mi> <mn>0</mn> </msub> <mo>+</mo> <mrow> <mo>(</mo> <msub> <mi>H</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>H</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mi>t</mi> <mo>/</mo> <msub> <mi>t</mi> <mn>0</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mrow> <mo>(</mo> <mn>0</mn> <mo>≤</mo> <mi>s</mi> <mo>≤</mo> <msub> <mi>s</mi> <mn>0</mn> </msub> <mo>,</mo> <mn>0</mn> <mo>≤</mo> <mi>t</mi> <mo>≤</mo> <msub> <mi>t</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>式中,θ0与θ1为旋转面片的起始与终止角度,H0与H1为旋转面的最小与最大轴向坐标;旋转面映射为长为s0,高为t0的矩形区域;(c)根据式(5)确定孔心在s-t平面的位置,设置为s-t平面局部坐标系的原点;(d)建立孔周曲线的参数方程: <mrow> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mi>s</mi> <mo>=</mo> <mi>s</mi> <mrow> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mi>t</mi> <mo>=</mo> <mi>t</mi> <mrow> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mn>0</mn> <mo>≤</mo> <mi>u</mi> <mo>≤</mo> <mn>1</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>采用拟合平面控制点方式确定的孔周曲线参数方程: <mrow> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mi>s</mi> <mo>=</mo> <munderover> <mi>Σ</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>N</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> <mo>·</mo> <msub> <mi>s</mi> <mi>i</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mi>t</mi> <mo>=</mo> <munderover> <mi>Σ</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>N</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> <mo>·</mo> <msub> <mi>t</mi> <mi>i</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mn>0</mn> <mo>≤</mo> <mi>u</mi> <mo>≤</mo> <mn>1</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>其中,Ni(u)为拟合方程的基函数;则空间曲面结构上的孔周曲线的参数方程为: <mrow> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mi>x</mi> <mo>=</mo> <mi>R</mi> <mrow> <mo>(</mo> <mi>z</mi> <mrow> <mo>(</mo> <mi>t</mi> <mrow> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mi>cos</mi> <mrow> <mo>(</mo> <msub> <mi>θ</mi> <mn>0</mn> </msub> <mo>+</mo> <mrow> <mo>(</mo> <msub> <mi>θ</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>θ</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mi>s</mi> <mrow> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> <mo>/</mo> <msub> <mi>s</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mi>y</mi> <mo>=</mo> <mi>R</mi> <mrow> <mo>(</mo> <mi>z</mi> <mrow> <mo>(</mo> <mi>t</mi> <mrow> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mi>sin</mi> <mrow> <mo>(</mo> <msub> <mi>θ</mi> <mn>0</mn> </msub> <mo>+</mo> <mrow> <mo>(</mo> <msub> <mi>θ</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>θ</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mi>s</mi> <mrow> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> <mo>/</mo> <msub> <mi>s</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mi>z</mi> <mo>=</mo> <msub> <mi>H</mi> <mn>0</mn> </msub> <mo>+</mo> <mrow> <mo>(</mo> <msub> <mi>H</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>H</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mi>t</mi> <mrow> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> <mo>/</mo> <msub> <mi>t</mi> <mn>0</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mn>0</mn> <mo>≤</mo> <mi>u</mi> <mo>≤</mo> <mn>1</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>(e)通过s-t平面与空间坐标系之间的映射关系,采用壳单元在有限元分析软件中建立三维带孔薄壁旋转曲面结构的有限元模型;(f)在有限元模型的基础上施加边界条件与载荷,建立带孔薄壁旋转曲面结构的力学模型;(g)根据结构特点及受载形式,确定设计变量分布,对称结构采用二分之一孔形设定设计变量,双对称结构采用四分之一孔形设定设计变量;(h)选取孔周最大等效应力最小为优化目标,结构体积作为约束函数,设定设计变量初始值与变化范围,建立带孔薄壁旋转曲面结构孔洞优化设计问题的优化模型;(i)选取优化算法进行优化设计。 |