CN111597641A - 一种适用于热塑性编织复合材料航天承力构件的材料-结构-制造一体化优化设计方法 - Google Patents
一种适用于热塑性编织复合材料航天承力构件的材料-结构-制造一体化优化设计方法 Download PDFInfo
- Publication number
- CN111597641A CN111597641A CN202010444328.8A CN202010444328A CN111597641A CN 111597641 A CN111597641 A CN 111597641A CN 202010444328 A CN202010444328 A CN 202010444328A CN 111597641 A CN111597641 A CN 111597641A
- Authority
- CN
- China
- Prior art keywords
- optimization
- composite material
- model
- design
- optimization design
- 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.)
- Granted
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
- G06F2113/00—Details relating to the application field
- G06F2113/22—Moulding
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Moulding By Coating Moulds (AREA)
- Casting Or Compression Moulding Of Plastics Or The Like (AREA)
Abstract
本发明提供一种适用于热塑性编织复合材料航天承力构件的材料‑结构‑制造一体化优化设计方法,包括:根据纤维体分比、纤维和基体性能参数,利用细观力学方法预测纤维增强复合材料的等效宏观弹性模量;进行复合材料标准拉伸实验测试,基于弹性常数和泊松比修正预测的等效宏观弹性模量;建立复合材料结构件的有限模型并对结构件外形尺寸和厚度进行参数化;构建优化模型,采用高效的自适应加点代理模型优化方法进行结构件的优化设计,直至优化结果收敛,得到构件的最优设计;结合成型工艺修改优化的构件构型,模压成型制备结构件。本发明优化设计效率高,缩短复合材料结构的设计周期,减少结构设计成本,为复合材料优化设计提供一种快速有效的方法。
Description
技术领域
本发明涉及热塑性复合材料构件优化设计领域,尤其涉及一种适用于热塑性编织复合材料航天承力构件的材料-结构-制造一体化优化设计方法。
背景技术
纤维增强复合材料具有优异的比强度、刚度、材料的各向异性、耐腐蚀性、抗疲劳、减震、吸能等特性,使其在现代航空航天、汽车领域的应用迅速增加。纤维增强复合材料表现出来的突出特性不仅被用来减轻结构的重量,也被用来提高结构的整体性能。在纤维复合材料之中,编织复合材料被认为比单向复合材料更有竞争力,因为它们在单层内具有良好的四面八方的控制力。而热塑性纤维增强复合材料与热固性复合材料相比密度小、比强度大、抗冲击性能优异、耐腐蚀、湿热性能好,热塑性复合材料储存方便、成型效率高、可熔融焊接和减少制造周期,可以多次加热熔融,重复使用,节约成本,些特性使热塑性机织物复合材料对许多结构应用具有吸引力。为适应当前社会各种领域对轻量化和整体性能的要求,并且随着热塑性纤维编织复合材料成型技术的不断发展,纤维编织复合材料将会在主承力结构中有更为广泛的应用。
由于热塑性纤维增强复合材料结构设计的灵活性,设计者可以根据不同的载荷工况建立优化模型来优化结构以达到更高的性能。然而现有的大多数轻量化研究中只是根据工程经验调整纤维分布、纤维含量、层合板铺层来对金属结构进行替换。但这种优化方式轻量化效果不明显,没有充分发挥复合材料灵活的可设计性。热塑性编织复合材料结构件优化设计不仅要实现对金属材料的简单替换,而且要针对具体工况建立优化模型,实现材料和结构件外形的优化设计。同时复合材料构件显式动力学仿真计算时间较长,启发式优化算法需要计算大量样本导致设计周期过长,高效的梯度类优化算法难以获得全局最优化解,都不适用于工程需求。在此应用背景下,本发明提出一种适用于热塑性编织复合材料航天承力构件的材料-结构-制造一体化快速高效的优化设计方法。
发明内容
根据上述提出的现有结构优化技术存在的热塑性复合材料轻量化不明显和结构计算优化效率较低的问题,而提供一种适用于热塑性编织复合材料航天承力构件的材料-结构-制造一体化优化设计方法。本发明采用的技术手段如下:
一种适用于热塑性编织复合材料航天承力构件的材料-结构-制造一体化优化设计方法,包括如下步骤:
步骤一:根据纤维体分比、纤维和基体性能参数,利用细观力学方法预测纤维增强复合材料的等效宏观弹性模量;
步骤二:进行复合材料标准拉伸实验测试,基于实验得到弹性常数和泊松比修正预测的等效宏观弹性模量;
步骤三:建立复合材料结构件的有限模型并对结构件外形尺寸和厚度进行参数化,实现结构件自动化建模分析和结果提取;
步骤四:构建优化模型,采用高效的自适应加点代理模型优化方法进行结构件的优化设计,直至优化结果收敛,得到结构件的最佳设计。
进一步地,所述步骤一中,预测的宏观弹性模量是通过将缎纹复合材料简化为正交层合板,先根据三维桥联模型公式,计算单向板的等效宏观弹性模量,根据层合板理论计算纤维增强复合材料的等效宏观弹性模量。
进一步地,所述步骤三中,自动化建模分析和结果提取是基于ABAQUS二次开发结构件的参数化程序。
进一步地,所述步骤四中,优化模型基于工程实际工况中的载荷条件进行设计,其中优化目标和约束包括结构件的强度、刚度及质量。
进一步地,所述的步骤四中自适应加点代理模型优化方法包括如下步骤:
通过拉丁超立方抽样方法生成预设量的初始样本点,选用克里金代理模型拟合设计变量-目标和设计变量-约束的数学关系;
采用综合目标最优和代理模型误差最大的期望最大加点准则确定新样本点,并计算目标响应值和约束值,添加到样本集中并更新代理模型;
重复优化过程直至满足收敛准则之一,所述收敛准则包括:新样本点与已有样本点距离达到阈值、最优解与上一代最优解的差达到阈值、迭代次数达到最大次数。
进一步地,所述步骤四还设有如下步骤:
步骤五:结合成型工艺修改优化的构件构型,模压成型制备结构件。
与现有技术相比,本发明具有以下优点:
本发明提出了一种适用于编织热塑性复合材料航天承力构件的材料-结构-制造一体化设计方法,基于理论公式进行等效弹性模量的预测并根据实验修正;针对具体工况建立优化模型,实现材料和结构件外形的优化设计,对复合材料结构件外形进行全体尺寸参数化,极大的提高了轻量化设计的自由度;同时使用的自适应加点代理模型优化方法采用样本点进行动态更新,减少了优化过程所需要的样本个数,极大的提高了计算效率。本发明为热塑性复合材料工程制造提供了一种高效的优化设计方法,更加符合实际工程制造实际要求,极大的缩短了预研周期,产生更高的经济效益。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图做以简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明一体化优化设计方法流程图。
图2为本发明实施例中桁条加固筒壳整体模型示意图。
图3为本发明实施例中3种复合材料桁条截面参数示意图,其中(a)为T型、(b)为L型和(c)为Ω型。
图4为本发明实施例中自适应加点代理模型优化方法优化框架示意图。
图5为本发明实施例中3种桁条最优的截面尺寸示意图,其中(a)为T型、(b)为L型和(c)为Ω型。
图6为本发明实施例中Ω型桁条考虑模压制备工艺修饰的截面示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,本发明公开了一种适用于热塑性编织复合材料航天承力构件的材料-结构-制造一体化优化设计方法,包括如下步骤:
步骤一:根据纤维体分比、纤维和基体性能参数,利用细观力学方法预测纤维增强复合材料的等效宏观弹性模量;
步骤二:进行复合材料标准拉伸实验测试,基于实验得到弹性常数和泊松比修正预测的等效宏观弹性模量;
步骤三:建立复合材料结构件的有限模型并对结构件外形尺寸和厚度进行参数化,实现结构件自动化建模分析和结果提取;
步骤四:构建优化模型,采用高效的自适应加点代理模型优化方法进行结构件的优化设计,直至优化结果收敛,得到结构件的最佳设计。
步骤五:结合成型工艺修改优化的构件构型,模压成型制备结构件。
所述步骤一中,本方法所预测的宏观弹性模量是通过将缎纹复合材料简化为正交层合板,先根据三维桥联模型公式,计算单向板的等效宏观弹性模量(如吕毅,吕国志,吕胜利.细观力学方法预测单向复合材料的宏观弹性模量[J].西北工业大学学报(6):787-790);再根据层合板理论,计算正交层合板面内柔度系数;根据工程常数的定义式1.11和1.12计算纤维增强复合材料的等效宏观弹性模量。
桥联矩阵Cij和单向板的等效宏观弹性模量如下:
c23=0(0.4)
对称层合板的面内应力和面内里的关系为
根据工程常数的定义,层合板用工程弹性常数表示的应变-面内力关系为:
所述步骤三中,自动化建模分析和结果提取是基于ABAQUS二次开发结构件的参数化程序。
所述步骤四中,优化模型基于工程实际工况中的载荷条件进行设计,其中优化目标和约束包括结构件的强度、刚度及质量。
所述的步骤四中自适应加点代理模型优化方法包括如下步骤:
通过拉丁超立方抽样方法生成预设量的初始样本点,选用克里金代理模型拟合设计变量-目标和设计变量-约束的数学关系;
采用综合目标最优和代理模型误差最大的期望最大加点准则确定新样本点,并计算目标响应值和约束值,添加到样本集中并更新代理模型;
重复优化过程直至满足收敛准则之一,所述收敛准则包括:新样本点与已有样本点距离达到阈值、最优解与上一代最优解的差达到阈值、迭代次数达到最大次数,动态代理模型可以达到减少优化过程所需要的样本个数,增加优化效率的目的。
实施例1
本实施例以为热塑性缎纹复合材料桁条金属筒壳的桁条结构设计,步骤如下:
1)根据预设的50%纤维体分比和表1的T300碳纤维和聚醚醚酮(PEEK)性能参数,利用式1.1到1.12预测纤维增强复合材料的等效宏观弹性模量,结果为拉伸模量E11=62.74GPa,E22=62.74GPa,剪切模量G12=3.46GPa,G23=3.3GPa,G13=3.3GPa,泊松比μ12=0.0419。
表1热塑性复合材料组分性能参数
2)进行标准拉伸实验进行复合材料基本力学性能测试得到弹性常数和泊松比,复合材料的拉伸强度为674.27MPa,拉伸弹性模量61.80GPa,泊松比为0.037;基于实验数据修正等效宏观弹性模量,结果为E11=61.8GPa,E22=61.8GPa,G12=3.46GPa,G23=3.3GPa,G13=3.3GPa,μ12=0.0375。
3)建立热塑性复合材料蒙皮桁条筒壳有限元模型,筒壳蒙皮材料为铝合金,桁条为碳纤维增强聚醚醚酮复合材料,如图2。
基于ABAQUS二次开发结构件的参数化分析程序,对L型、T型、Ω型3种桁条外形尺寸进行参数化,采用显式动力学对筒壳结构进行非线性屈曲分析,计算桁条筒壳结构的轴压极限承载力,实现结构件自动化建模、分析和后处理。如图3所示,复合材料T型(a)、L型(b)和Ω型(b)截面参数化示意图,保证三种桁条的初始质量一样,为1.925kg。
4)针对3种蒙皮桁条结构类型,以桁条高度、桁条宽度、桁条厚度、桁条数目为设计变量,整体桁条质量不超过初始质量为约束条件,开展金属蒙皮复合材料桁条筒壳结构轴压承载力最大的优化设计。并建立优化列式(13),设计变量为桁条高度H、桁条宽度L、桁条厚度t、桁条数目N,桁条总质量1.925kg为质量上限约束条件。
采用高效的自适应加点代理模型优化方法进行结构件的优化设计,优化框架如图4所示,直至优化结果收敛,得到3种桁条截面的最佳设计,如图5(a)(b)(c)。T型桁条的最优极限承载力为1503.59kN,L型桁条的最优极限承载力为1466.86kN,Ω型桁条的最优承载力为1563.28kN,三种型式的桁条承载力:Ω型桁条>T型桁条>L型桁条。在三种截面的桁条中,Ω型桁条筒壳的承载力最大,Ω型桁条的抗弯刚度最大,使得蒙皮桁条筒壳的整体刚度最大。
5)结合成型工艺修改优化的构件构型如图6所示,最后模压成型制备结构件。
最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (6)
1.一种适用于热塑性编织复合材料航天承力构件的材料-结构-制造一体化优化设计方法,其特征在于,包括如下步骤:
步骤一:根据纤维体分比、纤维和基体性能参数,利用细观力学方法预测纤维增强复合材料的等效宏观弹性模量;
步骤二:进行复合材料标准拉伸实验测试,基于实验得到弹性常数和泊松比修正预测的等效宏观弹性模量;
步骤三:建立复合材料结构件的有限模型并对结构件外形尺寸和厚度进行参数化,实现结构件自动化建模分析和结果提取;
步骤四:构建优化模型,采用高效的自适应加点代理模型优化方法进行结构件的优化设计,直至优化结果收敛,得到结构件的最佳设计。
2.根据权利要求1所述的适用于热塑性编织复合材料航天承力构件的材料-结构-制造一体化优化设计方法,其特征在于,所述步骤一中,预测的宏观弹性模量是通过将缎纹复合材料简化为正交层合板,先根据三维桥联模型公式,计算单向板的等效宏观弹性模量,根据层合板理论计算纤维增强复合材料的等效宏观弹性模量。
3.根据权利要求1所述的适用于热塑性编织复合材料航天承力构件的材料-结构-制造一体化优化设计方法,其特征在于,所述步骤三中,自动化建模分析和结果提取是基于ABAQUS二次开发结构件的参数化程序。
4.根据权利要求1所述的适用于热塑性编织复合材料航天承力构件的材料-结构-制造一体化优化设计方法,其特征在于,所述步骤四中,优化模型基于工程实际工况中的载荷条件进行设计,其中优化目标和约束包括结构件的强度、刚度及质量。
5.根据权利要求1或4所述的适用于热塑性编织复合材料航天承力构件的材料-结构-制造一体化优化设计方法,其特征在于,所述的步骤四中自适应加点代理模型优化方法包括如下步骤:
通过拉丁超立方抽样方法生成预设量的初始样本点,选用克里金代理模型拟合设计变量-目标和设计变量-约束的数学关系;
采用综合目标最优和代理模型误差最大的期望最大加点准则确定新样本点,并计算目标响应值和约束值,添加到样本集中并更新代理模型;
重复优化过程直至满足收敛准则之一,所述收敛准则包括:新样本点与已有样本点距离达到阈值、最优解与上一代最优解的差达到阈值、迭代次数达到最大次数。
6.根据权利要求1所述的适用于热塑性编织复合材料航天承力构件的材料-结构-制造一体化优化设计方法,其特征在于,所述步骤四还设有如下步骤:
步骤五:结合成型工艺修改优化的构件构型,模压成型制备结构件。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010444328.8A CN111597641B (zh) | 2020-05-22 | 2020-05-22 | 一种适用于热塑性编织复合材料航天承力构件的材料-结构-制造一体化优化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010444328.8A CN111597641B (zh) | 2020-05-22 | 2020-05-22 | 一种适用于热塑性编织复合材料航天承力构件的材料-结构-制造一体化优化设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111597641A true CN111597641A (zh) | 2020-08-28 |
CN111597641B CN111597641B (zh) | 2023-07-04 |
Family
ID=72189236
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010444328.8A Active CN111597641B (zh) | 2020-05-22 | 2020-05-22 | 一种适用于热塑性编织复合材料航天承力构件的材料-结构-制造一体化优化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111597641B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116644619A (zh) * | 2023-07-27 | 2023-08-25 | 北京理工大学 | 纤维预制体成型工艺的压力优化方法、装置、设备及介质 |
CN117786467A (zh) * | 2024-02-28 | 2024-03-29 | 上海交通大学四川研究院 | 基于自适应加点的飞机着陆风险预测的分类模型构建方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107742005A (zh) * | 2017-09-01 | 2018-02-27 | 杭州健途科技有限公司 | 一种纤维增强复合材料结构力学性能预测和控制方法 |
CN108710736A (zh) * | 2018-05-10 | 2018-10-26 | 北京航空航天大学 | 一种考虑应力约束的宏微一体结构拓扑优化方法 |
-
2020
- 2020-05-22 CN CN202010444328.8A patent/CN111597641B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107742005A (zh) * | 2017-09-01 | 2018-02-27 | 杭州健途科技有限公司 | 一种纤维增强复合材料结构力学性能预测和控制方法 |
CN108710736A (zh) * | 2018-05-10 | 2018-10-26 | 北京航空航天大学 | 一种考虑应力约束的宏微一体结构拓扑优化方法 |
Non-Patent Citations (2)
Title |
---|
吕毅 等: "细观力学方法预测单向复合材料的宏观弹性模量", 西北工业大学学报, vol. 24, no. 6, pages 787 - 790 * |
张春春 等: "横观各向同性基体复合材料的等效弹性常数", 应用数学和力学, vol. 39, no. 7, pages 750 - 765 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116644619A (zh) * | 2023-07-27 | 2023-08-25 | 北京理工大学 | 纤维预制体成型工艺的压力优化方法、装置、设备及介质 |
CN116644619B (zh) * | 2023-07-27 | 2023-09-22 | 北京理工大学 | 纤维预制体成型工艺的压力优化方法、装置、设备及介质 |
CN117786467A (zh) * | 2024-02-28 | 2024-03-29 | 上海交通大学四川研究院 | 基于自适应加点的飞机着陆风险预测的分类模型构建方法 |
CN117786467B (zh) * | 2024-02-28 | 2024-04-30 | 上海交通大学四川研究院 | 基于自适应加点的飞机着陆风险预测的分类模型构建方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111597641B (zh) | 2023-07-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107103138B (zh) | 一种激光喷丸变刚度轻量化方法 | |
CN109766656B (zh) | 一种基于拓扑优化的梯度点阵结构设计方法 | |
CN110222356B (zh) | 综合考虑稳定性与振动特性的板/壳结构轻量化拓扑优化设计方法 | |
CN111597641A (zh) | 一种适用于热塑性编织复合材料航天承力构件的材料-结构-制造一体化优化设计方法 | |
CN108829914B (zh) | 一种frp结构件的结构与工艺一体化设计方法 | |
Yan et al. | Concurrent multi-scale design optimization of composite frames with manufacturing constraints | |
CN113191077B (zh) | 一种基于连续纤维复材3d打印的变纤维含量拓扑优化方法 | |
Kanou et al. | Numerical modeling of stresses and buckling loads of isogrid lattice composite structure cylinders | |
Liu et al. | The compressive response of octet lattice structures with carbon fiber composite hollow struts | |
Qi et al. | A hybrid joining insert for sandwich panels with pyramidal lattice truss cores | |
CN112131770B (zh) | 一种考虑可靠性的功能梯度连续体结构轻量化设计方法 | |
CN101612996A (zh) | 一种板筋结构设计方法 | |
CN102354324B (zh) | 跨音速颤振模型复合材料单梁结构设计及其刚度计算方法 | |
Mahdi et al. | Stiffness Tailoring for Improved Buckling of Variable Angle Tow Composite Sandwich Plates | |
CN113515850A (zh) | 一种考虑纤维连续的纤维增强复合材料结构布局优化设计方法 | |
Liu et al. | Optimization of composite stiffened panels subject to compression and lateral pressure using a bi-level approach | |
Sun et al. | A discrete material optimization method with a patch strategy based on stiffness matrix interpolation | |
Aiello et al. | Buckling load design of sandwich panels made with hybrid laminated faces and transversely flexible core | |
Beemaraj et al. | Inverse Multi-Scale Robust Design of Composite Structures Using Design Capability Indices | |
Visnjic et al. | Shear-lag influence on maximum specific bending stiffness and strength of composite I-beam wing spar | |
CN111353229B (zh) | 一种实体结构光滑粒子动力学建模方法 | |
Coburn et al. | Local buckling of blade stiffened variable angle tow panels | |
CN110704954A (zh) | 变体飞行器柔性蒙皮胞状支撑体的拓扑优化方法及系统 | |
Wang et al. | Structure Analysis and Aeroelastic Optimization Method for Composite Wing-box Structure with Curved Fiber | |
Zhou | Quasi-static mechanical properties of optimized SC lattice structure based on periodic boundary theory |
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 |