碳纤维复合材料雷击损伤的数值模拟与防护研究综述

碳纤维复合材料雷击损伤的数值模拟与防护研究综述碳纤维复合材料凭借比强度高、比模量大、耐腐蚀等优异性能，已成为航空航天、风电、轨道交通等高端装备领域的核心结构材料。在民用航空领域，新一代客机中碳纤维复合材料的用量占结构总重量的比例已超过50%；在风电行业，单机容量10MW以上的叶片几乎全部采用碳纤维主梁结构。但碳纤维属于导电性能远低于金属的半导体材料，复合材料的层状叠层结构进一步导致其电导率、热导率呈现显著的各向异性，当装备在运行过程中遭遇雷击时，高达200kA的峰值电流会在材料内部产生集中的焦耳热效应与冲击应力，引发基体热解、纤维氧化、层间分层甚至贯穿性断裂等损伤，严重威胁结构的服役安全。因此，明确碳纤维复合材料雷击损伤的演化机制、构建高精度数值模拟方法、发展高效防护技术，已成为复合材料工程应用领域亟待解决的关键问题。碳纤维复合材料雷击损伤的多物理场耦合机制雷击作用下碳纤维复合材料的损伤演化是一个典型的强瞬态、多物理场耦合过程，涉及电流传导、热量传递、材料热化学分解、结构力学响应等多个环节的相互作用。相较于金属材料的雷击损伤仅以表面熔化、气化为主，复合材料的损伤呈现出显著的复杂性和隐蔽性。电流在复合材料内部的传导路径由纤维铺设方向与层间界面特性共同决定。碳纤维的轴向电导率通常可达10^4~10^5S/m，而横向电导率仅为轴向的1/100~1/1000，树脂基体的电导率更是低至10^-12S/m量级。当雷电流注入复合材料结构时，电流优先沿纤维轴向传导，若遇到层间界面或垂直于纤维方向的区域，会在局部形成电流集中点，导致该区域的焦耳热生成速率急剧提升。研究表明，当雷电流分量A（峰值200kA，作用时间10/350μs）注入碳纤维层合板时，电流注入点周边5mm范围内的电流密度可超过10^9A/m²，对应的焦耳热功率密度足以在微秒级时间内将树脂基体加热至分解温度。热场演化与材料热化学分解是损伤扩展的核心驱动力。树脂基体的热解过程通常始于300℃左右，在600℃以上基本完成，分解过程中会吸收大量潜热，同时生成CO、CO₂等气体产物和固体碳渣。热解气体的积聚导致层间孔隙压力快速上升，当压力超过层间剪切强度时，会引发层间分层损伤；而碳渣的产生会显著提升基体的导电性能，进一步改变电流的传导路径，形成“电流集中-焦耳热升温-基体热解-导电通道扩展”的正反馈效应。实验观测显示，200kA雷电流作用下，碳纤维复合材料层合板的内部最高温度可超过3000℃，不仅会导致基体完全分解，还会引发碳纤维的升华与氧化，纤维的强度损失率可达80%以上。冲击力学响应是雷击损伤的重要组成部分。雷电流的瞬态放电过程会产生两种力学效应：一是热冲击效应，材料局部快速升温导致的热膨胀受周围冷区域约束，产生高达数百MPa的热应力；二是等离子体冲击效应，雷击通道产生的高温等离子体与材料表面相互作用，形成峰值可达10MPa的冲击波，向材料内部传播。两种力学效应叠加，会在层合板内部产生复杂的应力波反射与折射，进一步加剧层间分层和基体开裂。对于厚度超过10mm的复合材料层合板，雷击造成的分层损伤面积通常是表面烧蚀面积的5~10倍，且大量损伤位于结构内部，常规无损检测手段难以完全识别，成为结构安全的隐形隐患。雷击损伤数值模拟方法的研究进展数值模拟是揭示复合材料雷击损伤演化规律、降低实验成本的重要手段，目前主流的模拟方法可分为多物理场耦合建模、等效损伤模型、多尺度模拟三大类，不同方法的适用场景和精度各有差异。多物理场直接耦合模型是当前应用最广泛的模拟方法，通过求解麦克斯韦方程、热传导方程、力学平衡方程以及材料热解动力学方程的耦合方程组，实现对雷击全过程的动态模拟。在电导率建模方面，现有研究通常将其视为温度与热解程度的函数：当温度低于基体热解起始温度时，电导率由纤维和基体的固有属性以及纤维铺设角度决定；当温度处于热解区间时，引入热解转化率变量，采用混合法则计算电导率的动态变化；当温度超过热解完成温度后，考虑碳渣的导电特性，将电导率提升2~3个数量级。热场计算中需要考虑焦耳热、热解潜热、热辐射、热对流等多种热源与散热项，力学场则引入温度相关的材料力学性能退化模型，以及基于内聚力单元的层间分层模拟方法。已有研究通过多物理场耦合模型，成功预测了不同雷电流参数下复合材料层合板的表面烧蚀面积、内部温度分布和分层损伤范围，模拟结果与实验测试的误差可控制在15%以内。但该类方法的计算量极大，通常只能对厘米级尺寸的试样进行模拟，难以应用于大型整体结构的雷击损伤分析。等效损伤模型通过简化物理过程，大幅提升了模拟效率，适用于工程结构的快速评估。这类模型的核心思路是将雷击的热-电-力耦合效应转化为等效的热输入或冲击载荷，避免直接求解复杂的多物理场方程。其中，“能量等效法”是最具代表性的方法，通过实验标定不同雷电流分量对应的总能量，将其以面热源或体热源的形式施加到结构表面，再通过热传导分析计算温度场分布，结合材料的热损伤阈值（基体热解温度、纤维升华温度）判断损伤区域。还有研究提出了“损伤阈值法”，基于大量实验数据建立雷电流参数（峰值、电荷量、作用时间）与损伤特征（烧蚀深度、分层面积、强度保留率）之间的经验公式，直接通过电流参数快速估算损伤程度。等效损伤模型的计算效率比多物理场耦合模型高两个数量级以上，但精度依赖于实验数据的积累，对于新型复合材料体系或复杂结构形式的适用性较差，通常需要针对具体材料进行参数标定。多尺度模拟方法是当前数值模拟领域的研究热点，旨在从微观、细观、宏观多个尺度揭示损伤演化的本质机制。微观尺度模拟聚焦于单根碳纤维与基体界面的电流传导与热损伤过程，通过分子动力学方法模拟树脂分子在高温下的断键、热解过程，以及碳纤维与基体界面的脱粘行为；细观尺度模拟以纤维束和基体组成的RepresentativeVolumeElement（代表性体积单元）为研究对象，考虑纤维的随机分布和接触状态，计算细观结构对电流传导路径和热场分布的影响，揭示层间电导率、热导率的各向异性来源；宏观尺度模拟则将细观尺度获得的等效材料参数输入到结构级模型中，实现对大型结构的雷击损伤分析。多尺度模拟方法能够建立材料微观结构与宏观雷击性能之间的关联，为新型抗雷击复合材料的设计提供理论指导，但目前不同尺度之间的参数传递和模型耦合仍存在诸多技术难点，距离工程化应用还有一定距离。雷击防护技术的研究现状与发展趋势为提升碳纤维复合材料的抗雷击性能，目前已发展出表面防护、结构改性、新型材料体系三大类防护技术，不同技术的防护效果、重量代价和应用场景各有不同。表面金属化防护是当前航空航天领域应用最成熟的防护方案，其核心思路是在复合材料表面铺设一层高导电的金属层，将雷电流快速疏导至结构其他部位，避免电流直接注入复合材料内部。常用的金属防护材料包括铜箔、铝箔、金属网、金属喷涂涂层等，其中膨胀铜网因重量轻、柔韧性好，能够适应复杂曲面结构，是民用客机复合材料外表面的主流防护方案。研究表明，厚度为0.1mm的膨胀铜网可将复合材料表面的面电导率提升至10^6S/m以上，能够承受200kA雷电流的冲击，表面烧蚀深度可降低90%以上。但金属防护层也存在显著缺点：一是密度大，会增加结构重量，对于民航客机而言，每增加1kg重量每年会增加约1000美元的燃油成本；二是金属与复合材料的热膨胀系数不匹配，在高低温循环环境下容易出现涂层开裂、脱落等问题；三是金属层容易发生腐蚀，需要额外的防腐处理，进一步提升了维护成本。结构改性防护技术通过调整复合材料自身的结构设计或组分，提升其本征抗雷击性能，避免额外的防护层重量。这类技术主要包括两种路径：一是在树脂基体中添加导电纳米填料，如碳纳米管、石墨烯、纳米金属颗粒等，构建贯穿基体的导电网络，提升复合材料的横向电导率和热导率，使电流能够更均匀地在层间传导，避免局部电流集中。实验结果显示，当基体中添加5wt%的碳纳米管时，复合材料的横向电导率可提升3个数量级，雷击后的分层损伤面积可降低40%以上。二是采用三维编织或针刺结构，通过Z向导电纤维增强层间的导电性能，同时提升层间剪切强度，抑制分层损伤的扩展。三维五向编织复合材料的雷击实验表明，其Z向电导率可比常规层合板高2个数量级，在相同雷电流作用下，分层损伤面积仅为层合板的30%左右。结构改性防护技术的优势是不会显著增加结构重量，且不存在界面匹配和腐蚀问题，但目前导电填料的添加会在一定程度上降低树脂基体的力学性能，三维编织结构的制备成本较高，还无法大规模应用于主承力结构。新型防护材料体系是近年来的研究热点，旨在兼顾优异的防护性能、低重量和多功能特性。其中，纳米金属线/聚合物复合导电薄膜具有重量轻、导电性好、柔韧性高等优点，面密度仅为传统铜网的1/5~1/3，面电导率可达10^5S/m以上，能够有效疏导雷电流，同时具备优异的抗腐蚀性能和耐疲劳性能。还有研究开发了具备焦耳热自愈功能的防护涂层，当雷击造成涂层局部损伤时，通过外部施加低电压产生焦耳热，使涂层中的形状记忆聚合物或微胶囊修复剂发生反应，实现损伤的自主修复，大幅提升结构的服役寿命。此外，随着智能材料技术的发展，集雷击电流监测、损伤自诊断、防护功能于一体的智能防护系统也开始得到探索，通过在结构中集成薄膜传感器，实时监测雷击电流的参数和损伤位置，为结构的维护决策提供数据支持。这些新型防护材料和技术目前仍处于实验室研究阶段，未来需要解决规模化制备、成本控制、与基体材料界面兼容性等问题，才能逐步实现工程应用。当前碳纤维复合材料雷击防护研究仍存在多个亟待突破的方向：一是需要建立更完善的多场耦合损伤模型，考虑极端高温下碳纤维的相变、热解气体的流动与排出等过程，进一步提