基于大质量法的随机加速度激励下结构拓扑优化设计方法
【专利摘要】本发明公开了一种基于大质量法的随机加速度激励下结构拓扑优化设计方法，用于解决现有随机载荷下的结构拓扑优化方法实用性差的技术问题。技术方案是采用大质量法将多点加速度激励转化为力激励施加到结构上，采用虚拟激励法结合模态加速度法计算随机激励下位移响应均方根，然后以结构指定位置位移响应均方根最小为目标，以结构质量为约束进行设计。相比【背景技术】的设计方法，本发明方法可以实现多点加速度随机激励，不局限于对结构进行单点加载。最终能够设计得到清晰有效的结构构型，从而能满足工程实际中考虑多点加速度加载的设计需求，实用性强。
【专利说明】
基于大质量法的随机加速度激励下结构拓扑优化设计方法
技术领域
[0001] 本发明涉及一种随机载荷下的结构拓扑优化方法，特别涉及一种基于大质量法的 随机加速度激励下结构拓扑优化设计方法。
【背景技术】
[0002] 实际工程中的结构经常受到各种峰值和谷值大小及序列随机出现的变频变幅载 荷，即随机载荷。如自然界中存在的风激励，航空航天飞行器在服役时受到的气动激励，汽 车在行驶过程中由于路面不平整受到的路面激励，地震激励等。随机激励对结构的正常工 作产生很大影响，甚至使结构破坏。因此在结构构型设计时考虑结构在随机激励下的性能 表现非常重要。
[0003] 文南犬l"Zhang ff .H. ,Liu H,Gao T. Topology optimization of large-scale structures subjected to stationary random exc i tat ion : An efficient optimization procedure integrating pseudo excitation method and mode acceleration method[J] ·Computers&Structures,2015,158:61-70 ·" 公开了一种随机载 荷下的结构拓扑优化方法。该方法应用最为广泛的完全二次结合法(CQC)进行随机响应分 析，当结构的自由度规模非常大时，其计算量会非常庞大，这给完全二次结合法在实际工程 中的应用造成极大的困扰。虚拟激励法(PEM)极大地提高了随机振动问题的求解效率，很好 地解决了完全二次结合法在计算效率上的问题。文献中对比了传统的虚拟激励法和将虚拟 激励法与模态加速度法进行结合得到的改进的虚拟激励法。虚拟激励法将随机响应功率谱 密度的求解转换为虚拟简谐响应的求解，传统的虚拟激励法用模态位移法(MDM)进行虚拟 简谐响应分析，在处理实际中的大规模自由度结构时，面临模态位移法的模态截断误差明 显增加，导致优化失败的问题。改进的虚拟激励法用计算效率非常接近，但计算精度远远高 于模态位移法的模态加速度(MAM)来进行虚拟简谐响应分析，因此可以大大提高随机响应 的精度。
[0004] 文献1公开的方法虽然能够实现随机激励结构拓扑优化，但是由于其激励是施加 在单个点的力类型随机激励，不能施加类似于地震激励的多点加速度随机激励。
【发明内容】
[0005] 为了克服现有随机载荷下的结构拓扑优化方法实用性差的不足，本发明提供一种 基于大质量法的随机加速度激励下结构拓扑优化设计方法。该方法采用大质量法将多点加 速度激励转化为力激励施加到结构上，采用虚拟激励法结合模态加速度法计算随机激励下 位移响应均方根，然后以结构指定位置位移响应均方根最小为目标，以结构质量为约束进 行设计。相比【背景技术】的设计方法，本发明方法可以实现多点加速度随机激励，不局限于对 结构进行单点加载。最终能够设计得到清晰有效的结构构型，从而能满足工程实际中考虑 多点加速度加载的设计需求，实用性强。
[0006] 本发明解决其技术问题所采用的技术方案:一种基于大质量法的随机加速度激励 下结构拓扑优化设计方法，其特点是包括以下步骤：
[0007] 步骤一、建立有限元模型，在拟施加激励位置外建立一大质量点，大质量点取结构 重量的1〇6倍，作为载荷的实际施加处，在大质量点与拟施加激励处的所有节点之间建立刚 性连接，设置拓扑设计变量nh初始值，h是正整数表示单元编号，Kh<Nh，Nh表示结构单元 总数量，给定材料密度P和杨氏模量E，给定质量约束上限反'
[0008] 步骤二、设置激励载荷，给出随机激励f(t)的功率谱密度矩阵Sf(co)，f(t)为p维 列向量，P为载荷中力的个数，t表示时间，Sf(c〇)为p维方阵，其下标f表示其为激励f(t)的 功率谱矩阵。ω为激励角频率，载荷的激励频段为[逛，表示激励角频率的下限，表 示激励角频率的上限。根据矩阵LDLT分解，存在下式。
[0010]其中Q为矩阵Sf(co)的秩，γΑρ维列向量表示第q个虚拟简谐激励，1彡q彡Q，上标 T表示向量或矩阵的转置。
[0011]根据大质量法原理存在下式：
[0013]其中，M、C、K分别表示质量矩阵、阻尼矩阵、刚度矩阵，:￡>、p、U分别表示加速度、 速度、位移，下标s表示结构非支撑处的自由度，下标g表示结构支撑处的自由度，Mgg是基底 大质量矩阵，为地震激励，y q为虚拟简谐激励向量。将式(2)的第2行展开，得下式：
[0015] 将式(3)左右两端左乘Mgf1，由于Mgf1中对角元素趋近于零，得到基础激励处的加 速度：
[0016] iig^ud (4)
[0017] 步骤三、根据当前设计变量值，采用以下材料插值模型分别计算每一个有限元单 元的材料密度Ph和杨氏模量Eh
[0018] Ph=P% (5)
[0020] 更新结构有限元模型中的相应材料属性并进行结构有限元分析。
[0021] 步骤四、从有限元分析结果中提取每个单元的刚度矩阵Kh和质量矩阵Mh，结构的前 1阶模态频率值模态振型供，炉为η行1列矩阵，η为结构总自由度数目，设结构 前1阶阻尼比为ξ?，ξ?为Rayleigh阻尼，按下式计算：
[0023] α与β为Rayleigh阻尼系数。
[0024]采用虚拟激励法结合模态加速度法计算结构自由度r的随机位移响应均方根0V 公式为
[0026]式中U表示位移，I I (gq(t))r| I表示复数(gq(t))r的模，gq(t)为η维列向量表示结构 在第q个虚拟简谐激励Yq下的位移响应，其第r项的计算公式为
[0028]式中不考虑结构刚体模态，计算所得位移为自由度r相对基础点的相对位移。式中 a为η维列向量，只有第r项为1，其它项均为0。内:为识的第i列。b为η行p列由0、1组成的载荷 分布矩阵，假如f(t)中第d个力施加在第ζ个自由度上，则b的第d列中只有第ζ个元素值是1， d列中其它元素值均为0。一4表示以自然常数e为底数的指数函数，ω为激励频率，j2 = -l。 式(10)中，
[0031]式中K为结构有限元整体刚度矩阵，Xq是第q个静力载荷by q下的相对位移向量，利 用惯性释放分析计算。
[0032]步骤五、定义拓扑优化模型：
[0034] 式中^为设计变量下限值，取0.001。°^表示结构自由度4勺随机位移响应均方值。 Μ表示结构质量。
[0035] 步骤六、将模型进行一次有限元分析;通过优化灵敏度分析，求得目标函数和约束 条件的灵敏度，选取梯度优化算法进行优化设计，得到优化结果。
[0036] 本发明的有益效果是:采用大质量法将多点加速度激励转化为力激励施加到结构 上，采用虚拟激励法结合模态加速度法计算随机激励下位移响应均方根，然后以结构指定 位置位移响应均方根最小为目标，以结构质量为约束进行设计。相比【背景技术】的设计方法， 本发明方法可以实现多点加速度随机激励，不局限于对结构进行单点加载。最终能够设计 得到清晰有效的结构构型，从而能满足工程实际中考虑多点加速度加载的设计需求。
[0037]本发明方法经过实施例225步迭代后得到设计结果。初始结构指定自由度r = 204 的随机位移响应均方根为1.409mm，设计得到结构的指定自由度r = 204的随机位移响 应均方根0??为〇.〇93mm，随机位移响应均方根降幅达到93.4%。由图1可见设计得到的结构 构型清晰有效，易于在实际工程中使用，实用性强。
[0038]下面结合附图和【具体实施方式】对本发明作详细说明。
【附图说明】
[0039] 图1是本发明方法实施例的设计结果图。
【具体实施方式】
[0040] 参照图1。本发明基于大质量法的随机加速度激励下结构拓扑优化设计方法具体 步骤如下：
[0041 ] (a)建立设计空间有限元模型:将长宽厚分别为0.8m，0.4m和0.005m的矩形平面结 构划分为80X40的正方形网格，在大质量点所在节点上施加质量为3 X107kg的大质量点， 约束大质量点处除竖直方向以外的所有自由度。设置拓扑设计变量nh初始值均为0.5。结构 单元总数目Nh = 3200。给定材料密度P = 7800kg/m3，杨氏模量E = 200GPa，质量约束上限 Λ/ =6.241(2。 〇
[0042] (b)设置激励载荷为作用在施加大质量点处节点竖直向上的随机载荷f(t)，结构 只有P=1个力载荷，因此随机载荷的功率谱密度矩阵Sf(?)为1维矩阵，其值设为2500N2/ (rad/s)，载荷的激励频段为[0,5027]rad/s。根据矩阵LDLT分解，存在下式。
[0044]此处Q=1为矩阵Sf( ω )的秩，γ 1 = 50表示第1个虚拟简谐激励。
[0045] (c)根据大质量法原理存在下式：
[0047]其中，M、C、K分别表示质量矩阵、阻尼矩阵、刚度矩阵，、办、U分别表示加速度、 速度、位移，下标s表示结构非支撑处的自由度，下标g表示结构支撑处的自由度，Mgg是基底 大质量矩阵，为地震激励(地震动加速度），y q为虚拟简谐激励向量。将式(2)的第2行展 开，得下式：
[0049] 将式(3)左右两端左乘Mgf1，由于Mgf1中对角元素趋近于零，可得基础激励处的加 速度：
[0050] Ug ?lid (4)
[0051] 这样就将加速度载荷等效为虚拟简谐力载荷施加到结构上了。
[0052] (d)根据当前设计变量值，采用以下材料插值模型分别计算每一个有限元单元的 材料密度Ph和杨氏模量Eh [0053] Ph=pnh (5)
[0055] 更新结构有限元模型中的相应材料属性并进行结构有限元分析。[0056] (e)从有限元分析结果中提取每个单元的刚度矩阵Kh和质量矩阵Mh，结构的前1 = 30阶模态频率值c〇i，模态振型为n = 6642行1 = 30列矩阵。设结构前1 = 30阶模态的阻 尼比为ξ?，ξ?为Rayleigh阻尼，按下式计算：
[0058] Rayleigh 阻尼系数α = 〇 · 〇1与β = 〇 · 00001。
[0059]采用虚拟激励法结合模态加速度法计算结构自由度r = 204(右端中点竖直方向对 应的自由度)的随机位移响应均方根的公式为
[0061]式中U表不位移，I I (gq(t))2Q4| I表不复数（gq(t))2Q4的模，gq(t)为6642维列向量表 示结构在第q个虚拟简谐激励Yq下的位移响应，其第r = 204项的计算公式为
[0063] 式中不考虑结构刚体模态，计算所得位移为自由度r = 204相对基础点的相对位 移。式中a为6642维列向量，只有第r = 204项元素值为1，其它项均为(Lb为6642行1列由0，1 组成的载荷分布矩阵，f(t)中只有1个力施加在第z = 204个自由度上，因此b的第1列中只有 第204个元素值是1，其它元素值均为Ο。-#表示以自然常数e为底数的指数函数，j 2 = _l。式 (9)中
[0064] Hi=( 〇2i-o2+2j|i〇i〇 )-1 (10)
[0065] Xq = K-Hbyq) (11)
[0066] 式中K为结构有限元整体刚度矩阵，Xq是第q个静力载荷by q下的相对位移向量，利 用惯性释放分析计算。
[0067] (f)定义拓扑优化模型：
[0069] 式中&为设计变量下限值，取0.001。σ?M表示自由度r = 204的随机位移响应均方 根。Μ表示结构质量。
[0070] (g)将模型进行一次有限元分析;通过优化灵敏度分析，求得目标函数和约束条件 的灵敏度，选取梯度优化算法进行优化设计，得到优化结果。
[0071] 本发明方法经过实施例225步迭代后得到设计结果。初始结构指定自由度r = 204 的随机位移响应均方根^^为1.409mm，设计得到结构的指定自由度r = 204的随机位移响应 均方根^7^?为〇.〇93mm，随机位移响应均方根降幅达到93.4%。由图1可见设计得到的结构构 型有效。本实施例中将载荷由背景文献中的单点随机力加载变化为更加符合实际情况的多 点随机加速度载荷，能更好地模拟工程实际中的受载情况。本实施例表明了本发明方法在 处理多点随机加速度激励下结构设计问题上的有效性。
【主权项】
1. 一种基于大质量法的随机加速度激励下结构拓扑优化设计方法，其特征在于包括W 下步骤： 步骤一、建立有限元模型，在拟施加激励位置外建立一大质量点，大质量点取结构重量 的1〇6倍，作为载荷的实际施加处，在大质量点与拟施加激励处的所有节点之间建立刚性连 接，设置拓扑设计变量nh初始值，h是正整数表示单元编号，l《h《化，化表示结构单元总数 量，给定材料密度P和杨氏模量E，给定质量约束上限态； 步骤二、设置激励载荷，给出随机激励f(t)的功率谱密度矩阵Sf(w)，f(t)为P维列向 量，P为载荷中力的个数，t表示时间，Sf(w)为P维方阵，其下标f表示其为激励f(t)的功率 谱矩阵；ω为激励角频率，载荷的激励频段为「ω，? ]，·兰表示激励角频率的下限，?表示激 励角频率的上限;根据矩阵LDLT分解，存在下式；(1) 其中Q为矩阵Sf(co)的秩，丫 q为P维列向量表示第q个虚拟简谐激励，l《q《Q，上标T表 示向量或矩阵的转置； 根据大质量法原理存在下式：其中，M、C、K分别表示质量矩阵、阻尼矩阵、刚度矩阵/&、&、U分别表示加速度、速度、 位移，下标S表示结构非支撑处的自由度，下标g表示结构支撑处的自由度，Mgg是基底大质量 矩阵，为地震激励，丫 q为虚拟简谐激励向量;将式(2)的第2行展开，得下式：将式(3)左右两端左乘Mgg-i，由于Mgg-i中对角元素趋近于零，得到基础激励处的加速度：(4) 步骤Ξ、根据当前设计变量值，采用W下材料插值模型分别计算每一个有限元单元的 材料密度Ph和杨氏模量Eh更新结构有限元模型中的相应材料属性并进行结构有限元分析； 步骤四、从有限元分析结果中提取每个单元的刚度矩阵Kh和质量矩阵Mh，结构的前1阶 模态频率值ω i，1《i《1，模态振型妒为η行1列矩阵，η为结构总自由度数目，设结构前1 阶阻尼比为Ci，Ci为Raylei曲阻尼，按下式计算：(巧 曰与β为Ra^ei曲阻尼系数； 采用虚拟激励法结合模态加速度法计算结构自由度r的随机位移响应均方根巧,f公式为(8) 式中U表示位移，II (gq(t))r||表示复数(gq(t))r的模，gq(t)为η维列向量表示结构在第 q个虚拟简谐激励γ q下的位移响应，其第r项的计算公式为巧) 式中不考虑结构刚体模态，计算所得位移为自由度r相对基础点的相对位移;式中a为η 维列向量，只有第r项为1，其它项均为0;斯'为口的第i列;b为η行Ρ列由0、1组成的载荷分布矩 阵，假如f(t)中第d个力施加在第Ζ个自由度上，贝化的第d列中只有第Ζ个元素值是l，d列中 其它元素值均为〇;e^t表示W自然常数e为底数的指数函数，ω为激励频率，j2 = -l;式（10) 中，式中K为结构有限元整体刚度矩阵，Xq是第q个静力载荷by q下的相对位移向量，利用惯 性释放分析计算； 步骤五、定义拓扑优化模型：式中3为设计变量下限值，取0.001;巧V表示结构自由度r的随机位移响应均方值;Μ表示 结构质量； 步骤六、将模型进行一次有限元分析;通过优化灵敏度分析，求得目标函数和约束条件 的灵敏度，选取梯度优化算法进行优化设计，得到优化结果。
【文档编号】G06F17/50GK106096119SQ201610398316
【公开日】2016年11月9日
【申请日】2016年6月7日
【发明人】朱继宏, 何飞, 张卫红, 郭文杰, 杨开科
【申请人】西北工业大学