CN109783848A - 单一缓上升雷电流分量作用下碳纤维复合材料损伤计算方法 - Google Patents

Abstract

本发明公开一种单一缓上升雷电流分量作用下碳纤维复合材料损伤计算方法，利用非破坏雷电流作用下得到的碳纤维复合材料的动态阻抗曲线，外推得到预先设计的雷电损伤仿真条件下碳纤维复合材料在规定单一缓上升速率雷电流分量作用下的各向异性电导率，作为碳纤维复合材料热电耦合模型中材料电导率的初始条件，更好地模拟碳纤维复合材料的真实雷电效应，更为准确地得到碳纤维复合材料的雷击损伤，特别是雷电损伤面积、深度与雷电流分量参数之间的相关关系，为碳纤维复合材料层合板配方、工艺的研究及其性能提升和工程应用奠定理论基础。

Description

单一缓上升雷电流分量作用下碳纤维复合材料损伤计算方法

技术领域

本发明属于碳纤维复合材料雷电损伤的仿真计算方法，涉及一种单一缓上升雷电流分量作用下碳纤维复合材料损伤计算方法。

背景技术

碳纤维复合材料既具有低密度、高强度、高模量、耐高温、耐化学腐蚀等特性，又具有纺织纤维的柔软可加工性，广泛应用于航空航天、军事及民用工业等各个领域。随着飞机设计的改进和碳纤维复合材料技术的进步，碳纤维增强型聚合物复合材料CFRP(CarbonFiber Reinforced Polymers)在大型民用飞机、军用飞机、无人机及隐形飞机上的用量不断增长，从1960年麦道公司DC-9机型上CFRP的用量不足1％，到2011年末波音B787的主翼、尾翼、机体、地板等结构的50％用的是CFRP材料，空客A350XWA上CFRP材料所占的比例达到53％。

相比较飞机中传统使用的铝、钢和钛合金材料，CFRP的电传导性能差。一般来讲，CFRP层合板经向方向的电阻率为10^-5Ω·m量级、横向平面方向的电阻率为10^-1Ω·m量级、深度/厚度方向的电阻率更大。这就使得CFRP层合板在雷击情况下无法像金属材料那样具有短时间使积累的电荷迅速转移或扩散的能力，这部分积聚的能量以焦耳热的形式使得CFRP温度急剧升高，从而导致CFRP的纤维断裂、树脂热解、深度分层等严重损伤。

欧盟和美军标规定了航空器雷电直接的试验要求和雷电流分量，其中雷电分量包括分量A(首次雷电回击分量)或Ah(首次雷电回击的过渡分量)、B(中间电流分量)、C/C*(持续电流分量)和D(后续回击分量)电流波，其中雷电流分量A、Ah和D波均为峰值高(分别为200kA、150kA、100kA)、上升速率快的雷电流分量；雷电流分量B可以是平均电流2kA、上升时间短、持续时间为几个毫秒的双指数波，也可以是上升较为缓慢的方波电流；雷电流分量C为上升时间缓慢、持续时间几百毫秒的电流波。

自CFRP问世以来，诸多学者就将研究热点集中在其机械特性的研究方面，得出了机械冲击参量与CFRP的抗拉伸强度、抗压缩强度及损伤区域和损伤深度之间的关系规律。目前而言，关于碳纤维复合材料雷电损伤的研究获得了越来越广泛的关注，许多研究者在实验手段缺乏的情况下，通过建立单一雷电流A分量作用下碳纤维复合材料雷电损伤的热电耦合模型，通过仿真计算初步获得了碳纤维复合材料雷电损伤面积、损伤深度的影响规律，其他雷电流分量，比如雷电流Ah分量、B分量、C分量、D分量或者H分量的雷电损伤研究并文献报道甚少。除此之外，现有研究文献中碳纤维雷电损伤的仿真模型一般以温度的临界值或者阻抗随热解度变化的规律而建立的，特别指出的是热电耦合模型的初始条件均是静态(直流)小电流情况下碳纤维复合材料的电导率，完全忽略了雷电流作用下碳纤维复合材料电导率与静态直流小电流的显著差异，以及碳纤维复合材料雷电流作用下的非线性特征，使得仿真计算的结果与实际雷击下的雷电损伤存在较大差异。发明专利ZL2015104538855公开了“碳纤维复合材料非破坏性雷电流作用下阻抗特性测量方法及测量装置”，相关文献的研究结果也表明：碳纤维复合材料由于其结构、工艺特性，在雷电流作用下其电导率呈现明显的非线性特性。

发明内容

本发明的目的在于提供一种单一缓上升雷电流分量作用下碳纤维复合材料损伤计算方法，准确地得到碳纤维复合材料的雷击损伤与雷电流分量参数之间的关系，为碳纤维复合材料配方、工艺的研究及应用提供理论依据。

为实现上述目的本发明采用如下方案：

单一缓上升雷电流分量作用下碳纤维复合材料损伤计算方法，包括如下步骤：

(1)、建立缓上升速率的雷电流分量试验平台，获得碳纤维复合材料缓上升速率、长持续时间雷电流分量作用下的动态伏安特性曲线；

(2)、对得到的碳纤维复合材料非破坏雷击作用下的伏安特性曲线进行数值拟合，获得碳纤维复合材料阻抗或电导率与波形参数之间的数学表达式；

(3)、按照预先设计的碳纤维复合材料雷电损伤仿真条件，将碳纤维复合材料的电导率按照步骤(2)获得的数学表达式进行外推，计算碳纤维复合材料在规定雷电流峰值的单一缓上升速率雷电流B分量或C分量作用下的各向异性电导率，雷电流峰值外推到100A至几千A，作为碳纤维复合材料缓上升速率雷电流热电耦合模型中材料电导率的初始条件；

(4)、设置拟建模仿真的碳纤维复合材料层合板的铺层结构，材料的密度、比热、热导率及机械强度参数；

(5)、设置碳纤维复合材料雷电损伤的仿真模型边界条件，包括环境温度、临界温度，雷击作用过程中碳纤维复合材料与周围环境的热传导及辐射系数；

(6)、划分碳纤维复合材料热电耦合模型的仿真计算网格，设定注入的单一缓上升雷电流分量参数，计算雷电流与碳纤维复合材料作用过程中存在的热电效应；

(7)、当碳纤维复合材料的温度升高到临界值之后，碳纤维复合材料内部树脂出现热解，随着材料热解度的增加，碳纤维复合材料的导电特性、导热特性及机械性能均会发生巨大变化，其导电特性由原来的绝缘或高阻状态转变为良导体；

(8)、根据仿真计算得到的碳纤维复合材料在单一雷电流分量作用下的温度及热解度分布，对其雷电损伤面积和损伤深度进行分析。

进一步，雷电流分量试验平台包括雷电流分量方波发生器，雷电流分量方波发生器高压端与被测碳纤维试样上表面电气连接，雷电流分量方波发生器低压端与被测碳纤维试样的下表面电气连接并接地；

还包括获取被试碳纤维复合材料上表面电压的脉冲电压取样单元和获取被试碳纤维复合材料下表面电流的雷电流取样单元，脉冲电压取样单元和雷电流取样单元与计算机测控分析单元相连。

进一步，所述雷电流分量方波发生器由RLC电路、CROWBA回路或L-C多链网络回路构成。

进一步，所述雷电流分量方波发生器包括可控直流充电电源及并联在可控直流充电上的储能电感-电容LC链，储能电感-电容LC链高压端串联放电开关，放电开关与被测碳纤维试样上表面电气连接，储能电感-电容LC链低压端与被测碳纤维试样的下表面电气连接并接地，通过控制储能电感-电容LC链的电感和电容容量获得缓上升速率的方波雷电流分量。

进一步，雷电流分量试验平台雷电流峰值范围为几A到几十A，上升时间从几百微秒至毫秒量级；或者上升时间为几十微秒幅值为几千A。

本发明具有以下有益效果：

本发明单一缓上升雷电流分量作用下碳纤维复合材料损伤计算方法，加入了碳纤维复合材料动态电导率的边界条件,即利用非破坏雷电流作用下得到的碳纤维复合材料的动态阻抗曲线，外推得到预先设计的雷电损伤仿真条件下碳纤维复合材料在规定单一缓上升速率雷电流分量作用下的各向异性电导率，作为碳纤维复合材料热电耦合模型中材料电导率的初始条件，更好地模拟碳纤维复合材料的真实雷电效应，更为准确地得到碳纤维复合材料的雷击损伤，特别是雷电损伤面积、深度与雷电流分量参数之间的相关关系，为碳纤维复合材料层合板配方、工艺的研究及其性能提升和工程应用奠定理论基础。

附图说明

图1是本发明雷电流分量试验平台原理框图；

图2是雷电流分量L-C多链网络发生回路；

图3是本发明缓上升速率非破坏雷电流分量作用下碳纤维复合材料动态电导率的测试流程图；

图4是本发明碳纤维复合材料在单一缓上升速率雷电流分量作用下雷电损伤的仿真计算流程图。

具体实施方式

下面结合附图和具体实施例对本发明作进一步详细描述，但不作为对本发明的限定。

参见图1，本发明的雷电流分量试验平台包括可控直流充电电源1、储能电感-电容LC链2、放电开关3、被测碳纤维试样4、雷电流取样单元5、脉冲电压取样单元6和计算机测控分析单元7。储能电感-电容LC链2并联在可控直流充电电源1上，储能电感-电容LC链2高压端串联放电开关3，可控直流充电电源1、储能电感-电容LC链2、放电开关3组成雷电流分量方波发生器，雷电流分量方波发生器高压端与被测碳纤维试样4上表面电气连接，雷电流分量方波发生器低压端与被测碳纤维试样4的下表面电气连接并接地。

还包括获取被试碳纤维复合材料4上表面电压的脉冲电压取样单元6和获取被试碳纤维复合材料4下表面电流的雷电流取样单元5，脉冲电压取样单元6和雷电流取样单元5与计算机测控分析单元7相连。通过控制储能电感-电容LC链2的电感和电容容量，可以获得缓上升速率的方波雷电流B分量。

参见图2，L-C多链网络回路可以产生缓上升速率的方波雷电流B分量波形，为了确保方波电流持续时间内峰值电流的稳定性，一般L-C多链方波回路的链数应大于8，图2中就是LC链数为12的方波电流发生回路。

按照标准要求，方波电流的持续时间为5ms，其LC多链网络的回路参数选择按照如下步骤：

依据下面的公式(1)：

公式(1)中，L_Σ、C_Σ分别为电感电容LC链中各链电感L₁、L₂……、L₁₁、L₁₂和各链电容C₁、C₂……、C₁₁、C₁₂之和。因此，只要电感电容LC多链回路中电感的电感量和电容的电容量满足方程(1)，就可以产生满足标准规定要求的方波雷电流B分量，可以按照表1选择参数：

表1：雷电流B分量回路参数的选择

图2中，可控直流充电电源由调压器Tr、变压器Tt、整流二极管D、充电电阻R₁组成；开关S与电阻R₂组成储能电容安全泄放回路。

参见图3，按照表1参数和图1电路建立试验平台，将被试碳纤维复合材料4串接在图1回路中，被试碳纤维复合材料4的上表面与雷电流发生回路的高压端电气连接，被试碳纤维复合材料4的下表面与雷电流发生回路的低压端电气连接。按照发明专利ZL2015104538855的测量方法，获得碳纤维复合材料雷电流作用下的准动态伏安特性曲线。对得到的碳纤维复合材料非破坏雷击作用下的伏安特性曲线进行数值拟合，获得雷电流B分量作用下碳纤维复合材料阻抗或电导率与波形参数之间的数学表达式。

同理，可以获得雷电流C分量作用下碳纤维复合材料阻抗或电导率与波形参数之间的数学表达式。

参见图4，本发明的碳纤维复合材料雷电损伤的仿真计算方法如下：

1)按照预先设计的碳纤维复合材料雷电损伤仿真条件，计算碳纤维复合材料在规定参数(电流峰值)的单一雷电流B分量(或C分量)作用下的各向异性电导率，作为碳纤维复合材料热电耦合模型中材料电导率的初始条件。

2)根据实际情况设置碳纤维复合材料层合板的铺层结构，材料的密度、比热、热导率及机械强度等参数；

3)设置碳纤维复合材料雷电损伤的仿真模型边界条件，包括环境温度，雷击作用过程中碳纤维复合材料与周围环境的热传导及辐射系数等。

4)划分碳纤维复合材料热电耦合模型的仿真计算网格，设定注入的雷电流分量参数，仿真计算雷电流与碳纤维复合材料作用过程中存在的热电效应。

5)当碳纤维复合材料的温度升高到临界值之后，碳纤维复合材料内部树脂出现热解。随着材料热解度的增加，碳纤维复合材料的导电特性、导热特性及机械性能均会发生巨大变化，其导电特性由原来的绝缘或高阻状态转变为良导体。

根据仿真计算得到的碳纤维复合材料在单一雷电流B分量(或C分量)作用下的温度及热解度分布，对其雷电损伤面积和损伤深度进行分析。

最后应该说明的是：以上实施例仅用于说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细说明，所属领域的普通技术人员应当理解：依然可以对本发明的具体实施方式进行修改或者等同替换，而未脱离本发明精神和范围的任何修改或者等同替换，其均应涵盖在本权利要求范围当中。

Claims (5)

1.单一缓上升雷电流分量作用下碳纤维复合材料损伤计算方法，其特征在于包括如下步骤：

(1)、建立缓上升速率的雷电流分量试验平台，获得碳纤维复合材料缓上升速率、长持续时间雷电流分量作用下的动态伏安特性曲线；

(2)、对得到的碳纤维复合材料非破坏雷击作用下的伏安特性曲线进行数值拟合，获得碳纤维复合材料阻抗或电导率与波形参数之间的数学表达式；

(3)、按照预先设计的碳纤维复合材料雷电损伤仿真条件，将碳纤维复合材料的电导率按照步骤(2)获得的数学表达式进行外推，计算碳纤维复合材料在规定雷电流峰值的单一缓上升速率雷电流B分量或C分量作用下的各向异性电导率，雷电流峰值外推到100A至几千A，作为碳纤维复合材料缓上升速率雷电流热电耦合模型中材料电导率的初始条件；

(4)、设置拟建模仿真的碳纤维复合材料层合板的铺层结构，材料的密度、比热、热导率及机械强度参数；

(5)、设置碳纤维复合材料雷电损伤的仿真模型边界条件，包括环境温度、临界温度，雷击作用过程中碳纤维复合材料与周围环境的热传导及辐射系数；

(6)、划分碳纤维复合材料热电耦合模型的仿真计算网格，设定注入的单一缓上升雷电流分量参数，计算雷电流与碳纤维复合材料作用过程中存在的热电效应；

(7)、当碳纤维复合材料的温度升高到临界值之后，碳纤维复合材料内部树脂出现热解，随着材料热解度的增加，碳纤维复合材料的导电特性、导热特性及机械性能均会发生巨大变化，其导电特性由原来的绝缘或高阻状态转变为良导体；

(8)、根据仿真计算得到的碳纤维复合材料在单一雷电流分量作用下的温度及热解度分布，对其雷电损伤面积和损伤深度进行分析。

2.根据权利要求1所述的方法，其特征在于：雷电流分量试验平台包括雷电流分量方波发生器，雷电流分量方波发生器高压端与被测碳纤维试样(4)上表面电气连接，雷电流分量方波发生器低压端与被测碳纤维试样(4)的下表面电气连接并接地；

还包括获取被试碳纤维复合材料(4)上表面电压的脉冲电压取样单元(6)和获取被试碳纤维复合材料(4)下表面电流的雷电流取样单元(5)，脉冲电压取样单元(6)和雷电流取样单元(5)与计算机测控分析单元(7)相连。

3.根据权利要求2所述的方法，其特征在于：所述雷电流分量方波发生器由RLC电路、CROWBA回路或L-C多链网络回路构成。

4.根据权利要求3所述的方法，其特征在于：所述雷电流分量方波发生器包括可控直流充电电源(1)及并联在可控直流充电(1)上的储能电感-电容LC链(2)，储能电感-电容LC链(2)高压端串联放电开关(3)，放电开关(3)与被测碳纤维试样(4)上表面电气连接，储能电感-电容LC链(2)低压端与被测碳纤维试样(4)的下表面电气连接并接地，通过控制储能电感-电容LC链(2)的电感和电容容量获得缓上升速率的方波雷电流分量。

5.根据权利要求1-4任一项所述的方法，其特征在于：雷电流分量试验平台雷电流峰值范围为几A到几十A，上升时间从几百微秒至毫秒量级；或者上升时间为几十微秒幅值为几千A。