CN107451308B - 一种复杂复合材料结构等效热传导系数多尺度计算方法 - Google Patents
一种复杂复合材料结构等效热传导系数多尺度计算方法 Download PDFInfo
- Publication number
- CN107451308B CN107451308B CN201610373839.9A CN201610373839A CN107451308B CN 107451308 B CN107451308 B CN 107451308B CN 201610373839 A CN201610373839 A CN 201610373839A CN 107451308 B CN107451308 B CN 107451308B
- Authority
- CN
- China
- Prior art keywords
- microscopic
- scale
- finite element
- heat conduction
- model
- 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.)
- Expired - Fee Related
Links
Images
Classifications
-
- 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]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明提出一种复杂复合材料结构等效热传导系数多尺度计算方法,采用尺度分离的方法,将宏观、细观、微观三尺度结构分离,根据不同尺度模型的几何特征,分别建立各个尺度分析模型;将三尺度问题转化为两个多尺度问题:宏观‑细观多尺度问题、细观‑微观多尺度问题,依次对着两个多尺度问题进行分析,将微观多尺度问题得到的等效模量最终返回给宏观多尺度问题。克服了传统结构分析方法计算效率低、精度差的缺点,有效提升了复合材料结构性能预测的效率和精度,使其可以用于指导复合材料的生产、研发等工作。本发明可应用于航空航天领域复杂复合材料结构设计、分析,以及其他复合材料工程领域的结构设计热、力学分析问题。
Description
技术领域
本发明涉及复合材料设计领域,是一种复杂复合材料结构分析设计方法,具体为一种复杂复合材料结构等效热传导系数多尺度计算方法。
背景技术
复合材料由于质量轻、强度高、具有较强的可设计性等特点,广泛应用与航空航天的结构器件中。而由于复合材料结构复杂,为了研究复合材料的性能,提升复合材料结构件的使用效率,近百年来,国内外学者提出了大量用以预测复合材料行为的理论。其核心是通过求解控制方程,从而确定结构内部位移、温度等物理量的分布,从而完成对材料性能的预测。
目前,复合材料材料性能预测方法主要分为四类:
第一类是解析法,其代表的方法有:稀疏法、Mori-Tanaka方法、自洽方法、广义自洽法。该类方法通过求解无限大基体内单夹杂或多夹杂问题,得到远场应变与单个夹杂平均应变之间关系,从而得到材料的有效模量。这种方法理论较为简单,但由于实际复合材料存在一定的边界,边界效应会使得计算结果产生一定的误差,此外,部分解析方法只适用于结构简单,体积分数较低的复合材料,这使得这类方法在复合材料性能预测上存在一定的局限。
第二类是半解析法,其代表的方法为变换场分析法。该方法在细观采用显式的本构关系来联系宏观与细观场,该方法需要给定均匀化与局部化的规则,对于多相材料以及非线性的非均匀材料,该方法等效本构关系所需要的内部变量会非常多,限制了该方法的应用。
第三类是数值方法,其代表的方法为数值均匀化方法,该方法将复合材料转化为一个多尺度分析问题,通过局部化和均匀化的方法,建立宏观积分点与细观体积代表单元之间的联系,从而完成材料性能的预测,与解析法相比,该方法的计算量较小,并且由于在计算时考虑到材料的细观形貌,因此计算精度较高。
然而现有的多尺度方法仅考虑了两尺度的信息,由于大多数复合材料采用铺层的形式,材料的细观尺度并非纤维与基体简单的组合形式,而是多种纤维与基体的组合,纤维的铺层角,每次纤维的排布形式都会很大程度的影响宏观材料的性能。此外,大多数多尺度分析软件用户开发由国外航空航天科研机构开发,由于种种原因这些软件都没有对外公开,并且大多数用于学术研究的程序因精度、计算代价等问题,限制了其在工程领域内的应用。
发明内容
为了避免现有技术的不足之处,本发明提出了一种复杂复合材料结构等效热传导系数多尺度计算方法,在方法中采用了复合材料结构分析的三尺度模型,由于该方法考虑了细观、微观结构,从而使得宏观结果计算精度得以提升;此外,该方法通过商业有限元软件ABAQUS的二次开发实现,从而增加了其通用性,使其能够更好地解决大规模的工程问题。
本发明的技术方案为:
所述一种复杂复合材料结构等效热传导系数多尺度计算方法,其特征在于:包括以下步骤:
步骤1:按照复合材料实际尺度建立宏观有限元分析模型,宏观有限元分析模型材料坐标系为(X1,X2,X3);通过显微CT扫描实验,得到复合材料细观结构的物理模型,根据复合材料细观结构物理模型的体积分数、增强相与基体相的几何特征以及排布形式、缺陷位置、铺层数量和铺层角信息,建立细观有限元模型,细观有限元模型材料坐标系记为(Y1,Y2,Y3);通过电子显微镜实验,得到复合材料微观单胞的物理模型,根据复合材料微观单胞物理模型增强相的体积分数、形状、以及缺陷位置,建立微观有限元模型,微观有限元模型材料坐标系记为(Z1,Z2,Z3);其中Yi=Xi/ξ,Zi=Yi/η,i=1、2、3,ξ、η分别为宏观-细观,细观-微观尺度间的桥接系数,且满足ξ<<1,η<<1;
步骤2:根据需要计算的复合材料,赋予微观有限元模型材料属性;
步骤3:将多尺度分析分为两步,首先通过细观-微观两尺度分析,得到细观尺度的等效热传导系数;根据细观尺度的等效热传导系数,通过宏观-细观两尺度分析,得到宏观结构的等效热传导系数:
步骤3.1:在周期性假设的条件下,将微观有限元模型的温度渐进展开式带入稳态热传导控制方程
中,得到微观等效的热传导系数表达式:
步骤3.2:采用等效热应变加载,将步骤3.1中的微观等效的热传导系数表达式转化为:
步骤3.3:得到步骤3.2中等效的微观有限元模型热传导系数矩阵后,根据细观有限元模型内每个铺层的铺层角,依据经典层合板理论,得到细观有限元模型每层铺层的等效热传导系数矩阵,并依此对细观有限元模型的热传导系数矩阵进行组装,形成总热传导系数矩阵:
步骤3.4:将步骤3.3得到的细观有限元模型总热传导系数矩阵赋予宏观有限元分析模型中,并对宏观有限元分析模型施加载荷,得到宏观有限元分析模型的响应。
有益效果
本发明提出的三尺度复合材料分析方法,有益效果是:
1、利用了多尺度方法,在复合材料分析过程中,充分考虑了细观、微观结构几何形貌对于宏观结构的影响,与传统的复合材料分析手段相比有更好地精度,此外,对于损伤、失效判断,可以通过观测细观结构应力分布的变化,明确损伤机理。
2、通过建立三尺度模型,使得在复合材料等效性质的计算时,考虑了纤维方向、铺层厚度因素对于整体热传导系数的影响,而传统的两尺度方法忽略了细观尺度铺层的厚度、纤维的方向因素。因此与传统多尺度方法相比有更好的精度。
3、三尺度方法能够通过基于ABAQUS平台的二次开发得以实现,具有较好的适用性,从而促进了多尺度方法在工程材料计算领域的应用。
本发明的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1:本发明的计算流程图;
图2:实施例中某型压力容器几何模型;
图3:压力容器横截面图;
图4:CT扫描下预制体的微观形貌;
图5:简化后的细观有限元模型;
图6:CT扫描下无纬布的微观形貌;
图7:针刺、网胎微观计算模型;
图8:无纬布微观计算模型;
图9:压力容器边界条件。
具体实施方式
下面详细描述本发明的实施例,所述实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
本实施例以某型压力容器等效材料性能计算为例,按照本发明技术方案进行实施,给出了详细的实施过程。
步骤1:压力容器由碳/碳复合材料构成,细观结构的预制体由无纬布、±45°铺层、0°、90°铺层组合而成。按照算例实际尺寸,如图2和图3所示,圆柱体长20mm。在商用有限元软件ABAQUS中建立压力容器宏观有限元分析模型,宏观有限元分析模型材料坐标系为(X1,X2,X3)。通过CT扫描和电镜扫描,分别得到压力容器结构的真实细观、微观模型。如图4和图6所示。
根据针刺碳/碳复合材料预制件微结构显微照片分析,可以确定微结构细观单胞的基本形式。细观单胞由若干不同铺层角的无纬布和复合网胎叠层铺设而成,在厚度方向通过针刺纤维束进行增强。网胎纤维在面内是杂乱分布的,因此是一种面内准各向同性材料。针刺和网胎类似,也属于各项同性材料。0°无纬布、环向90°无纬布及斜向无纬布的纤维排列较为紧凑,纤维体积分数较大。据此建立细观有限元模型,如图5。细观有限元模型材料坐标系记为(Y1,Y2,Y3)。
根据显微扫描得到材料的两种体积分数,在ABAQUS有限元软件中建立两种微观单胞,单胞1纤维体积分数为50%,用来模拟纤维含量较少的网胎以及针刺微观有限元模型,如图7。单胞2体积分数为81%,用来模拟纤维较为紧凑的无纬布微观有限元模型,如图8。微观有限元模型材料坐标系记为(Z1,Z2,Z3)。
Yi=Xi/ξ,Zi=Yi/η,i=1、2、3,ξ、η分别为宏观-细观,细观-微观尺度间的桥接系数,且满足ξ<<1,η<<1。
步骤2:根据需要计算的复合材料,赋予微观有限元模型材料属性。
步骤3:将多尺度分析分为两步,首先通过细观-微观两尺度分析,得到细观尺度的等效材料属性;根据细观尺度的等效材料属性,通过宏观-细观两尺度分析,得到宏观结构的等效材料属性。
当计算的等效材料性能为等效刚度矩阵,步骤3的具体步骤为:
步骤3.1:在ABAQUS中施加Tie约束,从而实现周期性边界条件的施加,使得对应面位移相同。在周期性假设的条件下,将微观有限元模型的位移渐进展开式带入弹性力学控制方程
中,得到微观等效的刚度表达式:
其中,上角标代表微观有限元模型,下角标代表6个不同应力的方向,下角标k,l代表3个不同位移的方向,的上角标代表微观有限元模型均匀化,下角标代表刚度矩阵中6个不同的方向,Y代表单胞体积,为微观位移特征函数,与位移对应k,l代表3个不同的位移特征函数的方向,Cijkl为单一组分材料的弹性模量,δmk为Kronecker张量,且满足:
这里实施例中在ABAQUS中设置6个线性扰动分析步,从而完成热载荷在不同方向(11、22、33、12、13、23)的加载。
步骤3.2:采用等效热应力加载,将步骤3.1中的微观等效的刚度表达式转化为:
将微观单胞的各个方向分析结果均匀化后,由于各个等效热载荷均为单位1载荷,因此均匀化后,得到微观结构等效材料属性,如表1,表2,至此微观分析结束。
表1:多尺度方法得到的针刺和网胎微观模型等效刚度矩阵
表2:微观多尺度方法得到的无纬布的等效刚度矩阵
步骤3.3:得到步骤3.2中等效的微观有限元模型刚度矩阵后,根据细观有限元模型内每个铺层的铺层角,依据经典层合板理论,得到细观有限元模型每层铺层的等效刚度矩阵,并依此对细观有限元模型的刚度矩阵进行组装,形成总刚度矩阵:
本实施例中,45°度铺层材料属性如表4所示,-45°度铺层材料属性如表5所示,90°铺层材料属性如表3所示,细观喷管压力容器等效材料属性如表6所示。
表3:经过坐标转换后得到的细观90°无纬布等效模量
表4:经过坐标转换后得到的细观45°无纬布等效模量
表5:经过坐标转换后得到的细观-45°无纬布等效模量
表6:多尺度方法得到的细观模型等效刚度矩阵
步骤3.4:将步骤3.3得到的细观有限元模型总刚度矩阵赋予宏观有限元分析模型中,并对宏观有限元分析模型施加载荷,本实施例中,在宏观模型左端面施加2方向载荷50MPa,约束右端面11、22、33三方向位移以及转角如图9所示,得到宏观有限元分析模型的响应。
当计算的等效材料性能为等效热传导系数,步骤3的具体步骤为:
步骤3.1:在周期性假设的条件下,将微观有限元模型的温度渐进展开式带入稳态热传导控制方程
中,得到微观等效的热传导系数表达式:
步骤3.2:采用等效热应变加载,将步骤3.1中的微观等效的热传导系数表达式转化为:
步骤3.3:得到步骤3.2中等效的微观有限元模型热传导系数矩阵后,根据细观有限元模型内每个铺层的铺层角,依据经典层合板理论,得到细观有限元模型每层铺层的等效热传导系数矩阵,并依此对细观有限元模型的热传导系数矩阵进行组装,形成总热传导系数矩阵:
步骤3.4:将步骤3.3得到的细观有限元模型总热传导系数矩阵赋予宏观有限元分析模型中,并对宏观有限元分析模型施加载荷,得到宏观有限元分析模型的响应。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (1)
1.一种复杂复合材料结构等效热传导系数多尺度计算方法,其特征在于:包括以下步骤:
步骤1:按照复合材料实际尺度建立宏观有限元分析模型,宏观有限元分析模型材料坐标系为(X1,X2,X3);通过显微CT扫描实验,得到复合材料细观结构的物理模型,根据复合材料细观结构物理模型的体积分数、增强相与基体相的几何特征以及排布形式、缺陷位置、铺层数量和铺层角信息,建立细观有限元模型,细观有限元模型材料坐标系记为(Y1,Y2,Y3);通过电子显微镜实验,得到复合材料微观单胞的物理模型,根据复合材料微观单胞物理模型增强相的体积分数、形状、以及缺陷位置,建立微观有限元模型,微观有限元模型材料坐标系记为(Z1,Z2,Z3);其中Yi=Xi/ξ,Zi=Yi/η,i=1、2、3,ξ、η分别为宏观-细观,细观-微观尺度间的桥接系数,且满足ξ<<1,η<<1;
步骤2:根据需要计算的复合材料,赋予微观有限元模型材料属性;
步骤3:将多尺度分析分为两步,首先通过细观-微观两尺度分析,得到细观尺度的等效热传导系数;根据细观尺度的等效热传导系数,通过宏观-细观两尺度分析,得到宏观结构的等效热传导系数:
步骤3.1:在周期性假设的条件下,将微观有限元模型的温度渐进展开式带入稳态热传导控制方程
中,得到微观等效的热传导系数表达式:
步骤3.2:采用等效热应变加载,将步骤3.1中的微观等效的热传导系数表达式转化为:
步骤3.3:得到步骤3.2中等效的微观有限元模型热传导系数矩阵后,根据细观有限元模型内每个铺层的铺层角,依据经典层合板理论,得到细观有限元模型每层铺层的等效热传导系数矩阵,并依此对细观有限元模型的热传导系数矩阵进行组装,形成总热传导系数矩阵:
步骤3.4:将步骤3.3得到的细观有限元模型总热传导系数矩阵赋予宏观有限元分析模型中,并对宏观有限元分析模型施加载荷,得到宏观有限元分析模型的响应。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610373839.9A CN107451308B (zh) | 2016-05-31 | 2016-05-31 | 一种复杂复合材料结构等效热传导系数多尺度计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610373839.9A CN107451308B (zh) | 2016-05-31 | 2016-05-31 | 一种复杂复合材料结构等效热传导系数多尺度计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107451308A CN107451308A (zh) | 2017-12-08 |
CN107451308B true CN107451308B (zh) | 2020-06-16 |
Family
ID=60484556
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610373839.9A Expired - Fee Related CN107451308B (zh) | 2016-05-31 | 2016-05-31 | 一种复杂复合材料结构等效热传导系数多尺度计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107451308B (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108875260B (zh) * | 2018-07-05 | 2022-11-25 | 中国人民解放军火箭军工程大学 | 一种热喷涂涂层性能参数计算方法 |
US11783100B2 (en) | 2018-09-14 | 2023-10-10 | Northwestern University | Integrated process-structure-property modeling frameworks and methods for design optimization and/or performance prediction of material systems and applications of same |
CN111027236B (zh) * | 2019-10-30 | 2021-09-14 | 武汉大学 | 岩土体的等效导热系数的细观尺度数值研究方法 |
CN112613162B (zh) * | 2020-12-02 | 2022-06-21 | 浙江大学 | 碳纳米管纤维增强复材有效热传导系数的多尺度计算方法 |
CN112560254B (zh) * | 2020-12-09 | 2022-09-06 | 西安电子科技大学 | 基于ct扫描的短纤维增强复合材料力学性能预测方法 |
CN112668113B (zh) * | 2020-12-22 | 2022-11-01 | 空气动力学国家重点实验室 | 一种复合材料多尺度防热优化方法 |
CN112949125B (zh) * | 2021-02-22 | 2024-05-24 | 南京航空航天大学 | 基于时域渐近理论的复合材料热传导双尺度分析方法 |
CN112949153B (zh) * | 2021-03-09 | 2022-10-28 | 西安交通大学 | 一种用于周期性结构复合材料高温下传热特性的快速预测方法 |
CN113761763B (zh) * | 2021-08-06 | 2023-05-30 | 上海索辰信息科技股份有限公司 | 微观和细观尺度结构rve多尺度宏观材料性质分析方法 |
CN114441590B (zh) * | 2021-12-22 | 2024-05-14 | 中国航天空气动力技术研究院 | 一种针对梯度防热材料传热及力学特性的确定方法及系统 |
CN114386286B (zh) * | 2022-01-18 | 2023-08-08 | 上海交通大学 | 基于高导热云母带制作的主绝缘导热系数计算方法及系统 |
CN114491874B (zh) * | 2022-04-07 | 2022-07-26 | 中国航发四川燃气涡轮研究院 | 一种二元喷管进口压力及载荷的快速分析方法 |
CN116822167B (zh) * | 2023-06-08 | 2024-04-16 | 西安交通大学 | 换热器热力耦合性能多尺度分析方法、系统、介质及设备 |
CN117034856B (zh) * | 2023-10-10 | 2024-01-30 | 深圳鸿芯微纳技术有限公司 | 等效传热系数确定方法、装置、电子设备及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105044146A (zh) * | 2015-07-02 | 2015-11-11 | 西安电子科技大学 | 一种两相复合材料的随机热均化分析方法 |
CN105183990A (zh) * | 2015-09-07 | 2015-12-23 | 哈尔滨工业大学 | 一种复合材料及其结构的多尺度响应分析方法 |
-
2016
- 2016-05-31 CN CN201610373839.9A patent/CN107451308B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105044146A (zh) * | 2015-07-02 | 2015-11-11 | 西安电子科技大学 | 一种两相复合材料的随机热均化分析方法 |
CN105183990A (zh) * | 2015-09-07 | 2015-12-23 | 哈尔滨工业大学 | 一种复合材料及其结构的多尺度响应分析方法 |
Non-Patent Citations (1)
Title |
---|
复合材料热传导系数均匀化计算的实现方法;张锐等;《复合材料学报》;20141231;1581-1587 * |
Also Published As
Publication number | Publication date |
---|---|
CN107451308A (zh) | 2017-12-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107451308B (zh) | 一种复杂复合材料结构等效热传导系数多尺度计算方法 | |
CN107451307B (zh) | 一种多尺度计算复杂复合材料结构等效刚度矩阵的方法 | |
CN107451309B (zh) | 一种多尺度计算复杂复合材料结构等效热膨胀系数的方法 | |
CN106066913B (zh) | 复杂复合材料结构等效材料性能多尺度计算方法 | |
WO2020211012A1 (zh) | 一种面向混杂纤维复合材料板壳结构的快速协同优化方法 | |
Koricho et al. | An experimental and finite element study of the transverse bending behaviour of CFRP composite T-joints in vehicle structures | |
Patni et al. | On the accuracy of localised 3D stress fields in tow-steered laminated composite structures | |
Waimer et al. | Integrative numerical techniques for fibre reinforced polymers-forming process and analysis of differentiated anisotropy | |
CN111950091A (zh) | 一种基于区间不确定性的复合材料层合板结构非概率可靠性拓扑优化方法 | |
CN106295024B (zh) | 一种考虑间隙及摩擦影响的复合材料螺栓连接载荷分布计算方法 | |
Kumar et al. | Investigations on the multistability of series-connected unsymmetric laminates | |
Ma et al. | Study on structural optimization of braid‐and‐lay integrated carbon fiber reinforced polymer for metro bogie | |
Yang et al. | Design and optimization of composite automotive hatchback using integrated material-structure-process-performance method | |
Dong et al. | Dimension variation prediction for composites with finite element analysis and regression modeling | |
CN105468826A (zh) | 复合材料的设计方法 | |
CN112528537A (zh) | 一种面向压缩稳定性的变刚度复合材料结构分析方法 | |
Cong et al. | Methodology for evaluating manufacturability of composite materials | |
Jia et al. | Validation and development of trace-based approach for composite laminates | |
Sun et al. | An integrated approach for rapidly and precisely predicting the spring-in of U-shaped composite parts with ply drop-offs | |
Hou et al. | Numerical Simulation and Multi-objective Optimization for Curing Process of Thermosetting Prepreg | |
Ma et al. | Parallel optimization of design and manufacturing—Carbon fiber battery pack for electric vehicles | |
Zhang et al. | Modal characteristics and manufacturing method of variable stiffness facesheets of sandwich plates made by robotic fiber placement | |
Onal et al. | Modeling of elastic, thermal, and strength/failure analysis of two-dimensional woven composites—a review | |
Sen | Estimation of elastoplastic thermal and residual stresses in a thermoplastic composite disc under uniform temperature effect | |
Tian et al. | Buckling optimization of curvilinear fiber-reinforced composite structures using a parametric level set method |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200616 Termination date: 20210531 |
|
CF01 | Termination of patent right due to non-payment of annual fee |