Wrinkles in fiber-reinforced resin composites: Micro-stress non-destructive testing
-
摘要: 纤维增强树脂复合材料在制造和服役过程中会不可避免地产生各类缺陷，且缺陷种类多、特征尺寸分散，检测难度大。缺陷检测的最终目的是对含缺陷构件的机械力学性能给出适应性评价，为此提出了一种特别适合于纤维增强树脂复合材料的集缺陷检测和性能评价为一体的微应力无损检测方法。该方法对构件施加一定的载荷使其处于微应力状态，结合全场位移的光学测量技术，捕捉缺陷导致的异常响应，以褶皱缺陷为例给出了具体的实施过程。首先基于褶皱缺陷特征响应的预测结果设计了特定的检测方案，并基于光栅投影测量技术创新性地提出了一种测量离面位移的新方法。试验结果表明，在轴向微应力加载下，利用提出的光栅投影测量方法可以探测到褶皱缺陷导致的离面位移畸变，畸变的位置和大小反映了缺陷的位置和严重程度。同时由于使用了光-力学的综合检测方法，可以跨越对缺陷具体形貌尺寸的探查，直接获得含缺陷构件在给定工况下的力学行为响应，为构件适应性评价提供参考依据。
-
关键词:
- 纤维增强树脂复合材料 /
- 褶皱缺陷 /
- 特征响应 /
- 光栅投影 /
- 无损检测
Abstract: Various defects with dispersive feature sizes can be accumulated inevitably in fiber-reinforced resin composites during their manufacturing and service processes, which are also difficult to be detected. The final objective of defects testing is to obtain the accuracy assessment of performance of defective structures. A micro-stress non-destructive defect detection and evaluation method is proposed for fiber-reinforced resin composites, which is especially suitable for composites measurement. Combing with the optical measurement technology for full-field displacement, the abnormal responses which caused by the defects of the structure under low stress level is able to be captured. The wrinkle defect detection is taken as an example to show the measurement processes. First, a specific detection scheme is designed based on the theoretical prediction of characteristic responses of wrinkles. Then, a new method of full-field displacement measurement is innovatively proposed based on the grating projection technology. Results show that under the axial tensile loading, the distorted out-of-plane displacements caused by wrinkles can be detected by the improved grating projection technology. The distorted displacements revealed the spatial distribution and severity of defects, and the influence on structural performance degradation can be evaluated based on the degree of displacement distortion. Further, by applying the optical-mechanical detection method, the mechanical responses of the defective component under a given working condition can be obtained directly, which can provide a reference for the adaptability evaluation of the component. -
纤维增强树脂复合材料具有显著和丰富的细观结构特征，在制造成型过程中不可避免或无法完全消除内部含有的缺陷，如夹杂、孔隙、纤维褶皱、分层等[1-3]。复合材料结构内部缺陷种类多、特征尺寸分散[4-6]、危害程度不一，检测难度大。缺陷检测的最终目的是对含缺陷构件的机械力学性能给出适应性评价。现有的无损检测技术(包括射线检测技术[7-9]、超声检测技术[10-14]、热成像技术[15, 16]等)在面对具有紧密多层细观结构的纤维复合材料时，缺陷识别存在一定的困难或障碍。
在复合材料界面缺陷检测方面，Li等[17]利用微CT技术对复合材料在冲击剪切载荷作用下的裂纹、脱粘和层间分层等损伤进行了检测。Chen [18]利用CT技术对复合材料桨叶蒙皮分层缺陷进行了试验研究和分析。虽然断层扫描CT技术对纤维褶皱及脱粘、分层缺陷的检测能力较好，但该方法受测试材料几何形状的限制，采集数据规模大，检测非常耗时，此外射线检测过程中存在辐射污染。超声波信号在遇到分层、孔隙这类缺陷时明显衰减，传统的超声检测方法可以较好地识别这类缺陷，而在识别褶皱缺陷方面有很大的困难。Larrañaga-Valsero等[11]比较了三种不同的超声技术在检测和表征用于风力涡轮机叶片的典型玻纤-碳纤混杂层合板内褶皱缺陷的有效性：基于全聚焦方法的全矩阵捕捉技术、相控阵检测技术和超声波水浸聚焦探测。试验结果表明，全矩阵捕捉方法得到了较好的试验结果，而相控阵技术和水浸聚焦测试应用于褶皱检测时误差较大，面临更多的挑战。值得一提的是激光剪切散斑技术在航空工业的复合结构无损检测方面已经有较成熟的应用，对复合材料层间的剥离可以进行可靠的评定[19]，但对于其他缺陷类型的识别，并未形成公认可靠有效的方法。
另一方面，运用上述缺陷检测方法，即使在一定程度上探查到局部缺陷位置和形貌，仍需借助大量数据经验的积累或复杂的仿真分析才能获得含缺陷的构件在给定工况下的行为响应 [20-22]。因此，从传统缺陷检测方法到结构性能评价之间，仍存在巨大的鸿沟需要跨越。为此，本文提出了一种特别适合于纤维增强树脂复合材料的集缺陷检测和性能评价为一体的微应力无损检测方法。
首先，该方法将缺陷视为对复合材料刚度矩阵的一种扰动，考虑缺陷的种类、严重程度和空间概率分布建立扰动的数学力学模型，并将该模型植入自主研发的三维复合材料有限元程序，获得以概率分布和统计信息表征的缺陷分布及其对复合材料结构的机械力学响应影响规律，以此奠定微应力无损检测的理论基础。
理论预测发现，复合材料中不同类型的缺陷对不同的载荷具有特殊的响应模式，这种特殊响应模式往往表现为位移场的畸变和扭曲。以褶皱缺陷为例，对含褶皱缺陷的层合板施加面内拉伸载荷，离面位移会产生明显的特殊响应如位移场的畸变或扭曲，而当施加横向压力载荷时，褶皱缺陷对板的挠度响应几乎没有影响[22]。这种高度敏感的响应和完全不敏感的响应，可分别称之为“敏感”响应和“钝感”响应，同时这种“敏感”响应和“钝感”响应也为我们提供了一种缺陷检测和性能评价的可能途径。
根据理论预测，可以为复合材料构件设计特定的加载方式并确定合适的载荷大小，使构件处于微应力水平，且对应的全场变形处于光学非接触测量方法可测范围内，通过对变形信息中蕴含的畸变信息的判断和提取，获知缺陷的位置。同时由于使用了光-力学的综合检测方法，可以跨越对缺陷具体形貌尺寸的探查，而直接获得含缺陷构件在给定工况下的力学行为响应，为构件适应性评价提供参考依据。
1. 纤维增强树脂复合材料的褶皱缺陷特征响应
1.1 纤维增强树脂复合材料中均匀型褶皱的细观力学模型
从几何描述上有均匀型和线性梯度型褶皱。由于本文中考虑的含褶皱试样整体厚度较薄，沿厚度方向每层褶皱的高度变化可近似忽略，因此文中的褶皱模型简化为均匀型褶皱，如图1所示。
![]()
图 1 纤维增强树脂复合材料均匀型褶皱Figure 1. Uniform wrinkle in fiber-reinforced resin compositesA—Fiber waviness; λ—Wavelength褶皱的形状由余弦公式描述[23]，其具体表达式如下：
z(x)=z0+Acos(2πx/λ),|x|⩽ (1) 式中：x 是平行于纤维的方向；z是沿厚度方向；z0是对应于无褶皱时层合板每层的高度；A和λ分别是褶皱的波幅和波长。
当一个代表性体积单元(RVE)包含一个完整的均匀型褶皱时，可以通过两步均匀化方法计算RVE的等效刚度矩阵，如图2所示。首先假设层合板不存在褶皱缺陷，进行沿厚度方向的均匀化计算，得到一个等效的层合板的刚度矩阵
{C^*} ；其次，令该等效层合板的纤维弯曲形成褶皱，在水平方向离散化并进行水平方向的均匀化计算得到最终的等效刚度矩阵{C^{**}} 。垂直方向均匀化和水平方向均匀化具体计算过程可作者论文参考文献[24]附录A。![]()
图 2 纤维增强树脂复合材料中含均匀型褶皱的RVE等效刚度计算模型Figure 2. Calculation model of equivalent stiffness of RVE with uniform wrinkle in fiber-reinforced resin compositesCij—Ply stiffness matrix; C*—Equivalent stiffness matrix of wrinkle-free laminates; C**—Equivalent stiffness matrix of laminates with uniform wrinkle; θk—Orientation angle of kth ply; φ—Out-of-plane misalignment1.2 纤维增强树脂复合材料中褶皱缺陷的分散性
在实际的复合材料构件中，褶皱的物理形态及其空间分布具有分散性，事先一般不能预知缺陷所在的位置。如Mukhopadhyay等[25]对弯曲复合材料结构内部的褶皱缺陷检测发现其呈现随机性特点，Lemanski等[26]利用多场影像分析技术对局部小范围的复合材料层合板截面处的纤维偏转角的分布特点进行统计学研究，纤维偏转角呈现近似正态分布的特点。此外，Riddle等[27]对已经完成制造或运行中的叶片缺陷进行了大量的统计与分析，发现褶皱波长、波幅及其偏转角度均近似呈现正态分布或威布尔分布。
基于上述理由，褶皱缺陷检测所采取的方法应该是一全场测量方法，同时还需考虑对构件施加合适的载荷类型及大小，使缺陷所造成的异常能够被显示出来并处于测量所覆盖的范围内。
为此本文建立了含分散性褶皱的纤维增强树脂复合材料薄板模型，如图3所示，该层合板长250 mm，宽25 mm，采用对称铺层，铺层顺序为[0/90/0/90]s，每层的厚度为0.1 mm，总厚度为0.8 mm。
![]()
图 3 含分散性褶皱的纤维增强树脂复合材料薄板模型Figure 3. Thin plate model for fiber-reinforced resin composites with dispersed wrinkles在下文的算例与实验中，组成该层合板的单向板采用碳纤维(CF)增强环氧树脂(EP)复合材料，其性能参数如表1所示。其中褶皱的几何尺寸服从下式所示的正态分布，而褶皱在结构上为随机分布。
表 1 碳纤维(CF)/环氧树脂(EP)复合材料弹性参数Table 1. Elastic parameters of carbon fiber (CF)/epoxy (EP) compositesE11 /GPa (E22 /E33) /GPa G23 /GPa G31 /GPa G12 /GPa ν21 ν32 ν31 133.3 9.09 3.16 7.24 7.23 0.261 0.436 0.261 Notes: E11, E22, E33—Elastic modulus (direction 11, 22, 33); G12, G23, G31—Shear modulus (direction 12, 23, 31); v21, v32, v31—Poisson’s ratio (direction 21, 32, 31). F(\gamma)=\frac{1}{\sqrt{2 {\text{π}}} S} \exp \left[-\frac{\left(\gamma-\gamma_{0}\right)^{2}}{2 S^{2}}\right] (2) 式中：γ表示褶皱的波幅与波长之比(A/λ称为波纹比)；γ0表示波纹比的均值。当γ为正数时，表示纤维在厚度方向的偏转沿z轴正向。相反地，当
\gamma 为负数时，意味着纤维在厚度方向的偏转沿z轴负方向。由于褶皱的方向有可能朝z轴的正向或负向，故此处波纹比的均值为0。1.3 CF/EP复材层合板特征响应分析
对含分散性褶皱缺陷的CF/EP复材层合板模型考虑不同的加载方式，如图4所示。
![]()
图 4 CF/EP复材层合板的不同加载方式Figure 4. Schematic of different loading modes for CF/EP composite laminatesσx—Axial tensile stress; My—Bending moment; τxy—Shearing stress基于作者自主开发的含分散性褶皱缺陷的有限元计算程序，图5给出了拉伸载荷下含分散性褶皱缺陷的层合板面内和离面位移响应。其中面内位移u、v分别表示沿x和y方向的位移，用位移等值线表示，离面位移w则采用三维曲面显示。作为对比，图6是不含缺陷的层合板在拉伸载荷下的位移。
![]()
图 5 含随机褶皱缺陷的CF/EP复材层合板拉伸位移响应Figure 5. Displacement responses of CF/EP composite laminates with randomly dispersed wrinkle defects under tensile loadu—In-plane displacement in x-direction; v—In-plane displacement in y-direction; w—Out-of-plane displacement in z-direction![]()
图 6 无缺陷的CF/EP复材层合板拉伸位移响应Figure 6. Displacement responses of pristine CF/EP composite laminates under tensile loadu—In-plane displacement in x-direction; v—In-plane displacement in y-direction; w—Out-of-plane displacement in z-direction比较图5和图6可见当复合材料层合板内部存在褶皱缺陷时，在单轴拉伸载荷作用下，其面内位移和离面位移均会产生特殊响应如位移场的畸变或扭曲，且离面位移的畸变程度尤为显著。
图7、图8分别为含随机褶皱缺陷及不含缺陷层合板的在弯曲载荷下的位移响应。对比可知，含随机褶皱缺陷层合板的沿y方向的面内位移响应对弯曲载荷敏感，而沿x方向的面内位移及离面位移不受褶皱缺陷的影响，其响应与无缺陷的理想材料一致。
![]()
图 7 含随机褶皱缺陷的CF/EP复材层合板弯曲位移响应Figure 7. Displacement responses of CF/EP composite laminates with randomly dispersed wrinkle defects under bending loadu—In-plane displacement in x-direction; v—In-plane displacement in y-direction; w—Out-of-plane displacement in z-direction![]()
图 8 无缺陷的CF/EP复材层合板弯曲位移响应Figure 8. Displacement responses of pristine CF/EP composite laminates under bending loadu—In-plane displacement in x-direction; v—In-plane displacement in y-direction; w—Out-of-plane displacement in z-direction图9、图10分别为含随机褶皱缺陷及不含缺陷层合板的剪切位移响应。对比可知，含随机褶皱缺陷层合板的面内位移响应对剪切载荷作用不敏感，无论是否含有缺陷，两种情况下面内位移响应非常接近。而含褶皱缺陷的层合板离面位移在靠近根部约束区域的两翼出现了独特的边缘扭曲现象，这与褶皱缺陷及各向异性复合材料的变形耦合有关。
![]()
图 9 含随机褶皱缺陷的CF/EP复材层合板剪切位移响应Figure 9. Displacement responses of CF/EP composite laminates with randomly dispersed wrinkle defects under shearing load![]()
图 10 无缺陷的CF/EP复材层合板剪切位移响应Figure 10. Displacement responses of pristine CF/EP composite laminates under shearing load上述结果表明，CF/EP复材层合板中的褶皱缺陷导致结构不同方向的位移对不同的加载方式具有完全不同的响应，在给定载荷下，某些方向的位移在分布模式和数量级上均出现显著的改变，而其他方向的位移可能完全不受缺陷影响。
这种高度敏感的响应和完全不敏感的响应，我们分别称之为“敏感”响应和“钝感”响应，同时这种“敏感”响应和“钝感”响应也为我们提供了一种用于缺陷检测和识别的工具。
根据上述三种载荷响应的计算，对于褶皱缺陷，其不同载荷下的位移响应敏感性如表2所示。
表 2 CF/EP复材层合板中褶皱缺陷对不同载荷的响应敏感性Table 2. Response sensitivity of wrinkle defects to different loading modes in CF/EP composite laminatesLoading mode Uniaxial tensile Bending Shearing Response sensitivity In-plane displacement u 
v 
Out-of-plane displacement w Note: *Sensitivity: 位移响应敏感度反映了位移畸变的程度，畸变程度越大，越容易被检出。从现实应用角度，应选择位移响应最敏感的方式进行加载荷和测量。因此对于褶皱缺陷可采用单轴拉伸或剪切加载方式，通过检测层合板的离面位移即可凸显褶皱缺陷引起的位移畸变。
课题组至今已针对褶皱和层间弱粘结缺陷开展了大量的理论分析工作[28-30]，作为比较，表3为含有层间弱粘结缺陷时，层合板在不同加载方式下的响应敏感性。
表 3 CF/EP复材层合板中层间弱粘结缺陷对不同载荷的响应敏感性[30]Table 3. Response sensitivity of weak bonding defects to different loading modes in CF/EP composite laminatesLoading mode Uniaxial tensile Bending Shearing Response sensitivity In-plane displacement u v Out-of-plane displacement w 综合表2与表3结论可见，单轴拉伸载荷下，褶皱或者弱粘结缺陷均导致离面位移的畸变，虽然此处不能根据位移的畸变反演缺陷的类型和具体形状尺寸，但反之，通过测量拉伸载荷下的离面位移，可以揭示畸变区域存在多种缺陷的可能性。
因此，本文所提出的方法，不是传统意义上的缺陷检测方法，而是为大型复材结构提供了一种全场测量快速识别异常区域以及越过缺陷的具体形式识别而直接根据位移畸变程度量化评价材料结构性能退化程度的一种方法。
2. CF/EP复材层合板检测方案的数值验证
尽管在实际纤维增强树脂复合材料结构中，褶皱的分布呈分散性，但是为了便于对上述检测理论进行分析计算和试验验证，本文制作了含单个确定性褶皱的试样进行相应的分析。
根据表2结论，考虑到单轴拉伸可以方便地在万能试验机上实现，因此本文选择拉伸加载方式和离面位移测量的组合作为检测的手段。
2.1 CF/EP复材层合板中褶皱缺陷制样
为形成不同程度的褶皱缺陷，本文采用钢棒法来制造缺陷，通过使用不同直径钢棒置入不同预浸料板之间，将褶皱缺陷波纹比引入CF/EP试样中。试样从底部往上共8层，试样I的钢棒放置在第2层与第3层的中间，试样II的钢棒放置在最底层。同时该钢棒使用聚四氟乙烯布包覆，固化过程完成后，由于聚四氟乙烯材料摩擦系数极小，因此可轻松将钢棒从层合板中取出。在制造过程中，由于树脂流动性较强，会聚集到由钢棒产生的凹槽中，待取出钢棒之后，可进一步将环氧树脂添加到由于钢棒拔出而留下的空腔中。钢棒法制作褶皱缺陷的示意图和制作过程如图11所示。
![]()
图 11 钢棒法制作含褶皱缺陷的CF/EP复材层合板试样Figure 11. Wrinkled specimen of CF/EP composite laminates made with steel bar method2.2 CF/EP试样结构与参数
含缺陷的CF/EP试样如图12所示，褶皱缺陷位于试样中部位置，缺陷的设计参数见表4。试样长250 mm，宽25 mm，采用对称铺层，铺层顺序为[0/90/0/90]s，共8层，厚度为0.8 mm。预浸料采用威海光威复合材料股份有限公司生产的USN 10000型号，单向纤维预浸料厚度为0.1 mm，纤维单位面积质量为100 g/m2，树脂含量为20%~40%。该单向板的力学性能与表1所示一致。
表 4 CF/EP试样中的缺陷参数Table 4. Defect parameters of the CF/EP specimensSpecimen number Layup sequences Wrinkle parameters Wavelength/mm Amplitude/mm Wrinkle ratio I [0/90/0/90]s 6.0 1.2 0.200 II [0/90/0/90]s 12.0 0.8 0.067 2.3 检测方案的有限元验证
对于图12所示的含褶皱缺陷试样，图13(a)、图13(b)给出了单轴拉伸下的离面位移有限元分析结果。对于试样Ⅰ，图13(a)表明当拉伸载荷大小为600 N时，对应的最大离面位移为0.169 mm，图13(b)表明当拉伸载荷大小为1200 N时，试样Ⅰ对应的最大离面位移为0.228 mm。作为对比，图13(c)是不含缺陷的理想试样在1200 N时的离面位移，由图可见，当结构中存在褶皱缺陷时，不但离面位移分布模式发生了显著改变，且离面位移的最大值增加了104的数量级，离面位移的最大值与褶皱缺陷所在的区域一致。
有限元分析结果也说明，在600~1200 N的拉伸载荷下，含褶皱缺陷试样离面位移的大小与下文中光学变形检测的范围相适应，因此根据有限元分析结果可合理设计拉伸试验的检测方案，包括加载方式、载荷大小和测试值的预估计。
![]()
图 13 CF/EP试样Ⅰ与无缺陷CF/EP试样离面位移有限元结果对比Figure 13. Comparison of finite element results on out-of-plane displacement between CF/EP specimen Ⅰ and wrinkle-free CF/EP specimen3. 双目结构光检测
光栅投影向来被用于三维表面的测量，其测量范围较大，测量精度可调，对环境要求不高。近年来，有文献将DIC(数字图像相关)方法与光栅投影方法相结合用以测量三维变形，其中光栅投影用于变形前后的三维曲面测量，DIC用于变形前后的图像数据匹配[31-32]。考虑到缺陷检测的应用场景往往面临大尺寸、大曲率结构及现场检测等实际需求，本文创新性地提出了独立利用双目结构光测量离面位移的新方法，即利用变形前后点云数据本身携带的信息进行数据重建和匹配。
根据上文分析结果，含褶皱缺陷的层合板试样，应通过检测其拉伸过程中的离面位移畸变情况以判断缺陷所在的位置。双目结构光检测直接得到的是表面高度信息，因此本文分别对加载前后的试样采集其表面高度信息z0和z1，从测试原理上可得离面位移d= z1− z0。但在实际测量中，试件拉伸后，在光测系统中物面相对位置发生变化，不仅试件本身的三维形貌发生变化，而且采集到的表面三维点云的数据量也发生变化，变形前后的点云不再存在一一对应的关系，本文创新性地提出了三维点云重构算法提取离面位移场，从而实现了基于光栅投影位移场测量的微应力缺陷检测方法。图14是点云重构算法的基本框架，该算法吸收了无网格法中形函数的构造方法，规格化节点上的点云数据取决于其邻域的半径、邻域内原始点云的数值及原始点云到该节点的距离作为权重，经最小二乘法拟合得到。
![]()
图 14 三维点云重构算法提取离面位移Figure 14. Schematic diagram of point cloud data reconstruction algorithm.3.1 实验装置
对褶皱缺陷试样进行轴向拉伸试验，拉伸试验采用双立柱台式Instron 3369型电子万能试验机进行，试验机负荷范围为0~50 kN。光学测试系统采用多频相移结构光投影结合三角测量法(见图15)获取试样表面的空间三维信息。
![]()
图 15 含褶皱缺陷的CF/EP复材层合板试样测试装置示意图Figure 15. Schematic diagram of the test device of wrinkled specimen of CF/EP composite3.2 实验结果
图16(a)为I号试样未加载时表面形貌点云，图16(b)是该试样加载到1200 N并保持稳定时采集的试样表面形貌点云。
![]()
图 16 CF/EP试样I的加载前后三维点云Figure 16. 3D point cloud of CF/EP specimen I before and after loading由于点云边缘数据误差较大，故提取图16中线框内有效区域数据进行分析。图17(a)和图17(b)将加载前后所检测到的试样表面有效点云信息分别以等高线的形式给出，由于复合材料试样为手工制样，表面平整度较差，在未施加载荷的情况下，其表面等高线并非理想的光滑曲线。
图17(c)是经点云数据处理并相减后，得到试样表面的离面位移分布。由于手工制造的褶皱缺陷在形式上也并非理想的余弦曲线，因此试样I实测的离面位移分布与图13所示的有限元分析结果有一定的偏差，有限元给出的理想褶皱缺陷所引起的离面位移在试样上分布较为均匀，但实测结果说明实际褶皱所造成的影响集中分布在试样中部的有限区域，这一分布与缺陷所在的区域吻合较好，在无缺陷区域，等值线显示其实际位移数值为0，说明远离缺陷区域，其位移的变化相对缺陷区域的变化是可以忽略的。
图18(a)、图18(b)是试样II加载前后的表面等高线，图18(c)是基于加载前后的点云信息所提取的离面位移，由图18可见：通过加载前后表面高度的测量并提取离面位移，可以很好地捕捉到缺陷所在的位置，离面位移的分布与缺陷位置具有对应一致的关系，且与计算预测结果符合较好。
表5为从200 N到1200 N的拉伸载荷下实测位移响应与有限元分析得到的缺陷处最大离面位移结果的对比。可见，两件试样在1200 N下的离面位移误差分别约为10.0%和18.3%，试样II误差相对较大可能与手工制样的缺陷形状误差及试样夹持端表面平整度较差造成加载过程中引起的附加载荷等有关。
表 5 CF/EP试样离面位移实测结果与有限元结果比较Table 5. Comparison of out-of-plane displacement of CF/EP specimens between measurement results and finite element resultsSpecimen number Load/N Measurement results/cm Finite element results/cm Error I 200 0.01078 0.0089 17.7% 500 0.01806 0.0154 14.7% 800 0.02093 0.0193 7.8% 1 000 0.02359 0.0212 10.1% 1 200 0.02508 0.0228 10.0% Ⅱ 200 0.00218 0.0027 20.1% 500 0.00918 0.0101 8.6% 1 200 0.01361 0.0115 18.3% 上述结果说明：由于褶皱缺陷导致的离面位移异常现象，可以由本文提出的双目结构光方法检测出来，且实测值与有限元结果的误差也在可接受范围内。
本文研究表明：针对大型纤维增强树脂复合材料构件，建立合理的缺陷细观力学模型并展开深入的力学响应分析可能是其缺陷检测理论中重要的一环，而在该理论的指导下通过光学全场测量方法，快速捕获响应异常区域，为后续在该异常区域进行精细化的缺陷识别和检测提供了一条可行的路径。
4. 量化评价
缺陷导致的位移畸变程度可以由缺陷区域的离面位移与无缺陷区域的离面位移之比(位移比)表示，以试样I为例，由于波纹比A/λ可以表征褶皱的严重程度，固定褶皱区范围，即取λ值为6 mm，调整波幅A作为变量来控制褶皱参数，使得A/λ取值范围为0.04~0.2。
图19为位移畸变程度随波纹比的变化趋势。虚线表示远离缺陷处的远场归一化位移(离面位移值与板厚比值)，实线表示缺陷区域的归一化位移。褶皱区域与远场非褶皱区域的位移比值表征了位移畸变程度。可见，随着波纹比的增加，远场位移几乎保持水平不变，而缺陷处的位移显著增大，因此位移畸变程度随波纹比增大而增加，同时说明畸变主要集中在缺陷所在的局部区域，对远场位移几乎没有影响。因此将位移比作为量化评价指标，可以实现不同严重程度的缺陷的比较与评价。
![]()
图 19 CF/EP复材层合板离面位移及位移比随波纹比的变化Figure 19. Normalized out-of-plane displacement and displacement ratio various with wrinkle ratio of CF/EP composite laminates此外，图20为CF/EP复材层合板离面位移及位移比随载荷的变化。可知，当波纹比为0.1875时，随着拉伸载荷在弹性范围内的增加，褶皱区域和非褶皱区域的归一化离面位移均呈线性增大的趋势，但其比值并不随加载而变化，因此若采用光-力学方法对试样进行加载，只要在弹性范围内，具体所施加的载荷大小并不会影响结论。
![]()
图 20 CF/EP复材层合板离面位移及位移比随载荷的变化Figure 20. Normalized out-of-plane displacement and displacement ratio various with applied node force of CF/EP composite laminates5. 结 论
本文建立了以概率分布和统计信息表征的缺陷力学模型，开展了含分散性缺陷的纤维增强树脂复合材料结构的力学响应规律研究，提出了一种集缺陷检测与性能评价为一体的方法，并以褶皱缺陷为例进行了具体的实施。研究结果表明：
(1) 含缺陷层合板在不同加载下对位移具有“敏感”或“顿感”响应特征，“敏感”响应反映了较大的位移畸变程度，可以为缺陷检测和识别提供理论指导。由本文数值预测结果可知，含褶皱缺陷的层合板在拉伸载荷下的离面位移响应最敏感。
(2) 针对含褶皱缺陷构件的“敏感”响应特征，创新性地提出了基于光栅投影测量技术的全场离面位移测量方法，所探测到的位移场畸变直观地反映了缺陷在试样上的位置。与剪切散斑等光测方法对大型结构的变形而导致过度的刚体运动使得散斑之间相关性丧失相比，本文基于双目结构光所提出的间接位移测量方法，测量视野大，不受刚体位移影响。
(3) 含褶皱层合板在拉伸加载下的光栅投影测量结果显示，褶皱缺陷处的离面位移异常增大，且实测的褶皱区域离面位移值印证了有限元预测结果。实测位移响应的畸变程度也为含缺陷复材结构性能评价提供了直接的参考依据。
(4) 当前纤维增强树脂复合材料的缺陷检测仍具有一定的难度和挑战性，一些复杂的复合材料构件仍依靠手工铺层裁剪，因此其表面平整度、尺寸精度与机加工产品具有很大的差距。与其他直接检测位移的方法相比，本文所提出的检测方法，可以消除表面初始误差对位移的影响，非常适合复合材料结构的变形检测。利用该方法首先在大型结构上快速发现位移畸变区域，然后再结合超声、CT等手段对该区域进行细致的扫描探查，为大型复材结构的缺陷检测工业化提供了一种可行的方法途径。
-
![]()
图 1 纤维增强树脂复合材料均匀型褶皱
Figure 1. Uniform wrinkle in fiber-reinforced resin composites
A—Fiber waviness; λ—Wavelength
![]()
图 2 纤维增强树脂复合材料中含均匀型褶皱的RVE等效刚度计算模型
Figure 2. Calculation model of equivalent stiffness of RVE with uniform wrinkle in fiber-reinforced resin composites
Cij—Ply stiffness matrix; C*—Equivalent stiffness matrix of wrinkle-free laminates; C**—Equivalent stiffness matrix of laminates with uniform wrinkle; θk—Orientation angle of kth ply; φ—Out-of-plane misalignment
![]()
图 3 含分散性褶皱的纤维增强树脂复合材料薄板模型
Figure 3. Thin plate model for fiber-reinforced resin composites with dispersed wrinkles
![]()
图 4 CF/EP复材层合板的不同加载方式
Figure 4. Schematic of different loading modes for CF/EP composite laminates
σx—Axial tensile stress; My—Bending moment; τxy—Shearing stress
![]()
图 5 含随机褶皱缺陷的CF/EP复材层合板拉伸位移响应
Figure 5. Displacement responses of CF/EP composite laminates with randomly dispersed wrinkle defects under tensile load
u—In-plane displacement in x-direction; v—In-plane displacement in y-direction; w—Out-of-plane displacement in z-direction
![]()
图 6 无缺陷的CF/EP复材层合板拉伸位移响应
Figure 6. Displacement responses of pristine CF/EP composite laminates under tensile load
u—In-plane displacement in x-direction; v—In-plane displacement in y-direction; w—Out-of-plane displacement in z-direction
![]()
图 7 含随机褶皱缺陷的CF/EP复材层合板弯曲位移响应
Figure 7. Displacement responses of CF/EP composite laminates with randomly dispersed wrinkle defects under bending load
u—In-plane displacement in x-direction; v—In-plane displacement in y-direction; w—Out-of-plane displacement in z-direction
![]()
图 8 无缺陷的CF/EP复材层合板弯曲位移响应
Figure 8. Displacement responses of pristine CF/EP composite laminates under bending load
u—In-plane displacement in x-direction; v—In-plane displacement in y-direction; w—Out-of-plane displacement in z-direction
![]()
图 9 含随机褶皱缺陷的CF/EP复材层合板剪切位移响应
Figure 9. Displacement responses of CF/EP composite laminates with randomly dispersed wrinkle defects under shearing load
![]()
图 10 无缺陷的CF/EP复材层合板剪切位移响应
Figure 10. Displacement responses of pristine CF/EP composite laminates under shearing load
![]()
图 11 钢棒法制作含褶皱缺陷的CF/EP复材层合板试样
Figure 11. Wrinkled specimen of CF/EP composite laminates made with steel bar method
![]()
图 13 CF/EP试样Ⅰ与无缺陷CF/EP试样离面位移有限元结果对比
Figure 13. Comparison of finite element results on out-of-plane displacement between CF/EP specimen Ⅰ and wrinkle-free CF/EP specimen
![]()
图 14 三维点云重构算法提取离面位移
Figure 14. Schematic diagram of point cloud data reconstruction algorithm.
![]()
图 15 含褶皱缺陷的CF/EP复材层合板试样测试装置示意图
Figure 15. Schematic diagram of the test device of wrinkled specimen of CF/EP composite
![]()
图 16 CF/EP试样I的加载前后三维点云
Figure 16. 3D point cloud of CF/EP specimen I before and after loading
![]()
图 19 CF/EP复材层合板离面位移及位移比随波纹比的变化
Figure 19. Normalized out-of-plane displacement and displacement ratio various with wrinkle ratio of CF/EP composite laminates
![]()
图 20 CF/EP复材层合板离面位移及位移比随载荷的变化
Figure 20. Normalized out-of-plane displacement and displacement ratio various with applied node force of CF/EP composite laminates
表 1 碳纤维(CF)/环氧树脂(EP)复合材料弹性参数
Table 1 Elastic parameters of carbon fiber (CF)/epoxy (EP) composites
E11 /GPa (E22 /E33) /GPa G23 /GPa G31 /GPa G12 /GPa ν21 ν32 ν31 133.3 9.09 3.16 7.24 7.23 0.261 0.436 0.261 Notes: E11, E22, E33—Elastic modulus (direction 11, 22, 33); G12, G23, G31—Shear modulus (direction 12, 23, 31); v21, v32, v31—Poisson’s ratio (direction 21, 32, 31). 
下载: 导出CSV
表 2 CF/EP复材层合板中褶皱缺陷对不同载荷的响应敏感性
Table 2 Response sensitivity of wrinkle defects to different loading modes in CF/EP composite laminates
Loading mode Uniaxial tensile Bending Shearing Response sensitivity In-plane displacement u v Out-of-plane displacement w Note: *Sensitivity: 表 3 CF/EP复材层合板中层间弱粘结缺陷对不同载荷的响应敏感性[30]
Table 3 Response sensitivity of weak bonding defects to different loading modes in CF/EP composite laminates
Loading mode Uniaxial tensile Bending Shearing Response sensitivity In-plane displacement u v Out-of-plane displacement w
下载: 导出CSV
表 4 CF/EP试样中的缺陷参数
Table 4 Defect parameters of the CF/EP specimens
Specimen number Layup sequences Wrinkle parameters Wavelength/mm Amplitude/mm Wrinkle ratio I [0/90/0/90]s 6.0 1.2 0.200 II [0/90/0/90]s 12.0 0.8 0.067
下载: 导出CSV
表 5 CF/EP试样离面位移实测结果与有限元结果比较
Table 5 Comparison of out-of-plane displacement of CF/EP specimens between measurement results and finite element results
Specimen number Load/N Measurement results/cm Finite element results/cm Error I 200 0.01078 0.0089 17.7% 500 0.01806 0.0154 14.7% 800 0.02093 0.0193 7.8% 1 000 0.02359 0.0212 10.1% 1 200 0.02508 0.0228 10.0% Ⅱ 200 0.00218 0.0027 20.1% 500 0.00918 0.0101 8.6% 1 200 0.01361 0.0115 18.3%
下载: 导出CSV
-
[1] 欧阳佳斯. 带波纹/褶皱类缺陷的纤维增强树脂基复合材料压缩性能研究[D]. 武汉: 武汉理工大学, 2016. OUYANG Jiasi. Study on the compressive mechanical pro-perty of fiberreinforced composites with waviness defect[D]. Wuhan: Wuhan University of Technology, 2016(in Chinese).
[2] TAY T E, SHEN F. Analysis of delamination growth in lami-nated composites with consideration for residual thermal stress effects[J]. Journal of Composite Materials,2002,36(11):1299-1320. DOI: 10.1177/0021998302036011592
[3] 丁珊珊. 考虑孔隙形貌的CFRP复合材料超声散射机理及孔隙率检测方法研究[D]. 大连: 大连理工大学, 2017. DING Shanshan. Ultrasonic scattering mechanism and porosity detection in CFRP composite materials considering void morphology[D]. Dalian: Dalian University of Technology, 2017(in Chinese).
[4] SAENZ-CASTILLO D, MARTÍN M I, CALVO S, et al. Effect of processing parameters and void content on mechanical properties and NDI of thermoplastic composites[J]. Composites Part A: Applied Science and Manufacturing,2019,121:308-320. DOI: 10.1016/j.compositesa.2019.03.035
[5] MUKHOPADHYAY S, JONES M I, HALLETT S R. Tensile failure of laminates containing an embedded wrinkle; numerical and experimental study[J]. Composites Part A: Applied Science and Manufacturing,2015,77:219-228. DOI: 10.1016/j.compositesa.2015.07.007
[6] NELSON L J, SMITH R A. Fibre direction and stacking sequence measurement in carbon fibre composites using Radon transforms of ultrasonic data[J]. Composites Part A: Applied Science and Manufacturing,2019,118:1-8. DOI: 10.1016/j.compositesa.2018.12.009
[7] SUTCLIFFE M P F, LEMANSKI S L, SCOTT A E. Measurement of fibre waviness in industrial composite components[J]. Composites Science and Technology,2012,72(16):2016-2023. DOI: 10.1016/j.compscitech.2012.09.001
[8] CREIGHTON C J, SUTCLIFFE M, CLYNE T W. A multiple field image analysis procedure for characterisation of fibre alignment in composites[J]. Composites Part A: Applied Science and Manufacturing,2001,32(2):221-229. DOI: 10.1016/S1359-835X(00)00115-9
[9] REVOL V, PLANK B, KAUFMANN R, et al. Laminate fibre structure characterisation of carbon fibre-reinforced polymers by X-ray scatter dark field imaging with a grating interferometer[J]. NDT & E International,2013,58:64-71.
[10] ZARDAN J P, GUEUDRE C, GORNELOUP G. Study of induced ultrasonic deviation for the detection and identification of ply waviness in carbon fibre reinforced polymer[J]. NDT & E International,2013,56(jun.):1-9.
[11] LARRAÑAGA-VALSERO B, SMITH R A, TAYONG R B, et al. Wrinkle measurement in glass-carbon hybrid laminates comparing ultrasonic techniques: A case study[J]. Composites Part A: Applied Science and Manufacturing,2018,114:225-240. DOI: 10.1016/j.compositesa.2018.08.014
[12] ZHANG Z, LIU M, LI Q, et al. Visualized characterization of diversified defects in thick aerospace composites using ultrasonic B-scan[J]. Composites Communications,2020,22:100435. DOI: 10.1016/j.coco.2020.100435
[13] PARK B, AN Y, SOHN H. Visualization of hidden delamination and debonding in composites through noncontact laser ultrasonic scanning[J]. Composites Science and Technology,2014,100:10-18. DOI: 10.1016/j.compscitech.2014.05.029
[14] IBRAHIM M E, SMITH R A, WANG C H. Ultrasonic detection and sizing of compressed cracks in glass-and carbon-fibre reinforced plastic composites[J]. NDT & E International,2017,92:111-121.
[15] ELHAJJAR R, HAJ-ALI R, WEI B. An infrared thermoelastic stress analysis investigation for detecting fiber waviness in composite structures[J]. Polymer-Plastics Technology and Engineering,2014,53(12):1251-1258. DOI: 10.1080/03602559.2014.886116
[16] KATUNIN A, DRAGAN K, DZIENDZIKOWSKI M. Damage identification in aircraft composite structures: A case study using various non-destructive testing techniques[J]. Composite Structures,2015,127:1-9.
[17] LI Y, SUN B, GU B. Impact shear damage characterizations of 3D braided composite with X-ray micro-computed tomography and numerical methodologies[J]. Composite Structures,2017,176(9):43-54.
[18] CHEN X. Fractographic analysis of sandwich panels in a composite wind turbine blade using optical microscopy and X-ray computed tomography[J]. Engineering Failure Analysis,2020,111:104475. DOI: 10.1016/j.engfailanal.2020.104475
[19] 洪友仁, 何浩培, 何小元. 剪切散斑: 一种光学测量技术及其应用[J]. 实验力学, 2006, 21(6):667-688. DOI: 10.3969/j.issn.1001-4888.2006.06.001 HONG Youren, HE Haopei, HE Xiaoyuan. Shearography: An optical measurement technique and applications[J]. Journal of Experimental Mechanics,2006,21(6):667-688(in Chinese). DOI: 10.3969/j.issn.1001-4888.2006.06.001
[20] FERREIRA L M, GRACIANI E, PARÍS F. Three dimensional finite element study of the behaviour and failure mechanism of non-crimp fabric composites under in-plane compression[J]. Composite Structures,2016,149:106-113. DOI: 10.1016/j.compstruct.2016.04.022
[21] LEONG M, HVEJSEL C F, THOMSEN O T, et al. Fatigue failure of sandwich beams with face sheet wrinkle defects[J]. Composites Science and Technology,2012,72(13):1539-1547.
[22] AYMERICH F, DORE F, PRIOLO P. Prediction of impact-induced delamination in cross-ply composite laminates using cohesive interface elements[J]. Composites Science and Technology,2008,68(12):2383-2390. DOI: 10.1016/j.compscitech.2007.06.015
[23] TAKEDA T. Micromechanics model for three-dimensional effective elastic properties of composite laminates with ply wrinkles[J]. Composite Structures,2018,189:419-427. DOI: 10.1016/j.compstruct.2017.10.086
[24] 申川川, 马利, 文安戈, 等. 纤维增强树脂复合材料中的褶皱缺陷: 分散性与虚拟测试[J]. 复合材料学报, 2022, 39(3):1332-1342. DOI: 10.13801/j.cnki.fhclxb.20210518.007 SHEN Chuanchuan, MA Li, WEN Ange, et al. Wrinkles in fiber-reinforced resin composites (Part I): Heterogeneity and virtual test[J]. Acta Materiae Composites Sinica,2022,39(3):1332-1342(in Chinese). DOI: 10.13801/j.cnki.fhclxb.20210518.007
[25] MUKHOPADHYAY S, JONES M I, HALLETT S R. Compressive failure of laminates containing an embedded wrinkle; experimental and numerical study[J]. Composites Part A: Applied Science and Manufacturing,2015,73:132-142. DOI: 10.1016/j.compositesa.2015.03.012
[26] LEMANSKI S L, SUTCLIFFE M P F. Compressive failure of finite size unidirectional composite laminates with a region of fibre waviness[J]. Composites Part A: Applied Science and Manufacturing,2012,43(3):435-444.
[27] RIDDLE T, CAIRNS D, NELSON J. Characterization of manu-facturing defects common to composite wind turbine blades: Flaw characterization[C]. 52nd AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics and Materials Conference. Denver, 2011.
[28] SHEN C, MA L, XU P, et al. Virtual testing of mechanical response of composite plates with normally distributed wrinkles[J]. Composite Structures,2019,229:111440. DOI: 10.1016/j.compstruct.2019.111440
[29] MA L, SHEN C C, WEN A G, et al. Dependence between displacement distortion and heterogeneity of vegetable fiber composites[J]. Mechanics of Advanced Materials and Structures, 2022, 29: 1889079 .
[30] SHEN C C, MA L, WEN A G, et al. Composite plates with randomly distributed weak bonding: Heterogeneity and virtual testing[J]. Mechanics of Advanced Materials and Structures,2022,29:1928344.
[31] FELIPE-SESÉ L, LÓPEZ-ALBA E, SIEGMANN P, et al. Integration of fringe projection and two-dimensional digital image correlation for three-dimensional displacements measurements[J]. Optical Engineering,2016,55(12):121711. DOI: 10.1117/1.OE.55.12.121711
[32] WU Z J, GUO W B, PAN B, et al. A DIC-assisted fringe projection profilometry for high-speed 3D shape, displacement and deformation measurement of textured surfaces[J]. Optics and Lasers in Engineering, 2021, 142: 106614.
-
期刊类型引用(3)
1. 凌圣博，姚启，王乐蓬，梅杰，吴振. 褶皱缺陷对复合材料拉伸性能的影响规律:试验及数值模型. 力学学报. 2024(06): 1740-1751 . 
百度学术
2. 陈正东，张启灿，吴周杰. 基于条纹投影的复杂结构多维度信息传感技术(内封面文章·特邀). 红外与激光工程. 2024(09): 81-99 . 百度学术
3. 王风全，梁森，吴召琦，田程. 基于分子动力学模拟异氰酸酯接枝改性热塑性树脂/纤维复合材料界面结合机理研究. 化工新型材料. 2023(12): 154-159 . 百度学术
其他类型引用(0)
-
计量
- 文章访问数: 1312
- HTML全文浏览量: 463
- PDF下载量: 118
- 被引次数: 3
