随机激励下大规模结构设计方法技术

本发明专利技术公开了一种随机激励下大规模结构设计方法，用于解决现有的结构设计方法实用性差的技术问题。技术方案是采用虚拟激励法结合模态加速度法计算随机激励下的位移响应均方根，然后以结构指定位置的位移响应均方根最小为目标，以结构质量为约束对结构进行设计。相比背景技术的设计方法，本发明专利技术方法明显提高了随机位移响应均方根的分析精度，最终设计得到清晰有效的结构构型，从而满足工程实际中大规模结构的设计需求。本发明专利技术方法实施例中结构的自由度数量为14762，是背景技术中结构自由度数量3782的3.9倍，设计得到的结构构型清晰有效，易于在实际工程中使用。

【技术实现步骤摘要】
随机激励下大规模结构设计方法
本专利技术涉及一种大规模结构设计方法，特别是涉及一种随机激励下大规模结构设计方法。
技术介绍
工程实际中的结构常常会承受各种随机振动激励，例如自然界中存在的风激励，航空航天飞行器在服役时受到的气动激励。结构在随机激励作用下破坏的事件时有发生，因此在结构构型设计时，考虑结构在随机激励下的性能表现非常重要。文献“ZhangQ,ZhangWH,ZhuJH,GaoT.Layoutoptimizationofmulti-componentstructuresunderstaticloadsandrandomexcitations.EngineeringStructures.2012.43:120-128”公开了一种随机激励下以结构指定位置的位移响应均方根最小化为目标，结构体积为约束的结构构型优化设计方法，该方法采用传统的CompleteQuadraticCombination(CQC)方法计算随机激励下的位移响应均方根。CQC方法计算随机响应时计算量巨大，并且在实际应用中存在由截断模态引起的分析误差，更值得注意的是分析误差会随着结构规模的增加而变大，因此文献中的方法只适用于设计结构规模较小的问题。文献所设计的结构中，最大的自由度数量为3782，规模很小，可以得到清晰的结构构型。但是当使用文献中的方法设计大规模结构时，设计得到的结构构型边界不清晰，构型复杂、结构上空洞较多，难以应用于工程实际(见图3)。
技术实现思路
为了克服现有的结构设计方法实用性差的不足，本专利技术提供一种随机激励下大规模结构设计方法。该方法采用虚拟激励法结合模态加速度法计算随机激励下的位移响应均方根，然后以结构指定位置的位移响应均方根最小为目标，以结构质量为约束对结构进行设计。相比
技术介绍
的设计方法，本专利技术方法可以明显提高随机位移响应均方根的分析精度，最终能够设计得到清晰有效的结构构型，从而能满足工程实际中大规模结构的设计需求。本专利技术解决其技术问题所采用的技术方案是：一种随机激励下大规模结构设计方法，其特点是采用以下步骤：(a)建立设计空间有限元模型，设置拓扑设计变量ηh初始值，h是正整数表示单元编号，1≤h≤Nh，Nh表示结构单元总数量。给定材料密度ρ和杨氏模量E。给定质量约束上限(b)设置激励载荷，给出随机激励f(t)的功率谱密度矩阵Sf(ω)，f(t)为p维列向量，p为载荷中力的个数，t表示时间，Sf(ω)为p维方阵，其下标f表示其为激励f(t)的功率谱矩阵。ω为激励角频率，载荷的激励频段为[ω，]。ω表示激励角频率的下限，表示激励角频率的上限。根据矩阵LDLT分解，存在下式成立其中Q为矩阵Sf(ω)的秩，γq为p维列向量表示第q个虚拟简谐激励，1≤q≤Q，上标T表示向量或矩阵的转置。(c)根据当前设计变量值，采用以下公式分别计算每一个有限元单元的材料密度ρh和杨氏模量Ehρh＝ηhρ更新结构有限元模型中的相应材料属性并进行结构有限元分析。(d)从有限元分析结果中提取每个单元的刚度矩阵kh和质量矩阵mh，结构的前l阶模态频率值ωi，1≤i≤l，模态振型为n行l列矩阵，n为结构总自由度数目。设置结构前l阶模态的阻尼比ζi。采用虚拟激励法结合模态加速度法计算结构自由度r的随机位移响应均方根的公式为式中u表示位移，||(gq(t))r||表示复数(gq(t))r的模，gq(t)为n维列向量表示结构在第q个虚拟简谐激励γq下的位移响应，其第r项的计算公式为式中a为n维列向量，只有第r项为1，其它项均为0。为的第i列。b为n行p列由0、1组成的载荷分布矩阵，假如f(t)中第d个力施加在第z个自由度上，则b的第d列中只有第z个元素值是1，d列中其它元素值均为0。ejωt表示以自然常数e为底数的指数函数，j2＝-1。上式中，xq＝k-1(bγq)式中k为结构有限元整体刚度矩阵，xq是第q个静力载荷bγq下的位移向量。结构自由度r的随机位移响应均方根对每一设计变量的灵敏度的计算公式为式中为偏微分符号，上式中式中式中Re表示复数的实部，Im表示复数的虚部，而式中式中(xq)h表示结构在第q个静力载荷bγq下单元h的位移向量，(ΛT)h表示结构在静力载荷向量a下单元h的位移向量的转置。g表示模态阶数，表示单元h的第i阶模态振型向量。上式中Λ＝k-1a读取每个单元的体积Vh，计算结构整体质量M及其对每一设计变量的灵敏度计算式分别为(e)根据当前设计变量值和灵敏度值，以结构自由度r的随机位移响应均方根为目标函数，结构整体质量M为设计约束，采用数学规划法对优化问题进行求解得到新的设计变量值；(f)重复步骤(c)至步骤(e)，直至最近两次迭代计算得到目标函数相对误差小于1％或达到预设的最大迭代次数，得到设计结果。本专利技术的有益效果是：该方法采用虚拟激励法结合模态加速度法计算随机激励下的位移响应均方根，然后以结构指定位置的位移响应均方根最小为目标，以结构质量为约束对结构进行设计。相比
技术介绍
的设计方法，本专利技术方法明显提高了随机位移响应均方根的分析精度，最终设计得到清晰有效的结构构型，从而满足工程实际中大规模结构的设计需求。本专利技术方法经过实施例160步迭代后得到设计结果。初始结构指定自由度r＝304的随机位移响应均方根为0.06m，设计得到结构的指定自由度r＝304的随机位移响应均方根为0.0015m，随机位移响应均方根降幅达到99.75％。由图2可见设计得到的结构构型清晰有效，易于在实际工程中使用。本专利技术方法实施例中结构的自由度数量为14762，是
技术介绍
中结构自由度数量3782的3.9倍。下面结合附图和具体实施方式对本专利技术作详细说明。附图说明图1是本专利技术方法实施例中几何模型与边界条件示意图。图2是本专利技术方法实施例的设计结果图。图3是
技术介绍
方法的设计结果图。具体实施方式参照图1-2。本专利技术随机激励下大规模结构设计方法具体步骤如下：(a)建立设计空间有限元模型：将长宽厚分别为0.8m，0.4m和0.001m的矩形平面结构划分为120×60的正方形网格，结构左边界固定。设置拓扑设计变量ηh初始值均为0.5，结构单元总数量Nh＝7200。给定材料密度ρ＝7800kg/m3和杨氏模量E＝200GPa。给定质量约束上限(b)设置激励载荷为作用在结构右边界中点竖直向下的随机载荷f(t)，结构只有p＝1个力载荷，因此随机载荷的功率谱密度矩阵Sf(ω)为1维矩阵，其值设为2500N2/(rad/s)，载荷的激励频段为[0，1884]rad/s。根据矩阵LDLT分解，存在下式成立。此处Q＝1为矩阵Sf(ω)的秩，γ1＝50表示第1个虚拟简谐激励。(c)根据当前设计变量值，采用以下公式分别计算每一个有限元单元的材料密度ρh和杨氏模量Ehρh＝ηhρ更新结构有限元模型中的相应材料属性并进行结构有限元分析。(d)从有限元分析结果中提取每个单元的刚度矩阵kh和质量矩阵mh，前l＝30阶模态频率值ωi，模态振型为n＝14762行l＝30列矩阵。设置结构前l＝30阶模态的阻尼比ζi均为0.03。采用虚拟激励法结合模态加速度法计算结构自由度r＝304的随机位移响应均方根的公式为式中u表示位移，||(gq(t))304||表示复数(gq(t))304的模，gq(t)本文档来自技高网...

【技术保护点】
一种随机激励下大规模结构设计方法，其特征在于包括以下步骤：(a)建立设计空间有限元模型，设置拓扑设计变量ηh初始值，h是正整数表示单元编号，1≤h≤Nh，Nh表示结构单元总数量；给定材料密度ρ和杨氏模量E；给定质量约束上限(b)设置激励载荷，给出随机激励f(t)的功率谱密度矩阵Sf(ω)，f(t)为p维列向量，p为载荷中力的个数，t表示时间，Sf(ω)为p维方阵，其下标f表示其为激励f(t)的功率谱矩阵；ω为激励角频率，载荷的激励频段为表示激励角频率的下限，表示激励角频率的上限；根据矩阵LDLT分解，存在下式成立Sf(ω)=Σq=1Q(γq)*(γq)T]]>其中Q为矩阵Sf(ω)的秩，γq为p维列向量表示第q个虚拟简谐激励，1≤q≤Q，上标T表示向量或矩阵的转置；(c)根据当前设计变量值，采用以下公式分别计算每一个有限元单元的材料密度ρh和杨氏模量Ehρh＝ηhρEh=15ηh5+ηh16E]]>更新结构有限元模型中的相应材料属性并进行结构有限元分析；(d)从有限元分析结果中提取每个单元的刚度矩阵kh和质量矩阵mh，结构的前l阶模态频率值ωi，1≤i≤l，模态振型为n行l列矩阵，n为结构总自由度数目；设置结构前l阶模态的阻尼比ζi；采用虚拟激励法结合模态加速度法计算结构自由度r的随机位移响应均方根的公式为σur=(∫ω‾ω‾Σq=1Q||(gq(t))r||2dω)1/2]]>式中u表示位移，||(gq(t))r||表示复数(gq(t))r的模，gq(t)为n维列向量表示结构在第q个虚拟简谐激励γq下的位移响应，其第r项的计算公式为式中a为n维列向量，只有第r项为1，其它项均为0；为的第i列；b为n行p列由0、1组成的载荷分布矩阵，假如f(t)中第d个力施加在第z个自由度上，则b的第d列中只有第z个元素值是1，d列中其它元素值均为0；ejωt表示以自然常数e为底数的指数函数，j2＝‑1；上式中，Hi=(ωi2-ω2+2jζiωiω)-1]]>xq＝k‑1(bγq)式中k为结构有限元整体刚度矩阵，xq是第q个静力载荷bγq下的位移向量；结构自由度r的随机位移响应均方根对每一设计变量的灵敏度的计算公式为∂σur∂ηh=12σur∫ω‾ω‾Σq=1Q||(gq(t))r||2∂ηhdω]]>式中为偏微分符号，上式中∂||(gq(t))r||2∂ηh=2||(gq(t))r||∂||(gq(t))r||∂ηh]]>式中∂||(gq(t))r||∂ηh=(Re((gq(t))r)Re(∂(gq(t))r∂ηh)+Im((gq(t))r)Im(∂(gq(t))r∂ηh))||(gq(t))r||-1]]>式中Re表示复数的实部，Im表示复数的虚部，而式中aT∂xq∂ηh=-(ΛT)h∂kh∂ηh(xq)h]]>∂Hi∂ηh=-Hi2(2ωi∂ωiηh+2jζiω∂ωiηh)]]>式中(xq)h表示结构在第q个静力载荷bγq下单元h的位移向量，(ΛT)h表示结构在静力载荷向量a下单元h的位移向量的转置；g表示模态阶数，1≤g≤l，表示单元h的第i阶模态振型向量；上式中Λ＝k‑1a∂kh∂ηh=75ηh4+115ηh5+ηhkh]]>∂mh∂ηh=1ηhmh]]>读取每个单元的体积Vh，计算结构整体质量M及其对每一设计变量的灵敏度计算式分别为M=Σh=1NhρhVh]]>∂M∂ηh=ρVh]]>(e)根据当前设计变量值和灵敏度值，以结构自由度r的随机位移响应均方根为目标...

【技术特征摘要】
1.一种随机激励下大规模结构设计方法，其特征在于包括以下步骤：(a)建立结构有限元模型，设置设计变量ηh初始值，h是正整数表示单元编号，1≤h≤Nh，Nh表示结构单元总数量；给定材料密度ρ和杨氏模量E；给定质量约束上限(b)设置激励载荷，给出随机激励f(t)的功率谱密度矩阵Sf(ω)，f(t)为p维列向量，p为载荷中力的个数，t表示时间，Sf(ω)为p维方阵，其下标f表示其为激励f(t)的功率谱矩阵；ω为激励角频率，载荷的激励频段为ω表示激励角频率的下限，表示激励角频率的上限；根据矩阵LDLT分解，存在下式成立其中Q为矩阵Sf(ω)的秩，γq为p维列向量表示第q个虚拟简谐激励，1≤q≤Q，上标T表示向量或矩阵的转置；(c)根据当前设计变量值，采用以下公式分别计算每一个有限元单元的材料密度ρh和杨氏模量Ehρh＝ηhρ更新结构有限元模型中的相应材料属性并进行结构有限元分析；(d)从有限元分析结果中提取每个单元的刚度矩阵kh和质量矩阵mh，结构的前l阶模态频率值ωi，1≤i≤l，模态振型为n行l列矩阵，n为结构总自由度数目；设置结构前l阶模态的阻尼比ζi；采用虚拟激励法结合模态加速度法计算结构自由度r的随机位移响应均方根的公式为式中u表示位移，||(gq(t))r||表示复数(gq(t))r的模，gq(t)为n维列向量表示结构在第q个虚拟简谐激励γq下的位移响应，其第r项的计算公式为式中a为n维列向量，只有第r项为1，其它项均为0；为的第i列；b为n行p列由0、1组成的载荷分布矩阵，假如f(t)中第d个力施加在第z个自由度上，则b的第d列中只有第z个元素值是1，d列中其它元素值均为0；ejωt表示以自然常数e为底数的指数函数，j2＝-1；上式中，xq＝k-1(bγq)式中k为结构有限元整体刚度矩阵，xq是第q个静力载荷bγq下的位移向量；结构自由度r的随机位移响应均方根对每一设计变量的...

【专利技术属性】
技术研发人员：张卫红，刘虎，牛草，蔡守宇，高彤，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：陕西;61