基于HFGMC理论的复合材料协同多尺度损伤模拟_白连乙.pdf

信息技术白连乙，等基于 HFGMC 理论的复合材料协同多尺度损伤模拟基金项目:国家自然科学基金项目(11872205)第一作者简介:白连乙(1997)，男，辽宁大连人，硕士研究生，研究方向为复合材料多尺度分析。DOI:1019344/j cnki issn16715276202301022基于 HFGMC 理论的复合材料协同多尺度损伤模拟白连乙，周储伟，马伟杰(南京航空航天大学 航空学院，江苏 南京 210016)摘要:为进行复合材料协同多尺度渐进损伤的分析，提出一种基于高保真通用胞元法和有限元相结合的方法。推导高保真通用胞元的计算格式，并通过一个典型周期性单胞的分析结果与有限元结果的对比，验证所建模型的精确性。采用所建立的方法对一个单向复合材料缺口试验件的渐进损伤过程进行协同多尺度模拟，揭示试验件宏观性能衰减与微观尺度损伤演化之间的内在联系。结果说明:该方法适用于复合材料多尺度分析。关键词:高精度通用胞元;复合材料;多尺度分析;渐进损伤中图分类号:TP3919文献标志码:A文章编号:1671-5276(2023)01-0091-04Synergistic Multiscale Damage Simulationof Composite Materials Based on HFGMC TheoryBAI Lianyi，ZHOU Chuwei，MA Weijie(College of Astronautics，Nanjing University of Aeronautice and Astronautice，Nanjing 210016，China)Abstract:A method combining HFGMC and FEM was developed to simulte the progressive damage of composite material at multiscales synergistically The computational formulars of HFGMC were deduced，and their precision was verified by the comparison ofthe reserch results between unicellulars and finite elements in a typical periodic unit cell The proposed combination method was usedto simulate the progressive damage of a unidirectional composite notched panel at both macroscale and microscale synergistically，which reveals the intrinsic connections between macroscale mechanical degradation and microscale damage evolution of thecomposite The results prove the validity of the proposed method in multiscale analysis of composite materialKeywords:high precision universal cell;composite materials;multiscale analysis;progressive damage0引言复合材料是指由增/韧强相(纤维、颗粒等)和基体相两种或两种以上不同性质的材料用物理或化学方法在宏观尺度上组合成的新材料，其综合性能优于组分材料。随着复合材料在重要工程领域中应用的不断拓展，对性能的要求也愈来愈高。因此，须更深入了解其破坏机理，以便充分发挥其性能潜力，优化材料设计。复合材料在微观层面是不均匀的，其结构方面具有明显的多尺度特征，其性能与微观结构、组分材料密切相关1。多尺度研究分析方法能够充分利用宏观的高效性及微观尺度的精确性，是复合材料性能研究的重要方法之一。目前多尺度分析方法可以分为 3 类:层级多尺度、协同多尺度和并发多尺度2。其中协同多尺度既可实现宏/微观之间的耦合分析，将复合材料结构的损伤破坏等分析建立在具体的损伤机制之上，且计算规模又不至于过大，是一种很具有吸引力的方法。高精度通用胞元法(high fidelity general method cell，HFGMC)是一种兼顾性较高的求解精度与效率的微观力学模型，在复合材料多尺度分析中应用广泛。该方法在ABOUDI J 提出的通用胞元法(GMC)3 基础上发展而来，已被广泛应用于复合材料性能的分析。如 ABOUDI J 等4 采用 HFGMC 研究了复合材料非在力和温度载荷下的非线性性能，王明路等5 基于修正的 HFGMC 对复合材料细观尺度进行了分析。美国 NASA 采用 HFGMC 对复合材料各种性能进行了大量的研究，促进了该理论的发展。本文基于 HFGMC 提出了一种复合材料协同多尺度损伤分析方法，通过对比有限元结果，验证了本文方法的精度，进一步对复合材料缺口件进行渐进损伤分析，证明了本文方法在复合材料多尺度分析中的可行性。1HFGMC 理论假设复合材料结构是由一系列周期性单胞(repeatedunit cell，UC)拼装而成，如图 1(a)所示，此复合材料的性能可通过分析 UC 获得。UC 采用 HFGMC 分析，一个 UC 划分成 NN个方形子胞，如图 1(b)所示。19信息技术白连乙，等基于 HFGMC 理论的复合材料协同多尺度损伤模拟图 1HFGMC 模型根据 UC 的非均质结构，赋予每个子胞材料属性，如或是纤维的或基体的、或是纤维/基体界面的。假设 VE中每个子胞的位移函数，如式(1)所示。u(，)i=ijxj+W(，)i(00)+y()2W(，)i(10)+y()3W(，)i(01)+123y(，)22h24()W(，)i(20)+123y(，)23l24()W(，)i(02)(1)式中:i=2，3 为平面的两个方向;ijxj表示 UC 的宏观位移分量;(，)为子胞在 UC 中的编号位置;h、l为子胞的边长;y()2、y()3是子胞局部坐标系;W(，)i(mn)是子胞位移函数中的待定参数。采用有限体积直接平均细观力学中的思想6，将子胞的界面平均量作为 HFGMC 模型的未知量，如图 2 所示。图 2子胞边界的平均位移和平均应变以 F1 面 2 方向为例，位移与应变的界面平均量可定义成如下形式:u(，)2F1=1ll/2l/2u(，)2h2，y()3()dy()3(，)22F1=1ll/2lr/2(，)22h2，y()3()dy()3(2)将子胞界面平均位移代入几何方程式(3)与本构方程式(4)中，得到子胞边界平均位移与平均应力之间的关系，并用矩阵形式表达(即子胞的“单元刚度矩阵”)，如式(5)所示。()ij=12(iu()j+ju()i)(3)(，)ij=S(，)ijkl(，)kl(4)22F122B123F123B132F232B233F233B2=K11K12K15K16K17K18K21K22K25K26K27K28K33K34K35K36K37K38K43K44K45K46K47K48K51K52K53K54K55K56K61K62K63K64K65K66K71K72K73K74K75K76K81K82K83K84K85K86u2F1u2B1u3F1u3B1u2F2u2B2u3F2u3B2+C22C23C22C232C442C442C442C44C23C33C23C33223323(5)相邻子胞界边界应该满足平均意义下的应力连续条件，此外在 UC 边界上的子胞还应满足周期性边界条件。联立位移连续条件可建立四组方程，以其中一组为例，如式(6)所示。K(，)12u(，)2B1+(K(，)11+K(+1，)22)u(+1，)2B1+K(+1，)21u(+2，)2B1+K(，)16u(，)2B2+K(，)15u(，+1)2B2+K(+1，)26u(+1，)2B2+K(+1，)25u(+1，+1)2B2+K(，)18u(，)3B2+K(，)17u(，+1)3B2+K(+1，)28u(+1，)3B2+K(+1，)27u(+1，+1)3B2=(C(+1，)22C(，)22)22+(C(+1，)23C(，)23)33(6)求解式(6)即可得到子胞边界上的平均位移。将上述方程组写为矩阵形式，如式(7)所示。K U=Coe U=K1Coe(7)式中:K 为 4NN阶方阵，包含子胞材料和几何参数信息;Coe为 4NN3 阶系数矩阵，包含子胞的材料参数;U和 分别表示子胞待求的边界平均位移和平均应变。将假设的位移表达式带入几何方程并进行平均积分，得到子胞平均应变与边界位移之间的关系，可以将子胞的平均应变以宏观应变的形式表示，如式(8)所示。(，)ij=A(，)(8)即可得到 UC 的宏观平均应力，如式(9)所示。=1HLN=1N=1hl(，)ij=1HLN=1N=1hlC(，)A(，)=C*T(9)式中 C*=1HLN=1N=1hlC(，)A(，)，即为当前 UC 所代表复合材料的宏观弹性矩阵。29信息技术白连乙，等基于 HFGMC 理论的复合材料协同多尺度损伤模拟2算例与分析21HFGMC 精度验证本文选取单向纤维增强复合材料理想化的 UC 为例，分别建立了 HFGMC 模型和 FEM 模型，两者网格剖分一致，如图 3 所示。有限元分析采用 ABAQUS 软件，组分材料参数如表 1 所示。图 3单向复合材料 UC 子胞划分示意图表 1组分材料力学性能参数名称弹性模量/MPa泊松比强度/MPa基体3 50003870纤维72 0000223 400界面3 500038200HFGMC 和 FEM 两种模型均施加平均 1%的横向拉伸应变，并且对有限元模型施加周期性边界条件(HFGMC自动满足周期性边界条件)。两种模型得到的应变、应力场分布如图 4 所示。从图中可以看出，两者的结果具有很好的一致性，对于受拉方向应变 E11 的最大值，HFGMC模型比 FEM 的误差小 23%，且发生的位置也相同，但FEM 模型中 E11 分布区域更宽一些。对于剪应变 E12，两个模型预测的分布规律也是一致的，FEM 模型(绝对值)最大值的计算比 HFGMC 模型高出 6%，误差在工程应用的可接受范围之内。两种模型应力分析结果与应变分析类似，FEM 模型正应力 S11 的最大值比 HFGMC 模型高54%，正应力 S12 的最大值则高65%。图 5 为上述两种模型对于单向复合材料对应不同纤维体积含量下横向拉伸模量、泊松比和切变模量预测的对比，可以看出在所研究的纤维体积含量范围内，两者对上述弹性常数的预测均比较接近，横向拉伸模量、切变模量和泊松比，FEM 的预测值比 HFGMC 分别高出 12%、32%、25%。22基于 HFGMC 的复合材料多尺度分析本节所采用的复合材料缺口件是 90单铺层的玻璃纤维层合板，具体尺寸为 200mm60mm。考虑到应力集中现象，在缺口处进行网格细化，并选取上节的 UC 模型描述该复合材料的微观结构。试验件级的宏观分析采用FEM 模型，在 ABAQUS 平台上计算，模型网格如图 6(a)所示。假设 FEM 模型的每个单元(积分点)中包含了无数的微观 UC，所以一个单元(积分点)的损伤分析就对应于一个 HFGMC 的 UC，见图 6(b)。基体子胞材料的破坏判断选用改进的 TsaiHill 失效准则，纤维与界面子胞材料采用最大应力准则，当子胞判定失效时，将该子胞的弹性矩阵中元素的数值折减为原来的 01。图 4HFGMC 和 FEM 模型得到应变、应力场图5不同纤维体积含量下两种模型的平均弹性常数预测值图 6单向复合材料协同多尺度模型39信息