复合材料带曲率C型梁热隔膜预成型试验与仿真分析.pdf

第2 期2023年6 月纤维复合材料FIBERCOMPOSITESNo.242Jun.2023复合材料带曲率C型梁热隔膜预成型试验与仿真分析亚纪轩，叶亮，丁蕊，刘传军（中建材（上海）航空技术有限公司，上海2 0 0 12 6）摘要随着复合材料自动铺带和铺丝技术的成熟，并且广泛应用于复合材料结构的生产与制造。通过与铺带或铺丝技术相结合使得复合材料热隔膜预成型技术成为制作L型、C型，Z型等细长零件的重要手段，此种方式可以减少大量的人工操作，并且能够用于实现自动化工业生产。复合材料热隔膜预成型是一个使预浸料叠层变形的复杂过程，本文对带曲率C型梁的预成型过程进行仿真分析，首先通过试验得到预浸料铺层间摩擦系数，然后对湿态预浸料的力学性能进行测试。利用测试得到参数，采用有限元方法对C型梁预成型过程中铺层间的滑移进行模拟得到了铺层变形后的褶皱缺陷。通过专用的热隔膜成型设备对C型梁预成型进行试验研究，将试验结果与仿真结果进行对比分析，验证了预成型仿真手段的可行性，为后续热隔膜预成型过程中的参数优化提供指导。关键词复合材料；热隔膜预成型；有限元仿真分析；褶皱；层间滑移Experiment and Simulation Analysis of HotDiaphragm Forming about CompositeC-Beam with CurvatureYA Jixuan,YE Liang,DING Suo,LIU Chuanjun(CNBM(Shanghai)Aviation Technology Co.,Ltd.,Shanghai 200126)ABSTRACT With the maturity of composite material automatic tape laying（A T L)a n d a u t o ma t i c f i b e r p l a c e me n t(AFP),it is widely used in the production and manufacture of composite structure.Combining with ATL or AFP method,hot diaphragm forming(HDF)has become an important means of making L,C,Z and other slender structures.This meth-od can reduce lots of manual operations and can be used for automatic industrial production.The HDF is a complex deform-ation process.In this paper,the simulation of the forming process of the C-beam with curvature is carried out.Firstly,thefriction coefficient between the plies is obtained through the test,and then the mechanical properties of the wet prepreg aretested.The interlayer slipping during the forming process is simulated to obtain the wrinkle defects after the deformation ofthe composite plies.Experimental research on C-beam preforming is carried out with special thermal diaphragm forming e-quipment,and the experiment results were compared with simulation results to verify the feasibility of the preforming simula-tion method and provide guidance for parameter optimization in the subsequent HDF process.KEYWORDS composite material;hot diaphragm forming;finite element simulation;wrinkle;interlayer slipping通讯作者：亚纪轩，硕士，工程师。研究方向为复合材料结构设计。Em a i l：y a j x c n b m a.c n2期1引言复合材料具有高刚度比，高强度比和耐腐蚀性等优异性能，在航空航天等领域得到了广泛的应用 。然而，复合材料部件的制造过程通常复杂且耗时，为应对这一挑战，热隔膜预成型（HDF）技术逐渐得以推广应用，该技术能够以高精度和高效率生产复杂形状的复合材料零件2 。HDF工艺首先加热预浸料叠层，然后通过隔膜对预浸料叠层施加压力将其成型为所需的形状3。王立冬4-6 通过试验研究了碳纤维预浸料层间滑移特性，并且建立了热隔膜预成型过程的有限元模型并与试验结果进行了对比分析。靳晓博7 对碳纤维复合材料平直C型梁的热隔膜成型过程进行试验与仿真模拟。陈超等8 对热隔膜预成型的原理进行分析，并探讨了自动铺放与热隔膜预成型组合制造工艺的工程化应用特性。元振毅9 等对制件热隔膜预成型过程中的层间滑移情况及固化后的回弹变形进行了相关的研6%纤维方向（a）预浸料湿态力学性能测试试样示意图图1预浸料湿态力学性能测试试样方案试验机设备型号为Instron5966，将样品安装在测试夹具上，如图2 所示。DIC的测试条件为相机像素：40 0 0 pixel3000pixel，视场：36.6 mm49.9mm，光源：白光光源，标距段：35mm。由于热隔膜预成型过程中预浸料叠层的温度为60，所以在加载前将保温箱的温度稳定在6 0，然后加载。DIC设备由高分辨率相机和图像处理软件组成，在试样变形过程中捕获试样的图像，跟踪样品的运动，样品的变形由DIC设备监测和记录。通过对DIC设备捕获的图像进行后处理，确定样品的应变，进一步获得预浸料的力学性能，试验流程如图3所示，试验过程中，DIC设备拍摄的应变全场图如图4所示。Q通过试验机载荷与试验件的截面积相除得到，根据公式（1）得到1，T，T L根据试验DIC应变场能够直接得到&x，8 y，x y，根据公式（2），复合材料带曲率C型梁热隔膜预成型试验与仿真分析43究。赵月青10-12 对热隔膜预成型过程中层间的摩擦以及面内剪切变形特性进行了分析。在本文中，对带曲率C型复合材料梁结构的热隔膜预成型的模拟和试验进行研究。首先通过测试获得热隔膜预成型仿真模型的相关参数。然后，采用有限元方法建模对HDF过程进行仿真模拟，并与试验结果进行对比分析，验证了热隔膜预成型仿真模型的准确性。2热隔膜预成型参数试验2.13预浸料湿态力学性能测试本文中使用的碳纤维预浸料牌号为Hexcel的M21C，使用数字图像相关（DIC）设备对湿态预浸料的力学性能进行测试。本试验通过偏轴拉伸的试验方法获得预浸料的纤维方向的模量E1，横向模量E2，剪切模量G12。参考GB/T32788.4-2016标准，将预浸料按照标准切割成所需的尺寸。纤维与加载方向夹角10，喷涂黑白漆制备散斑，如图1所示。5-2057-20可以得出8 1，8 T，r 从而得出Ei，E2，G i 2。OL=cos0 0,0r=sin QxTr=-cose sing Ox（b）预浸料湿态力学性能试验试样图2 预浸料湿态力学性能试验过程(1)44纤维复合材料2023年将试样夹持在Instron5966试验机上设置图像采集调节2 D-DIC拍顿率，采集参摄系统考图像试验机与软件同步开启一运行计算，得代入计算得到到应变场ex，Ej,E2,G12eex选定子区，设定步长，DIC计算采集并保存图片擦滑移特性。图3试验流程图eyy 1 Lagrange0.0180.0163250.014650.0129750.01130.0096250.007950.0062750.00460.0029250.00125-0.0042-0.0021-0.00377!0.00545-0.00712!-0.0088图4DIC设备拍摄应变场全场图Lr=-2cos sing 8x+2coso sine 8,+(cos-sin0):xy8L=8,cos+8,sin-x sinacos8r=（8 x-8,)s i n 2+xycos2共进行三组试验，对试验结果处理后，湿态预浸料力学性能的平均值如表1所示。表1预浸料湿态力学性能参数Ei(GPa)平均值15.562.2预浸料层间摩擦系数测试预浸料层间摩擦系数的测试目的是为了得到预浸料在热隔膜预成型过程中，预浸料层间滑移的特性。本文通过特制的试验夹具对预浸料层间的摩擦系数进行测量，试验夹具示意图如图5所示。利用弹簧压缩对预浸料施加载荷，同时利用加热片对预浸料进行加热。硅橡胶加热片通过控温箱来调节温度，巡检仪用来监测预浸料实时温度。在测试压力为0.0 9 3MPa的条件下，测试预浸料在6 0 下的摩上铺层下铺层，中间铺层-（a）测试夹具示意图图5摩擦系数测试方案将待测试预浸料按纤维0 方向裁切，尺寸分别为37 0 mm35mm、19 0 m m 35m m。以37 0 mm35mm试样为U型层，19 0 mm35mm试样为内层，按照图6 所示进行贴合；将试样夹持在模具当中，以螺栓固定试样下底部分，上端使用拉力机夹头夹紧；通过弹簧向模具施加固定压力，压力为(2)0.09MPa，如图6 所示；温度设置为6 0，以控温箱控制试样温度，巡检仪监测试样温度；待温度稳定后开始拉伸测试，测试速率为1mm/min。E2(MPa)G122.(MPa)111.76843.95加热片4铝板硅橡胶1均压板3硅橡胶保温板加强片A区域0铺层B区域0铺层(a)试验试样图6 摩擦系数测试试样与加持试验结果载荷位移曲线如图7 所示，曲线分为试样夹具加压弹簧一控温箱一巡检仪(b)试样测试橡胶加热垫均压板试样压缩弹簧(b)试样夹持2期三个阶段，第一阶段为线性阶段，载荷随着位移的变化线性增加；第二阶段为“屈服”阶段，载荷随位移变化增加的斜率逐渐减小；第三阶段为稳定阶段，载荷随着位移的变化缓慢的增加。454035-30.25201510+03型梁热隔膜预成型试验本文中选用的C型梁长度为1m，曲率半径不复合材料带曲率C型梁热隔膜预成型试验与仿真分析24Displacement(mm)图7 试验载荷位移曲线45是常数，最小的曲率半径为7 39 9.4mm，逐渐增大到10 6 0 5.1mm，如图8 所示。试验件腹板面宽度140mm，两侧缘条高度为30 mm，腹板与缘间的R角半径为7 mm。C型梁热隔膜预成型试验使用的是Fill厂家的热隔膜预成型设备，首先通过铺丝机完成扇形预浸料叠层的制作，铺层的顺序为4545/0 2/9 0/0 2 45/-45，共计18 层，每层厚度为0.18 7 mm；然后再利用激光投影对叠层进行定位孔的加工，防止叠层在预成型过程中发生整体的滑动；将叠层转移到模具上，并使用销钉将叠层上的定位孔与模具上的定位孔进行定位，再把无孔隔离膜放在叠层上，降低隔膜与叠层之间的摩擦，过程如图9 所示。最大值-10 6 0 5.1358 9 868最小值-7 39 9.40 1434图8 C型梁曲率半径测量（a）铺丝完成扇形预浸料叠层（b）预浸料叠层投影打定位孔（c）预浸料叠层安装于模具上将预浸料叠层与模具放在热隔膜设备的平台上，然后对热隔膜设备中设置抽真空的步骤以及加温程序，如图10（a）所示，首先是升温阶段，以5/m i n 的速率将隔膜上温度增加到6 5，保温40min（此时料片处热电偶监测得到温度为6 0），保温到2 5min时，先以0.0 2 KPa/s的速度慢速抽真（d）预浸料叠层上铺放隔离膜图9 C型梁热隔膜预成型试验准备过程空至17 KPa，然后以10 KPa速率快速抽真空至-93KPa，然后保持真空，保温40 min结束后，对隔膜降温，降温过程中保持真空，温度降低至40 时，真空逐渐进行释放，试验过程中隔膜的变形如图10 （b)所示。试验设备设定的程序结束后，将模具从试验台46Temp.(C)膜上温度6540Initial Temp.Vac.初始温度/真空-17-93Vac.(kPa)真空（a）热隔膜预成型设备温度与真空控制曲线图（b）热隔膜预成型过程图10 C型梁热隔膜预成型试验上推出，试验件如图11（a）所示，将隔离膜从变形后的预浸料片上拿下，C型梁外侧缘条（曲率半径大）处发现多处褶皱，如图11（b）中缘条处的纤维复合材料40min白色线标记；在C型梁内侧缘条（曲率半径小）处25min无褶皱，C型梁腹板平面上无褶皱。由于扇形预浸uruu料叠层在向下弯折的过程中，内外两侧缘条的受力u/.S状态是不一样的，内侧缘条处预浸料叠层在弯折的.02kPaxs过程中是受拉伸，外侧缘条处预浸料叠层在弯折的Time(min)-10kPa/s时间2023年过程中是受压缩，导致材料的堆积，从而形成褶皱。4C型梁热隔膜预成型仿真C型梁的预成型仿真根据试验件模型以及试验参数进行模型的建模。预浸料的本构模型采用正交各向异性材料本构，使用工程常数进行定义，使用的参数如2.1节中描述；预浸料层间的滑移模型采用库伦摩擦，摩擦系数等于载荷除以接触面积与压力的乘积，根据2.2 中数据处理得到不同滑移距离下的摩擦系数，然后作为仿真模型的输人；热隔膜采用超弹性材料本构；模具与热隔膜平台采用刚体属性。模具与平台的网格大小选用4mm，预浸料的网格大小为2 mm，如图12 所示。模具与平台采用固支约束，约束所有自由度，在预成型料片对应销钉定位位置处的节点采用固支约束，如图13所示，仿真过程中采用多载荷步加载，隔膜上的压力与升温过程与试验保证一致。（a）试验件预成型完成后变形结果（b）试样外侧缘条褶皱缺陷图11C型梁热隔膜预成型试验结果隔膜具图12 子预成型仿真模型仿真计算完成后，预浸料的变形结果图14所结果进行对比，试验件的褶皱位置并不能与仿真结示，C型梁腹板与内侧缘条处没有产生褶皱缺陷，果完全对应上，试验结果中包含8 个褶皱，仿真结C型梁外侧缘条处会生成褶皱。将仿真结果与试验果中包含12 个褶皱，仿真结果中的4个褶皱位置2期能够与试验结果对应上，如图15所示。图14仿真变形结果云图复合材料带曲率C型梁热隔膜预成型试验与仿真分析图13模型中销钉处固支约束47预浸料叠层在滑移的过程中，在缘条的边缘会产生一个斜坡，如图16 中所示。取仿真结果的一个剖面对结果进行处理，得到在内缘条处的斜坡的角度为32.0，外缘条处的斜坡角度为33.6，如图17 所示。由于试验件是由铺丝扇形板直接进行的热隔膜试验，没有切割到净尺寸，所以无法得到对应的斜坡角度的试验结果。此角度的理论结果可以通过在三维软件中作图得出，首先做平面与C型梁内型面做交线，然后将交线向厚度方向偏置3.366（0.18 7 18）m m，然后将偏置线两端修剪到与交线的长度，然后测量线端点连线与缘条直线的夹角，内缘条侧夹角为31.6，外缘条侧夹角为33.3，如图18 所示，与仿真结果的偏差分别为1.2%,0.9%。图15仿真与试验结果对比图16 叠层滑移斜坡示意图33.63图17 仿真结果斜坡角度预测的完全对应，由于预浸料叠层滑移产生的斜5结语坡，仿真与试验结果对比角度偏差在2%以内。通本文首先进行了预浸料力学性能试验与预浸料过热隔膜预成型仿真手段用于在零件设计阶段使用铺层之间的摩擦系数试验；然后将试验得到力学性热隔膜预成型方式的可制造性评估，从而对零件的能参数与摩擦系数应用于C型梁热隔膜预成型仿真可制造性提供指导是可行的。后续将对热隔膜预成模型的建模；最后将C型梁热隔膜预成型仿真结果型试验参数进行进一步的优化，探索预成型参数对与试验结果进行对比，仿真结果模拟出在C型梁外于零件预成型褶皱缺陷的影响。侧缘条处形成褶皱，但是褶皱的具体位置并不能够48纤维复合材料2023年31.63153633.331920图18 理论斜坡角度8陈超，唐中华，张帅，郭长龙，倪敏轩，刘秀，王希，孙萌参考文献1杜善义，先进复合材料与航空航天J复合材料学报2 0 0 7(01):1-12.2吴志恩复合材料热隔膜成型【J复合材料制造技术.2018.3郑广强，吴文贵自动铺带技术在航空复合材料制造领域的应用J航空制造技术，2 0 13，435（15）：40-43.4王立冬碳纤维预浸料层间滑移性能研究及在热隔膜成型中的应用【D上海交通大学，2 0 2 0.5 Lidong Wang,Xiongqi Peng.Characterization of inter-ply slippingbehaviors in hot diaphragm preforming:Experiments and modelling.Compos Part A.2019,121:28-35.6 Lidong Wang,Xiongqi Peng.Development and application of a hy-perelastic model for diaphragm considering the influence of tempera-ture.Int.J.Computational Materials Science.2019,1950010.7】靳晓博碳纤维复合材料C型梁热隔膜预成型方法研究【D.大连理工大学，2 0 2 0.萌复合材料热隔膜预成型技术的工程化应用研究J纤维复合材料，2 0 2 1，38（0 3）：32-35.9元振毅，王永军，杨凯，刘天骄，王俊彪，魏生民，赵国强，魏俊鹏：热固性树脂基复合材料热隔膜成型过程数值仿真J复合材料学报，2 0 16，33（0 7）：1339-1350.10Zhao Y,Gu Y,Zhang T,et al.Characterization of intra-ply shearbehaviors of unidirectional prepregs during hot diaphragm formingprocess J.Polymer Composites,2021,42(2):1008-1020.11Zhao Y,Zhang T,Li H,et al.Characterization of prepreg-prepregand prepreg-tool friction for unidirectional carbon fiber/epoxyprepreg during hot diaphragm forming process J.Polymer Tes-ting,2020,84:106440.12Fei C,Ping C,Yueqing Z.DEFORMATION ANALYSIS OF C-SPAR COMPOSITE STRUCTURE DURING THE PREFORMINGPROCESSING J.