碳纤维复合材料雷击损伤的电-热-力耦合多尺度模拟研究综述

碳纤维复合材料雷击损伤的电-热-力耦合多尺度模拟研究综述一、碳纤维复合材料雷击损伤机理与多尺度特性碳纤维复合材料（CFRP）由高强度碳纤维作为增强相、树脂基体作为粘结相复合而成，其各向异性的电学、热学与力学性能使其雷击损伤演化过程呈现显著的多尺度耦合特征。从微观尺度看，单根碳纤维直径约为5-7μm，电导率沿纤维轴向可达10^4-10^5S/m，而垂直于纤维方向的电导率仅为10^-2-10^0S/m，树脂基体的电导率更是低至10^-12S/m，这种电导率的巨大差异导致雷击电流注入时电流优先沿碳纤维轴向传导，在纤维-基体界面处容易形成电流集中。热性能方面，碳纤维的轴向热导率可达800W/(m·K)，而树脂基体的热导率仅为0.2W/(m·K)，且树脂的玻璃化转变温度通常在120-200℃之间，分解温度约为300-400℃，当雷击产生的焦耳热无法及时导出时，局部温度会迅速升高引发基体热解、纤维氧化甚至升华。力学性能上，碳纤维轴向拉伸强度可达3-7GPa，但横向强度仅为轴向的1/10-1/20，热应力与基体分解产生的气体压力共同作用下，会引发界面脱粘、层间分层、纤维断裂等多种损伤形式。介观尺度下，复合材料铺层结构是雷击损伤演化的核心载体，单层厚度通常为0.125-0.25mm，不同铺层方向的纤维排列方向差异导致层间电导率、热导率突变。当雷击电流从注入点向内部扩散时，在相邻铺层的界面处会产生电流折射现象，电流密度在垂直于纤维方向的分量显著增大，进一步加剧局部焦耳热的产生。同时，层间树脂基体的热解会产生大量CO、CO2和小分子烃类气体，内部气压可达1-10MPa，当气压超过层间断裂韧性时就会引发分层损伤，分层面积通常可达雷击附着点面积的10-100倍。层内的纤维束编织结构（如平纹、斜纹、缎纹）还会形成电流传导的“优先路径”，使得损伤区域呈现与编织纹路一致的分布特征。宏观尺度下，复合材料结构的整体构型、厚度、铺层顺序以及表面防护体系直接决定雷击损伤的最终形态。对于飞机机翼、机身等大型复合材料结构，雷击电流通常需要从注入点传导到接闪点，电流传导路径上的结构不连续处（如铆钉孔、结构转角、装配间隙）都会形成电流集中，导致局部损伤加剧。采用金属网、导电涂层等表面防护措施时，防护层的电导率、厚度以及与复合材料基体的结合强度会影响电流的分流效果，当防护层无法承受大电流冲击时会发生汽化、剥落，失去对下层复合材料的保护作用。二、电-热-力耦合多尺度模拟的理论框架当前主流的多尺度模拟框架通常采用“微观-介观-宏观”三级嵌套结构，不同尺度之间通过等效参数传递实现耦合。微观尺度模拟的核心是建立单根碳纤维与基体的RepresentativeVolumeElement（RVE，代表体积单元）模型，考虑纤维随机分布、界面相性能等细观特征，通过施加电场、温度场和应力场边界条件，计算微观尺度下的电流密度分布、温度场演化和应力应变响应。该尺度的控制方程包括电场的麦克斯韦方程、热传导的傅里叶方程以及力学平衡方程，三者通过焦耳热作为电-热耦合项，热膨胀系数、温度相关的材料性能参数作为热-力耦合项，电流导致的电致应力作为电-力耦合项，实现三场的全耦合求解。介观尺度模拟通常以单层复合材料或层合板的单元作为研究对象，将微观尺度计算得到的等效电导率、热导率、力学性能参数作为介观模型的输入，同时引入层间界面单元来模拟分层损伤的萌生与扩展。介观尺度的损伤演化模型通常考虑三种损伤模式：一是基体热解损伤，通过引入温度相关的基体性能退化模型，当温度超过热解温度时，基体的电导率、热导率和力学性能按一定的损伤变量线性或非线性下降；二是纤维损伤，当纤维温度超过升华温度（约3000℃）或者轴向拉应力超过拉伸强度时，纤维发生断裂，其承载能力和导电能力完全丧失；三是界面损伤，采用内聚力单元（CohesiveZoneModel）描述纤维-基体界面和层间界面的损伤，当界面处的切应力和正应力的组合达到断裂准则时，界面发生脱粘或分层，同时损伤后的界面电导率、热导率会出现数量级的下降，进一步影响后续的电流和热传导。宏观尺度模拟则将介观尺度得到的层合板等效性能参数应用于整体结构的雷击分析，重点关注电流在整个结构中的传导路径、温度场的分布范围以及整体结构的力学响应变化。对于大型复杂结构，通常采用多尺度子模型法，先通过宏观模拟得到雷击区域的边界条件，再将高应力、高电流密度区域作为子模型进行介观甚至微观尺度的精细化模拟，在保证计算精度的同时降低计算成本。时间尺度上，雷击电流的作用时间通常为几微秒到几百微秒，而损伤演化过程则会延续到毫秒甚至秒级，因此模拟过程需要采用多时间步长技术，电流脉冲作用阶段采用极小的时间步长（纳秒到微秒级）捕捉瞬时的电流和温度变化，电流消失后采用较大的时间步长（微秒到毫秒级）模拟热扩散和应力释放过程。三、多尺度模拟的关键技术进展在等效参数homogenization（均匀化）方法方面，传统的Mori-Tanaka方法、自洽方法虽然计算效率高，但难以考虑复合材料的随机细观结构和非线性损伤演化。近年来，基于计算homogenization的数值均匀化方法得到广泛应用，该方法直接对微观RVE模型进行多场耦合求解，通过体积平均得到等效的电、热、力性能参数，能够准确反映损伤演化过程中材料性能的非线性变化。对于编织复合材料，还发展了基于拓扑结构重构的均匀化方法，通过micro-CT扫描获取真实的纤维束编织结构，建立与实际结构一致的RVE模型，使得等效参数的误差控制在5%以内。多尺度耦合的界面处理技术是模拟准确性的关键。传统的节点耦合方法在不同尺度网格不匹配时容易出现应力集中和计算误差，近年来发展的网格重叠法（Arlequin方法）通过在不同尺度模型的重叠区域引入能量加权函数，实现宏观和介观模型的力平衡和位移协调，避免了边界条件施加带来的误差。此外，基于域分解的并行计算技术使得大规模多尺度模拟成为可能，通过将整个结构划分为多个子域，每个子域采用不同尺度的模型进行并行求解，子域之间通过边界条件进行数据交换，计算效率相比传统串行方法提升了1-2个数量级，能够实现平方米级复合材料结构的雷击多尺度模拟。多场耦合的损伤本构模型不断完善。早期的模拟通常忽略损伤对电、热性能的影响，或者只考虑单向的性能退化，近年来发展的fullycoupled（全耦合）损伤本构模型将损伤变量同时引入电导率、热导率和力学本构关系中，能够模拟“电流集中-焦耳热升高-材料损伤-电流和热传导路径改变-进一步电流集中”的正反馈演化过程。对于基体热解过程，还引入了化学动力学模型，考虑热解反应的速率、气体生成量和孔隙率变化，将热解产生的气体压力作为额外的力学载荷引入应力平衡方程，能够更准确地预测分层损伤的萌生和扩展。针对雷击过程中的等离子体-材料相互作用，部分研究还将外部等离子体场与材料内部的电-热-力场进行耦合，考虑电弧附着点的移动、等离子体对材料的烧蚀作用，进一步提高了模拟的真实性。四、典型应用场景与模拟结果验证在航空复合材料结构的雷击防护设计中，多尺度模拟技术已经得到初步应用。针对飞机机翼前缘的CFRP结构，通过多尺度模拟可以预测不同雷击电流分量（A分量、B分量、C分量、D分量）作用下的损伤深度、分层面积和剩余强度，为防护层的厚度、铺层顺序优化提供依据。研究表明，当采用0.1mm厚的铝网作为表面防护层时，能够将复合材料表面的电流密度降低80%以上，基体热解深度减少60%，但铝网与复合材料界面处的热应力容易导致界面脱粘，通过在铝网与复合材料之间引入导电胶层可以降低界面热应力，减少脱粘损伤的发生。针对复合材料风力发电机叶片的雷击防护，多尺度模拟能够模拟不同雷击位置、不同电流幅值下的损伤演化过程。叶片的后缘区域由于厚度较薄、铺层角度复杂，雷击损伤通常更为严重，模拟结果显示，当100kA峰值电流作用于叶片后缘时，内部最高温度可达2500℃以上，分层面积可达0.5m^2，通过在叶片前缘和后缘铺设连续的导电接闪系统，能够将90%以上的雷击电流引导到接闪器，大幅降低复合材料主体结构的损伤风险。模拟结果的验证通常采用实验室雷击试验结合无损检测的方法。通过脉冲电流发生器模拟标准雷击电流波形，对复合材料试样进行冲击试验，采用超声波C扫描、micro-CT扫描、红外热成像等技术检测损伤的深度、面积和形态，与模拟结果进行对比。大量验证试验表明，当前先进的电-热-力耦合多尺度模拟对损伤面积的预测误差可以控制在15%以内，对损伤深度的预测误差可以控制在10%以内，能够满足工程设计的需求。部分研究还通过实时的高速摄影、光纤传感技术监测雷击过程中的温度和应变变化，进一步验证了多尺度模拟对瞬态演化过程预测的准确性。五、现存挑战与未来发展方向尽管碳纤维复合材料雷击损伤的多尺度模拟已经取得显著进展，但仍存在诸多亟待解决的挑战。首先，多尺度、多场耦合的计算复杂度极高，当前的模拟通常需要数小时甚至数天的计算时间，难以应用于大型复杂结构的全尺寸分析和多方案优化，发展高效的模型降阶方法、机器学习加速的多尺度模拟方法是未来的重要发展方向。其次，极端条件下的材料性能参数匮乏，雷击过程中材料处于数千摄氏度、数十兆帕的极端环境，目前对高温下碳纤维与基体的电导率、热导率、力学性能的实验测量数据十分有限，特别是热解过程中多相混合状态的性能变化规律仍不明确，制约了模拟精度的进一步提升。此外，多损伤模式的耦合演化机理仍需深入研究，当前的模型通常分别考虑基体热解、纤维断裂、界面分层等损伤模式，但对于不同损伤模式之间的相互作用机制（如分层导致的电流路径改变进一步加剧基体热解，纤维断裂产生的裂纹引发应力集中加速界面脱粘）的研究仍不够深入，缺乏统一的多损伤耦合准则。雷击过程中的动态效应，如电流脉冲的动态传播、应力波的反射与叠加、热冲击效应等，在现有模拟中通常被简化或忽略，而这些动态效应对损伤的萌生和扩展有显著影响，尤其是在纳秒到微秒级的初始阶段，动态效应的作用甚至超过准静态效应。未来的研究方向将集中在四个方面：一是发展数据驱动的多尺度模拟方法，通过机器学习算法建立微观结构与宏观性能之间的映射关系，将多尺度模拟的计算效率提升2-3个数量级，实现大型结构的实时雷击损伤预测；二是构建极端条件下的材料性能数据库，结合原位实验测试和分子动力学模拟，获取宽温度范围（室温到3000℃）、宽应变率范围（准静态到10^6s^-1）下的复合材料各组分性能参