发明名称 一种加筋板结构的优化设计方法
摘要 本发明一种加筋板结构的优化设计方法,先建立加筋板结构的自由振动方程,然后标识加筋板结构的各项参数,若是不变参数,则标识为已知参数,并给定其具体参数值,将余下的参数设定为待优化参数,并给定其优化范围,并设定加筋板结构固有频率的优化目标,最后选择合适的优化算法,结合自由振动方程,对设定的待优化参数进行最优解搜索,寻找出使加筋板结构获得最接近优化目标的参数值,根据参数优化结果形成加筋板;本发明可指导任意边界条件下加筋板结构的优化设计,通过适当调整参数来获得期望的结构固有振动频率,可选取各项参数中的一种或多种作为优化参数进行优化,可自动快速寻找到各项参数的最优值,最终实现趋于用户期望值的结构优化目标。
申请公布号 CN102930137B 申请公布日期 2016.05.18
申请号 CN201210376675.7 申请日期 2012.09.29
申请人 华侨大学 发明人 欧达毅
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 泉州市文华专利代理有限公司 35205 代理人 车世伟
主权项 一种加筋板结构的优化设计方法,其特征在于具体包括如下步骤:步骤1、建立加筋板结构的自由振动方程,用于模拟任意的边界条件以及任意位置的加筋处理:使用公式<img file="FDA0000890070070000011.GIF" wi="549" he="95" />计算固有频率其中,ω为自振角频率,<img file="FDA0000890070070000012.GIF" wi="84" he="79" />为对应模态;{M}为加筋板结构的整体质量矩阵,由母板单元等效质量矩阵{M<sub>p</sub>}<sub>e</sub>和加筋条单元等效质量矩阵{M<sub>s</sub>}<sub>e</sub>所组成,{K}为加筋板结构的整体刚度矩阵,由母板单元等效刚度矩阵{K<sub>p</sub>}<sub>e</sub>、加筋条单元等效刚度矩阵{K<sub>s</sub>}<sub>e</sub>和支撑边界单元等效刚度矩阵{K<sub>b</sub>}<sub>e</sub>所组成,其中:{M<sub>p</sub>}<sub>e</sub>=ρ<sub>p</sub>h∫∫{N}<sup>T</sup>{H<sub>p</sub>}{N}dxdy,{M<sub>s</sub>}<sub>e</sub>=ρ<sub>s</sub>∫{N}<sup>T</sup>{T<sub>s</sub>'}<sup>T</sup>{H<sub>s</sub>}{T<sub>s</sub>'}{N}dl,{K<sub>p</sub>}<sub>e</sub>=∫∫{B<sub>p</sub>}<sup>T</sup>{D<sub>p</sub>}{B<sub>p</sub>}dxdy,{K<sub>s</sub>}<sub>e</sub>=∫({T<sub>s</sub>}{B<sub>p</sub>})<sup>T</sup>{D<sub>s</sub>}({T<sub>s</sub>}{B<sub>p</sub>})dl,<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mrow><mo>{</mo><msub><mi>K</mi><mi>b</mi></msub><mo>}</mo></mrow><mi>e</mi></msub><mo>=</mo><mo>&Integral;</mo><mfenced open = "(" close = ")"><mtable><mtr><mtd><mrow><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>k</mi><mrow><mi>t</mi><mi>b</mi></mrow></msub><msup><mrow><mo>{</mo><msub><mi>N</mi><mi>w</mi></msub><mo>}</mo></mrow><mi>T</mi></msup><mo>{</mo><msub><mi>N</mi><mi>w</mi></msub><mo>}</mo><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>k</mi><mrow><mi>r</mi><mi>b</mi></mrow></msub><msup><mrow><mo>{</mo><mfrac><mrow><mo>&part;</mo><msub><mi>N</mi><mi>w</mi></msub></mrow><mrow><mo>&part;</mo><msub><mover><mi>n</mi><mo>&RightArrow;</mo></mover><mi>b</mi></msub></mrow></mfrac><mo>}</mo></mrow><mi>T</mi></msup><mo>{</mo><mfrac><mrow><mo>&part;</mo><msub><mi>N</mi><mi>w</mi></msub></mrow><mrow><mo>&part;</mo><msub><mover><mi>n</mi><mo>&RightArrow;</mo></mover><mi>b</mi></msub></mrow></mfrac><mo>}</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>k</mi><mrow><mi>n</mi><mi>b</mi></mrow></msub><mo>&lsqb;</mo><msup><mrow><mo>{</mo><msub><mi>N</mi><mi>u</mi></msub><mo>}</mo></mrow><mi>T</mi></msup><mo>{</mo><msub><mi>N</mi><mi>u</mi></msub><mo>}</mo><msup><mrow><mo>(</mo><mover><mi>x</mi><mo>&RightArrow;</mo></mover><mo>&CenterDot;</mo><msub><mover><mi>n</mi><mo>&RightArrow;</mo></mover><mi>b</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>{</mo><msub><mi>N</mi><mi>v</mi></msub><mo>}</mo></mrow><mi>T</mi></msup><mo>{</mo><msub><mi>N</mi><mi>v</mi></msub><mo>}</mo><msup><mrow><mo>(</mo><mover><mi>y</mi><mo>&RightArrow;</mo></mover><mo>&CenterDot;</mo><msub><mover><mi>n</mi><mo>&RightArrow;</mo></mover><mi>b</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>&rsqb;</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>k</mi><mrow><mi>m</mi><mi>b</mi></mrow></msub><mo>&lsqb;</mo><msup><mrow><mo>{</mo><msub><mi>N</mi><mi>u</mi></msub><mo>}</mo></mrow><mi>T</mi></msup><mo>{</mo><msub><mi>N</mi><mi>u</mi></msub><mo>}</mo><msup><mrow><mo>(</mo><mover><mi>x</mi><mo>&RightArrow;</mo></mover><mo>&CenterDot;</mo><msub><mover><mi>n</mi><mo>&RightArrow;</mo></mover><mi>&Gamma;</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>{</mo><msub><mi>N</mi><mi>v</mi></msub><mo>}</mo></mrow><mi>T</mi></msup><mo>{</mo><msub><mi>N</mi><mi>v</mi></msub><mo>}</mo><msup><mrow><mo>(</mo><mover><mi>y</mi><mo>&RightArrow;</mo></mover><mo>&CenterDot;</mo><msub><mover><mi>n</mi><mo>&RightArrow;</mo></mover><mi>&Gamma;</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>&rsqb;</mo></mrow></mtd></mtr></mtable></mfenced><msub><mi>d&Gamma;</mi><mi>b</mi></msub><mo>,</mo></mrow>]]></math><img file="FDA0000890070070000013.GIF" wi="1389" he="486" /></maths>其中,下标p,b,s分别对应于母板、支撑边界和加筋条,下标u,v,w分别对应面内法向位移、面内切向位移和板中面扰度,ρ为材料密度,h为母板厚度,{N}为单元形函数,{N}<sup>T</sup>为{N}的转置运算矩阵,{B}为应变矩阵,{D}为抗弯刚度矩阵,{H}为关于h的对角矩阵函数,{T<sub>s</sub>}和{T<sub>s</sub>'}为变换矩阵,是加筋条对应局部坐标轴与全局坐标轴的夹角函数,通过这两个变换矩阵可实现对任意位置加筋条的模拟;l为加筋条单元轴,k<sub>tb</sub>,k<sub>rb</sub>,k<sub>nb</sub>和k<sub>mb</sub>为支撑边界的参数,分别代表加筋板结构支撑边界的横向刚度、转动刚度、面内法向刚度和面内切向刚度,这些参数可为复数形式,实部表示边界支撑的弹性特性,虚部表示边界支撑的阻尼特性,模拟弹性边界支撑时,虚部值取0,通过改变这些参数可实现对任意边界条件的模拟,这些参数可以是常数也可以是位置的函数;<img file="FDA0000890070070000021.GIF" wi="60" he="70" />和<img file="FDA0000890070070000022.GIF" wi="68" he="71" />为单元边界轮廓Γ<sub>b</sub>的单位法向量和单位切向量;上述母板单元等效质量矩阵{M<sub>p</sub>}<sub>e</sub>、加筋条单元等效质量矩阵{M<sub>s</sub>}<sub>e</sub>、母板单元等效刚度矩阵{K<sub>p</sub>}<sub>e</sub>、加筋条单元等效刚度矩阵{K<sub>s</sub>}<sub>e</sub>、支撑边界单元等效刚度矩阵{K<sub>b</sub>}<sub>e</sub>由加筋板结构的各项参数推导获得,该参数包括:母板形状、母板尺寸、母板材料、边界条件、加筋条形状、加筋条数量、加筋条材料、加筋条尺寸、加筋位置;步骤2、标识加筋板结构的各项参数,其中参数包括:母板形状、母板尺寸、母板材料、边界条件、加筋条形状、加筋条数量、加筋条材料、加筋条尺寸、加筋位置,若是不变参数,则标识为已知参数,并给定其具体参数值,将余下的加筋板结构参数设定为待优化参数,并给定其优化范围;步骤3、设定加筋板结构固有频率的优化目标,设定优化目标可以采用以下两种方式:方式一:直接设定固有频率的具体优化目标值;方式二:不设定具体优化目标值,而是要求加筋板结构获得最大固有频率;步骤4、选择合适的优化算法,结合步骤1的自由振动方程,对步骤2设定的待优化参数进行最优解搜索,寻找出使加筋板结构获得最接近优化目标的参数值;步骤5、根据步骤4的参数优化结果形成加筋板。
地址 362000 福建省泉州市丰泽区城东华侨大学