《万方数据-数字化期刊群》全文上网期刊
CNKI《中国学术期刊(网络版)》全文收录期刊
《中文科技期刊数据库》(维普网)全文收录期刊
超星期刊域出版平台、博看网全文收录期刊
日本科学技术振兴机构数据库收录

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

考虑局部模态约束的卫星结构优化

刘嫣洁 杨苗苑 李子轩

刘嫣洁, 杨苗苑, 李子轩. 考虑局部模态约束的卫星结构优化[J]. 航天器环境工程, 2023, 40(6): 611-616 doi:  10.12126/see.2023077
引用本文: 刘嫣洁, 杨苗苑, 李子轩. 考虑局部模态约束的卫星结构优化[J]. 航天器环境工程, 2023, 40(6): 611-616 doi:  10.12126/see.2023077
LIU Y J, YANG M Y, LI Z X. Structural optimization for satellites with a consideration of local modal constraints[J]. Spacecraft Environment Engineering, 2023, 40(6): 611-616 doi:  10.12126/see.2023077
Citation: LIU Y J, YANG M Y, LI Z X. Structural optimization for satellites with a consideration of local modal constraints[J]. Spacecraft Environment Engineering, 2023, 40(6): 611-616 doi:  10.12126/see.2023077

考虑局部模态约束的卫星结构优化

doi: 10.12126/see.2023077
详细信息
    作者简介:

    刘嫣洁,工程师,主要从事卫星机械总体设计与结构设计的研究

  • 中图分类号: V415.2; O224

Structural optimization for satellites with a consideration of local modal constraints

  • 摘要: 在卫星结构设计中,基于传统的全局频率约束问题的结构优化方法不能直接应用于模式置换现象的局部频率约束问题。为此提出一种近似方法:首先假设在优化过程中结构整体质量矩阵的变化可以忽略不计,且指定局部频率的模态振型在优化前后是一致的;然后基于这2个假设,将原始局部频率约束的优化问题转化为局部区域节点位移约束的优化问题,节点位移约束由当前局部频率和期望局部频率之间的比值确定。数值算例和实际工程优化实例验证了该方法的合理、有效。
  • 图  1  近似优化模型优化过程

    Figure  1.  Optimization process of the approximation optimization model

    图  2  结构板的初始有限元模型

    Figure  2.  Initial FEM of the structural plate

    图  3  初始结构板的局部振型

    Figure  3.  Local modal shape of the initial structural plate

    图  4  静力学分析得到的初始结构板节点位移

    Figure  4.  Nodal displacements of the initial structural plate obtained from static analysis

    图  5  加强筋加固方式及加强筋横截面尺寸

    Figure  5.  Reinforcement method and cross section dimension of the reinforcing stiffener

    图  6  结构板近似模型迭优化代历程

    Figure  6.  Iteration history for approximation model of the structural plate

    图  7  优化前后结构板的局部模态振型

    Figure  7.  Local modal shapes of the structural plate before and after optimization

    图  8  卫星的初始局部模态振型

    Figure  8.  Initial local modal shape of the satellite

    图  9  静力学分析得到的卫星节点位移

    Figure  9.  Nodal displacements of the satellite obtained by static analysis

    图  10  卫星结构近似模型优化迭代历程

    Figure  10.  Iteration history for approximation model of the satellite structure

    图  11  优化前后卫星的局部模态振型

    Figure  11.  Local modal shapes of the satellite before and after optimization

    表  1  结构板近似模型优化结果

    Table  1.   Optimization results of approximation model of the structural plate

    参数初始设计近似优化结果
    设计变量H/m0.010.041 4
    W/m0.010.011 7
    t1/m0.0020.001
    t2/m0.000 20.001
    局部振动模态频率/Hz22.50459.136
    局部振动模态的阶数33
    节点最大位移/m4.140.47
    局部频率相对误差/%1.44
    结构总质量/kg125.1125.3
    下载: 导出CSV

    表  2  卫星结构近似模型优化结果

    Table  2.   Optimization results of the satellite structural approximation model

    参数初始设计近似优化结果
    设计变量H/m0.010.051 8
    W/m0.010.002 6
    局部振动模态频率/Hz42.23459.667
    局部振动模态的阶数4669
    节点最大位移/m0.0590.022
    局部频率相对误差/%0.54
    结构总质量/kg8 188.38 194.8
    下载: 导出CSV
  • [1] 朱继宏, 张卫红, 邱克鹏. 结构动力学拓扑优化局部模态现象分析[J]. 航空学报, 2006, 27(4): 619-623

    ZHU J H, ZHANG W H, QIU K P, Local modal analysis of structural dynamic topology optimization[J]. Acta Aeronautica ET Astronautica Sinica, 2006, 27(4): 619-623
    [2] RAMANA G. Structural optimization with frequency constraints: a review[J]. AIAA Journal, 1993, 31(12): 2296-2303 doi:  10.2514/3.11928
    [3] KHAN M R, WILLMERT K D. An efficient optimality criterion method for natural frequency constrained structures[J]. Computers and Structures, 1981, 14(5): 501-507
    [4] VENKAYYA V B, TISCHLER V A. Optimization of structures with frequency constraints[J]. American Society of Mechanical Engineers, 1983, 54: 239-259
    [5] VANDERPLAATS G N, SALAJEGHEH E. An effective approximation technique for frequency constraints in frame optimization[J]. Intercalation Journal for Numerical Methods in Engineering, 1988, 26(5): 1057-1069 doi:  10.1002/nme.1620260505
    [6] BENDSOE M P, SIGMUDE O. Topology optimization: theory, methods and applications[M]. Berlin: Springer, 2004: 1213-1214
    [7] LI J H Y, CHEN S Y, HUANG H. Topology optimization of continuum structure with dynamic constraints using mode identification[J]. Journal of Mechanical Science and Technology, 2015, 29(4): 1407-1412 doi:  10.1007/s12206-015-0311-y
    [8] 叶红玲, 沈静娴, 隋允康. 频率约束的三维连续体结构动力拓扑优化设计[J]. 力学学报, 2012, 44(6): 1037-1045

    YE H L, SHEN J X, SUI Y K. Research on dynamic topology optimization of three dimension continuum structures with frequency constraints[J]. Acta Mechanical Sinica, 2012, 44(6): 1037-1045
    [9] DEATON J D, GRANDHI R V. A survey of structural and multidisciplinary continuum topology optimization[J]. Structural and Multidisciplinary Optimization, 2014, 49: 1-38 doi:  10.1007/s00158-013-0956-z
    [10] CHAUR M C, RICHARD P S. Identification of local modes using modal mass distributions[C]//International Modal Analysis Conference. Las Vegas, NV, USA, 1989: 772-776
    [11] TENEK L, HAGIWARA I. Eigen-frequency maximization of plates by optimization of topology using homogenization and mathematical programming[J]. Journal of Materials Science: Materials in Electronics, 1994, 37(4): 667-677
    [12] PEDERSEN N L. Maximization of eigenvalues using topology optimization[J]. Structural and Multidisciplinary Optimization, 2000, 20(1): 2-11 doi:  10.1007/s001580050130
    [13] CHENG G B, WANG B. Constraint continuity analysis approach to structural topology optimization with frequency objective constraint[C]//Proceedings of 7th Word Congress on Structural and Multidisciplinary Optimization. Seoul, 2007
    [14] XU B, JIANG J S, TONG W H, et al. Topology group concept for truss topology optimization with frequency constraints[J]. Journal of Sound and Vibration, 2003, 261(5): 911-925 doi:  10.1016/S0022-460X(02)01021-0
  • [1] 金玮安, 张磊.  大热流背景下蜂窝板结构参数的多目标优化 . 航天器环境工程, 2023, 40(2): 128-133. doi: 10.12126/see.2022139
    [2] 张少渤, 李圣山, 田晓景, 韩毅, 张志坚, 陈倬, 任清峰.  面向批产的小卫星总装流程优化 . 航天器环境工程, 2023, 40(5): 575-580. doi: 10.12126/see.2023053
    [3] 朱泽昊, 周徐斌, 刘兴天, 杜冬, 冯彦军.  控制力矩陀螺微振动时的固有频率波动机理研究 . 航天器环境工程, 2023, 40(3): 233-240. doi: 10.12126/see.2022117
    [4] 裴彦伟, 谷松, 赵春娟, 赵相禹.  空间多载荷高精度拼接支撑结构的优化设计 . 航天器环境工程, 2022, 39(1): 76-82. doi: 10.12126/see.2022.01.011
    [5] 王聪, 叶田园, 赵越阳, 常正实.  电压和频率对CO2介质阻挡放电特性影响的实验研究 . 航天器环境工程, 2020, 37(6): 610-616. doi: 10.12126/see.2020.06.012
    [6] 钱志英, 韩世泽, 马为佳, 余快, 王晓姝, 高行素.  航天器振动试验中的频率漂移现象研究 . 航天器环境工程, 2018, 35(4): 342-347. doi: 10.12126/see.2018.04.006
    [7] 周兴广, 田振强, 王冰, 秦亚明, 贺智国.  三轴振动试验系统工作平台的谐振频率分析 . 航天器环境工程, 2017, 34(1): 35-39. doi: 10.12126/see.2017.01.006
    [8] 李太平1,2, 翁海宽1,2, 江浩1,2, 洪岩1, 齐晓军1,2预紧力对系统频率漂移的影响 . 航天器环境工程, 2017, 34(6): 636-641. doi: 10.12126/see.2017.06.011
    [9] 王东升1, 任万发1, 刘青林1, 姚海艳2振动试验夹具共振频率设计要求研究 . 航天器环境工程, 2014, 31(1): 37-41. doi: 10.12126/see.2014.01.007
    [10] 孙毅, 唐民, 于庆奎.  典型卫星轨道的位移损伤剂量计算与分析 . 航天器环境工程, 2013, 30(5): 487-492.
    [11] 陈 靖1, 张 翔2, 陈卫东2, 史瑾文3.  微小卫星解锁分离装置主结构设计分析及优化 . 航天器环境工程, 2012, 29(6): 681-686.
    [12] 陈 磊, 洪建忠, 杨永生, 蒋春梅, 李心耀.  相对变形约束下某离心机吊篮拓扑优化设计 . 航天器环境工程, 2012, 29(1): 100-103.
    [13] 卫洪涛, 孔宪仁, 王本利, 张相盟.  非线性连接结构对一个典型卫星频率漂移的影响 . 航天器环境工程, 2012, 29(3): 297-303.
    [14] 韩 潇, 祁 妍, 刘波涛.  真空容器大门法兰结构设计及优化 . 航天器环境工程, 2010, 27(4): 493-495.
    [15] 付仕明, 徐小平, 裴一飞.  空间站集成全局热数学模型的建模和分析 . 航天器环境工程, 2010, 27(1): 75-79.
    [16] 张逸波.  卫星振动压环优化设计 . 航天器环境工程, 2009, 26(1): 47-50.
    [17] 代 福, 熊胜明, 龚自正.  重复频率脉冲激光辐照下光学薄膜元件 温升的有限元分析 . 航天器环境工程, 2009, 26(6): 510-513. doi: 10.12126/see.2009.06.003
    [18] 代 福, 龚自正, 牛锦超, 等.  高重复频率DPL激光辐照光学薄膜元件温升规律实验研究 . 航天器环境工程, 2009, 26(3): 232-235.
    [19] 空间碎片防护结构设计优化技术研究进展 . 航天器环境工程, 2005, 22(6): 330-334.
    [20] 支架结构随机振动响应优化研究 . 航天器环境工程, 2005, 22(01): 38-40.
  • 加载中
图(11) / 表ll (2)
计量
  • PDF下载量:  10
  • 文章访问数:  61
  • HTML全文浏览量:  7
文章相关
  • 中图分类号:  V415.2; O224
  • 收稿日期:  2023-05-22
  • 修回日期:  2023-11-27
  • 网络出版日期:  2023-12-28
  • 刊出日期:  2023-12-25

考虑局部模态约束的卫星结构优化

doi: 10.12126/see.2023077
    作者简介:

    刘嫣洁,工程师,主要从事卫星机械总体设计与结构设计的研究

  • 中图分类号: V415.2; O224

摘要: 在卫星结构设计中,基于传统的全局频率约束问题的结构优化方法不能直接应用于模式置换现象的局部频率约束问题。为此提出一种近似方法:首先假设在优化过程中结构整体质量矩阵的变化可以忽略不计,且指定局部频率的模态振型在优化前后是一致的;然后基于这2个假设,将原始局部频率约束的优化问题转化为局部区域节点位移约束的优化问题,节点位移约束由当前局部频率和期望局部频率之间的比值确定。数值算例和实际工程优化实例验证了该方法的合理、有效。

English Abstract

刘嫣洁, 杨苗苑, 李子轩. 考虑局部模态约束的卫星结构优化[J]. 航天器环境工程, 2023, 40(6): 611-616 doi:  10.12126/see.2023077
引用本文: 刘嫣洁, 杨苗苑, 李子轩. 考虑局部模态约束的卫星结构优化[J]. 航天器环境工程, 2023, 40(6): 611-616 doi:  10.12126/see.2023077
LIU Y J, YANG M Y, LI Z X. Structural optimization for satellites with a consideration of local modal constraints[J]. Spacecraft Environment Engineering, 2023, 40(6): 611-616 doi:  10.12126/see.2023077
Citation: LIU Y J, YANG M Y, LI Z X. Structural optimization for satellites with a consideration of local modal constraints[J]. Spacecraft Environment Engineering, 2023, 40(6): 611-616 doi:  10.12126/see.2023077
    • 随着机械复杂性和动力学性能要求的显著提高,动力学优化已成为结构设计中的关键问题,尤其是对于航天器结构而言[1]。如果卫星结构的固有频率与运载火箭发动机的激励频率相同,则在发射阶段会发生共振,而通过卫星结构传递的高振幅振动将会对星内有效载荷造成破坏。因此,卫星结构的固有频率必须高于目标值以避免发生共振,从而提高结构自身和有效载荷的安全性[2]。结构优化方法已被应用于解决实际卫星结构设计中的频率约束问题[3-6]。然而,大多数研究考虑的是全局频率约束,例如一阶横向和一阶纵向固有频率,很少有研究关注到对于改善关键有效载荷安装区域的局部结构刚度同样非常重要的局部频率约束。

      局部频率约束和全局频率约束的优化问题之间没有本质区别[7]。然而,由于结构尺寸的改变而引起的模态置换是局部频率约束优化问题的主要难点[8]:自动识别某一确定的局部频率难度很大,且其阶数在结构优化过程中通常不是恒定的,从而导致优化程序收敛困难[9-10],使得确定局部频率约束的优化问题的最优值变得极为复杂。因此,全局频率约束的结构优化方法不能直接应用于局部频率约束的结构优化问题。近几十年中,关于局部模态领域的研究较多:Tenek和Hagiwara[11]将单元密度限制在较大的阈值之上以避免局部模态,这就是动态结构频率最大化设计中材料布局的优化结果中出现明显灰色区域的原因。Pedersen[12]提出一种改进的可变惩罚因子方法来避免SIMP(solid isotropic microstructure with penalization)的局部模态。Cheng和Wang[13]提出一种利用多项式惩罚低密度区域单元刚度矩阵的方法。然而,目前的研究大多集中在连续体结构的拓扑优化方面,对于尺寸优化问题的研究较少。此外,目前关于局部模态识别的参考方法很少。

      本文提出一种解决局部频率约束的结构优化问题的近似方法,其基本思想是将具有局部频率约束的原始结构动力学优化问题转化为具有节点位移约束的近似结构静力学优化问题,以回避在优化过程中使用局部模态阶数和局部模态识别。方法仍以最小化结构质量为目标建立近似优化模型,并将横截面尺寸作为设计变量,将局部区域在一定局部频率下的节点位移作为约束;节点位移约束由局部区域的初始位移以及初始局部频率与期望局部频率之间的比值导出。最后给出数值算例和实际的卫星结构优化问题验证所提出方法的合理、有效性。

    • 具有频率约束的传统动力学结构优化问题可以表述为非线性数学规划问题,即

      $$ \left\{\begin{array}{l}{\mathrm{find}} :{\boldsymbol{X}}={\left\{{x}_{1},{x}_{2},\cdot\cdot\cdot ,{x}_{n}\right\}}^{\mathrm{T}}\\ {\mathrm{min}} :F\left(\boldsymbol{X}\right)={\text {Weight}}\\ {\mathrm{s}}.{\mathrm{t}}.\;\;g\left(\boldsymbol{X}\right)={\tilde{\omega }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}-{\omega }_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}{\text{≤}} 0\\ {x}_{i}^{\mathrm{l}}{\text{≤}} {x}_{i}{\text{≤}} {x}_{i}^{\mathrm{u}}\left(i=\mathrm{1,2},\cdot\cdot\cdot ,n\right)\end{array}\right. 。 $$ (1)

      式中:xi (i=1, 2, ···, n)为设计变量;F(X) 为结构质量(Weight);g(X) 为频率约束;ωlocal为当前设计结构的局部频率;$ {\tilde{\omega }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}} $为局部频率的期望值;$ {x}_{i}^{\mathrm{l}} $$ {x}_{i}^{\mathrm{u}} $分别为xi的下限和上限。

    • 为了提高局部结构刚度,将在局部区域增加加强筋以提高局部模态频率。为了近似转换原始优化问题,提出以下2个必不可少的假设:

      1)由于加强筋引起的质量增加与结构的总质量相比非常小,所以在优化过程中结构整体质量矩阵的变化可以忽略不计;

      2)指定局部频率的模态振型在优化前后是一致的,即$ {\varnothing }_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}^{0}=r\cdot {\varnothing }_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}} $,其中,$ {\mathit{\varnothing }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}^{0} $$ {\mathit{\varnothing }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}} $分别为优化前(初始结构)和优化后(当前结构)的模态振型。

      结构修改对单元刚度影响重大,其变化不应忽视[14]。因此,假设刚度矩阵是优化问题中设计变量的函数,即

      $$ {\boldsymbol{K}}={\boldsymbol{K}}({\boldsymbol{X}})。$$ (2)

      考虑到可以将初始结构的局部模态振型视为在一定静载荷下的节点位移,本文根据模态振型定义了静载荷$ \tilde{F} $local

      $$ {\tilde{F}}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}=\boldsymbol{K}\left({\boldsymbol{X}}_{0}\right)\cdot {\boldsymbol{\varnothing }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}^{0}=\boldsymbol{K}\left({\boldsymbol{X}}_{0}\right)\cdot {\boldsymbol{u}}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}^{0} \text{,} $$ (3)

      式中:K(X0)为初始结构的单元刚度矩阵;$ {\boldsymbol{u}}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}^{0} $为初始结构的节点位移。

      基于假设2),静载荷也可以表示为

      $$ {\tilde{F}}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}=\boldsymbol{K}\left({\boldsymbol{X}}_{0}\right)\cdot r\cdot {\boldsymbol{\varnothing }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}=\boldsymbol{K}\left(\boldsymbol{X}\right)\cdot {\boldsymbol{u}}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}} \text{,} $$ (4)

      式中:K(X)为当前结构的单元刚度矩阵;$ \boldsymbol{u} $local为当前结构的节点位移。

      无阻尼结构系统的特征值问题可以表述为

      $$ (\boldsymbol{K}-{{\omega }_{i}}^{2}\cdot \boldsymbol{M})\cdot {\boldsymbol{\varnothing }}_{i} = 0\text{,} $$ (5)

      式中:K为刚度矩阵;M为质量矩阵;$ {\mathit{\varnothing }}_{i} $为第i阶特征向量或固有模态;$ {\omega }_{i} $为相应的第i阶特征值或固有频率。通过求解该特征值问题,可以计算出结构的第i阶固有频率,

      $$ {{\omega }_{i}}^{2}=\frac{{\boldsymbol{\varnothing }}_{i}^{\mathrm{T}}\cdot \boldsymbol{K}{\cdot \boldsymbol{\varnothing }}_{i}}{{\boldsymbol{\varnothing }}_{i}^{\mathrm{T}}\cdot \boldsymbol{M}{\cdot \boldsymbol{\varnothing }}_{i}} 。 $$ (6)

      定义放大倍数P为局部频率的期望值$ {\tilde{\omega }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}} $与初始结构的局部频率$ {\omega }_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}^{0} $的比值,即

      $$ P=\frac{{\tilde{\omega }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}}{{\omega }_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}^{0}} 。 $$ (7)

      可以将当前结构的局部频率表述为

      $$ {\left({\omega }_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}\right)}^{2}=\frac{{\left({\boldsymbol{\varnothing }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}\right)}^{\mathrm{T}}\cdot \boldsymbol{K}\left(\boldsymbol{X}\right)\cdot {\boldsymbol{\varnothing }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}}{{\left({\boldsymbol{\varnothing }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}\right)}^{\mathrm{T}}\cdot \boldsymbol{M}\cdot {\boldsymbol{\varnothing }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}} 。 $$ (8)

      同样基于假设2),当前结构的局部频率也可表述为

      $$ {\left({\omega }_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}\right)}^{2}=\frac{{{(\boldsymbol{\varnothing }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}^{0})}^{\mathrm{T}}\cdot \boldsymbol{K}\left(\boldsymbol{X}\right)\cdot {\boldsymbol{\varnothing }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}^{0}}{{{(\boldsymbol{\varnothing }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}^{0})}^{\mathrm{T}}\cdot \boldsymbol{M}\cdot {\boldsymbol{\varnothing }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}^{0}} 。 $$ (9)

      将式(4)代入式(8),得到

      $$ {\left({\omega }_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}\right)}^{2}=\frac{{{(\mathrm{\varnothing }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}^{0})}^{\mathrm{T}}\cdot {\tilde{F}}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}\cdot {\boldsymbol{\varnothing }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}^{0}}{{{(\mathit{\varnothing }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}^{0})}^{\mathrm{T}}\cdot \boldsymbol{M}\cdot {\boldsymbol{\varnothing }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}^{0}}\cdot \frac{1}{{\boldsymbol{u}}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}} 。 $$ (10)

      将式(3)代入式(10),得到

      $$\begin{split} {\left({\omega }_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}\right)}^{2}= & \frac{{{(\boldsymbol{\varnothing }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}^{0})}^{\mathrm{T}}\cdot \boldsymbol{K}\left({\boldsymbol{X}}_{0}\right)\cdot {\boldsymbol{\varnothing }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}^{0}}{{{(\boldsymbol{\varnothing }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}^{0})}^{\mathrm{T}}\cdot \boldsymbol{M}\cdot {\boldsymbol{\varnothing }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}^{0}}\cdot \frac{{\boldsymbol{u}}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}^{0}}{{\boldsymbol{u}}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}}=\\ & {\left({\omega }_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}^{0}\right)}^{2}\cdot \frac{{\boldsymbol{u}}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}^{0}}{{\boldsymbol{u}}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}} 。\end{split} $$ (11)

      将式(7)代入式(11),可以进一步将期望的局部频率$ {\tilde{\omega }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}} $与当前结构的局部频率$ {\omega }_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}} $联系起来,得到

      $$ {\left({\tilde{\omega }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}\right)}^{2}={P}^{2}\cdot {\left({\omega }_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}\right)}^{2}\cdot \frac{{\boldsymbol{u}}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}}{{\boldsymbol{u}}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}^{0}} 。 $$ (12)

      至此,原始问题中的局部频率约束可表述为

      $$\begin{split} g\left(\boldsymbol{X}\right)= & {\tilde{\omega }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}-{\omega }_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}=\frac{{\left({\tilde{\omega }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}\right)}^{2}-{\left({\omega }_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}\right)}^{2}}{{\tilde{\omega }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}+{\omega }_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}}=\\ & \frac{{\left({\omega }_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}\right)}^{2}}{{\tilde{\omega }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}+{\omega }_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}}\cdot ({P}^{2}\cdot \frac{{\boldsymbol{u}}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}}{{\boldsymbol{u}}_{\mathbf{l}\mathbf{o}\mathbf{c}\mathbf{a}\mathbf{l}}^{0}}-1){\text{≤}} 0 ;\end{split} $$ (13)

      由于结构的频率肯定为非负,所以可以用节点位移约束代替频率约束,则式(13)可以表述为

      $$ {g}^{*}\left(\boldsymbol{X}\right)={P}^{2}\cdot \frac{{\boldsymbol{u}}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}}{{\boldsymbol{u}}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}^{0}}-1{\text{≤}} 0 。 $$ (14)

      从最终的节点位移约束表达式g*(X)可以看出,最优化工程中需要约束结构中每个节点的位移,这给近似优化模型的建立带来了极大的工作量,同时也增加了优化问题求解的复杂性。然而,结构的节点位移与模态振型相同。也就是说,对于局部振动模态,相对明显的节点位移几乎都分布于局部振动区域。因此,在建立近似优化模型时仅约束局部振动区域的节点位移是一种合理的简化方式。

    • 当前的近似优化模型可以表示为非线性规划数学问题,

      $$ \left\{\begin{array}{l}{\mathrm{find}} :{\boldsymbol{X}}={\left\{{x}_{1},{x}_{2},\cdot\cdot\cdot ,{x}_{n}\right\}}^{\mathrm{T}}\\ {\mathrm{min}} :f\left(\boldsymbol{X}\right)\\ s.t.\;\;{g}^{*}\left(\boldsymbol{X}\right)={\boldsymbol{u}}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}-{\dfrac{1}{{P}^{2}}}\cdot {\boldsymbol{\varnothing }}_{\mathrm{l}\mathrm{o}\mathrm{c}\mathrm{a}\mathrm{l}}^{0}{\text{≤}} 0\\ {x}_{i}^{\mathrm{l}}{\text{≤}} {x}_{i}{\text{≤}} {x}_{i}^{\mathrm{u}}\left(i=\mathrm{1,2},\cdot\cdot\cdot ,n\right)\end{array}\right. 。 $$ (15)

      此近似优化模型不再受局部模态的影响,局部模态置换现象不会改变优化模型中的约束条件。

      定义相对误差$ e=\left|{P}_{\mathrm{T}}-{P}_{\mathrm{P}}\right|/{P}_{\mathrm{T}} $,其中,PTPP分别为局部频率的理论放大倍数和实际放大倍数。可用相对误差的大小来度量近似方法的精度。

    • 应用近似方法的优化过程(见图1)概述如下:

      图  1  近似优化模型优化过程

      Figure 1.  Optimization process of the approximation optimization model

      1)在Patran中建立结构的初始有限元模型,并进行详细的模态分析。

      2)从模态分析结果中提取局部振动模态的节点平衡力。

      3)根据步骤2)获得的节点平衡力建立静载荷边界条件。

      4)使用上述信息创建近似优化模型,并进行优化分析。

      5)更新分析模型并检查优化结果是否收敛到最佳值,如是,优化结束;如否,则返回步骤2)。

      6)重复步骤2)~步骤5),直到满足收敛标准,结束优化过程。

      值得注意的是,由于局部模态的阶数等必要的数据需要手动获取,故需手动建模以获得具有更高近似程度的优化模型。如何获得局部模态的节点平衡力并建立有限元模型中的静态载荷条件是创建近似优化模型时的2个难点,本文进行算例验证时将使用Msc.Patran/Nastran软件完成这2项操作。

    • 本章将通过2个数值算例验证本文所提出的近似优化方法的合理性和有效性。第1个简单数值算例的参数来自于Pedersen[12]论文中的一个例子;第2个算例是某大型卫星的实际结构优化问题。

    • 该算例是简单的结构板优化问题,以局部频率不大于60 Hz为优化约束,以结构的质量最小为优化目标,初始有限元模型(FEM)如图2所示,共包含20×20个有限元网格(网格尺寸50 mm×50 mm),图中白色和黑色区域的材料分别为M1M2,其中:M1的弹性模量E1=2.0×1011 N/m2,泊松比v1=0.3,质量密度ρ1=7800 kg/m3M2E2=2.0×105 N/m2v2=0.3,ρ2=78 kg/m3。结构板的尺寸为L1×L2×T=1 m×1 m×0.002 m,T为板厚。局部模态存在高振动密度区域的设计特征,因此结构板的黑色区域出现局部振动,如图3所示。

      图  2  结构板的初始有限元模型

      Figure 2.  Initial FEM of the structural plate

      图  3  初始结构板的局部振型

      Figure 3.  Local modal shape of the initial structural plate

      局部振动模态的节点平衡力可直接从结构模态分析结果中提取,并据此建立静载荷条件。根据新建立的工况进行静力学分析,可以看到,图3所示初始结构板的模态振型与图4所示静力分析的节点位移分布保持一致,与第2章提到的近似优化思路完全相符。

      图  4  静力学分析得到的初始结构板节点位移

      Figure 4.  Nodal displacements of the initial structural plate obtained from static analysis

      在最大位移发生处设置1个加强筋以增大局部频率。在近似模型中,设计变量为加强筋的横截面尺寸(参见图5),通过约束局部区域的节点位移来确定最小结构质量。初始设计对所有4个设计变量施加下限和上限。原优化问题中以局部频率不大于60 Hz作为优化约束,转化至近似优化问题中即以加强筋上的最大节点位移不大于0.557 m作为优化约束。

      图  5  加强筋加固方式及加强筋横截面尺寸

      Figure 5.  Reinforcement method and cross section dimension of the reinforcing stiffener

      近似模型优化结果如表1所示,迭代历程如图6所示。可以看到,经过4次优化迭代后得到结构板的最小质量为125.3 kg,仅比初始质量增加了0.2 kg,符合2.1节提出的假设1)。

      表 1  结构板近似模型优化结果

      Table 1.  Optimization results of approximation model of the structural plate

      参数初始设计近似优化结果
      设计变量H/m0.010.041 4
      W/m0.010.011 7
      t1/m0.0020.001
      t2/m0.000 20.001
      局部振动模态频率/Hz22.50459.136
      局部振动模态的阶数33
      节点最大位移/m4.140.47
      局部频率相对误差/%1.44
      结构总质量/kg125.1125.3

      图  6  结构板近似模型迭优化代历程

      Figure 6.  Iteration history for approximation model of the structural plate

      为了验证近似方法的准确性,使用当前参数对优化后的结构板进行模态分析并提取局部模态振型。如图7所示,结构板优化前后局部模态振型保持一致,符合2.1节提出的假设2)。经进一步观察,局部频率仅存在1.44%的相对误差,符合设计要求。

      图  7  优化前后结构板的局部模态振型

      Figure 7.  Local modal shapes of the structural plate before and after optimization

    • 本算例中,某大型卫星被用于验证本文近似方法在大型卫星结构优化中的效率。首先建立由壳单元和梁单元组成的原始结构有限元模型并进行模态分析,结果表明卫星结构舱内某水平舱板上存在局部振动模态振型(如图8所示),且局部振动频率小于火箭激励频率(60 Hz)。此时卫星及星内设备的安全存在重大共振隐患,因此,需针对水平舱板进行卫星结构优化设计,优化设计约束为局部频率不大于60 Hz,优化目标为卫星结构质量最小。

      图  8  卫星的初始局部模态振型

      Figure 8.  Initial local modal shape of the satellite

      与算例1的优化过程相同,通过静力学分析得到整星的节点位移结果,如图9所示。可以看到,图8的局部模态振型与图9的节点位移分布保持一致。

      图  9  静力学分析得到的卫星节点位移

      Figure 9.  Nodal displacements of the satellite obtained by static analysis

      通过在出现节点最大位移的区域设置矩形横截面加强筋来增大对应的局部频率。以加强筋的横截面尺寸作为优化设计变量(矩形截面的长为H、宽为W),以节点位移不大于0.029 m作为优化约束,以卫星结构质量最小为优化设计目标建立近似优化模型。

      近似模型的优化结果如表2所示,迭代历程如图10所示。可以看到,经过5次优化迭代后得到卫星的总质量为8 194.4 kg,比卫星初始质量仅增加6.5 kg,符合2.1节提出的假设1)。此外,在最后2次迭代中,优化目标的值变化不大。

      表 2  卫星结构近似模型优化结果

      Table 2.  Optimization results of the satellite structural approximation model

      参数初始设计近似优化结果
      设计变量H/m0.010.051 8
      W/m0.010.002 6
      局部振动模态频率/Hz42.23459.667
      局部振动模态的阶数4669
      节点最大位移/m0.0590.022
      局部频率相对误差/%0.54
      结构总质量/kg8 188.38 194.8

      对优化后的卫星进行模态分析并提取局部模态振型。对比优化前后的局部模态振型(图11)可见,优化前后卫星的局部模态振型保持相似,符合2.1节提出的假设2)。经进一步观察,局部频率仅存在0.54%的相对误差,符合设计要求。

      图  10  卫星结构近似模型优化迭代历程

      Figure 10.  Iteration history for approximation model of the satellite structure

      图  11  优化前后卫星的局部模态振型

      Figure 11.  Local modal shapes of the satellite before and after optimization

    • 本文提出一种处理局部频率约束的结构优化问题的近似方法。其核心思想是将具有局部频率约束的动力学优化问题转化为具有节点位移约束的静力学优化问题,则可以不受局部频率非恒定阶数的影响。在节点位移约束的推导过程中,分别对优化前后的质量矩阵和局部振型进行了近似。本文中的数值算例和实际工程优化算例均验证了该方法合理可行,可以有效处理具有局部频率约束的结构优化问题。

参考文献 (14)

目录

    /

    返回文章
    返回