CN112699484B - 一种大载荷、高应力、多路径传载结构的优化方法 - Google Patents
一种大载荷、高应力、多路径传载结构的优化方法 Download PDFInfo
- Publication number
- CN112699484B CN112699484B CN202011612688.0A CN202011612688A CN112699484B CN 112699484 B CN112699484 B CN 112699484B CN 202011612688 A CN202011612688 A CN 202011612688A CN 112699484 B CN112699484 B CN 112699484B
- Authority
- CN
- China
- Prior art keywords
- optimization
- load
- multipath
- force transmission
- optimal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及结构强度分析技术领域,特别涉及一种大载荷、高应力、多路径传载结构的优化方法。该方法包括:根据多路径传载结构的已知设计空间,建立其有限元模型;对该有限元模型进行进行变密度法的拓扑优化,得到拓扑优化最佳传力路线及最佳传力模型;其中,以结构能量最小化为目标函数,以第一结构体积比、集中载荷施加点的变形作为约束条件,以结构对称平面作为几何约束限制,冻结边界施加区域;基于最佳传力模型,以结构重量最小化为目标函数,以第二结构体积比、集中载荷施加点的变形作为约束条件,并考虑加工工艺设置最大/最小原件尺寸限制,以结构对称平面作为几何约束限制,进行多路径传载的尺寸优化,得到多路径传载结构的尺寸优化参数。
Description
技术领域
本发明涉及结构强度分析技术领域,特别涉及一种大载荷、高应力、多路径传载结构的优化方法。
背景技术
大型飞机受力形式日趋复杂化,大载荷、较高结构应力水平亟需规划多路径传载结构形式传递飞机飞行载荷。现有多路径传载结构传力形式分析主要基于工程理论分析,分析精度有限,本文采用多约束优化方法,精确对大载荷、高应力传载结构进行传载表征。
发明内容
发明目的:针对大载荷、高应力水平的飞机传力结构,规划多路径传载结构形式,形成较优的结构传力方案和结构参数,最终形成一种大载荷、高应力、多路径传载结构的优化方法。
技术方案:
一种大载荷、高应力、多路径传载结构的优化方法,其特征在于,包括:
根据多路径传载结构的已知设计空间,建立其有限元模型;
对该有限元模型进行进行变密度法的拓扑优化,得到拓扑优化最佳传力路线及最佳传力模型;其中,以结构能量最小化为目标函数,以第一结构体积比、集中载荷施加点的变形作为约束条件,以结构对称平面作为几何约束限制,冻结边界施加区域;
基于最佳传力模型,以结构重量最小化为目标函数,以第二结构体积比、集中载荷施加点的变形作为约束条件,并考虑加工工艺设置最大/最小原件尺寸限制,以结构对称平面作为几何约束限制,进行多路径传载的尺寸优化,得到多路径传载结构的尺寸优化参数。
得到多路径传载结构的尺寸优化参数之后,所述方法还包括:
根据多路径传构的载结尺寸优化参数,设计多路径传载结构的具体结构;
根据具体结构,建立静态求解有限元模型,并静力求解,得到结构应力应变分布图;
基于结构应力应变分布图,对具体结构进行缘条稳定性、支柱刚度验证;
若某个结构不满足任一验证条件,则对该结构根据该验证条件进行局部修正;最终修正的参数和无需修正的参数作为多路径传载结构的最佳结构参数。
对该有限元模型进行进行变密度法的拓扑优化之后,所述方法还包括:
根据多路径传载结构的已知设计空间,对集中载荷施加点的单传力路径和双传力路径进行工程分析,选出最优传力路径;
对比拓扑优化最佳传力路线和最优传力路径;
若两者符合,则认为拓扑优化最佳传力路线合理;若两者不符合,则修改拓扑优化的约束参数再次进行优化,直到两者符合。
变密度法将结构拓扑优化问题转换为以单元材料密度为优化变量的非线性优化问题;非线性问题能转换为二次规划问题,即序列规划方法,序列规划方法能通过灵敏度分析,寻找最优方向来更新单元密度分布,最终获得优化收敛结果。
静态求解有限元模型的网格尺寸是3mm,单元类型采用2次修正的四面体单元。
有限元模型的网格尺寸是6mm,单元类型采用2次修正的四面体单元。
缘条稳定性约束:
式中,E、μe为材料的弹性模量、泊松比;b、δ为加载边的宽度、加载边的宽度,Kc为压缩临界应力系数,η为塑性修正系数。
支柱刚度约束:
a)支柱截面积大于腹板面积40%;
Amin=0.4×A;
A、Amin是指隔波构件截面积、壁板构件截面积;
b)支柱构件满足最小惯性矩要求;
Iu、Imin指隔波构件惯性矩、最小惯性矩;d、t、he、ks是指支柱间距、腹板厚度、上下缘条形心间的距离、剪切屈曲系数。
有益效果:创新性地对大载荷、高应力、多路径传载进行传力路线拓扑优化;对较优传力方案进行尺寸优化,得到较优的结构尺寸参数;优化结构进行精确结构参数表征,有效减轻结构重量。
附图说明
图1为多路径传载结构的结构示意图;
图2为两种传力路径的示意图;
图3为拓扑优化有限元模型的示意图;
图4为拓扑优化流程图;
图5为拓扑优化材料分布图;
图6为尺寸优化结构尺寸分布图;
图7为设计出的多路径传载结构的结构图;
图8为设计出的多路径传载结构的优化结构静力求解结果的示意图;
图9为多路径传载结构的最佳结构参数图。
具体实施方式
某大型飞机多路径传载结构[4],此结构通过集中点[1]传递集中载荷P,此结构两个约束位置为集中点[2]和集中点[3],见图1。
(1)基于图1结构进行工程受力分析,集中载荷施加点P有两种传力路径:一是集中剪力Q和对集中点[5]的弯矩M均通过集中点[5]进行传递,二是弯矩M通过集中点[4]剪力Q1、集中点[5]剪力Q2形成的一对力偶进行平衡,剪力通过两点的剪力差进行平衡,见图2;
Q=P,M=PL1
式中L1、L2是指集中载荷施加点到集中点[5]、集中点[4]的水平距离。
(2)对两种不同的传力路径进行传力路线的拓扑优化分析,形成较优传力路线:以结构能量最小化为目标函数,以结构体积比(设置优化前后体积比为0.2)、集中点[1]的变形(点[1]的变形值为25m)作为约束条件,以结构XZ对称平面作为几何约束限制,区域1、区域2、区域3作为冻结边界区域,运用变密度法进行拓扑优化分析,模型离散采用2次修正的四面体单元,见图3,此数值拓扑优化模型的数学描述如下:
式中:C(x)指能量函数,V(x)、V0、f指优化前体积、优化后体积、优化前后体积比,Xc、Xmax是指优化变形、最大变形数值。
(3)应用变密度法将结构拓扑优化问题转换为以单元材料密度为优化变量的非线性优化问题。将非线性问题转换为一系列线性或二次规划问题,即序列规划方法。序列规划方法通过灵敏度分析,寻找最优方向来更新单元密度分布,最终获得优化收敛结果,具体拓扑优化流程见图4;
(4)基于步骤(3),得到满足约束条件的拓扑材料优化分布结果,见图5;
(5)图5结构,进行尺寸优化分析,以结构重量最小化为目标函数,以体积比(体积比0.1)和集中点[1]变形(最大变形25mm)作为约束条件,并设置最大/最小原件尺寸限制(单元最小尺寸1.5mm、最大尺寸3mm)及XZ平面对称结构限制,进行尺寸优化,数学表述如下,优化流程同(3),较优的结构参数方案见图6;
式中,wij(x1,…,xk)、V0、Vf、uf、δ1、δ、δ2是指重量函数、优化前体积、优化后体积、优化后变形数值、优化最小厚度、优化厚度、优化最大厚度。
(6)依据拓扑优化结果、尺寸优化结果进行结构设计,并对图7结构进行静力应力求解,见图8,得到满足缘条稳定性约束和支柱刚度约束最终结构参数,见图9,静力稳定性约束和支柱刚度约束数学描述如下:
缘条稳定性约束:
式中,E、μe为材料的弹性模量、泊松比;b、δ为加载边的宽度、加载边的宽度,Kc为压缩临界应力系数,η为塑性修正系数。
支柱刚度约束:
c)支柱截面积大于腹板(两支柱构件之间)面积40%;
Amin=0.4×A
A、Amin是指隔波构件截面积、壁板构件截面积。
d)支柱构件满足最小惯性矩要求;
Iu、Imin指隔波构件惯性矩、最小惯性矩;d、t、he、ks是指支柱间距、腹板厚度、上下缘条形心间的距离、剪切屈曲系数。
此结果已经得到试验良好验证。
某大型民机的增生装置,通过支臂结构将30吨载荷传递至盒段结构,支臂结构的形式及重量直接决定了其方案的可实施性。
Claims (5)
1.一种大载荷、高应力、多路径传载结构的优化方法,其特征在于,包括:
根据多路径传载结构的已知设计空间,建立其有限元模型;
对该有限元模型进行变密度法的拓扑优化,得到拓扑优化最佳传力路线及最佳传力模型;其中,以结构能量最小化为目标函数,以第一结构体积比、集中载荷施加点的变形作为约束条件,以结构对称平面作为几何约束限制,冻结边界施加区域;
基于最佳传力模型,以结构重量最小化为目标函数,以第二结构体积比、集中载荷施加点的变形作为约束条件,并考虑加工工艺设置最大/最小原件尺寸限制,以结构对称平面作为几何约束限制,进行多路径传载的尺寸优化,得到多路径传载结构的尺寸优化参数;
得到多路径传载结构的尺寸优化参数之后,所述方法还包括:
根据多路径传构的载结尺寸优化参数,设计多路径传载结构的具体结构;
根据具体结构,建立静态求解有限元模型,并静力求解,得到结构应力应变分布图;
基于结构应力应变分布图,对具体结构进行缘条稳定性、支柱刚度验证;
若某个结构不满足任一验证条件,则对该结构根据该验证条件进行局部修正;最终修正的参数和无需修正的参数作为多路径传载结构的最佳结构参数;
缘条稳定性约束:
式中,E、μe为材料的弹性模量、泊松比;b、δ为加载边的宽度、加载边的宽度,Kc为压缩临界应力系数,η为塑性修正系数;
支柱刚度约束:
a)支柱截面积大于腹板面积40%;
Amin=0.4×A;
A、Amin是指隔波构件截面积、壁板构件截面积;
b)支柱构件满足最小惯性矩要求;
Iu、Imin指隔波构件惯性矩、最小惯性矩;d、t、he、ks是指支柱间距、腹板厚度、上下缘条形心间的距离、剪切屈曲系数。
2.根据权利要求1所述的方法,其特征在于,对该有限元模型进行进行变密度法的拓扑优化之后,所述方法还包括:
根据多路径传载结构的已知设计空间,对集中载荷施加点的单传力路径和双传力路径进行工程分析,选出最优传力路径;
对比拓扑优化最佳传力路线和最优传力路径;
若两者符合,则认为拓扑优化最佳传力路线合理;若两者不符合,则修改拓扑优化的约束参数再次进行优化,直到两者符合。
3.根据权利要求1所述的方法,其特征在于,变密度法将结构拓扑优化问题转换为以单元材料密度为优化变量的非线性优化问题;非线性问题能转换为二次规划问题,即序列规划方法,序列规划方法能通过灵敏度分析,寻找最优方向来更新单元密度分布,最终获得优化收敛结果。
4.根据权利要求1所述的方法,其特征在于,静态求解有限元模型的网格尺寸是3mm,单元类型采用2次修正的四面体单元。
5.根据权利要求1所述的方法,其特征在于,有限元模型的网格尺寸是6mm,单元类型采用2次修正的四面体单元。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011612688.0A CN112699484B (zh) | 2020-12-29 | 2020-12-29 | 一种大载荷、高应力、多路径传载结构的优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011612688.0A CN112699484B (zh) | 2020-12-29 | 2020-12-29 | 一种大载荷、高应力、多路径传载结构的优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112699484A CN112699484A (zh) | 2021-04-23 |
CN112699484B true CN112699484B (zh) | 2022-11-22 |
Family
ID=75512669
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011612688.0A Active CN112699484B (zh) | 2020-12-29 | 2020-12-29 | 一种大载荷、高应力、多路径传载结构的优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112699484B (zh) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018126465A1 (zh) * | 2017-01-09 | 2018-07-12 | 大连理工大学 | 一种用于消除薄膜结构拉伸褶皱的优化设计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100989190B1 (ko) * | 2008-08-29 | 2010-10-20 | 한양대학교 산학협력단 | 등가정하중을 이용한 위상최적설계방법 |
-
2020
- 2020-12-29 CN CN202011612688.0A patent/CN112699484B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018126465A1 (zh) * | 2017-01-09 | 2018-07-12 | 大连理工大学 | 一种用于消除薄膜结构拉伸褶皱的优化设计方法 |
Non-Patent Citations (2)
Title |
---|
基于变密度法的大型天线系统支撑结构优化设计;袁镭等;《机械设计与研究》;20130220(第01期);全文 * |
某型号离心机吊篮拓扑优化设计;王东升等;《航天器环境工程》;20090615(第03期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112699484A (zh) | 2021-04-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Cavagna et al. | NeoCASS: an integrated tool for structural sizing, aeroelastic analysis and MDO at conceptual design level | |
CN105701297B (zh) | 一种基于多点自适应代理模型的反射面天线机电耦合设计方法 | |
CN111027110B (zh) | 一种连续体结构拓扑与形状尺寸综合优化方法 | |
CN112100774A (zh) | 一种基于变密度法的应力和应变能双约束的拓扑优化方法 | |
CN107341316B (zh) | 设计相关压力载荷作用下的结构形状-拓扑联合优化方法 | |
CN109344524B (zh) | 一种薄板结构加强筋分布优化方法 | |
CN110222356B (zh) | 综合考虑稳定性与振动特性的板/壳结构轻量化拓扑优化设计方法 | |
CN110569519B (zh) | 考虑非设计域的三维连续体结构动静力学性能拓扑优化设计方法 | |
CN103646131A (zh) | 一种考虑气动弹性约束的复合材料机翼多目标优化设计方法 | |
CN103440378B (zh) | 基于应力约束的机翼翼梁结构拓扑优化方法 | |
CN109255144B (zh) | 一种起落架刹车构件轻量化设计方法 | |
CN110334427B (zh) | 一种基于工程梁理论的tbw布局飞机机翼结构分析方法 | |
CN105447269A (zh) | 一种飞机在定常过载状态下的非结构质量配平方法 | |
Gray et al. | Geometrically nonlinear high-fidelity aerostructural optimization for highly flexible wings | |
CN110704953B (zh) | 一种大展弦比机翼静气弹性能设计敏度的分析方法 | |
CN112784359B (zh) | 薄壁梁约束扭转极限承载力迭代计算方法 | |
CN107341309B (zh) | 一种基于垂尾载荷的机身与尾翼连接铰点载荷分配方法 | |
CN107526866B (zh) | 基于特征驱动的翼面结构拓扑优化方法 | |
Nejat et al. | Topology optimization of lightweight periodic lattices under stiffness and stability constraints | |
CN104978450A (zh) | 一种直升机振动主动控制位置优选方法 | |
CN109902350B (zh) | 对变截面梁的截面惯性矩进行模型修正中克服模态交换的方法 | |
Fugate et al. | Aero-Structural Modeling of the Truss-Braced Wing Aircraft Using Potential Method with Correction Methods for Transonic Viscous Flow and Wing-Strut Interference Aerodynamics | |
CN112699484B (zh) | 一种大载荷、高应力、多路径传载结构的优化方法 | |
Zhang et al. | A morphing wing with cellular structure of non-uniform density | |
CN101612996A (zh) | 一种板筋结构设计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |