第1期纤维复合材料㊀No.1㊀1192024年3月FIBER ㊀COMPOSITES ㊀Mar.2024复合材料纵横剪切模量的有限元计算与试验对比研究王天琦,魏喜龙,王威力,李㊀刚,魏㊀程,单鹏宇,田忠恩(哈尔滨玻璃钢研究院有限公司,哈尔滨150028)摘㊀要㊀层合板(ʃ45ʎ/90ʎ/0ʎ)是准各向同性层合板,其刚度在板面内各个方向上都相同㊂使用计算软件,编写程序代码,运用有限元的方法,划分有限元网格,计算了层合板(ʃ45ʎ/90ʎ/0ʎ)的纵横剪切模量[1]㊂将计算结果与复合材料纵横剪切模量试验结果进行对比,得到数值模拟和试验结果之间的误差为0.89%,一致性很好,说明该计算方法可用于工程计算模拟㊂关键词㊀复合材料层合板;准各向同性;有限元方法;纵横剪切模量Research on Finite Element Method and Testing forDetermination Interlaminar Shear Modulusof Composite MaterialsWANG Tianqi,WEI Xilong,WANG Weili,LI Gang,WEI Cheng,SHAN Pengyu,TIAN Zhongen(Harbin FRP Institute Co.,Ltd.,Harbin 150028)ABSTRACT ㊀Laminated plates (ʃ45ʎ/90ʎ/0ʎ)are quasi isotropic laminates with the same stiffness in all directions within the plate ing computational software,program code was written,and finite element method was used to di-vide the finite element mesh and calculate the longitudinal and transverse shear modulus of laminated plates (ʃ45ʎ/90ʎ/0ʎ)[1].The calculated results are compared with the transverse and longtitude shear modulus test results,it was found that the error between the numerical simulation and experimental results was 0.89%,indicating good consistency,indicating that this calculation method can be used for engineering calculation simulation.KEYWORDS ㊀composite laminates;quasi -isotropic;finite element method;longitudinal and transverse shear modulus通讯作者:王天琦,硕士研究生,助理工程师㊂研究方向为结构力学㊂E -mail:151****4822@1㊀引言纤维复合材料由于其强度高㊁耐腐蚀等优点被广泛应用于航空航天和海洋等领域中㊂在复合材料的研发过程及工程应用中,经常要对不同的复合材料试件进行多种形式的性能试验,以获取所需的材料和结构性能参数㊂纤维增强复合材料的剪切性能是其基本性能之一,是复合材料设计中必须考虑的重要问题㊂许多结构铺层的设计㊁优化都是以剪切强度㊁模量数据作为依据㊂本文采用有限元数值模拟和试验对比的方法,对聚合物基复合材料层合板纵横剪切性能进行分析[2]㊂研究有限元数据模拟方法可在无试验条件下快速为复合材料结构件提供试验数据的依据㊂纤维复合材料2024年㊀2㊀准各向同性层合板的力学试验2.1㊀复合材料层合板纵横剪切试件制备按照GB /T 3355-2014要求制备准各向同性层合板,试样铺层为[45ʎ/-45ʎ]5s ,试件尺寸如图1所示㊂图1中h 表示试样厚度,为2mm [3]㊂试验中所用到的材料规格及供应厂家如表1所示㊂图1㊀纵横剪切试件示意图表1㊀试验材料名称规格供应商单向预浸料SYT49S /8485-A200(S186)/环氧树脂双酚A 环氧树脂南通星辰合成材料有限公司2.2㊀试验结果按照GB /T 3355-2014进行力学试验,试验在标准环境下进行,记录试件的宽度㊁厚度㊁试件受到的破坏载荷,并计算剪切强度和剪切模量㊂纵横剪切试件参数及试验结果如表2所示㊂表2㊀纵横剪切试件参数及试验数据汇总表编号123456平均值宽度b /mm 24.9824.9524.9625.0325.0325.0425.00厚度h /mm 2.41 2.45 2.36 2.38 2.42 2.46 2.41破环载荷/N 3017309225792703310931162936剪切模量/GPa 4.09 4.05 4.04 3.97 3.93 4.04 4.00剪切强度/MPa28.9029.6525.0825.930.0028.4928.02㊀㊀45ʎ纵横剪切强度计算如公式(1)所示㊂S 45ʎ=P 45ʎmax 2bh(1)公式(1)中,P 45ʎmax 为最大破坏载荷,N;b 为试件宽度,mm;h 为试件厚度,mm㊂45ʎ纵横剪切模量计算如公式(2)所示㊂G 12=әP 45ʎ2bh (әε2-әε1)(2)公式(2)中,әP 45ʎ为载荷-应变曲线上初始直线段的载荷增量,kN;b 为试件宽度,mm;h 为试件厚度,mm;әε1为对应于әP 45ʎ的纵向应变;әε2为对应于әP 45ʎ的横向应变[4]㊂由表2中可看出,试验数据有一定的离散,但离散程度不大㊂3㊀准各向同性层合板纵横剪切模量计算方法㊀㊀采用经典层合板理论(CLT),四节点举行层合板有限元计算方法计算变形与应变的关系,得到纵横剪切模量的计算值㊂3.1㊀经典层合板理论与均质材料所组成的结构不同,复合材料层合板结构的分析必须立足于对每一单层的分析㊂由于存在不同的组分层,决定了层合结构的厚度方向具有宏观非均质性㊂为了得到层合结构的刚度特性,必须弄清楚各单层的刚度特性;为了对层合结构的强度做出判断,必须首先对各单层的强度做出判断㊂因此,单层的宏观力学分析是层合结构分析的基础[5]㊂经典层合板理论,基于Kirchhoff 假设,即直法线假设和法线长度保持不变,Z 向应力可以忽略假设,而建立的薄层合板中面变形方程㊂当受到外部载荷作用时,层合板将发生面内伸缩或者弯曲变形㊂层合板坐标如图2所示[6]㊂直线段变形前后关系如图3所示,图3中以XZ 平面内的变形为例,B 是中面上的一点,C 是截面上的任意点,β是层合板中面的转角,如公式(3)所示㊂u c =u 0-z c sinβʈu c -z c β(3)公式(3)中,β=∂w 0∂x,注:在小变形情况下,021㊀1期复合材料纵横剪切模量的有限元计算与试验对比研究图2㊀层合板坐标图图3㊀直线段变形前后关系才有βʈsinβ,根据上式,层合板上任意一点的位移u 可以表示,如公式(4)所示㊂u =u 0-z∂w 0∂x(4)同理在YZ 平面内,也可以得到任意一点位移v的表达式如公式(5)所示㊂v =v 0-z∂w 0∂y(5)另外,经典层合板理论中,任意一点的位移w 与其中性面上的面外位移w 0相等,如公式(6)所示㊂w =w 0(x ,y )(6)公式(6)中,u ㊁v ㊁w 分别为沿板厚范围内x,y,z 方向的位移,中面上的点x,y,z 方向的位移为u 0㊁v 0㊁w 0,其中称为板的挠度[6]㊂后续计算过程,依据上述理论进行㊂3.2㊀变形与应变的关系为了得到层合板的刚度矩阵和柔度矩阵,并且在有限元方法中运用,引入层合板的应力-应变关系㊂为简化问题,对所研究的层合板作如下假设[7]㊂(1)层合板各单层之间粘结良好,可作为一个整体结构板,并且粘结层很薄,其本身不发生形变,即各单层板之间变形连续㊂(2)层合板虽由多层单层板叠合而成,但其总厚度仍符合假设,即厚度t 与跨度L 之比为(150~1100)<t L <18~110(3)整个层合板是等厚度的㊂在上述假设前提下,基于经典层合板理论中位移的表达式如公式(7)所示㊂u (x ,y ,z )=u 0(x ,y )-z ∂w 0∂xv (x ,y ,z )=v 0(x ,y )-z∂w 0∂y(7)w (x ,y ,z )=w 0(x ,y )引入几何方程,应变与位移的关系式如公式(8)所示㊂εx =∂u ∂x ,εy =∂v ∂y ,εz =∂w ∂zγyz =∂w ∂y +∂v ∂z ,γzx =∂u ∂z +∂w ∂x ,γxy =∂y ∂x +∂u∂y(8)将层合板的变形带入上式可得如公式(9)所示㊂εx =∂u ∂x =∂u 0∂x -z ∂2w 0∂x 2=εx 0+zk xεy =∂v ∂y =∂v 0∂y -z ∂2w 0∂y 2=εy 0+zk yγxy =∂u ∂y +∂v ∂x =∂u 0∂y +∂v 0∂x -2z ∂2w 0∂x ∂y=γxy 0+zk xy(9)式中,k x ,k y 为中面的弯曲挠曲率,即曲率半径的倒数,k xy 为中面的扭曲率,εx0㊁εy0㊁γxy0为中面应变[8]㊂将总应变整合后形式如下,其中第一部分代表的是面内的应变,第二部分代表的是弯曲引起的应变,如公式(10)所示[9]㊂{ε}={εo }+z {k }{εo }∂u 0∂x ∂v 0∂y ∂u 0∂y +∂v 0∂x ìîíïïïïïïïïüþýïïïïïïïï(10)121纤维复合材料2024年㊀{k }=∂2w ∂x 2∂2w ∂y 2-2∂2w ∂x ∂y ìîíïïïïïïïïüþýïïïïïïïï将公式(8)代入单层板应力-应变关系式,可以得到层合板中第k 层的应力-应变关系式[10]:σx σy σz éëêêêêùûúúúú=Q 11㊀Q 12㊀Q 16Q 12㊀Q 22㊀Q 26Q 11㊀Q 26㊀Q 66éëêêêêùûúúúúεxo εyo γxyo ìîíïïïïüþýïïïï+z k x k y zk xy ìîíïïïïüþýïïïï(11)式(11)中Q 为刚度矩阵,由公式(9)可知,层合板应变由中缅应变和弯曲应变两部分组成,沿厚度线性分布;而应力除与应变有关外,还与各单层刚度特性有关,若各层刚度不相同,则各层应力不连续分布,但在每一层内是线性分布的[11]㊂3.3㊀四节点矩形层合板有限元计算纵横剪切模量有限元法的基本理念是假设将一个连续体分割成数目有限的小体(单元),彼此间只在数目有限的指定点(节点)处互相连结,组成一个单元的集合体以代替原来的连续体,再在节点上引进等效力以代替实际作用于单元上的外力㊂因此,可以选择一个简单的函数来近似地表示位移分量的分布规律,建立位移和节点力之间的关系,把有无限个自由度的连续体理想化为有限个自由度的单元集合体,使问题简化为适合于数值解法的结构型问题[12]㊂基于复合材料层合板理论,结合3.1及3.2中的理论计算结果,建立并划分的层合板有限元网格如图4所示,在节点1,2,3上施加了边界条件-即将所有自由度固定㊂图4㊀层合板试件有限元数值模拟图在节点13和15上,沿x 轴施加25N 的力,在节点14上,沿x 轴施加50N 的力,施加载荷为100N㊂层合板试件厚度如公式(12)所示㊂H =n ㊃δ=10㊃0.24=2.4mm(12)式中,δ为单层厚度0.24mm㊂将节点8和节点11之间沿X 轴的相对伸长量定义为ε1,沿Y 轴的相对伸长量定义为ε2㊂则沿X 轴的剪切强度和相对伸长量如公式(13)~(16)所示㊂τ=F2bH=0.83333MPa (13)ε1=u (11)-u (8)a =0.0001274(14)ε2=u (11)-u (8)a=0.000018(15)γ=ε1-ε2(16)即剪切模量如公式(17)所示㊂G 12ᶄ=әτәγ=4035.61MPa (17)试验结果与有限元数值模拟结果之间的误差如公式(18)所示㊂ψ=G 12ᶄ-G 12G 12=0.89%(18)数值模拟中用到的参数如表3所示㊂表3㊀纵横剪切数值模拟结果汇总表数值模拟参数结果数值两个剪应变点之间的剪应力差值/MPa 0.833333两个剪应变点之间的剪应变差值/(rad /rad)0.0002纵向应变算数平均值/(mm /mm)0.0001274环向应变算数平均值/(mm /mm)0.000018试验结果剪切模量/MPa 4000有限元数值模拟剪切模量/MPa4035.61误差0.89%由表2中数据可知,试验结果与数值模拟得出的层合板剪切模量相差不大,误差仅在1%以内,存在误差的主要原因是模拟使用的摩擦参数与试验实际情况有差别,并且试验使用的板料厚度也存在一定的偏差,这些因素的综合影响导致了有限元模拟与试验值之间的偏差[13]㊂4㊀结语基于经典层合板理论,运用计算软件并使用有限元分析方法对复合材料层合板的纵横剪切试验进行数值模拟,提出了一个含有前后处理的完整的四221㊀1期复合材料纵横剪切模量的有限元计算与试验对比研究节点有限元程序算法,该程序算法可以为模拟层合板的拉伸㊁压缩等力学性能提供计算基础[14]㊂本文研究分析得到了层合板的纵横剪切模量,根据数值模拟后的数据显示,有限元方法在一定程度上可以得到比较准确的预测复合材料层合板的力学性能,能够指导复合材料层合板铺层结构[15],有限元数据和试验数据之间的误差在1%左右,使用有限元数值模拟方法计算剪切模量的方法是可行的,可以满足工程需要㊂参考文献[1]杨涵,周仕刚,薛元德.(ʃ45ʎ/90ʎ/0ʎ)s准各向同性层合板的拉伸强度[J].纤维复合材料,2010(4):5.DOI:10. 3969/j.issn.1003-6423.2010.04.001.[2]鲍宏琛,刘广彦.准各向同性纤维增强复合材料层合板的开孔拉伸破坏模拟[J].复合材料学报,2016,33(5):7.DOI: 10.13801/ki.fhclxb.20160112.005.[3]刘宝良,夏军.基于Mindlin理论的复合材料层合板的有限元分析[J].黑龙江科技大学学报,2004,14(3):171-173. DOI:10.3969/j.issn.1671-0118.2004.03.010.[4]朝鲁,张鸿庆,唐立民.一个计算微分方程(组)对称群的Mathematica程序包及其应用[J].计算物理,1997,14 (003):375-379.DOI:10.1007/BF02947208.[5]邵明.有限元法在复合材料层合板力学性能中的应用研究[D].浙江大学,2013.[6]李真,陈秀华,汪海.基于重合网格法的含孔复合材料层合板有限元分析[C]//第17届全国复合材料学术会议.中国力学学会;中国宇航学会;中国航空学会;中国复合材料学会, 2012.DOI:ConferenceArticle/5af1b392c095d71bc8cb4150.[7]李斯华李良伟.复合材料层合板的层间应力有限元分析[J].四川建材,2013,039(001):63,66.[8]ЖилинП.А.Векторыитензорывторогорангавтрехмерномпространстве.СПб.:Нестор,2001.276с.[9]ЖилинП.А.Рациональнаямеханикасплошныхсред:учеб.пособие.СПб.:Изд-воПолитехн.ун-та,2012.584с.[10]СорокинФ.Д.Прямоетензорноепредставлениеуравненийбольшихперемещенийгибкогостержнясиспользованиемвектораконечногоповорота//ИзвестияРАН.Механикатвердоготела.1994.ɴ1.C.164-168.[11]КувыркинГН,ГоловинНН,ПетрикевичББ,et al.Разработкатеоретическихосновиметодовоценкинесущейспособностиипрочностнойнадежности(рискаразрушения)термонапряженныхэлементовконструкцийизкомпозитныхматериаловнаосновеуглерода[R].Российскийфондфундаментальныхисследований,1994.[12]ДемидовАС,ХомовскийЯН.Особенностирасчётаивыборкритериевнапряженно-деформированногосостояниядлятонкостенныхконструктивныхэлементовизуглерод-углеродныхкомпозитныхматериалов[J].Sciences of Europe, 2016(9-4(9)):80-84.[13]王言磊,郝庆多,欧进萍.复合材料层合板面内剪切实验方法的评价[J].玻璃钢/复合材料,2007(03):6-8.DOI:10. 3969/j.issn.1003-0999.2007.03.002.[14]张纪奎,郦正能,关志东,等.复合材料层合板固化压实过程有限元数值模拟及影响因素分析[J].复合材料学报,2007, 24(2):6.DOI:10.3321/j.issn:1000-3851.2007.02.022.[15]许良梁,矫桂琼,卢智先,等.含分层复合材料层合板剪切屈曲的实验研究[J].机械强度,2007,29(4):5.DOI:10. 3321/j.issn:1001-9669.2007.04.011.321。
