一种大载荷、高应力、多路径传载结构的优化方法技术

本发明专利技术涉及结构强度分析技术领域，特别涉及一种大载荷、高应力、多路径传载结构的优化方法。该方法包括：根据多路径传载结构的已知设计空间，建立其有限元模型；对该有限元模型进行进行变密度法的拓扑优化，得到拓扑优化最佳传力路线及最佳传力模型；其中，以结构能量最小化为目标函数，以第一结构体积比、集中载荷施加点的变形作为约束条件，以结构对称平面作为几何约束限制，冻结边界施加区域；基于最佳传力模型，以结构重量最小化为目标函数，以第二结构体积比、集中载荷施加点的变形作为约束条件，并考虑加工工艺设置最大/最小原件尺寸限制，以结构对称平面作为几何约束限制，进行多路径传载的尺寸优化，得到多路径传载结构的尺寸优化参数。的尺寸优化参数。的尺寸优化参数。

【技术实现步骤摘要】
一种大载荷、高应力、多路径传载结构的优化方法


[0001]本专利技术涉及结构强度分析
，特别涉及一种大载荷、高应力、多路径传载结构的优化方法。

技术介绍

[0002]大型飞机受力形式日趋复杂化，大载荷、较高结构应力水平亟需规划多路径传载结构形式传递飞机飞行载荷。现有多路径传载结构传力形式分析主要基于工程理论分析，分析精度有限，本文采用多约束优化方法，精确对大载荷、高应力传载结构进行传载表征。

技术实现思路

[0003]专利技术目的：针对大载荷、高应力水平的飞机传力结构，规划多路径传载结构形式，形成较优的结构传力方案和结构参数，最终形成一种大载荷、高应力、多路径传载结构的优化方法。
[0004]技术方案：
[0005]一种大载荷、高应力、多路径传载结构的优化方法，其特征在于，包括：
[0006]根据多路径传载结构的已知设计空间，建立其有限元模型；
[0007]对该有限元模型进行进行变密度法的拓扑优化，得到拓扑优化最佳传力路线及最佳传力模型；其中，以结构能量最小化为目标函数，以第一结构体积比、集中载荷施加点的变形作为约束条件，以结构对称平面作为几何约束限制，冻结边界施加区域；
[0008]基于最佳传力模型，以结构重量最小化为目标函数，以第二结构体积比、集中载荷施加点的变形作为约束条件，并考虑加工工艺设置最大/最小原件尺寸限制，以结构对称平面作为几何约束限制，进行多路径传载的尺寸优化，得到多路径传载结构的尺寸优化参数。
[0009]得到多路径传载结构的尺寸优化参数之后，所述方法还包括：
[0010]根据多路径传构的载结尺寸优化参数，设计多路径传载结构的具体结构；
[0011]根据具体结构，建立静态求解有限元模型，并静力求解，得到结构应力应变分布图；
[0012]基于结构应力应变分布图，对具体结构进行缘条稳定性、支柱刚度验证；
[0013]若某个结构不满足任一验证条件，则对该结构根据该验证条件进行局部修正；最终修正的参数和无需修正的参数作为多路径传载结构的最佳结构参数。
[0014]对该有限元模型进行进行变密度法的拓扑优化之后，所述方法还包括：
[0015]根据多路径传载结构的已知设计空间，对集中载荷施加点的单传力路径和双传力路径进行工程分析，选出最优传力路径；
[0016]对比拓扑优化最佳传力路线和最优传力路径；
[0017]若两者符合，则认为拓扑优化最佳传力路线合理；若两者不符合，则修改拓扑优化的约束参数再次进行优化，直到两者符合。
[0018]变密度法将结构拓扑优化问题转换为以单元材料密度为优化变量的非线性优化
问题；非线性问题能转换为二次规划问题，即序列规划方法，序列规划方法能通过灵敏度分析，寻找最优方向来更新单元密度分布，最终获得优化收敛结果。
[0019]静态求解有限元模型的网格尺寸是3mm，单元类型采用2次修正的四面体单元。
[0020]有限元模型的网格尺寸是6mm，单元类型采用2次修正的四面体单元。
[0021]缘条稳定性约束：
[0022][0023]式中，E、μ
e
为材料的弹性模量、泊松比；b、δ为加载边的宽度、加载边的宽度，K
c
为压缩临界应力系数，η为塑性修正系数。
[0024]支柱刚度约束：
[0025]a)支柱截面积大于腹板面积40％；
[0026]A
min
＝0.4
×
A；
[0027]A、A
min
是指隔波构件截面积、壁板构件截面积；
[0028]b)支柱构件满足最小惯性矩要求；
[0029][0030]I
u
、I
min
指隔波构件惯性矩、最小惯性矩；d、t、h
e
、k
s
是指支柱间距、腹板厚度、上下缘条形心间的距离、剪切屈曲系数。
[0031]有益效果：创新性地对大载荷、高应力、多路径传载进行传力路线拓扑优化；对较优传力方案进行尺寸优化，得到较优的结构尺寸参数；优化结构进行精确结构参数表征，有效减轻结构重量。
附图说明
[0032]图1为多路径传载结构的结构示意图；
[0033]图2为两种传力路径的示意图；
[0034]图3为拓扑优化有限元模型的示意图；
[0035]图4为拓扑优化流程图；
[0036]图5为拓扑优化材料分布图；
[0037]图6为尺寸优化结构尺寸分布图；
[0038]图7为设计出的多路径传载结构的结构图；
[0039]图8为设计出的多路径传载结构的优化结构静力求解结果的示意图；
[0040]图9为多路径传载结构的最佳结构参数图。
具体实施方式
[0041]某大型飞机多路径传载结构[4]，此结构通过集中点[1]传递集中载荷P，此结构两个约束位置为集中点[2]和集中点[3]，见图1。
[0042](1)基于图1结构进行工程受力分析，集中载荷施加点P有两种传力路径：一是集中
剪力Q和对集中点[5]的弯矩M均通过集中点[5]进行传递，二是弯矩M通过集中点[4]剪力Q1、集中点[5]剪力Q2形成的一对力偶进行平衡，剪力通过两点的剪力差进行平衡，见图2；
[0043]Q＝P,M＝PL1[0044][0045]式中L1、L2是指集中载荷施加点到集中点[5]、集中点[4]的水平距离。
[0046](2)对两种不同的传力路径进行传力路线的拓扑优化分析，形成较优传力路线：以结构能量最小化为目标函数，以结构体积比(设置优化前后体积比为0.2)、集中点[1]的变形(点[1]的变形值为25m)作为约束条件，以结构XZ对称平面作为几何约束限制，区域1、区域2、区域3作为冻结边界区域，运用变密度法进行拓扑优化分析，模型离散采用2次修正的四面体单元，见图3，此数值拓扑优化模型的数学描述如下：
[0047][0048]式中：C(x)指能量函数，V(x)、V0、f指优化前体积、优化后体积、优化前后体积比，X
c
、X
max
是指优化变形、最大变形数值。
[0049](3)应用变密度法将结构拓扑优化问题转换为以单元材料密度为优化变量的非线性优化问题。将非线性问题转换为一系列线性或二次规划问题，即序列规划方法。序列规划方法通过灵敏度分析，寻找最优方向来更新单元密度分布，最终获得优化收敛结果，具体拓扑优化流程见图4；
[0050](4)基于步骤(3)，得到满足约束条件的拓扑材料优化分布结果，见图5；
[0051](5)图5结构，进行尺寸优化分析，以结构重量最小化为目标函数，以体积比(体积比0.1)和集中点[1]变形(最大变形25mm)作为约束条件，并设置最大/最小原件尺寸限制(单元最小尺寸1.5mm、最大尺寸3mm)及XZ平面对称结构限制，进行尺寸优化，数学表述如下，优化本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种大载荷、高应力、多路径传载结构的优化方法，其特征在于，包括：根据多路径传载结构的已知设计空间，建立其有限元模型；对该有限元模型进行进行变密度法的拓扑优化，得到拓扑优化最佳传力路线及最佳传力模型；其中，以结构能量最小化为目标函数，以第一结构体积比、集中载荷施加点的变形作为约束条件，以结构对称平面作为几何约束限制，冻结边界施加区域；基于最佳传力模型，以结构重量最小化为目标函数，以第二结构体积比、集中载荷施加点的变形作为约束条件，并考虑加工工艺设置最大/最小原件尺寸限制，以结构对称平面作为几何约束限制，进行多路径传载的尺寸优化，得到多路径传载结构的尺寸优化参数。2.根据权利要求1所述的方法，其特征在于，得到多路径传载结构的尺寸优化参数之后，所述方法还包括：根据多路径传构的载结尺寸优化参数，设计多路径传载结构的具体结构；根据具体结构，建立静态求解有限元模型，并静力求解，得到结构应力应变分布图；基于结构应力应变分布图，对具体结构进行缘条稳定性、支柱刚度验证；若某个结构不满足任一验证条件，则对该结构根据该验证条件进行局部修正；最终修正的参数和无需修正的参数作为多路径传载结构的最佳结构参数。3.根据权利要求1所述的方法，其特征在于，对该有限元模型进行进行变密度法的拓扑优化之后，所述方法还包括：根据多路径传载结构的已知设计空间，对集中载荷施加点的单传力路径和双传力路径进行工程分析，选出最优传力路径；对比拓扑优化最佳传力路线和最优传力路径；若两者符合，则认为拓扑优化最佳传力路线合理；若两...

【专利技术属性】
技术研发人员：郑茂亮，赵占文，鬲钰焯，
申请(专利权)人：中国航空工业集团公司西安飞机设计研究所，
类型：发明
国别省市：