Aerospace Engineering

Trans-scale computation of cryogenic composite tank based on multi-fiber RVE and clustering analysis

  • Dan HE , 1 ,
  • Chenhui FENG 1 ,
  • Xin CHANG 2a, 2b
Expand
  • 1. Key Laboratory of Liaoning Province for Composite Structural Analysis of Aerocraft and Simulation,Shenyang Aerospace University,Shenyang 110136,China
  • 2a. a. Light Alloy Research Institute,Central South University,Changsha 410083,China
  • 2b. State Key Laboratory of Precision Manufacturing for Extreme Service Performance,Central South University,Changsha 410083,China

Received date: 2024-04-15

  Online published: 2024-09-10

Abstract

Using the carbon fiber-reinforced composite tank can remarkably reduce the weight of the launch vehicle.However,the analysis method for cryogenic tanks subjected to mechanical-thermal loads remains to be studied,especially to accurately consider the microthermal stress produced between fiber and matrix in the cryogenic environment.A representative volume model containing multiple fibers was adopted,combined with the matrix and fiber failure criteria,to establish a microscopic stress field and failure prediction model.The k-means clustering method was used for dimensionality reduction calculation,and an efficient and high-fidelity trans-scale analysis method for composite tanks was proposed.The results of illustrative examples show that the proposed method can accurately predict the elastic constants and failure strength of the composite single-layer plate according to the thermal and mechanical constants of fiber and matrix.The leakage failure process of a composite tank subjected to mechanical-thermal loads was simulated,and the critical load and failure state were given.

Cite this article

Dan HE , Chenhui FENG , Xin CHANG . Trans-scale computation of cryogenic composite tank based on multi-fiber RVE and clustering analysis[J]. Journal of Shenyang Aerospace University, 2024 , 41(4) : 1 -10 . DOI: 10.3969/j.issn.2095-1248.2024.04.001

推进剂贮箱是航天运载器动力系统与结构系统的关键部件,承担着传递载荷与盛装燃料的作用,其质量约占结构总质量的60%左右1。相比于铝合金贮箱,低温铝合金内衬/无内衬复合材料贮箱能够实现减重20%~40%的目标2,其中,无内衬贮箱轻量化效果更好,但设计难度也更大3。复合材料贮箱设计中不仅要关注结构的稳定性,更需关注贮箱是否发生功能失效。由于贮箱在服役过程中同时承受机械载荷和低温载荷,复合材料中纤维与基体的热不匹配性会引起不可忽略的细观应力4,可能直接导致复合材料树脂基体失效、产生微裂纹,进而引起燃料泄漏导致复合材料贮箱发生功能失效。如X-33飞行器所采用的多瓣蜂窝夹层结构液氢贮箱,测试时复合材料层压板中的热应力促进了微裂纹的产生,最终引起了燃料渗漏,导致试验失败25。因此,分析无内衬复合材料贮箱中组分间热应力对贮箱失效的影响十分关键,使用基于经典层合理论及宏观失效准则的传统宏观分析方法无法考虑纤维与基体之间的热应力,需对无内衬复合材料贮箱开展考虑宏细观特征的跨尺度计算研究。
近年来,跨尺度分析在工程领域得到了快速的发展。Mayes等6提出了多连续体理论(multicontinuum theory),基于三角形代表体元(representative volume element,RVE)得到了材料的本构,并对纤维基体分别进行了失效判断。李星等7提出了基于应力的复合材料跨尺度失效准则,通过对正方形RVE与六边形RVE选取参考点求得宏观至细观的应力放大系数,建立了损伤判定与损伤演化的方法,跨尺度计算与试验结果吻合较好。Huang等8针对低温复合材料贮箱力学性能评价问题提出了跨尺度分析方法。在细观尺度上,用基于六边形RVE来表征复合材料细观组分,对贮箱开展渐进失效分析并与Tsai-Wu准则进行了对比。任明法等9采用了基于六边形RVE的宏细观一体化分析方法,开展了复合材料跨尺度分析。结果表明,在温度载荷下,对基体采用细观最大应力准则可以有效地预测复合材料贮箱基体开裂。常鑫10提出了考虑制造缺陷的复合材料多尺度分析方法,并建立了相应的多尺度分析框架,对复合材料机翼盒段结构与复合材料贮箱结构开展研究,获得了两种结构的失效模式与失效载荷。Tian等11建立了六边形RVE并采用局部精确均匀化理论(locally exact homogenization theory)获得了不同破坏模式碳纤维复合材料的刚度退化系数,并基于Hashin准则与最大应力准则对低温贮箱进行了跨尺度渐进损伤分析。上述方法能够获得比较准确的预测结果,但因为他们所采用的RVE中仅包含单根纤维,所以未考虑温度作用下多纤维引起的基体局部应力集中现象。González等12证明了包含30根纤维足以表征单向复合材料细观特征。杨雷等4研究了温度的纤维分布形式对单向复合材料强度的影响,结果表明纤维随机分布的最大热应力比周期性分布的情况明显偏大,热应力的存在降低了复合材料的横向拉伸/压缩强度。多纤维RVE的效果虽然好于单纤维RVE,但对于贮箱结构而言,跨尺度渐进损伤分析需在多个尺度间开展参数传递,包含多根纤维的RVE分析效率较低,限制了其应用。
为了提高跨尺度计算效率,对多纤维RVE中单元采用聚类方法可减少RVE自由度。Liu等13提出了自洽聚类分析(self-consistent clustering analysis,SCA),在离线阶段,使用k-means聚类对具有相似力学行为的材料进行聚类以实现数据压缩;在在线阶段,宏观有限元模型与压缩后的材料信息进行实时传递计算,形成了FEM-SCA的跨尺度计算方法。随后Liu等14将SCA方法应用于应变软化问题,通过用户定义材料实现每一个积分点处包含微损伤算法的SCA模型,开展了双开孔板的拉伸跨尺度分析,获得了宏观模型中不同积分点的微观应力场和损伤场。Han等15基于FEM-SCA方法,提出了编织复合材料的多尺度分析方法,获得了T型结构在循环弯曲作用下宏观尺度与中尺度的应力、应变场及材料的非线性行为。张力等16将微观界面纳入多尺度失效分析中,采用缩聚后的放大因子进行多尺度信息传递,并对2D编织SiC/SiC复合材料开展了渐进损伤失效分析。王猛17基于改进的微观力学失效理论多尺度方法预测了编织复合材料拉伸性能,并开展了含孔编织材料拉伸失效分析与编织复合材料板的冲击失效分析研究,验证了其多尺度方法的有效性。上述研究均表明,使用k-means聚类方法对RVE中单元进行缩聚,能够有效地提高跨尺度计算效率。
为此,本文建立了多纤维的RVE,采用k-means聚类方法对RVE单元进行聚类分析,降低了多纤维RVE在跨尺度分析中的计算成本。通过对IM7/977-3复合材料力学性能参数及复合材料单向板强度的预测,证明了该方法的准确性,并建立了基于多纤维RVE和聚类分析的复合材料贮箱跨尺度分析框架。

1 细观力学模型

1.1 多纤维代表体元

在跨尺度分析中,细观尺度下通常以代表单元为研究对象,因此,合理地选取代表单元十分重要18。在建立单向复合材料RVE时,RVE中纤维的分布形式主要有两种:一种是认为纤维在空间上是规则排列的,建立的代表体元有正方形RVE、六边形对角分布RVE与六边形RVE等,如图1所示。学者们基于此类RVE开展了相关研究6-8,并取得了一定的成果;另一种是采用多纤维RVE,即单个RVE中包含多根纤维。显微照片下实际的纤维分布如图2a所示,对于高体分比的复合材料而言,一般认为其纤维分布存在一定的统计规律,纤维之间的距离与角度由制造过程确定,生成与实际纤维间距离、夹角分布具有相似统计规律的RVE,这种形式的RVE受载荷时能够将相邻纤维间存在的应力集中现象体现得更好。学者们开发出了如随机序列膨胀法(random sequential expansion,RSE)等生成算法19,用于生成包含多根纤维的RVE。本文针对单向复合材料,建立多纤维的正方形RVE用于表征复合材料细观尺度上的特征。RVE中纤维的体分比为70%20,包含40根直径为 4.5   μ m的纤维,边长 L = 30.15   μ m,厚度 t = L / 20。RVE及网格模型如图2b所示。
图1 典型的纤维规则分布的RVE
图2 多纤维代表体元

1.2 周期性边界条件

为确保RVE之间满足力的平衡与位移连续条件5,需对RVE施加周期性边界条件。RVE上相对面的网格应满足一一对应的关系,每个节点是相邻RVE节点的位移函数。周期性边界条件如式(1)所示22
u i ( l c / 2 , y , z ) = u i ( - l c / 2 , y , z ) + u i x l c u i ( x , w c / 2 , z ) = u i ( x , - w c / 2 , z ) + u i y w c u i ( x , b , h c / 2 ) = u i ( x , y , - h c / 2 ) + u i z h c i = x , y , z
式中: l c w c h c分别为RVE在xyz方向的长度; u i为位移。周期性边界条件的推导可参见文献[51222]。以ABAQUS为例,施加周期性边界条件可通过Python脚本快速实现。

1.3 等效参数计算方法

单向复合材料可视为正交各向异性材料,其热力学本构如式(2)所示。
σ ¯ = C ε ¯ - Δ T α T
首先,将1.2节中的周期性边界条件施加于多纤维RVE,确保RVE上节点力的平衡与位移连续。其次,不考虑温度的影响,在线弹性阶段对RVE分别施加6组单位应变载荷,可得到RVE细观应力场。此时,对1组细观应力场中单元应力进行体平均得到等效应力,由等效应力除以施加的单位应变载荷可得到等效刚度矩阵 C 的1列。最后,对其余细观应力场进行计算得到完整的刚度矩阵。由刚度矩阵求逆得到柔度矩阵,根据工程弹性常数与柔度系数的关系23可求得材料的等效参数。体平均公式如式(3)所示6
σ i = k = 1 N σ i k V k V
式中: σ i k V k分别为第k个单元中第i个应力分量与单元体积;N为RVE中总单元数;V为所有单元总体积。若不考虑应变载荷,对RVE施加单位温度载荷,由等效刚度矩阵 C式(4)可求解获得热膨胀系数。
α = - 1 Δ T C - 1 σ

1.4 强度准则与刚度折减模型

服役于低温环境的贮箱组分材料间存在的较大热应力会影响结构的极限载荷924,可能导致材料提前发生失效。因此,在复合材料贮箱跨尺度分析中,需在细观尺度下对组分材料失效情况进行判定。本文跨尺度分析采用细观最大应力准则判断组分材料失效,细观最大应力准则如式(5)所示。
σ i ± S = 1
式中: σ i为材料主应力,当组分材料中主应力超过强度极限时材料发生失效; + S - S分别为材料的拉伸与压缩强度。基于细观判断-宏观折减方法,并采用Camanho刚度折减模型对失效单元的宏观刚度性能进行折减。同时,组分发生失效需对材料刚度性能进行折减。Camanho刚度折减模型适用于跨尺度分析25,能够对材料中不同失效形式实现相应的折减,折减方式简单且不会造成矩阵奇异26,Camanho刚度折减模型如表1所示。
表1 Camanho刚度折减模型26
失效类型 折减系数
纤维拉伸失效 E 11 d = 0.07 E 11
纤维压缩失效 E 11 d = 0.07 E 11
基体拉伸失效 E 22 d = 0.2 E 22 , G 12 d = 0.2 G 12 , G 23 d = 0.2 G 23
基体压缩失效 E 22 d = 0.4 E 22 , G 12 d = 0.4 G 12 , G 23 d = 0.4 G 23

2 复合材料跨尺度分析方法与聚类分析

2.1 跨尺度分析方法

为了得到贮箱结构在宏观响应下细观组分场的应力,本文采用Huang等8提出的分层级跨尺度分析方法获得复合材料中纤维与基体的应力场。任一点的细观应力由该点在宏观尺度上的响应计算而来,叠加公式如式(6)所示。
σ u e = H e ε M + Δ T S e
式中: ε M为宏观结构中积分点处的应变; Δ T为温度变化量; H e S e分别为6组单位应变载荷与1组单位温度载荷下RVE中积分点的应力场; σ u e为当前积分点处温度与宏观应变作用下的细观应力场。
由于在细观尺度建立的RVE包含多根纤维,且RVE模型被精细离散化处理,若将所有单元信息用于跨尺度计算会严重降低计算效率。相关研究表明13-14,在细观尺度上,使用k-means聚类方法将相似力学行为的单元划分为一类,能够大幅降低计算成本。

2.2 k-means聚类分析

由于细观尺度下RVE被精细地离散,通过对RVE施加7组单位载荷,纤维与基体中分别含有大量力学行为相似的单元。将每组单位载荷下积分点的应力信息按 σ 11 , σ 22 , σ 33 , τ 12 , τ 13 , τ 23展开,形成N行42列的数据格式,如表2所示。其中: A i j中角标i表示第i组单位载荷; j表示单元中第j个应力分量; A 11 A 16通过施加宏观应变载荷为 ε = 1,0 , 0,0 , 0,0获得; A 71~ A 76则由施加单位温度载荷获得。采用k-means聚类方法对纤维基体单元分别进行聚类,即完成式(7)的最小二乘优化过程13
S = a r g m i n J = 1 k n S J A n - A ¯ J 2
表2 应力信息格式

A

积分点

A 11 A 66 A 71 A 76
Ele1,int1 281070 3476.0 0.1549 - 9.3 E - 5
Ele1,int2 281070 3476.0 0.1549 - 9.3 E - 5
Ele2,int1 281080 3199.5 0.1548 8.68 E - 4
EleN,int1 281090 3782.5 0.1547 - 1.56 E - 4
本文将RVE中包含的纤维基体单元进行聚类,其中每一类中的单元在空间上可以是不连续的。对每一类中的应力信息进行体平均,将获得的等效应力信息代替原RVE中的应力信息用于跨尺度细观应力计算。使用跨尺度分析方法,本文认为由宏观响应与等效应力信息计算的组分细观应力满足失效条件时,即认为该类中所有单元发生失效,从而导致RVE发生失效。根据1.4节强度准则与刚度折减模型,对发生失效的单元采用细观失效判断、宏观性能折减的方式实现结构的跨尺度失效分析。结合2.1节,最终形成基于聚类分析的复合材料贮箱跨尺度分析方法,如图3所示。
图3 基于聚类分析的复合材料贮箱跨尺度分析方法

3 复合材料贮箱跨尺度计算

3.1 等效参数预测

纤维与基体的力学性能如表34所示。RVE划分为约116 000个单元,由C3D8R六面体单元及少量C3D6楔形单元组成。
表3 IM7纤维力学性能620
材料参数 数值
E f   1 / G P a 277.1
E f   2 = E f   3 / G P a 12.9
μ 12 = μ 13 0.28
μ 23 0.45
G f   12 = G f 13 / G P a 10.72
G f   23 / G P a 4.45
α f   1 / ( 10 - 6 · K - 1 ) -0.9
α f   2 = α f   3 / ( 10 - 6 · K - 1 ) 7.2
X f t / M P a 4 082
X f c / M P a 2 266
表4 977-3基体力学性能620
材料参数 数值
E m / G P a 3.45 × ( 2 - e 0.718 ( T - 293 ) / 293 )
α m / ( 10 - 6 · K - 1 ) - 0.159 × 10 - 3 T 2 + 0.233 T + 1.673
μ m 0.41
X m t / M P a 81.3
X f t / M P a 350.2
根据IM7纤维与977-3基体的力学性能,结合式(2)、(3)开展IM7/977-3复合材料单层板的等效参数预测,所得结果如表5所示,最大误差为0.9%,单层板等效参数预测值与试验结果取得很好的一致性,验证了所建立RVE模型的合理性。
表5 单层板等效参数预测值与试验结果20对比
材料参数 试验 预测结果 误差/%
E 11 / G P a 193.39 194.6 0.6
E 22 = E 33 / G P a 8.85 8.81 0.4
G 12 = G 13 / G P a 4.67 4.658 0.2
μ 12 = μ 13 0.32 0.317 0.9

3.2 单向板跨尺度强度预测

从纤维与基体的力学参数出发,基于2.2节中所建立的包含多纤维和聚类分析的复合材料贮箱跨尺度分析方法,开展单向板跨尺度强度预测。建立的单向板有限元模型中长、宽、总厚度及单层厚度分别为160、40、2、0.125 mm,铺层为 0 16。首先,对建立的RVE施加7组单位载荷,根据2.2节与式(3)对RVE中纤维与基体分别开展聚类与体平均,所得的等效应力信息被用于细观失效的计算。随后,结合式(6)所示的跨尺度分析方法实现计算组分应力,并依据强度准则判断组分失效,对发生失效的单元进行相应的刚度折减,最终实现单向板跨尺度强度预测。分别对模型开展纵向拉伸、纵向压缩、横向拉伸与横向压缩等载荷下单向板强度性能计算。单向板有限元模型如图4所示,模型共包含约150 000个实体单元。
图4 单向板有限元模型
表6汇总对比了纵向拉伸时直接求解与聚类降阶求解的未知量数目及计算时间。跨尺度计算中,细观应力状态基于应力、应变进行求解,同时将温度的影响纳入考虑,每个类包含7个未知量。由表6可以看出聚类-32-32降阶求解的未知量数目是直接求解的 1 / 1   810,计算耗时约为 1 / 107,采用更低的聚类数量时计算效率甚至更高。对精细离散化RVE所有单元使用聚类降阶的方法,未知量数目得到了显著的降低,跨尺度计算所需的时间也得到了降低。
表6 直接求解与聚类降阶求解对比
不同求解方案 单元数/聚类数

未知

量数

计算

时间/s

直接求解 115 866 811 062 76 920
聚类-16-16 32 224 704
聚类-24-24 48 336 714
聚类-32-32 64 448 719
图5为跨尺度方法仿真得到的复合材料单向板应力应变曲线与试验数据的对比,表7将强度的试验值与预测值进行了对比。Godines等20对IM7/977-3复合材料开展了大量试验,由于试验中两个应变片所展现的非线性与不同的响应20,使用跨尺度分析方法预测的纵向压缩应力、应变曲线略有不同。基于聚类分析的单向板跨尺度强度预测值与试验值误差均在10%以内,仿真结果与试验数据较为吻合,从而验证了基于聚类分析的复合材料跨尺度计算方法的可行性。
图5 单向板跨尺度分析与试验应力应变曲线对比
表7 单向板强度的试验值与预测值20对比
强度 试验值/MPa 预测值/MPa 误差/%
X T 2 877.79 2 873.96 -0.13
X C 1 680 1 596.5 -4.9
Y T 53.6 58.48 9.1
Y C 227.5 205.76 -9.5

3.3 复合材料贮箱跨尺度分析

贮箱作为存放燃料的部件,服役于极低的温度环境下并承受一定的内压载荷,首先要确保结构的盛装性能24。组分材料中纤维强度远高于基体强度,导致基体比纤维早一步发生失效,由此形成微裂纹,可能引起燃料泄露。因此,应重点关注组分材料中基体失效情况。基于3.2节中的结果,本节采用前文建立的基于聚类的跨尺度方法与折减模型对低温复合材料贮箱开展强度预测。选取纤维基体聚类数均为16,采用基体穿孔24作为贮箱最终失效模式,即认为某位置所有单元均发生基体失效则该结构发生失效。假设纤维性能不随温度发生变化,基体性能与温度有关,如表4所示。根据1.3节内容对-253 ℃低温下复合材料参数进行预测,复合材料在低温下的性能与常温性能有所不同,如表8所示。贮箱几何参数如表9所示,铺层为 ± 90 / ± 5.74 / 90 / ± 5.74 / 90 / ± 5.74 S。建立的1/8贮箱有限元模型如图6所示,共包含大约220 000个实体单元。贮箱有限元模型边界上施加对称边界条件,定义初始温度为20 ℃,降至工作温度-253 ℃后,逐渐施加内压增量0.015 MPa,开展该贮箱跨尺度渐进失效分析。
表8 预测IM7/977-3复合材料在低温下热力学性能
材料参数 数值
E 11 / M P a 195.11
E 22 = E 33 / M P a 10.21
G 12 = G 13 / M P a 5.682
G 23 / M P a 3.33
μ 12 = μ 13 0.318
μ 23 0.537
α 11/(10-6 K-1 -0.841
α 22 = α 33/(10-6 K-1 -7.78
表9 复合材料贮箱模型参数 (mm)
模型参数 数值
直筒段长度L 270
直筒段外径R 202.5
封头区域长半轴a 200
封头区域短半轴b 150
极孔半径r 20
材料单层厚度t 0.125
图6 1/8贮箱有限元模型
复合材料贮箱的跨尺度渐进失效过程如图7所示,此处仅展示贮箱第1、3、15层分别在载荷0.86、0.90、0.95 MPa时的损伤情况。标记为深灰色的单元表示单层中发生的基体失效区域,标记为黑色的单元表示该区域基体未发生失效。结果显示,载荷为0.86 MPa时基体失效首先出现于第3层中直筒段与封头段交汇处,随着载荷增加失效区域向贮箱直筒段和厚度方向扩展,最后在载荷0.95 MPa时,复合材料贮箱中发生了基体贯穿失效,形成了渗漏通道导致此时贮箱发生功能失效。
图7 复合材料贮箱跨尺度渐进失效过程

4 结论

本文基于多纤维RVE,采用k-means聚类方法对RVE中的单元进行聚类开展降维计算,建立了基于聚类分析的复合材料贮箱跨尺度分析方法,给出了完整的实施流程,并与文献中的试验结果进行了对比验证,得到主要结论如下:
(1)本文所提方法能够基于纤维和基体的性能准确地预测复合材料单向板的弹性常数及破坏强度,误差在10%以内。
(2)本文研究基于一个受温度与内压载荷共同作用的贮箱,结果表明,在-253 ℃、0.86 MPa内压作用下,该贮箱中基体失效首先发生在封头与直筒段连接处。随着内压增加失效区域也逐渐扩展,最终在0.95 MPa时发生基体穿孔,形成渗漏通道,此时贮箱发生功能失效。
(3)相比于传统的多尺度分析,本文所提方法未知量个数与计算时间分别降为原来的 1 / 1   810 1 / 100,具有较好的计算效率,能够快速实现复合材料贮箱的跨尺度失效预测。
(4)由于结构级试验的难度,目前仅从单向板层面开展了强度验证,但未对复合材料贮箱开展试验验证,可在未来开展此工作,针对结构级试验的失效位置和失效载荷开展进一步的验证工作。
1
黄诚,刘德博,吴会强,等.我国航天运载器复合材料贮箱应用展望[J].沈阳航空航天大学学报201633(2):27-35.

2
湛利华,关成龙,黄诚,等.航天低温复合材料贮箱国内外研究现状分析[J].航空制造技术201962(16):79-87.

3
黄诚,雷勇军.大型运载火箭低温复合材料贮箱设计研究进展[J].宇航材料工艺201545(2):1-7.

4
杨雷,刘新,高东岳,等.考虑纤维随机分布的复合材料热残余应力分析及其对横向力学性能的影响[J].复合材料学报201633(3):525-534.

5
Yang L Li Z Xu H,et al.Prediction on residual stresses of carbon/epoxy composite at cryogenic temperature[J].Polymer Composites201940(9):3412-3420.

6
Mayes J S Hansen A C.Composite laminate fai-lure analysis using multicontinuum theory[J].Composites Science and Technology200464(3/4):379-394.

7
李星,关志东,刘璐,等.复合材料跨尺度失效准则及其损伤演化[J].复合材料学报201330(2):152-158.

8
Huang C Ren M Li T,et al.Trans‑scale modeling framework for failure analysis of cryogenic composite tanks[J].Composites Part B:Engineering201685:41-49.

9
任明法,张笑闻,常鑫,等.低温推进剂复合材料贮箱基体开裂预测方法[J].宇航学报201738(9):1005-1012.

10
常鑫.含孔隙和分层缺陷大型复合材料构件强度多尺度分析[D].大连:大连理工大学,2020.

11
Tian L Zhao H Wang G,et al.Progressive da-mage analysis for multiscale model of linerless composite cryotank and integrated design[J].AIAA journal.202260(3):1873-1882.

12
González C Llorca J.Mechanical behavior of unidirectional fiber‑reinforced polymers under transverse compression:microscopic mechanisms and modeling[J].Composites Science and Technology200767(13):2795-2806.

13
Liu Z Bessa M A Liu W K.Self‑consistent cluste-ring analysis:An efficient multi‑scale scheme for inelastic heterogeneous materials[J].Computer Methods in Applied Mechanics and Engineering2016306:319-341.

14
Liu Z Fleming M Liu W K.Microstructural material database for self‑consistent clustering analysis of elastoplastic strain softening materials[J].Computer Methods in Applied Mechanics and Engineering2018330:547-577.

15
Han X Gao J Fleming M,et al.Efficient multiscale modeling for woven composites based on self‑consistent clustering analysis[J].Computer Methods in Applied Mechanics and Engineering2020364:112929.

16
张力,王猛,陈强,等.考虑微观界面的2D编织SiC/SiC复合材料宏-细-微多尺度渐进损伤失效分析[J].工程力学202239(3):233-248.

17
王猛.碳纤维增强复合材料宏-细-微观损伤失效研究[D].南京:东南大学,2022.

18
陈玉丽,马勇,潘飞,等.多尺度复合材料力学研究进展[J].固体力学学报201839(1):1-68.

19
Yang L Yan Y Ran Z,et al.A new method for generating random fibre distributions for fibre reinforced composites[J].Composites Science and Technology201376:14-20.

20
Godines C Dormohammadi S Abdi F,et al.Da-mage tolerant composite design principles for aircraft components under static service loading u-sing multi‑scale progressive failure analysis[J].Journal of Composite Materials201751(10):1393-1419.

21
Naya F González C Lopes C S,et al.Computational micromechanics of the transverse and shear behavior of unidirectional fiber reinforced polymers including environmental effects[J].Compo-sites Part A:Applied Science and Manufacturing201792:146-157.

22
Huang H Talreja R.Effects of void geometry on elastic properties of unidirectional fiber reinforced composites[J].Composites Science and Techno-logy200565(13):1964-1981.

23
沈观林,胡更开.复合材料力学[M].北京:清华大学出版社,2006.

24
黄诚.航天运载器低温复合材料贮箱结构设计方法研究[D].长沙:国防科技大学,2017.

25
张笑闻.考虑界面相的低温复合材料贮箱结构多尺度分析方法研究[D].大连:大连理工大学,2020.

26
姜海林.基于等效位移法的低速冲击载荷下复合材料层合板损伤力学模型[D].大连:大连理工大学,2020.

Outlines

/