网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

含预填块复合材料帽型单筋板弯曲性能研究  PDF

  • 戴征征
  • 余章杰
  • 张琪
  • 陆方舟
  • 蔡登安
  • 周光明
南京航空航天大学机械结构力学及控制国家重点实验室,南京 210016

中图分类号: TB332

最近更新:2023-02-22

DOI:10.16356/j.1005-2615.2022.01.012

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
EN
目录contents

摘要

通过有限元分析对含预填块复合材料帽型单筋板弯曲承载能力和失效机制进行了预测,并完成试验验证。首先基于有限元软件ABAQUS,建立了含预填块帽型单筋板实体模型,完成了对筋条和蒙皮复合材料铺层及粘聚区界面的模拟。同时,本文根据Hashin准则和Tsai‑Wu准则编写了UMAT子程序,完成了对模型的有限元仿真计算,并比较了两种准则的优劣。结果表明:有限元结果中模型的承载能力、损伤机制和失效位置和试验结果吻合,载荷‑位移曲线大致一致。达到极限载荷时,共固化成型的含预填块复合材料帽型单筋板在筋条的短切纤维块与泡沫交界处发生“压断”破坏,导致结构承载能力下降,试验和仿真均未出现筋条与蒙皮脱粘现象。并且,增加筋条铺层层数和短切纤维块长度可有效提高结构承载能力。

复合材料以其卓越的力学性能及减重优势在飞机结构中得到广泛应

1。复合材料加筋板是一种典型的复合材料整体化结构,相比于单一的复合材料层合板,加筋板通过在层合板表面添加不同类型的筋条,可使结构的极限承载能力大幅度增2‑3。加强筋条的形式有很多种,如T形、L形、I形及帽4‑5,其中帽型加筋结构因其具有天然的损伤容限性能高等优良特点,而被大量应用于机身、后压力框等结构6。帽型筋条的切面尺寸较大,其两腹板通过凸缘与蒙皮相连形成一个闭合剖面,具有很高的受压稳定性,可以承受重载,能够满足多种载荷传递路径要7。但是,在服役过程中,机身桶段帽型加筋结构几乎一直承受复杂弯曲载8,相比于其轴向强度、刚度,帽型筋板结构弯曲性能较差,帽板和底板连接处易发生脱粘失9,同时为与其他结构相连接,帽型筋板在使用环境中经常需要在沿筋条打入销钉等连接固定件,导致筋条和蒙皮在此处更容易出现损伤。这不仅限制复合材料的力学优势的发挥,也对飞机安全构成极大的威10‑11

针对二次胶接复合材料帽型筋板结构,孙晶

12通过四点弯曲试验模拟了帽型筋板受后屈曲载荷时的蒙皮/筋条界面性能,建立了渐进损伤模型,分析了蒙皮/筋条界面的脱粘机理和失效过程;Turon13基于内聚力模型研究了复合材料加筋板脱粘失效问题,将有限元软件模拟结果与试验结果进行比较,两者结果相近;王力立14对复合材料帽型筋结构四点弯曲试验进行了有限元模拟,并对复材板损伤过程进行了模拟,筋条破坏方式均为加载点的压塌破坏。以上帽型筋结构数值仿真研究中,筋条和蒙皮的有限元模型均简化为壳单元模型,对复材板的初始失效判定均采用ABAQUS中内置的二维Hashin准则,模型与实际结构有所差别。

针对不同成型工艺对帽型筋板结构拉、压承载能力的影响,国内外也有许多研究。Kim

15对不同工艺成型帽型加筋板进行拉拔试验,对比分析了共固化、共胶接和二次胶接成型加筋板工艺稳定性以及拉拔高度对承载能力的影响;徐荣章16通过试验研究发现筋条与蒙皮共固化成型的三筋条复合材料帽型筋板承受轴压载荷时,结构破坏主要是由于蒙皮压缩分层断裂破坏而非筋条与蒙皮分层破坏导致。

总体而言,针对复合材料帽型加筋结构,现有的研究大多集中于筋条与复材板二次胶接成型的帽型筋,研究重点大多为筋条与底板的界面粘接强度问题。现有研究表明,二次胶接成型的帽型筋结构连接界面强度远远低于复合材料面内强度,特别是加筋结构承受弯曲、轴压载荷时,界面脱粘往往导致二次胶接帽型筋结构提前破

17。但截至目前,针对筋条与底板共固化成型复合材料帽型筋板结构弯曲破坏形式仍缺乏研究。

本文研究对象为含预填块复合材料帽型单筋板的抗弯性能,帽型筋板的筋条与蒙皮采用共固化一体成型,以避免界面脱胶导致的结构承载能力下降;同时通过在结构承压部位填充短切纤维预填块(即实际使用环境中连接销钉部位等有集中载荷作用的地方),避免承载点过早被压塌导致结构承载能力下降。本文首先建立含预填块复合材料帽型单筋板三维实体单元模型,并根据Tsai‑Wu准则、Hashin准则及Camanho线性退化准则,编写相关渐进损伤子程序。通过数值模拟,对含预填块复合材料帽型单筋板弯曲破坏失效模式及损伤机制进行预测,并比较了两种损伤判据的对计算结果的影响。同时本文还研究了预填块长度对复合材料帽型单筋板承载能力的影响,为结构优化及实际应用提供参考。

1 试验与方法

1.1 分析对象参数

本文所研究的复合材料帽型筋板由筋条、蒙皮及预填块3部分在热压罐中共固化成

18,蒙皮面板长420 mm,宽100 mm。共固化前,预填块四周提前裹有1层J272‑A胶膜,预填块由短切纤维块和泡沫组成,中间部分短切纤维块长度为100 mm,两端的泡沫长度均160 mm,试件下方粘接有加强片。帽型筋板试验件包含两种规格,分别编号为M1、M2,均由碳纤维预浸料CCF300/BA9916‑Ⅱ双向布逐层铺叠,固化后单层厚度为0.18 mm,蒙皮铺层数为8层,铺层方式为[(±45°)/(0°,90°)]2s。筋条为帽型,M1筋条铺层层数为8层,铺层方式为[(±45°)/(0°,90°)]2s;M2筋条铺层层数为12层,铺层方式为[(±45°)/(0°,90°)]3s,试件具体尺寸如图1所示。试件所选用双向布基本力学性能参数见表1,短切纤维和泡沫的力学性能参数分别见表23表1中,Eijij=1,2,3)为纤维方向弹性模量,XTYTZT分别为3个方向的拉伸强度,XCYCZC分别为3个方向的压缩强度,νijGijSijij=1,2,3)指3个方向的泊松比、切变模量和剪切强度。表2中,E1为短切纤维块弹性模量,G为短切纤维块切变模量,ν为泊松比,σ表示短切纤维强度极限。表3中,E2为泡沫弹性模量,ν为泊松比,σe表示泡沫屈服极限。E33G13G23ν13ν23ZTZCS13S23参考文献[10]中的数据,其余层合板、短切纤维块及泡沫力学性能参数均为同批成型的基础性能试验件试验测量所得数据。

图1  帽型筋板试件示意图

Fig.1  Schematic of hat‑stiffened panel

表1  CF3031/BA9916‑Ⅱ材料性能参数
Table 1  Material properities of C3031/BA9916‑Ⅱ
参数数值
E11E22/GPa 60
E33/GPa 7.7
ν121323 0.05
G12/GPa 6.4
G13=G23/GPa 6.4
XT/MPa 630
XC/MPa 520
YT/MPa 630
YC/MPa 520
ZT/MPa 50
ZC/MPa 155
S12/MPa 140
S13/MPa 105
S23/MPa 83
表2  短切纤维块材料性能参数
Table 2  Material properities of chopped fiber block
E1/GPaG/GPaνσ/MPa
30 5.2 0.3 180
表3  PMI泡沫材料性能参数
Table 3  Material properities of PMI foam
E2/MPaνσe/MPa
65 0.3 1.49

1.2 数值模拟方法

试验前先通过有限元分析软件ABAQUS建立含预填块复合材料帽型筋三维渐进损伤有限元模型,预测试件的承载能力和损伤情况。分析时将筋条帽板和蒙皮底板视为层合板,模型单层厚度为层合板固化后平均单层厚度0.18 mm,单层板本构关系参考经典层合板理论确

19。初始失效准则选用包含Yeh分层失效准则的三维Hashin失效准20,并与Tsai‑Wu失效判21的仿真计算结果进行对比,三维Hashin失效判据可区别经向拉伸失效、经向压缩失效、纬向拉伸失效、纬向压缩失效、拉伸分层失效和压缩分层失效6种失效模式。

经向拉伸失效(σ110

σ11XT2+τ12S122+τ13S132=1 (1)

经向压缩失效(σ11<0

σ11XC2=1 (2)

纬向拉伸失效(σ22+σ330

σ22+σ33YT2+1S232τ232-σ22σ33+τ12S122+τ13S132=1 (3)

纬向压缩失效(σ22+σ33<0

σ22+σ332S232+σ22+σ33YCYC2S232-1+1S232τ232-σ22σ33+τ12S122+τ13S132=1 (4)

拉伸分层失效(σ330

τ23S232+τ13S132+σ33ZT2=1 (5)

压缩分层失效(σ33<0

τ23S232+τ13S132=1 (6)

Tsai‑Wu准则是以张量形式提出的强度理论,通过增加方程中项数的方法以提高强度准则与试验结果的一致

22,但不能区别具体损伤形式,Tsai‑Wu准则认为应力空间中的破坏面存在下列形式

Fiσi+Fijσiσj=1      i,j=1,2,,6 (7)

式中:FiFij分别为一阶和二阶强度系数张量,对于复合材料单层板,可以通过破坏强度确定FiFij的取值,其中F12F13F23计算时系数取-1,即

F1=1XT-1XC,F2=1YT-1YCF3=1ZT-1ZC,F11=1XTXCF22=1YTYC,F33=1ZTZCF44=1S122,F55=1S232F66=1S132,F12=-(F11×F22)12F23=-(F22×F33)12,F13=-(F11×F33)12 (8)

复合材料层合板发生失效后损伤单元的材料性能发生退化,工程弹性常数通常要做刚度折减处理,对于材料性能退化方法,已有大量的研究。本文选择Camanho线性退化准

23,具体退化参数见表4

表4  层合板刚度退化方式
Table 4  Degradation mode of laminate
Failure modeStiffness degradation
Tensile fracture of fiber E11'=0.07E11
Compression fracture of fiber E11'=0.14E11
Tensile cracking of matrix

E22'=02E22G12'=0.2G12

G23'=0.2G23

Compression fracture of matrix

E22'=04E22G12'=0.4G12

G23'=0.4G23

Delamination E33'=G23'=G13'=ν23'=ν13'=0

本文根据本构方程、三维Hashin准则、Tsai‑Wu准则和完全退化模型编写UMAT子程序研究含预填块复合材料帽型加筋板渐进损伤过程及失效机制。

采用0厚度Cohesive单元模拟J272‑A胶膜,研究预填块与复材板的粘接界面胶膜损伤情况。具体通过ABAQUS软件内置的Traction‑separation法则定义Cohesive单元的属性,通过二次应力准则判断Cohesive单元损伤,采用基于能量的线性扩展准则来控制界面刚度的退化方式和退化速

24‑25,胶膜性能参数见表5

表5  J272‑A胶膜材料属性
Table 5  Material properities of J272‑A

E/

(MPa⋅mm-1

G1/

(MPa⋅mm-1

G2/

(MPa⋅mm-1

tn0/

MPa

ts0/

MPa

tt0/

MPa

Gc/

(J⋅m-2

3 000 1 500 1 500 15 20 20 1.8

预填块中的短切纤维视作各向同性材料,采用Mises应力准则进行损伤判定,泡沫定义为塑性材料,采用ABAQUS中内置的可压溃泡沫模型建模。

按照实际参数通过Catia建立帽型筋三维实体模型,如图2所示。实体模型装配完成后导入Hypermesh中,按实际铺层层数划分实体网并导入ABAQUS中进行材料属性赋予、装配及单元类型选择等,蒙皮、筋条和预填块的单元类型均采用8节点、6面体减缩积分单元,蒙皮和筋条接触部分单元采用共节点处理,蒙皮与预填块、筋条与预填块均设置0厚度粘聚力接触。

图2  帽型单筋板三维实体模型

Fig.2  Three-dimensional solid model of hat-stiffened panel

图3所示,建立与试验环境等直径的压头和支撑,对支撑施加位移约束,同时建立参考点与压头耦合,通过参考点施加位移载荷。

图3  有限元模型加载示意图

Fig.3  Loading diagram of finite element model

1.3 试验方法

三点弯曲试验在WDW‑E200D微机控制电子万能试验机上进行,该试验机加载精度为0.5%,试验的固定加载速率为1 mm/min。试验夹具跨距为350 mm,压头直径为30 mm,加载点在试验件中心,支撑圆柱直径为30 mm,试验加载图如图4所示。

图4  帽形筋板三点弯曲试验图

Fig.4  Three-point bending test diagram of hat-stiffened panel

为了研究含预填块复合材料帽型加筋板的变形行为,需要在特定位置粘贴应变片(花)获取应变数据并分析,尤其重点关注筋条和蒙皮中泡沫和短切纤维分界面变形情况,在这些位置增设应变片(花),应变片(花)布置位置如图5所示,应变数据采集通过JM3813应变测量仪,共计48个通道。

图5  帽型筋板应变片(花)位置示意图

Fig.5  Schematic of strain gauge locations of hat-stiffened panel

正式试验前进行预试,检查试验加载系统是否正常,并对试验机位移、载荷传感器进行标定。试验过程中每1 kN进行一次应变数据采集,并及时记录试验声响,辅助分析试验件结构损伤,试验后对断面部位进行重点观察,分析结构破坏机理,每种规格的试验件包含3件,分别编号为M1‑1、M1‑2、M1‑3和M2‑1、M2‑2、M2‑3。

2 结果与讨论

模型M1和M2的有限元仿真预测载荷‑位移曲线及试验载荷位移曲线如图6所示。Hashin准则和Tsai‑Wu准则计算得到的极限载荷见表6,试验结果得到的极限载荷见表7。有限元仿真中两种损伤准则的载荷‑位移曲线趋势大致相同,两种损伤准则给出了极限载荷的所在区间。其中Tsai‑Wu准则作为宏观失效准则,对模型的损伤预测更为保守,给出了偏安全的极限载荷,且计算量较小,工程应用上价值更高。Hashin准则给出了具体损伤形式,便于针对损伤形式对帽型筋板结构设计和复材铺层进行改进优化,以提高承载能力。

图6  复合材料帽型单筋板有限元计算载荷-位移曲线

Fig.6  Load-displacement curves of composite hat-stiffened panel by finite element calculation

表6  复合材料帽型筋板有限元仿真结果
Table 6  Finite element simulation results of composite hat‑stiffened panel
ModelDamage criterionUltimate load /kN
M1 Hashin 15.15
Tsai‑Wu 13.76
M2 Hashin 24.13
Tsai‑Wu 22.82
表7  帽型筋板三点弯曲试验结果
Table 7  Three‑point bending test results of hat‑stiffened panel
Specimen

Ultimate

load/kN

Average

load/kN

Coefficient of

variation/%

M1‑1 13.62 14.60 5.8
M1‑2 15.18
M1‑3 15.00
M2‑1 23.39 23.03 3.4
M2‑2 22.14
M2‑3 23.57

试验件M1‑1和M2‑1均加载至第1次掉载即停止加载,留待切剖进一步观察试验件内部损伤破坏情况。试验结果载荷‑位移曲线重合度较高,极限载荷离散系数较小,说明了试验具备较好的可重复性。试验件M1‑1、M2‑2极限载荷分别为13.62 kN,略小于同组其他试验件,可能是试验件加工误差导致。

同时,对比仿真和试验的载荷位移曲线和极限载荷,有限元仿真所预测的结构承载能力和结构整体刚度与试验结果大致相同,说明所建立的有限元模型的准确性。

2.1 有限元结果

图7所示,以模型M1的Hashin准则作为损伤判据的计算结果为例,重点关注达到极限载荷前的载荷‑位移曲线。

图7  M1模型Hashin准则有限元计算载荷-位移曲线

Fig.7  Load-displacement curve of M1 model (Hashin criterion)

该部分的载荷位移曲线可以分为2个阶段:第1阶段,载荷随位移呈线性关系,随着位移增加,载荷迅速升高;第2阶段,曲线斜率出现减小,载荷仍随位移增加但增速变缓,当位移达到5.1 mm,载荷增至极限载荷15.15 kN后,载荷迅速降低。

结合载荷位移曲线、损伤变量状态及应力应变云图,分析含预填块复合材料帽型单筋板的损伤过程和失效机制,第1阶段1方向应力S11图8所示,图中F指模型承载大小,U指压头向下位移量。由于短切纤维与泡沫刚度不匹配,应力主要集中在筋条短切纤维与泡沫交界处,泡沫与短切纤维应力水平较低,泡沫第1阶段最大Mises应力为0.69 MPa,短切纤维第1阶段最大Mises应力为120 MPa左右,此时结构完整,筋条、蒙皮与短切纤维均未出现损伤,各个部位共同承载。

图8  模型M1第1阶段应力云图

Fig.8  Stress contour of model M1 in the first stage

第2阶段分为两个过程:第1个过程对应载荷‑位移曲线斜率降低,载荷仍缓慢增加到极限载荷,对应的筋条和胶层损伤演化过程及对应的应力云图如图910所示。当位移达到4.1 mm,载荷达到13.49 kN,由于面内剪切应力S12、层间剪切应力S13S23的共同作用,筋条的±45°铺层在短切纤维与泡沫分界处开始出现出现分层失效,同时此处的胶膜开始出现损伤并向外扩展,导致结构刚度下降,载荷‑位移曲线斜率降低。其中SDV1~6分别对应经向拉伸失效、经向压缩失效、纬向拉伸失效、纬向压缩失效、拉伸分层失效和压缩分层失效6种失效模式。CSDMG为内聚力模型损伤变量,其值不为0时表示胶膜发生损伤;其值为1表示胶膜完全损伤。

图9  第2阶段第1过程损伤演化过程和对应应力云图

Fig.9  Damage evolution process and corresponding stress contour of the first process in the second stage

图10  胶膜损伤扩展云图

Fig.10  Contour of damage evolution for adhesive layer

第2个过程对应载荷‑位移曲线达到极限载荷并迅速掉载,该过程主要出现筋条的经向压缩失效和纬向压缩失效,对应的损伤演化过程如图11所示,筋条的应力应变云图如图12所示。当位移达到4.7 mm,载荷达到14.72 kN,筋条(0°/90°)层泡沫和短切纤维分界处开始出现经向压缩失效并迅速向周边扩展;当位移达到5.1 mm,载荷达到极限载荷15.15 kN,筋条±45°层泡沫和短切纤维分界处出现纤维压缩破坏和基体压缩破坏,载荷迅速降低。与第1个过程分层失效相比,经向压缩失效与纬向压缩失效扩展更为迅速,且出现失效后模型载荷‑位移曲线很快到达载荷最大值,说明纤维压缩对筋条结构破坏起到了决定性的作用。

图11  筋条第2阶段第2过程损伤演化过程

Fig.11  Damage evolution process of the second progresss in the second stage of the rib

图12  筋条第2阶段第2过程应力云图

Fig.12  Stress contour of the second process in the second stage of the rib

筋条“压断”第1次掉载后,随着加载位移继续增加,结构仍能继续承载,蒙皮成为主要承载结构,此过程蒙皮的应力云图如图13所示;当位移增至15.5 mm,载荷达到8.95 kN,出现第2次掉载,蒙皮出现纤维拉伸断裂,结构完全破坏,计算终止,具体损伤情况及应力云图如图14所示。

图13  蒙皮第3阶段应力云图

Fig.13  Stress contour of the skin in the third stage

图14  蒙皮破坏位置及对应应力云图

Fig.14  Damage location and corresponding stress cloud of skin

仿真全过程短切纤维预填块仅在与泡沫分界处出现细微损伤,承载状况良好,第4阶段短切纤维预填块损伤情况及应力云图如图15所示。

图15  短切纤维块损伤情况及应力云图

Fig.15  Damage situation and stress contour of chopped fiber block

2.2 试验结果

以M1‑3为例,分析由载荷‑位移曲线可以看出:起始阶段,载荷随位移呈线性增长,在达到极限载荷前未出现任何波动,试验件未发生损伤。当位移达到4.95 mm,载荷为15 kN时,筋条在靠近短切纤维与泡沫交界处发生断裂破坏,载荷剧烈下降至4.25 kN。但观察试验件,筋条与蒙皮未出现分层现象,蒙皮在此时未出现明显损伤,可以继续承载。继续下压压头,载荷在4 kN附近保持稳定,当位移达到13.95 mm,载荷出现突降,蒙皮短切纤维与泡沫交界处断裂破坏,筋条与蒙皮仍未出现分层现象。

利用高速切割机,沿短切纤维与泡沫分界处的截面对试验件进行切开,通过KH‑7700数字显微镜对切割面内部损伤情况进行观察,放大倍数为5倍。如图16所示,筋条靠内铺层在短切纤维与泡沫分界处出现分层损伤,胶膜出现局部脱粘现象,远离短切纤维与泡沫交界处的切面则未发现筋分层及泡沫脱粘现象,试验现象与仿真结果高度一致。

图16  帽型筋板试验件损伤位置剖面图

Fig.16  Sectional view of damage positions of hat-stiffened panel

试验件M1‑3损伤破坏情况如图17所示,试验件破坏形式与有限元仿真预测结果吻合。

图17  帽型筋板损伤破坏示意图

Fig.17  Schematic of damage positions of hat-stiffened panel

2.3 分析与拓展

由试验所得的应变‑载荷曲线可以看出,对应位置的应变数据大致保持一致,说明压头是加载在试验件中心对称面上的,试验加载和设计是合理的。M1试验件试验所得的极限载荷平均值为14.60 kN,对应的Hashin损伤判据仿真计算极限载荷为15.15 kN,与试验结果仅差3.6%;M2试验件试验所得的极限载荷平均值为23.03 kN,对应的Hashin损伤判据仿真计算极限载荷为24.13 kN,与试验结果仅差4.8%,误差均在5%以内,试验和有限元仿真的载荷‑位移曲线趋势也大致保持一致。

对比有限元模型和试验件损伤破坏的位置,两者结构损伤破坏的位置具有高度一致性:达到极限载荷后第1次掉载,损伤破坏的位置出现在筋条短切纤维与泡沫交界处;第2次掉载,是由于蒙皮在短切纤维与泡沫交界处发生破坏,说明了所建立的渐进损伤有限元模型的可靠性。

为了验证有限元计算得到的应变数据的准确性,取载荷为12 kN时,筋条和蒙皮短切纤维与泡沫分界处具有代表性应变测点的应变数据进行对比,即在有限元模型中取22、25、38、39、44、45测点的应变值与试验结果进行对比,如表8所示,39号测点应变试验与仿真差值较大,可能是该测点应变片损坏或粘接不够牢固而产生的测量误差,除去该测点外各个测点仿真和试验应变数据大致吻合,进一步说明有限元建模的可靠性。

表8  部分测点应变值与仿真值对比
Table 8  Comparison of strain values and simulated values

Strain gauge

number

Simulation

value/με

Test

value/με

Error/%
22 -1 469 -1 553 5.7
25 -1 454 -1 329 8.6
38 -6 523 -7 129 8.5
39 -6 392 -5 639 11.7
44 2 013 1 972 2.1
45 1 982 1 807 9.7

试验验证了渐进损伤有限元模型的准确性后,本文还通过有限元模型继续研究了短切纤维块长度对帽形筋承载能力的影响,不同纤维块长度的有限元模型及数值模拟破坏载荷如表9所示,极限载荷‑短切纤维块长度曲线见图18

表9  不同纤维块长度模型参数及极限载荷
Table 9  Model parameters and ultimate loads of different models
ModelM1‑0M1‑60M1‑100M1‑120M1‑140M1‑400
Length of chopped fiber block/mm 0 60 100 120 140 400
Ultimate load /kN 3.59 11.70 15.15 16.95 18.45 24.20

图18  极限载荷与短切纤维块长度关系曲线

Fig.18  Relationship between the ulimit load and the length of chopped fiber block

通过分析发现,短切纤维块长度增加可以显著提高帽形筋结构的承载能力。对于不含短切纤维块的帽型筋结构,承载能力较弱,极限载荷仅为3.59 kN,破坏方式为筋条在加载点处发生压塌破坏;当短切纤维长度为60 mm时,极限载荷提高至11.7 kN,结构整体刚度和承载能力均得到很大提高,极限载荷较不含短切纤维预填块模型提高了226%,同时,破坏位置也发生改变,含短切纤维预填块块的模型在筋条短切纤维与泡沫交界处“压断”。由图18可看出,结构弯曲承载能力随着短切纤维块长度增加,呈不断增加的趋势,但增速逐渐降低。

3 结 论

(1)本文建立了含预填块复合材料帽型单筋板实体模型,基于ABAQUS软件,根据复合材料单层板Hashin、Tsai‑Wu失效判据,编写UMAT子程序,结合粘聚区模型进行有限元分析。对含预填块复合材料帽型单筋板弯曲破坏形式、承载能力进行预测,得到含预填块复合材料帽型单筋板载荷‑位移曲线、弯曲失效渐进损伤机理及相应的应力云图,同时比较了两种失效判据对计算结果的影响。其中Tsai‑Wu准则作为宏观失效准则,对模型的损伤预测更为保守,给出了偏安全的极限载荷,且计算量较小,工程应用上价值更高。Hashin准则给出了具体损伤形式,便于针对损伤形式对帽型筋板结构设计和复材铺层进行改进优化,以提高承载能力。

(2)通过含预填块复合材料帽型单筋板三点弯曲试验,对有限元仿真结果进行了验证。试验结果显示:极限载荷试验值和仿真值相差比较小,破坏形式吻合较好,载荷‑位移曲线大致一致,验证了有限元仿真分析的准确性。仿真与试验结果共同表明:共固化成型的含预填块复合材料帽型单筋板弯曲破坏形式并非筋条与蒙皮脱粘或分层,而是筋条在靠近短切纤维与泡沫交界处发生“压断”破坏。

(3)分析M1与M2的试验结果和仿真结果,容易发现筋条铺层层数对含预填块复合材料帽型单筋板弯曲承载能力存在较大影响,M2比M1多4层铺层,极限载荷提升53.2%,提高筋条铺层层数可显著提高结构承载能力。

(4)本文根据有限元仿真对结构中短切纤维块长度进行了参数化研究,仿真结果显示:短切纤维块的存在改变了复合材料帽型单筋板弯曲破坏位置和形式,能较大的提高帽型筋结构承载能力,且随着短切纤维块长度增加,其极限载荷呈不断增加的趋势,但增速逐渐降低。实际应用中,可根据需要在帽型筋板压力集中部位填充长度合适的短切纤维预填块(即实际使用环境中连接销钉部位等有集中载荷作用的地方),避免承载点过早被压塌导致结构失效或承载能力下降。

参考文献

1

徐瑀童. 复合材料层合板低速冲击与剩余强度研究[D]. 南京南京航空航天大学2018. [百度学术] 

XU Yutong. Research on low-speed impact and residual strength of composite laminates[D]. NanjingNanjing University of Aeronautics and Astronautics2018. [百度学术] 

2

杨乃宾章怡宁.复合材料飞机结构设计[M]. 北京北京航空工业出版社20029-10. [百度学术] 

Yang NaibinZhang Yining. Composite aircraft structure design[M]. BeijingBeijing Aviation Industry Press20029-10. [百度学术] 

3

刘莹. 复合材料加筋壁板结构的承载能力研究[D]. 沈阳东北大学2008. [百度学术] 

Liu Ying. Study on the bearing capacity of composite reinforced siding structure[D]. ShenyangNortheastern University2008. [百度学术] 

4

徐建印春伟童明波.复合材料T型加筋结构连接界面性能研究[J].江苏航空2012S1):111-113. [百度学术] 

Xu JianYin ChunweiTong Mingbo. Study on the interface performance of composite T-shaped reinforced structure[J].Jiangsu Aviation2012S1):111-113. [百度学术] 

5

王震鸣.复合材料加筋板壳的承载能力问题[J].复合材料学报19841): 40-48. [百度学术] 

Wang Zhenming. Bearing capacity of composite stiffened plates and shells[J]. Acta Materiae Compositae Sinica19841): 40-48. [百度学术] 

6

沈真张晓晶.复合材料飞机结构强度设计与验证概论[M].上海上海交通大学出版社2011. [百度学术] 

Shen ZhenZhang Xiaojing. Introduction to strength design and verification of composite aircraft structure[M]. ShanghaiShanghai Jiaotong University Press2011. [百度学术] 

7

Ambur D R. Design and evaluation of a foam-filled hat-stiffened panel concept for aircraft primary structure applications: NASA Technical Memorandum 109175[R].[S.l.]: NASA1995. [百度学术] 

8

李吻李勇.Z-pin增强复合材料帽型单加筋板弯曲性能[J].航空学报20163712): 3843-3852. [百度学术] 

Li WenLi Yong. Bending behavior of Z-pin reinforced composite hat-shaped single stiffened plate[J]. Acta Aeronautica Sinica20163712): 3843-3852. [百度学术] 

9

孙晶晶张晓晶. 复合材料帽型加筋结构界面脱粘的试验和数值分析[C]//第17届全国复合材料学术会议. [S.l.]中国航空学会201285-91. [百度学术] 

Sun JingjingZhang Xiaojing. Experimental and numerical analysis of interfacial debonding of composite hat-shaped reinforced structures[C]//Proceedings of the 17th National Conference on Composite Materials. [S.l.]Chinese Aeronautical Society201285-91. [百度学术] 

10

郭丽君陆方舟周光明.碳纤维/环氧树脂复合材料缠绕接头拉伸失效机制[J].复合材料学报2020379): 1-11. [百度学术] 

Guo LijunLu FangzhouZhou Guangming. Tensile failure mechanism of carbon fiber/epoxy resin composite winding joints[J]. Acta Materiae Compositae Sinica2020379): 1-11. [百度学术] 

11

孙唯. 复合材料帽型加筋结构的稳定性分析[D]. 哈尔滨哈尔滨工程大学2015. [百度学术] 

Sun Wei. Stability analysis of composite hat-shaped reinforced structure[D]. HarbinHarbin Engineering University2015. [百度学术] 

12

孙晶晶. 复合材料帽型筋条脱粘的失效机理分析[D]. 上海上海交通大学2013. [百度学术] 

Sun Jingjing. Failure mechanism analysis of debonding of composite cap ribs[D]. ShanghaiShanghai Jiaotong University2013. [百度学术] 

13

TURON ADÁVILA C GCAMANHO P Pet al. An engineering solution for mesh size effects in the simulation of delamination using cohesive zone models[J]. Engineering Fracture Mechanics20077410): 1665-1682. [百度学术] 

14

王力立陈宏林国伟. 复合材料帽型单筋板弯曲破坏数值分析[J].兵器装备工程学报2019404):202-208. [百度学术] 

Wang LiliChen HongLin Guoweiet al. Numerical analysis of bending failure of composite cap-shaped single-ribbed plate[J]. Journal of Ordnance Equipment Engineering2019404): 202-208. [百度学术] 

15

KIM G HCHOI J HKWEON J H. Manufacture and performance evaluation of composite hat-stiffened panel[J]. Composite Structures2010929): 2276-2284. [百度学术] 

16

徐荣章凌晓涛贾利勇.成型工艺对复合材料帽型加筋板轴压特性的影响[J].材料工程2018464):152-157. [百度学术] 

Xu RongzhangLing XiaotaoJia Liyong. The effect of molding process on the axial compression characteristics of composite hat-shaped stiffened plate[J]. Materials Engineering2018464): 152-157. [百度学术] 

17

李文晓吴文平.PMI泡沫芯材填充GFRP帽型筋梁弯曲性能研究[J].工程塑料应用2009371): 56-59. [百度学术] 

Li WenxiaoWu Wenping. Research on bending properties of GFRP cap-shaped rib beam filled with PMI foam core material[J]. Engineering Plastics Application2009371): 56-59. [百度学术] 

18

沈观林. 复合材料力学[M]. 北京清华大学出版社1996. [百度学术] 

SHEN Guanlin. Mechanics of composite materials[M]. BeijingTsinghua University Press1996. [百度学术] 

19

王耀先.复合材料力学与结构设计[M].上海华东理工大学出版社2012. [百度学术] 

WANG Yaoxian. Mechanics and structural design of composite materials[M]. ShanghaiEast China University of Science and Technology Press2012. [百度学术] 

20

HASHIN Z. Fatigue failure criteria for unidirectional fiber composites[J]. Journal of Applied Mechanics1981484): 846. [百度学术] 

21

Tsai S WWu E M. A general theory of strength for anisotropic materials[J]. Journal of Composite Materials197151): 58-80. [百度学术] 

22

WANG S SCHIM S M. Fatigue damage and degradation in random short-fiber SMC composite[J]. Journal of Composite Materials1983172): 114-134. [百度学术] 

23

HAVAR TDRECHSLER K. Design and progressive failure analysis of 3D-reinforced composite force introduction loops[C]//Proceedings of AIAA/ASME/ASCE/AHS/ASC StructuresStructural Dynamics, and Materials Conference. [S.l.]AIAA2013. [百度学术] 

24

孙立许希武徐焜.含脱粘损伤的复合材料加筋板压缩破坏渐进损伤数值分析[J].材料科学与工程学报2009275): 740-744. [百度学术] 

Sun LiXu XiwuXu Kun. Numerical analysis of progressive damage of composite stiffened plates with debonding damage[J]. Journal of Materials Science and Engineering2009275): 740-744. [百度学术] 

25

雷良超.复合材料缠绕接头力学性能研究[D]. 南京南京航空航天大学2018. [百度学术] 

LEI Liangchao. Research on mechanical properities of the winding composite lug[D]. NanjingNanjing University of Aeronautics&Astronautics2018. [百度学术] 

您是第位访问者
网站版权 © 南京航空航天大学学报
技术支持:北京勤云科技发展有限公司
请使用 Firefox、Chrome、IE10、IE11、360极速模式、搜狗极速模式、QQ极速模式等浏览器,其他浏览器不建议使用!