短纤维增强热塑性塑料（SFRTs）工艺诱导各向异性热机械性能的有限元分析（FEA）框架

《Polymers》：FEA Framework for Short-Fiber-Reinforced Thermoplastics with Process-Induced Anisotropic Thermomechanical Properties

　　

本研究提出了一种用于注塑成型短纤维增强热塑性塑料（short-fiber-reinforced thermoplastics, SFRTs）的有限元分析（finite element analysis, FEA）框架，该框架考虑了工艺诱导各向异性（process-induced anisotropy）、温度及湿度因素。研究人员将完全取向材料模型（fully oriented material model, FOMM）与来自注塑成型仿真的纤维取向张量（fiber orientation tensors, FOTs）相结合，以确定PA66-GF30的局部正交各向异性热机械性能。该模型采用细观力学公式（micromechanical formulations）描述弹性、强度及热膨胀。研究人员通过逆标定（inverse calibration）利用标准化复合材料数据同时识别了基体杨氏模量（Young’s modulus）、泊松比（Poisson’s ratio）、拉伸强度（tensile strength）及线热膨胀系数（coefficient of linear thermal expansion）。标定前，预测精度在很大程度上取决于性能属性、温度及湿度条件，其中在120 °C及湿态（conditioned-state）拉伸强度条件下偏差最大。标定后，直接标定的纵向杨氏模量与拉伸强度与其目标值偏差均在±0.4%以内。非直接标定的横向模量与拉伸强度分别与其辅助估计值（auxiliary estimates）偏差在±5.4%及±1.9%以内。对于干态（dry-as-molded）材料，方向性热膨胀系数与其标定目标偏差在±5.2%以内。对于湿态材料，采用干态热膨胀值作为目标值，且偏差保持在±5.9%以内。这些结果证明了成功的标定与内部一致性，而独立验证需要未在识别过程中使用的方向性测量。该框架实现了结构有限元分析中环境依赖性各向异性特性的实用映射。

热塑性聚合物因其优异的力学与热性能、高制造效率及替代传统材料的潜力而日益重要。添加短纤维（如玻璃纤维、碳纤维）增强可显著提升其力学性能。短纤维增强热塑性塑料（short-fiber-reinforced thermoplastics, SFRTs，指以短纤维增强的热塑性复合材料）通常采用注塑成型，导致热机械性能呈现显著的局部各向异性。充模过程中，基体聚合物链在高剪切及快速冷却区域沿流动场取向，增强纤维与流动基体动态交互，形成依赖于流动历史的工艺相关取向状态，导致壁厚方向上纤维取向与基体分子取向的不均匀分布。这些制造诱导特征可能引起残余应力与模塑后变形，并因热塑性材料的粘弹性而进一步松弛。若将注塑件简单视为各向同性连续体，强度评估可能产生显著误差。目前工业常用的有限元分析（finite element analysis, FEA，一种数值计算方法）工作流往往忽略制造诱导的局部材料特征，导致初始状态与强度响应预测精度不足。

此外，热塑性塑料的力学性能高度依赖温度与变形速率，与聚合物链的热激活运动密切相关，尤其在玻璃化转变温度（glass transition temperature, T_g，聚合物从玻璃态向高弹态转变的特征温度）附近。低于T_g时材料刚度与强度较高；接近并超过T_g时，链段运动加剧，刚度与强度显著降低，粘弹性效应（如蠕变与应力松弛）愈发重要。聚酰胺类材料吸收的水分起到增塑剂作用，降低T_g，使材料依据其处于干态（dry-as-molded，成型后未吸湿状态）、湿态（conditioned，标准温湿度平衡吸湿后状态）或饱和态（saturated，最大吸湿平衡状态）而呈现显著不同的力学响应。吸湿通常导致刚度与拉伸强度降低，韧性增加，因此分析至少需考虑干态与湿态两种条件。

除温度与湿度外，SFRTs的数值建模还需考虑熔接线（weld line，注塑件中两股熔体汇合形成的界面）等局部缺陷，其常导致局部刚度与强度下降。同时，注塑成型后残余应力来源包括不均匀冷却、保压压力、近凝固壁面处的流动诱导效应以及基体与纤维间的相互作用，可导致尺寸偏差及表面应力开裂。综上所述，面向各向同性结构材料的传统工程方法不能直接移植到注塑成型SFRTs。现有强度评估程序常依赖简化的各向同性材料描述，忽略纤维取向诱导的局部各向异性，且常遗漏模塑后效应与缺陷，可能导致仿真与实验结果出现显著偏差，进而影响部件设计决策及金属替代的可行性。

针对上述问题，研究人员D.B.与M.B.提出了一种用于注塑成型短纤维增强热塑性塑料的集成有限元分析框架，该框架能够考虑工艺诱导各向异性、温度及湿度影响。研究通过结合细观力学公式、完全取向材料模型（fully oriented material model, FOMM，纤维完全单向对齐的理想化材料模型）、来自注塑成型仿真的纤维取向张量（fiber orientation tensors, FOTs，统计描述纤维取向状态的张量）以及取向平均方法，生成空间变化的正交各向异性热机械性能，用于结构FEA。对于完整温度与湿度相关基体性能数据不可得的情况，研究人员引入了逆标定（inverse calibration，通过匹配试验数据反识别材料参数的方法）路线。研究表明，所开发框架在完整各向异性材料数据集不可用时，为将工艺诱导的局部各向异性和环境因素纳入结构FEA提供了实用手段。

为开展上述研究，研究人员采用了以下关键技术方法：基于Halpin–Tsai模型的正交各向异性弹性性能计算、基于Van Hattum–Bernardo模型的正交各向异性强度预测、基于Rosen–Hashin模型的正交各向异性热膨胀计算等细观力学公式；利用二阶纤维取向张量（FOTs），采用Orthotropic Smooth Closure（OSC，正交各向异性平滑闭合近似）由二阶张量重构四阶张量，实现取向平均；以CAMPUS Plastics数据库中PA66-GF30的标准化复合材料数据为参考，通过逆标定同时识别基体杨氏模量（Young’s modulus）、泊松比（Poisson’s ratio）、拉伸强度（tensile strength）及线热膨胀系数（coefficient of linear thermal expansion）；使用ANSYS Parametric Design Language（APDL）编写程序，实现FOT数据从注塑成型仿真到结构FEA模型的映射与局部各向异性性能分配。

**FOMM制备与非标定性能预测**

研究人员基于初始组分性能，采用Halpin–Tsai、Van Hattum–Bernardo及Rosen–Hashin细观力学模型，计算了PA66-GF30完全取向材料模型（FOMM）的假设性正交各向异性热机械性能，并验证了正交各向异性刚度与柔度矩阵的热力学许可性。

**逆标定识别与精度评估**

研究人员建立了逆标定问题，针对每个温度与湿度条件，在允许参数域（Ω_T,h）内同步优化基体杨氏模量（E_m）、泊松比（ν_m）、拉伸强度（UTS_m）及线热膨胀系数（α_m）。通过FOT-based取向平均将FOMM性能转换为试样级有效性能，以最小化目标函数J_T,h。标定前，非标定FOMM经取向平均后的预测精度强烈依赖于性能类型与温湿条件：纵向杨氏模量在?20 °C与20 °C时与参考数据偏差在?4.3%至?3.1%以内，但在120 °C时干态与湿态分别欠预测约22.7%与23.6%；横向模量偏差范围达?31.0%至0.5%，最大偏差发生在120 °C；湿态拉伸强度过预测最大达51.4%，干态热膨胀系数在120 °C时纵向偏差达84.6%，湿态纵向热膨胀偏差在120 °C时达192.8%。经逆标定后，直接标定的纵向杨氏模量和拉伸强度在所有温度与湿度条件下与目标偏差均在±0.4%以内；非直接目标的横向模量与拉伸强度分别与辅助估计值偏差在±5.4%与±1.9%以内；干态方向性热膨胀系数与目标偏差在±5.2%以内；湿态以干态热膨胀值为替代目标，偏差在±5.9%以内，体现出内部一致性。

**局部FOT映射与试样制备路线的影响**

研究人员通过注塑成型仿真对比了直接注塑成型单侧边浇口试样与注塑成型板材机加工试样的壁厚方向FOT分布。直接注塑试样纵向分量a₁₁在全厚度范围内占优（约0.85–0.92），呈对称分布；板材外层a₁₁约0.83–0.84，而芯层a₂₂升高至约0.76，a₃₃仅约0.04，形成显著分层取向结构。体积平均后，直接注塑试样a₁₁=0.866、a₂₂=0.069、a₃₃=0.065；板材试样a₁₁=0.707、a₂₂=0.267、a₃₃=0.026。FEA映射结果显示，所开发软件能够有效保留从注塑成型仿真到结构FEA模型的工艺诱导FOT分布特征。

**正交各向异性刚度对加载角度的依赖性**

结构FEA演示表明，PA66-GF30的有效杨氏模量随加载角度β从0°（沿名义流动方向）增至90°（横向加载）而连续下降。干态下降幅约35%至85%，湿态下约33%至87%，且随暴露温度升高，角度相关刚度衰减更显著。角刚度敏感系数（angular stiffness-sensitivity measure, S_E(β)，量化刚度随加载角度变化率的指标）在β约30°至33°处达到峰值，最大范围约0.8%至1.9%/度。研究指出，仅最大化纤维取向并非充分的工艺优化目标，局部取向状态应与预期结构载荷路径相匹配。

**Closure近似敏感性评估**

研究人员采用Orthotropic Fitted Closure（OFC，正交各向异性拟合闭合近似）替代OSC进行补充评估，但未重新标定。结果表明，在干态与湿态下，拉伸强度与热膨胀系数对closure选择相对不敏感，最大差异分别为5.5%与3.9%。弹性模量对closure选择较为敏感，尤其对于取向集中度较低的板材机加工试样在120 °C时，横向模量E₂₂的对称绝对百分比差异在干态下达13.0%，湿态下达13.7%。尽管如此，替换closure未改变研究的主要定性结论，但在取向状态较弱且高温条件下，定量刚度预测应考虑closure相关不确定性。

**讨论与结论**

讨论部分强调，供应商提供的性能应视为特定工艺相关取向状态的有效响应，而非制造无关的材料特性。非标定FOMM的精度强烈依赖于性能类型与环境条件，当基体主导行为显著时，基于初始组分性能的非标定路线可靠性降低。逆标定显著改善了与选定复合材料级目标的一致性，但所识别的基体参数应视为整个细观力学与FOT-based取向平均框架的有效模型参数，而非未增强PA66的本征属性，且不应在未经证实的情况下迁移至其他细观力学模型。结构FEA结果进一步表明，刚度高度依赖于加载方向、温度及湿度条件。本框架实现了成功的标定与内部一致性验证，但独立验证完整三维各向异性响应仍需依赖识别过程中未使用的方向性测量。主要不确定性来源包括辅助横向参考值、湿态热膨胀替代目标、热膨胀试样的未知制造历史以及所采用的细观力学与代表性FOT假设。