综述：基于力学行为分析的复合材料有限元方法综述

《Mathematics》：Mechanical Behavior Analysis-Based Finite Element Method of Composites: A Review

　　

复合材料的力学行为由所考虑复合材料的工程常数值决定。如果研究人员研究均质各向同性复合材料，则需要两个工程弹性常数来表征材料；如果涉及横向各向同性复合材料，则需要五个弹性常数。对于更复杂的材料，可能需要确定更多的弹性常数以获得复合材料在实际应用中的行为。在本文中，作者介绍了用于计算纤维复合材料工程常数的主要经典方法，这些方法与有限元方法（FEM）并行使用，并强调了使用直接有限元方法（FEM）来实现这一目标的优点（和缺点）。相同纤维的排列提供了规律性，使得计算更容易，并且在某些情况下可以应用简单的方法。本研究批判性地介绍了已成为经典的结果、当前结果和不寻常的实例。所有这些发现都在使用有限元方法（FEM）的背景下进行了讨论，无论是作为主要计算方法还是作为应用经典方法的有用辅助手段。本文重点介绍了使用有限元方法（FEM）进行此目的的研究。针对所讨论的不同方法以及整体领域，强调了未来的研究方向。

**1. 引言**
复合材料广泛应用于航空航天、汽车、船舶和能源等工业领域，因其高强度-重量比、可定制力学性能和低制造成本而具有显著优势。确定纤维增强复合材料的弹性常数面临诸多挑战，包括材料复杂性、环境条件（如高温、湿度、载荷场景）、建模与测试困难、尺度效应以及跨学科方法的需求。这些常数决定材料的力学行为和特定应用中的性能选择。随着复合材料应用领域的快速发展，现代时期对生产方法、建模方法以及材料性能提升的需求日益增加。本文综述了与有限元方法（FEM）并行使用的经典计算方法，并强调了直接使用有限元方法（FEM）的优点和缺点。研究批判性地介绍了经典结果、当前结果和不寻常的实例，并讨论了这些发现与有限元方法（FEM）作为主要计算方法或经典方法辅助工具的关系。文献中的研究涵盖了从早期界限理论（如Hill和Hashin的工作）到具体复合材料模型（如芳纶纤维、Kevlar应用、天然纤维和纺织品）的发展。此外，纳米纤维（如石墨烯、碳纳米管）和混合复合材料（结合玻璃、碳、芳纶等不同纤维）的兴起，以及多尺度建模（连接宏观和微观尺度）和机器学习（ML）技术的应用也成为重要趋势。预测模型逐步集成损伤力学（如分层、基体开裂、纤维-基体界面退化）以更准确预测寿命和失效。实验测试（如加速测试和自动化实验）对于验证模型至关重要。复合材料在不同环境（温度、湿度、紫外线、热循环）下的力学性能理解是跨学科（材料科学、机械工程、计算建模等）合作的结果。本文旨在综述以聚合物为基体的复合材料力学特性的主要建模方法，其结果可谨慎推广至陶瓷或金属复合材料。通过提供有限元方法（FEM）单独或与其他经典方法结合的最新应用，本文完善了过去五十年的研究综述，提供了更全面和当前的观点。

**2. 模型与方法以确定弹性系数**
由于复合材料的材料组合、浓度和微观结构范围极广，通过实验确定所有组合的力学特性几乎不可能且成本高昂。因此，通过计算模型（即合适的力学模型）估算这些性能成为关键，模型需通过实验验证。以下介绍主要模型及其与有限元方法（FEM）的并行应用。

**2.1 微力学模型**
**2.1.1 一般考虑**
微力学分析是基于微观结构信息计算弹性常数的基本方法。它考虑纤维与基体间的力学相互作用，通过了解复合材料的微观几何（组分形状、材料、几何尺寸、取向、分布等）、应力-应变关系、刚度和断裂机制，在宏观尺度上获得材料性能（如弹性模量、强度、热导率）。评估可涉及损伤起始（开裂）的理解、微观结构优化（纤维体积分数、几何、排列等）以及不同加载条件下的行为确定。

**2.1.2 模型**
假设单向复合材料中纤维呈周期性排列（如方形排列），纤维为线性弹性且各向异性，聚合物基体为非线性粘弹性且各向同性。通过代表性单元细胞（RUC）或代表性体积元（RVE）进行分析，其中纤维截面为圆形，每个子细胞代表四分之一纤维及其基体。位移采用一阶线性展开表示，基于局部坐标系建立平均应力与平均应变的关系，从而确定整体本构定律的系数。通过界面牵引连续性和位移连续性，最终得到复合材料的宏观弹性常数。

**2.1.3 应用与未来发展**
微力学模型可用于单向复合材料在纵向和横向法向载荷下的模拟。有限元方法（FEM）广泛用于确定应力-应变场，例如结合内聚力模型研究纤维-基体界面效应（如金属-基体长纤维复合材料），或用于连续双向纤维增强复合材料（CBFRC）的宏观力学性能估计。模型可忽略泊松效应或采用简单方法（如DYNA3D中的显式单元）。未来发展方向包括：多物理场耦合（机械、热、环境）；时间相关界面行为模型；考虑微观结构变异性；多尺度方法；以及机器学习（ML）的应用，如物理信息神经网络（PINNs）将本构蠕变法则集成到机器学习（ML）模型中，或基于高保真微力学模拟的代理模型探索设计空间。

**2.1.4 有限元方法（FEM）在方法发展中的应用**
各向异性弹性在确定复合材料力学常数中至关重要，因材料性能依赖于纤维取向。例如正交各向异性材料需要多个弹性常数（3个杨氏模量E_x、E_y、E_z；3个剪切模量G_xy、G_yz、G_xz；6个泊松比）。经典分析方法难以满足所有设计和计算需求，有限元方法（FEM）因其特性和多功能性成为重要辅助工具。然而，该方法具有高数学复杂度（一般情况需21个独立弹性常数）、实验验证成本高、结果对输入数据精度敏感（需考虑纤维真实取向、孔隙率、层压缺陷等）以及难以物理分离相互作用等缺点。尽管如此，各向异性弹性分析能提供更精确的力学行为描述，有利于设计优化。

**2.2 分析方法**
**2.2.1 界限估计：模型与方法**
变分方法提供弹性常数的上下界（如体积模量、剪切模量、泊松比）。以方形排列的平行圆柱纤维增强聚合物为例，Hill和Hashin的理论给出了界限公式，如体积模量K、纵向杨氏模量E₁₁、泊松比ν的上下界，以及轴向剪切模量G_axial的界限。Voigt和Reuss界限（均匀变形和均匀应力）是更基本的界限。分析表明，当上下界差异很大时，工程实用性有限。通过使用附加形态信息（如广义Hashin-Shtrikman变分原理），可得到更精确的界限。有限元方法（FEM）可用于计算或细化界限，通过显式建模微观结构和各向异性弹性，获得更接近真实行为的界限。有限元方法（FEM）的优势包括高精度、可建模真实几何（纤维非均匀分布、界面特性等），但计算成本高且需更多输入数据。界限估计方法（Voigt、Reuss、Hashin-Shtrikman、Hill或Mori-Tanaka）常用于初步设计阶段的快速评估，但精度有限，且不适用于非线性行为。

**2.2.2 应用与未来发展**
界限方法的应用实例包括：利用Hashin界限估算复合材料物理特性；利用三维有限元模型预测TC33/环氧单向碳纤维增强聚合物（CFRP）层合板在双轴载荷下的损伤起始和演化，其中采用扩展Drucker-Prager塑性屈服模型定义基体弹塑性行为；利用改进的层合板理论分析薄QI-CFRP试样的面外模量。Mori-Tanaka方法是一种高效估计弹性常数的方法，应用于玻璃增强聚氯乙烯（PVC）的蠕变分析，并与三维数值模拟对比验证其在预测纤维长度和体积分数影响方面的适用性。该方法还用于纤维增强混凝土（FRC）的线性弹性行为研究，其刚度估计与有限元方法（FEM）结果高度一致。未来发展方向包括：将Mori-Tanaka形式与机器学习（ML）集成以更高效预测有效性能；多尺度多物理建模（如压电纤维复合材料、磁致活性聚合物）；纳入非线性时间相关行为；改进不规则、梯度或层次结构的界限技术；生物启发和智能复合材料（如自修复复合材料的均匀化方法）；以及数字孪生与实验验证（如X射线断层扫描（X-ray CT）结合有限元方法（FEM）进行实时验证）。

**2.3 均匀化粘弹性系数**
**2.3.1 均匀化方法**
均匀化理论通过建模复合材料复杂微结构（通常为重复结构）来创建具有可比宏观性能的均质材料。该方法基于周期性结构，利用平均化过程确定工程常数。应力张量需满足弹性理论方程（包括对称性和边界条件），本构关系（Hooke定律）中的系数可通过平均应变和应力场获得，或通过变分公式求解。由于精确应力-应变场难以解析确定，通常需借助数值方法（主要为有限元方法（FEM））获取均匀化系数。

**2.3.2 应用与未来发展**
均匀化方法广泛用于评估连续纤维增强单向复合材料的力学行为，例如碳纤维增强聚合物。有限元方法（FEM）用于确定应力-应变场并验证实验结果。该方法已应用于不同树脂和增强材料的多种情况，包括碳纤维增强聚合物（CFRP）的边界条件影响、长和短纤维增强粘弹性复合材料的蠕变行为、三维互锁编织复合材料的多尺度问题、短纤维和长纤维复合材料的各向异性力学行为，以及SiC纤维增强钛基复合材料等。均匀化方法和多尺度模型被证明有效，宏观行为受微观尺度组分行为、加载条件和环境因素影响。未来发展方向包括：多尺度或多物理建模；降阶模型（ROM）降低计算成本；基于机器学习（ML）的本构模型；高分辨率实验数据（扫描电子显微镜（SEM）、显微计算机断层扫描（micro-CT））直接集成到均匀化过程中，并通过数字体积相关（DVC）验证蠕变应变场；考虑微损伤演化（开裂、纤维断裂）的均匀化方法。这些进展在碳纤维增强聚合物（CFRP）、玻璃纤维增强聚合物（GFRP）、碳/碳（C/C）复合材料、金属基复合材料（MMC）和陶瓷基复合材料（CMC）的高温结构应用中有重要意义。

**3. 有限元方法（FEM）用于计算复合材料弹性常数**
**3.1 有限元方法（FEM）作为确定复合材料力学性能的有用工具**
有限元方法（FEM）主要用于确定复合材料结构的应力和应变场，是工程师和研究人员的常用工具，在大多数情况下提供优异结果。在断裂研究中，有限元方法（FEM）可用于处理裂纹和界面几何，采用多层级方法和混合单元近似校正技术。分层损伤（如复合材料在车削、钻孔、铣削过程中的诱导损伤）可通过有限元方法（FEM）和离散元方法模拟。对于纤维增强复合材料（FRC）的损伤行为，有限元方法（FEM）可用于微观和宏观层级，如玻璃纤维增强环氧树脂在库仑-摩尔准则下的损伤起始和发展。多尺度有限元方法（FEM）模型和边界积分方程方法与有限元方法（FEM）结合使用，可处理裂纹尖端的奇异行为。有限元方法（FEM）还用于预测非线性行为，例如改进的固定点迭代法确定应力场，避免大外载荷的收敛问题。增强有限元方法（A-FEM）可研究穿过单元内部的任意不连续面，且与商业有限元方法（FEM）软件兼容。

**3.2 有限元方法（FEM）作为确定力学性能的方法**
以随机取向短天然纤维增强聚合物复合材料为例，采用有限元方法（FEM）进行准静态数值模拟。准静态问题采用增量迭代算法（如牛顿-拉夫逊（Newton-Raphson）方法）处理非线性。变分公式通过时间有限差分得到增量公式，逐步求解小变形和小位移。通过共振法实验确定剪切模量G和弹性模量E，并与有限元方法（FEM）计算的共振频率进行验证和比较。

**3.2.1 随机取向短纤维增强复合材料问题的准静态公式**
考虑粘弹性复合材料，忽略惯性效应。采用Kelvin-Voigt本构关系，包括弹性和粘弹性张量。控制方程包括平衡方程、Dirichlet和Neumann边界条件以及初始条件。由于数据依赖时间，过程不可逆且依赖加载路径，需通过增量加载和迭代求解。

**3.2.2 准静态问题的变分公式**
引入函数空间（如V, H¹(Ω)），利用格林（Green）公式得到变分公式，为有限元方法（FEM）离散化提供基础。

**3.2.3 时间离散化的增量公式**
通过时间有限差分近似位移时间导数，将变分问题离散化为增量问题，逐增量计算小变形和小位移，最终叠加得到总响应。

**3.2.4 有限元近似和增量迭代算法描述**
在有限维空间内插值位移场，通过标准有限元过程得到矩阵形式的平衡方程（刚度矩阵K和载荷向量），通过牛顿-拉夫逊（Newton-Raphson）方法求解每一步的固定点问题。

**3.3 有限元方法（FEM）在振动分析中确定工程常数的应用**
以两端固定的碳纤维增强聚合物梁为例，通过有限元方法（FEM）计算固有频率，再基于经典梁理论公式反推弹性模量或剪切模量。该方法快速、成本低，适用于多种材料（包括合成纤维或混合纤维-基体系统），且与非破坏性测试（NDT）结合可提供准确实时行为建模。通过比较有限元方法（FEM）所得频率与经典解析公式，可得到均匀化材料的弹性常数。例如，扭转振动中剪切模量G可通过特征值计算（公式：G = (4 L² f_n² ρ I_p) / (n² I_t)）获得，其中f_n为特征频率，L为梁长度，ρ为密度，I_p为极惯性矩，I_t为扭转常数。通过平均多个特征值可提高估计精度。该方法简单易用，但受有限元方法（FEM）近似误差和经典梁理论假设影响。

**4. 讨论与结论**
尽管有限元方法（FEM）常用于研究纤维增强复合材料的力学行为，该领域仍有大量未实现的研究潜力。未来方向包括：多尺度建模，通过代表性体积元（RVE）连接微观行为（纤维、基体、界面）与宏观结构响应，可建立层级或并发多尺度有限元方法（FEM）模型。有限元平方（FE²）方法是一种计算均匀化方法，通过同时求解宏观和微观有限元模型，在宏观尺度获取整体结构行为，在微观尺度捕获纤维-基体相互作用，适用于时间相关问题（如蠕变），但计算成本极高。界面行为（如内聚和分层）的建模可通过渐进退化规律实现。非线性基体行为和渐进损伤（如纤维断裂）仍需进一步发展非线性粘弹性/粘塑性本构模型。实验验证（如数字图像相关（DIC）和原位显微镜）和基准数据集开发至关重要。复合材料性能的显著分散性要求耦合概率方法（蒙特卡洛、随机有限元方法（FEM）、结构可靠性）。增材制造（3D打印）引入的复杂纤维取向和孔隙可通过有限元方法（FEM）优化长期可靠性。人工智能（AI）和机器学习（ML）加速了有限元方法（FEM）的集成，提高预测效率和降低计算成本。有限元方法（FEM）在工业过程模拟（如树脂流动、残余应力）中的应用是新方向。天然纤维复合材料需考虑环境因素（湿度、生物降解）。数字孪生和加速测试可优化蠕变测试过程，通过实时数据输入预测损伤。有限元方法（FEM）向多物理场分析（热-机械、电-机械、流体-结构）的扩展为功能复合材料和智能结构提供广阔前景。随机复合材料的研究仍是新领域，需注意代表性体积元（RVE）必须为周期性平行四边形以确保有效材料常数的正确性。