Modified silicone alkyne hybrid resin for RTM and its composite
-
摘要: 树脂传递模塑(RTM)成型工艺是一种液体闭模成型工艺,成品表面光滑、安全环保、成本较低,极具发展潜力。硅炔杂化树脂是一种耐高温性能优异的有机无机杂化树脂,在航空航天领域均有广泛应用,但其与纤维的结合性能较差,复合材料力学性能不佳,成为制约其发展应用的主要因素。本文以苯并噁嗪树脂(PBZ)和氨基稀释剂(PAD)共混改性聚(间二乙炔基苯-甲基氢硅烷)树脂(PSA),制备了适用于RTM成型工艺的改性硅炔杂化树脂(PPA)。通过DSC、FTIR、流变仪、旋转黏度计和热重等方法对PPA树脂的RTM成型工艺、固化行为及耐热性能进行分析。结果表明:PPA树脂的加工窗口宽,可以实现RTM成型的目标,但固化温度高于PSA树脂;PPA树脂的耐热性能优异,其中PPA-1 (PSA∶PBZ∶PAD质量比为5∶1∶2)树脂在氮气和空气中质量损失5%的温度(Td5)分别为585.4℃和568.3℃,1000℃质量保留率分别为88.3% 和26.0%。复合材料力学性能测试表明,石英纤维增强PPA树脂复合材料(QF/PPA)力学性能随PBZ含量的增加而逐渐升高,其中QF/PPA-1常温下弯曲强度为346.2 MPa,层间剪切强度为21.4 MPa,较QF/PSA复合材料分别提高了120.6%和72.6%,400℃热老化2 h后的弯曲强度和层间剪切强度分别为256.5 MPa和17.1 MPa,400℃热老化2 h后的力学性能保留率超过70%。Abstract: Resin transfer molding (RTM) forming process is a liquid closed mold forming process, which has a smooth surface of the finished product, safety and environmentally friendly, low cost, and has great development potential. Silicone alkyne hybrid resin is an organic-inorganic hybrid resin with excellent high-temperature resistance, widely used in aerospace field, but its poor bonding performance with fibers and poor mechanical properties of composites have become the main factors restricting its development and application. Poly(m-diacetylbenzene-methylhydrosilane) resin (PSA) was modified by blending benzoxazine resin (PBZ) and amino diluent (PAD) to prepare modified silicone alkyne hybrid resin (PPA) suitable for RTM molding process. The RTM molding process, curing behavior, and heat resistance of PPA resin were analyzed by DSC, FTIR, rheometer, rotational viscometer, and thermogravimetry. The results show that the processing window of PPA resin is wide, which can achieve the goal of RTM molding, but the curing temperature is higher than that of PSA resin. The heat resistance of PPA resin is excellent, and the temperature of 5% mass loss (Td5) of PPA-1 resin (PSA∶PBZ∶PAD mass ratio 5∶1∶2) resin is 585.4℃ and 568.3℃, respectively. The mass retention rates at 1000℃ is 88.3% and 26.0%, respectively. The mechanical property tests show that the mechanical properties of quartz fiber reinforced PPA resin composites (QF/PPA) gradually increase with the increase of PBZ content. The flexural strength of QF/PPA-1 at room temperature is 346.2 MPa, and the interlaminar shear strength is 21.4 MPa, which is 120.6% and 72.6% higher than that of QF/PSA composites, respectively. After thermal aging at 400℃ for 2 h, the flexural strength and interlaminar shear strength are 256.5 MPa and 17.1 MPa, respectively, and the retention rate of mechanical properties after thermal aging at 400℃ for 2 h exceeds 70%.
-
复合材料以其比重小、比强度和比模量高等特点,在航空航天领域被广泛应用[1],常见于飞机机翼、尾翼和机身的加筋壁板结构。在压缩、剪切等面内载荷作用下,加筋板结构最常见的失效模式就是失稳(或称屈曲)[2]。
求解加筋板的屈曲问题常采用数值方法,大多通过ABAQUS、ANSYS等商业有限元软件对加筋板屈曲和后屈曲问题进行分析[3-7]。邵青等[8]讨论了侧边边界条件对轴压复合材料加筋板屈曲和后屈曲性能的影响。孙中雷等[9]指出在截面积相等的前提下,帽形加筋板比T型、工字型加筋板的屈曲载荷更高。相对于数值方法,工程简化分析方法求解速度更快,可通过设计手册/指南查取相应计算公式、设计曲线和修正系数,从而对加筋板的屈曲载荷进行快速估算[10-12],但计算精度不高,很大程度上依赖设计经验。Mo等[13]和葛东云等[14]在工程简化分析方法中考虑了蒙皮厚度、加筋板曲率和筋条间距的影响,给出了轴压复合材料帽形加筋曲板的初始屈曲修正公式。于振波[15]采用刚度等效的方法对理想边界条件下的复合材料平板压缩屈曲计算方法进行了修正。杨俊清等[16]考虑了帽形筋条底脚对蒙皮的支撑作用,研究了多种边界条件下按不同的厚度折算方式时工程简化分析方法对复合材料帽形加筋壁板轴压屈曲分析的适用性。汪厚冰等[17]采用几种工程分析方法对轴压复合材料帽形加筋板进行了屈曲载荷计算,通过与试验结果对比研究,得到一种简单有效的复合材料帽形加筋板轴压屈曲载荷计算方法。陈金睿等[18]考虑了筋条下缘条对蒙皮板元宽度的影响,利用加筋板横截面特性提出了一种里兹能量法的改进计算方法。
综上所述,目前对于复合材料帽形加筋板初始屈曲的工程修正分析主要围绕蒙皮板元宽度、刚度修正及边界支持条件修正等方面开展研究。本文首先对典型复合材料帽形加筋板进行了轴压稳定性试验,然后基于加筋板的弯曲刚度比提出了一种复合材料帽形加筋板蒙皮局部屈曲工程修正计算方法。该方法通过选择合适的蒙皮板元宽度及考虑帽形筋条下缘条的蒙皮板元刚度等效处理方法,采用加筋板的弯曲刚度比对等效蒙皮板元的工程简支解和固支解进行线性修正。随后采用有限元方法及本文修正方法对两类带筋条下缘条的帽形加筋板进行了对比分析,采用有限元方法、工程简化分析方法、能量法及本文修正方法对公开文献中的带横框帽形加筋板试验数据进行了对比分析,验证了本文方法的合理性和有效性,其结果满足工程精度要求,为复合材料帽形加筋板的结构设计与分析提供了一种新的快速分析思路。
1. 复合材料帽形加筋板轴压稳定性试验
1.1 试件概述
试件为带两根横框的五筋条复合材料帽形加筋板结构,加筋板采用X850碳纤维增韧树脂预浸料,筋条与蒙皮共胶接。横框为2024-T42铝合金结构,由L形型材和槽形型材连接而成且在帽形筋条交叉处打孔以保证帽形筋条连续。筋条铺层[45/0/0/−45/90]s,蒙皮铺层[45/−45/−45/90/45/0]s,定义铺层0°方向为筋条轴向,单层厚度0.185 mm。试件材料属性见表1。试件共两件,两端做灌胶处理,如图1所示。其中1#试件加筋板长度1005 mm (单侧灌胶端长度50 mm),宽度955 mm,两横框间距为626 mm,筋条间距为210 mm;2#试件加筋板尺寸长度890 mm (单侧灌胶端长度60 mm),宽度1155 mm,两横框间距为554 mm,筋条间距为260 mm。筋条截面尺寸参数见图2。
表 1 单层材料参数Table 1. Lamina material parametersMaterial E11/GPa E22/GPa G12/GPa ν12 X850 162.0 9.1 4.6 0.331 Notes:E11—Elastic modulus in fiber direction; E22—Elastic modulus transverse fiber direction; ν12—Poisson’s ratio; G12—Shear modulus. 1.2 试验概述
轴向压缩试验在2000 kN油压四柱拉压试验机(LY2000,长春新实验机有限责任公司)上进行。试验机加载构件与压缩夹具上部相连。试样放置在夹具中间,调整夹具前后位置,以调整试验机加载中心。为确保试验件承受沿筋条方向的纯压载荷,并保证失稳前对中度,将试验件垂直放置在试验机支持平台上,横框的左右两侧通过螺栓连接到搭接板上,搭接板可随试验件压缩在滑轨中上下移动,以模拟横框的面外扭转约束,支持与加载方式如图3所示。调正压心后,由试验机的上平台施加轴向压缩载荷,载荷合力作用点(试验机加载中心)与加筋板截面压心一致。
加筋板应变片位置如图4所示,图中L表示两横框间距,L1和L2分别表示灌胶端边缘到横框底边的距离。应变片按图4中位置从上到下、从左到右依次编号,编号共7位编码:第一位为试件编号,“1”为1#试件,“2”为2#试件;第二位表示零件类型,“1”为蒙皮,“2”为筋条,“3”为横框;第三位表示行数,从上到下蒙皮共5行、筋条共5行、横框共2行,分别计数;第四位表示列数,从左到右内侧蒙皮(有筋条一侧)共4列、外侧蒙皮(无筋条一侧)共9列、筋条共5列、横框共4列,分开计数;第五位表示应变片位于加筋板外侧或内侧,“1”为外侧,“2”为内侧;第六位表示应变片在筋条上的位置,如图4(c)所示;第七位为应变片方向,“1”、“2”、“3”分别为0°、45°、90°方向,其中0°方向为筋条轴线方向。编码方法如表2所示。
表 2 应变片编码方法Table 2. Coding method of strain gaugeCoding location Coding meaning Coding method First number Specimen number 1—Specimen 1#; 2—Specimen 2# Second number Part type 1—Skin; 2—Stiffener Third number Row number of strain gauge 5 rows on skin, 5 rows on stiffeners, 2 rows on 2 frames, counting from up to down respectively Fourth number Column number of strain gauge 4 columns on inner skin, 9 columns on outer skin, 5 columns on 5 stiffeners, 4 columns on frame, counting from left to right respectively Fifth number Strain gauge on inner skin or outer skin 1—Outer skin; 2—Inner skin Sixth number Strain gauge location on stiffener 1—Crown; 2—Left web; 3—Right web; 4—Left bottom flange; 5—Right bottom flange Seventh number Strain gauge direction 1—0° direction; 2—45° direction; 3—90° direction 1#试件加载到315 kN蒙皮开始屈曲(图5(a)),2#试件加载到185 kN蒙皮开始屈曲(图5(b)),此时蒙皮载荷应变曲线的斜率出现反号,筋条的载荷应变曲线斜率发生变化。
1.3 有限元模型验证
采用ABAQUS软件对上述试验进行线性特征值屈曲分析。单元类型采用四节点常规壳减缩积分(S4R)单元,网格尺寸经收敛性分析[18]取为10 mm×10 mm。为模拟两端灌胶区,边界条件一端固支,另一端只释放轴压方向的自由度,并施加单位压缩载荷。试验件网格划分及一阶屈曲模态见图6 (以1#试件为例)。计算得到1#试件屈曲载荷286.9 kN,2#试件屈曲载荷164.6 kN,与试验值误差分别为−8.9%和−11.0% (误差正值表示计算值大于试验值,误差负值表示计算值小于试验值,下同),表明该有限元建模方法较准确,结果可信。
2. 基于弯曲刚度比的复合材料帽形加筋板蒙皮局部屈曲工程修正算法
2.1 复合材料帽形加筋板蒙皮局部屈曲工程简化分析方法概述
由于复合材料加筋板筋条间蒙皮局部屈曲是最常见的屈曲模式,因此仅讨论该屈曲模式下的加筋板屈曲载荷计算。在工程简化分析方法中,复合材料帽形加筋板筋条间蒙皮局部屈曲分析可简化为侧边简支或固支的各向异性平板屈曲分析,其中板宽为b,一般取相邻帽腰之间的距离或者下缘条自由边缘之间的距离;板长为a,取加筋板两个横框之间的距离。加载边边界条件视加持情况取为简支或者固支。以加载边简支为例[12],四边简支正交各向异性矩形层压平板的屈曲载荷计算公式为
Nx=π2D22b2[D11D22(ba)2m2+2(D12+2D66D22)+(ab)21m2] (1) 式中:Nx为单位长度的轴压屈曲载荷;m为沿板方向的屈曲半波数;D11、D12、D22和D66为蒙皮弯曲刚度系数,下同。计算时取m为1,2,3···,计算相应的一组Nx,其中最小的Nx即为板的屈曲载荷Nxcr。
两加载边简支、两侧边固支的正交各向异性矩形层压平板屈曲载荷计算公式为
Nxcr=π2√D11D22b2[K−2.40(1−D12+2D66√D11D22)] (2) 式中,K为压缩屈曲系数,按正交各向异性矩形平板的轴压K-λ曲线(图7)查取,其中
λ=(a/b)4√D22/D11 。加载边固支的计算公式详见文献[12],此处不再赘述。2.2 基于弯曲刚度比的复合材料帽形加筋板蒙皮局部屈曲修正方法
由于筋条对于蒙皮的实际支持不能简单地简化为简支或者固支,这种支持的强弱可以通过加筋板筋条与蒙皮的面内刚度、扭转刚度、弯曲刚度的比值反映[18]。由于筋条的面内刚度、扭转刚度、弯曲刚度具有一定相关性,为简化计算量,仅采用弯曲刚度比对工程简化分析进行修正。该方法计算步骤如下:
(1) 根据工程简化分析方法,将加筋板的屈曲计算简化为加载边简支或者固支的蒙皮板元的屈曲计算。其中板宽为b,取相邻帽腰之间的距离;板长为a,取加筋板两个横框之间的距离。
(2) 采用工程简化分析方法的平板屈曲公式分别计算蒙皮板元侧边简支的屈曲载荷
NSSxcr 和侧边固支的屈曲载荷NCCxcr 。考虑筋条下缘条对蒙皮刚度的加强作用,根据刚度叠加原则,等效后的蒙皮板元的弯曲刚度为b {\boldsymbol{D}} + 2{b_{{\text{bf}}}} {\boldsymbol{D}}_{\text{f}}^{} ,经简化后采用下式对蒙皮板元的弯曲刚度系数进行等效处理:{\boldsymbol{D}}_{\text{R}}^{} = {\boldsymbol{D}}_{}^{} + {\boldsymbol{D}}_{\text{f}}^{}\frac{{2{b_{{\text{bf}}}}}}{b} (3) 式中:
{{\boldsymbol{D}}_{\text{R}}} 为蒙皮板元的修正弯曲刚度系数矩阵;{\boldsymbol{D}} 为未修正的蒙皮板元弯曲刚度系数矩阵;{\boldsymbol{D}}{}_{\text{f}} 为筋条下缘条的弯曲刚度系数矩阵;bbf为筋条下缘条宽度;b为蒙皮板元宽度。(3) 基于加筋板的弯曲刚度比对侧边简支解和侧边固支解进行修正,修正公式如下:
{N_{{\text{cr}}}} = N_{x{\text{cr}}}^{{\text{SS}}} + (N_{x{\text{cr}}}^{{\text{CC}}} - N_{x{\text{cr}}}^{{\text{SS}}})\frac{{(D - 3)}}{3} (4) 式中,D为加筋板的弯曲刚度比,即加筋板n根筋条的弯曲刚度之和与蒙皮的弯曲刚度的比值。加筋板截面的中性轴位置及筋条/蒙皮的弯曲刚度采用文献[2]中方法计算。将上述得到的Ncr除以蒙皮板元厚度即得到加筋板的局部屈曲应力。该方法适用于带筋条下缘条的复合材料帽形加筋板的筋间蒙皮初始屈曲分析。
对第一节的两个帽形加筋板试验件分别采用工程简化分析方法和修正方法进行计算。横框支持简化为简支,蒙皮板元宽度选为相邻帽腰之间的距离,分别为151 mm和201 mm。采用侧边简支公式(1)计算得到屈曲载荷分别为249 kN和154 kN。采用侧边固支公式(2)计算时,λ分别为4.4和2.9,K值分别取为7.0和7.2,计算得到屈曲载荷分别为383 kN和244 kN。采用修正方法计算时首先采用式(3)将下缘条刚度与蒙皮板元刚度进行等效处理,然后采用式(1)计算得到侧边简支载荷分别为285 kN和173 kN,采用式(2)计算得到侧边固支载荷分别为447 kN和275 kN,然后计算得到加筋板弯曲刚度比分别为3.501和3.456,最后采用修正式(4)计算得到修正后的屈曲载荷,分别为312 kN和189 kN,与试验值误差分别为1.0%和2.2%,均在3%以内。详细结果见表3。可以看出,传统的工程简化分析方法计算的屈曲载荷误差较大,而采用修正方法,其计算精度明显提高。
表 3 复合材料帽形加筋板轴压稳定性试验值与计算值对比Table 3. Comparison of calculated and experimental values of local buckling load of hat-stiffened composite panelNo. D Experimental
buckling load /kNEngineering simplified
method/kNUpdating
method/kNError/% Simply
supportedBuilt-in Simply
supportedBuilt-in Updating
method1# 3.501 315 249 383 312 −21.0 21.6 −1.0 2# 3.456 185 154 244 189 −16.8 31.9 2.2 Note: D—Flexural stiffness ratio. 3. 算例分析
以不带横框的三筋条复合材料帽形加筋板为研究对象,分别采用有限元方法、工程简化分析方法和本文修正方法进行对比分析,验证修正方法的合理性和有效性。随后引用文献试验数据,以带横框的三筋条复合材料帽形加筋板为研究对象,综合比较工程简化分析方法、能量法及本文修正方法,形成适用于帽形加筋板初步分析的工程估算方法。
3.1 不带横框的三筋条帽形加筋板屈曲分析
研究对象为两类含筋条下缘条的帽形筋条加筋板,如图8所示。两类加筋板的几何尺寸和材料属性如表4和表5所示。表4中a为加筋板长度,B为加筋板宽度,bc为筋条帽顶宽度,bw为筋条帽腰宽度,bf为筋条下缘条宽度,θ为帽腰与蒙皮夹角,bss为相邻筋条间蒙皮宽度,如图9所示。对于每一类加筋板,其几何尺寸不变,分别给定3种不同的加筋板铺层,共6种帽形加筋板,铺层顺序见表6。
有限元分析采用ABAQUS有限元软件进行线性特征值分析,模型采用四节点常规壳减缩积分(S4R)单元,并对网格尺寸进行收敛性分析[18];加筋板两侧边自由,两加载边采用参考点耦合约束,一参考点施加单位压缩载荷,另一参考点建立简支边界条件。
表 4 三筋条帽形加筋板有限元模型几何参数Table 4. Geometric parameters of finite element models of hat-stiffened panels with three stiffenersType a
/mmB
/mmbc
/mmbw
/mmbf
/mmθ/(°) bss
/mm1 300 278 10.4 10.7 16.5 82.9 70 2 590 532 25 33 25 60 103 Notes: a—Length of stiffened panel; B—Width of stiffened panel; bc—Width of stiffener’s crown; bw—Width of stiffener’s web; bf—Width of stiffener’s bottom flange; θ—Angle between stiffener’s web and skin; bss—Width of skin adjacent stiffener. 表 5 三筋条帽形加筋板单层材料参数Table 5. Lamina material parameters of hat-stiffened panels with three stiffenersType E11/GPa E22/GPa G12/ GPa {\nu _{12}} 1 98 10.8 5.2 0.31 2 154 8.5 4.5 0.35 图10为两类加筋板的一阶屈曲模态图。表7为分别采用有限元方法、工程简化分析方法和修正方法计算的6种加筋板的屈曲载荷。表中误差正值表示计算值大于有限元值,误差负值表示计算值小于有限元值。结果分析如下:
(1) 6种加筋板屈曲模式均为筋条间蒙皮局部屈曲,筋条无明显扭转、弯曲变形,符合本文研究前提;
(2) 采用工程简化分析方法,侧边简支解误差绝对值的算术平均值为48.7%,表明简支解比有限元值明显偏小。工程经验表明按简支和按固支计算的屈曲载荷相差可达一倍以上,对于帽形筋条工程上一般简化为固支[12]。因此按简支计算误差较大;
表 6 三筋条帽形加筋板有限元模型铺层顺序及铺层厚度Table 6. Stacking sequence and thickness of finite element models of hat-stiffened panels with three stiffenersNo. Skin Crown Web Flange Stacking
sequenceTotal
thickness/
mmStacking
sequenceTotal
thickness/
mmStacking
sequenceTotal
thickness/
mmStacking
sequenceTotal
thickness/
mm1.1 [45/0/0/−45/90]s 1 [09/−45/452/−452]s 2.8 [−45/452/−452]s 1 [03/−45/452/
−45/−45]s1.5 1.2 [45/0/−45/90/
45/−45/0]s1.4 [09/−45/452/−452]s 2.8 [−45/452/−452]s 1 [03/−45/452/
−45/−45]s1.5 1.3 [0/90]3s 1.2 [45/0/90/0/−45/
0/45/90]s1.5 [45/0/90/0/−45/
0/45/90]s1.5 [45/0/90/0/−45/
0/45/90]s1.5 2.1 [45/−45/−45/
90/45/0]s1.2 [45/0/0/−45/90]s 0.9 [45/0/0/−45/90]s 0.9 [45/0/0/−45/90]s 0.9 2.2 [45/−45/902/
45/0/90/0]s1.5 [0/45/0/0/−45/45]s 1.2 [0/45/0/0/−45/45]s 1.2 [0/45/0/0/−45/45]s 1.2 2.3 [45/−45/902/45/
0/90/0]s1.5 [45/0/0/−45/90]s 0.9 [45/0/0/−45/90]s 0.9 [45/0/0/−45/90]s 0.9 表 7 两类帽形加筋板轴压稳定性有限元值与计算值对比Table 7. Comparison of calculated and FEM values of local buckling load of two kinds of hat-stiffened composite panelsNo. D FEM
buckling
load/kNEngineering simplified method/kN Updating method /kN Error/% Simply supported Built-in Simply
supportedBuilt-in Updating method 1.1
1.2
1.3
2.1
2.2
2.32.750
4.713
4.842
4.014
4.019
3.81513.0
38.5
23.9
26.4
52.7
43.27.7
15.9
9.4
16.0
28.5
23.113.1
29.0
19.1
24.7
46.5
37.714.2
38.1
26.0
27.9
54.6
41.2−40.8
−58.7
−60.7
−39.4
−45.9
−46.50.8
−24.7
−20.1
−6.4
−11.8
−12.79.2
−1.0
8.8
5.7
3.6
−4.6(3) 采用工程简化分析方法,侧边固支解误差绝对值的算术平均值为12.8%,相较于侧边简支解误差有所减小,表明将帽形筋条简化为固支有一定合理性,但是由于筋条间蒙皮间距较小,筋条下缘条对于蒙皮屈曲的加强作用较明显,而工程简化分析方法没有考虑筋条下缘条的加强作用,因此计算结果仍有一定误差;
(4) 采用修正方法之后的屈曲载荷误差绝对值的算术平均值为5.5%,满足工程精度要求,相较于工程简化分析方法计算精度有所提高,表明修正方法合理有效。
3.2 带横框的三筋条帽形加筋板屈曲分析
采用试验方法、工程简化分析方法、能量法和本文修正方法对三筋条带横框的复合材料帽形加筋板进行轴压屈曲分析,加筋板外形见图11,详细尺寸参数、铺层信息、单层材料参数及稳定性试验结果详见文献[16]。
采用工程简化分析方法进行屈曲分析时,将加筋板简化为长度为500 mm、宽度为195 mm (帽腰之间的距离)的平板,采用式(1)计算得到侧边简支屈曲载荷143.3 kN。侧边固支屈曲载荷258 kN[16]。
采用修正方法计算得到弯曲刚度比3.641,分别采用式(1)和式(2)计算侧边简支和侧边固支屈曲载荷。计算固支载荷时λ=2.2,K值为7.2。计算得到简支载荷205.3 kN,固支载荷348.0 kN。采用式(4)计算得到修正后的屈曲载荷为235.8 kN。
采用能量法[18]进行屈曲分析时,将加筋板简化为长度为500 mm、宽度为195 mm的平板,帽形筋条的自由截面扭转刚度为7.48×108 N·mm2,代入到能量法计算模型中得到屈曲载荷为238.2 kN,计算结果见表8。表中误差正值表示计算值大于试验值,误差负值表示计算值小于试验值。结果分析如下:
(1) 对于带横框的帽形加筋板,工程简化分析方法、能量法和本文修正方法均将横框简化为简支。除工程简化分析方法简支解误差较大之外,工程简化分析方法固支解、能量法和本文修正方法计算精度满足工程要求,表明对帽形筋条简化为固支的一般工程假设是合理的,也表明对横框支持的简化是合理的;
(2) 从计算精度角度看,修正方法计算误差3.8%,能量法计算误差2.8%,两者计算精度相当,比工程简化分析方法精度稍高。其原因在于工程简化分析方法将帽形筋条简化为固支,对于某些带下缘条的帽形加筋板计算误差较大,能量法则简化为弹性支持,修正方法则基于弯曲刚度比对简支解和固支解综合修正,是对工程简化分析方法的修正和改善;
(3) 从时间成本上看,3种方法相比有限元方法均为快速计算方法,可编写计算程序,用时较短。
综上所述,工程简化分析方法算法简单,依赖设计员工程经验,适用于方案阶段结构参数打样;修正方法和能量法简化方式更合理,适用范围更广,适合于初步设计阶段结构参数化分析和稳定性校核,可相互对照使用。
表 8 带横框的三筋条帽形加筋板局部屈曲载荷计算值与试验值对比Table 8. Comparison of calculated and experimental values of local buckling load of hat-stiffened panels with three stiffeners and two framesD Experimental
buckling
load/kNEngineering simplified method/kN Energy method/kN Updating method/kN Error/% Simply supported Built-in Simply supported Built-in Energy method Updating method 3.641 245.0 143.3 258 238.2 235.8 −41.5 5.3 −2.8 −3.8 4. 结 论
(1) 对典型复合材料帽形加筋板进行了轴压稳定性试验,提出了一种基于加筋板弯曲刚度比的复合材料帽形加筋板蒙皮局部屈曲工程修正计算方法,所预测的屈曲载荷与稳定性试验值误差在3%以内。
(2) 对于带筋条下缘条的三筋条复合材料帽形加筋板,工程简化分析方法侧边简支解误差绝对值的算术平均值为48.7%,侧边固支解误差绝对值的算术平均值为12.8%,表明将帽形筋条简化为固支有一定合理性,但是由于没有考虑筋条下缘条对蒙皮的加强作用,因此计算结果仍有一定误差。修正方法误差绝对值的算术平均值为5.5%,验证了修正方法的合理性和有效性,对带筋条下缘条的复合材料帽形加筋板筋间初始屈曲载荷估算有一定的参考意义。
(3) 对于典型的带横框三筋条帽形加筋板,修正方法计算误差3.8%,与能量法计算精度相当,比工程简化分析方法精度稍高。工程简化分析方法算法简单,依赖设计员工程经验,适用于方案阶段结构参数打样;修正方法适用于初步设计阶段结构参数化分析和稳定性评估,可与能量法相互对照使用。
-
表 1 不同比例的改性硅炔杂化树脂(PPA)配方
Table 1 Modified silicone alkyne hybrid resin (PPA) formula with different proportions
Sample Resin mass ratio PSA PBZ PAD PPA-1 5 1 2 PPA-2 5 2 2 PPA-3 5 3 2 Note: PAD—Amino diluent. 表 2 PPA树脂的黏度
Table 2 Viscosity of PPA resin
Sample Viscosity/(mPa·s) 25℃ 30℃ 40℃ PSA 740 350 100 PPA-1 1000 610 230 PPA-2 1700 1000 400 PPA-3 4000 1800 640 表 3 PPA树脂的凝胶时间
Table 3 Gelation time of PPA resin
Sample Gel time/min 160℃ 180℃ 200℃ PPA-1 23.5 11.0 5.1 PPA-2 20.1 9.3 4.0 PPA-3 18.2 8.1 3.0 表 4 PPA树脂的TGA数据
Table 4 TGA data of PPA resin
Sample Td5/℃ Mass retention/% Air N2 Air N2 PSA 563.1 689.8 31.26 91.36 PPA-1 568.3 613.7 26.01 89.13 PPA-2 555.4 601.8 23.04 88.78 PPA-3 538.6 595.9 21.66 88.74 Note: Td5—Temperature of 5% mass loss. 表 5 QF/PPA复合材料热老化10 h质量损失
Table 5 Mass loss of QF/PPA composites thermal aging
Sample Mass loss/wt% QF/PSA 4.28 QF/PPA-1 2.06 QF/PPA-2 2.50 QF/PPA-3 2.80 -
[1] SHEN Y, YUAN Q, HUANG F, et al. Effect of neutral nickel catalyst on cure process of silicon-containing polyarylacetylene[J]. Thermochimica Acta,2014,590:66-72. DOI: 10.1016/j.tca.2014.06.002
[2] LANGE N, DIETRICH P M, LIPPITZ A, et al. New azidation methods for the functionalization of silicon nitride and application in copper-catalyzed azide-alkyne cycloaddition (CuAAC)[J]. Surface and Interface Analysis, 2016, 48(7):621-625.
[3] 陈元俊, 郭康康, 王帆, 等. 含氟硅芳炔树脂的合成与性能研究[J]. 华东理工大学学报(自然科学版), 2020, 46(3):368-375. CHEN Yuanjun, GUO Kangkang, WANG Fan, et al. Synthesis and performance of fluorine-containing silicon arylacetylene resins[J]. Journal of East China University of Science and Technology,2020,46(3):368-375(in Chinese).
[4] GAO G, ZHANG S, WANG L, et al. Developing highly tough, heat-resistant blend thermosets based on silicon-containing arylacetylene: A material genome approach[J]. ACS Applied Materials & Interfaces,2020,12(24):27587-27597.
[5] MA M, GONG C, LI C, et al.The synthesis and properties of silicon-containing arylacetylene resins with rigid-rod 2, 5-diphenyl-[1, 3, 4]-oxadiazole moieties[J]. European Polymer Journal, 2020, 143(4): 110192.
[6] YANG R, WANG Y, HAO B, et al. Synthesis of ortho-methyltetrahydrophthalimide functional benzoxazine containing phthalonitrile group: Thermally activated polymerization behaviors and properties of its polymer[J]. High Performance Polymers,2021,33(2):196-204. DOI: 10.1177/0954008320954519
[7] CORRIU R, GERBIER P, GUÉRIN C, et al. Poly[(silylene) diacetylene]/finemetal oxide powder dispersions: Use as precursors to silicon-based composite ceramics[J]. Jour-nal of Materials Chemistry,2000,10(9):2173-2182. DOI: 10.1039/b002788j
[8] ITOH M, INOUE K, HIRAYAMA N, et al. Fiber reinforced plastics using a new heat-resistant silicon based polymer[J]. Journal of Materials Science,2002,37(17):3795-3801. DOI: 10.1023/A:1016538014803
[9] ITOH M, INOUE K, IWATA K, et al. New highly heat-resistant polymers containing silicon: Poly(silyleneethynylene-phenyleneethynylene)s[J]. Macromolecules,1997,30(4):694-701. DOI: 10.1021/ma961081f
[10] ITOH M, IWATA K, ISHIKAWA J I, et al. Various silicon-containing polymers with Si(H)C≡C units[J]. Journal of Polymer Science, Part A: Polymer Chemistry,2001,39(15):2658-2669. DOI: 10.1002/pola.1242
[11] 袁航, 孟庆杰, 张昊, 等. 新型含硅聚芳炔树脂基透波复合材料的制备与性能[J]. 复合材料学报, 2021, 38(11):3629-3639. DOI: 10.13801/j.cnki.fhclxb.20210210.009 YUAN Hang, MENG Qingjie, ZHANG Hao, et al. Preparation and properties of novel silicon-containing polyarylacetylene resin based wave-transparent composite[J]. Acta Materiae Compositae Sinica,2021,38(11):3629-3639(in Chinese). DOI: 10.13801/j.cnki.fhclxb.20210210.009
[12] 王志德, 白小陶, 贾宇翔, 等. 硼硅炔杂化聚合物的合成及耐热性能[J]. 高分子材料科学与工程, 2022, 38(10):17-22. WANG Zhide, BAI Xiaotao, JIA Yuxiang, et al. Synthesis and heat resistance of borosilyne hybrid polymers[J]. Polymer Materials Science and Engineering,2022,38(10):17-22(in Chinese).
[13] 司书帅, 袁荞龙, 黄发荣. 含炔基酚醛树脂改性PSA及碳纤维布/PSA-EPAN复合材料性能[J]. 复合材料学报, 2018, 35(3):545-552. SI Shushuai, YUAN Qiaolong, HUANG Furong. Properties of PSA modified by alkynyl-containing phenolic resin and carbon fabric/PSA-EPAN composites[J]. Acta Materiae Compositae Sinica,2018,35(3):545-552(in Chinese).
[14] TONG Y, YUAN Q, HUANG F. Preparation and characterization of silicon-containing polyarylacetylene/montmorilloite nanocomposites[J]. Journal of Macromolecular Science Part B—Physics,2019,58(4):469-488. DOI: 10.1080/00222348.2019.1590982
[15] 王茂源, 束长朋, 贾宇翔, 等. 双酚A型邻苯二甲腈树脂改性硅炔杂化树脂及其复合材料性能[J]. 复合材料学报, 2021, 38(11):3652-3660. DOI: 10.13801/j.cnki.fhclxb.20210114.003 WANG Maoyuan, SHU Changpeng, JIA Yuxiang, et al. Bisphenol A type o-phthalonitrile resin modified silicoalkyne hybrid resin and its composite properties[J]. Acta Materiae Compositae Sinica,2021,38(11):3652-3660(in Chinese). DOI: 10.13801/j.cnki.fhclxb.20210114.003
[16] 王佳明, 贾明印, 董贤文, 等. 树脂传递模塑成型工艺研究进展[J]. 塑料工业, 2021, 49(11):9-14, 43. DOI: 10.3969/j.issn.1005-5770.2021.11.003 WANG Jiaming, JIA Mingyin, DONG Xianwen, et al. Research progress on the molding process of resin transfer mould[J]. China Plastics Industry,2021,49(11):9-14, 43(in Chinese). DOI: 10.3969/j.issn.1005-5770.2021.11.003
[17] 金东升, 张庆茂, 张朋, 等. RTM在耐高温复材结构中的应用问题与对策[C]//第二十一届全国复合材料学术会议(NCCM-21). 北京, 2020: 359-364. JIN Dongsheng, ZHANG Qingmao, ZHANG Peng, et al. Application problems and counter measures of RTM in high temperature resistant composite structures[C]//The 21st National Academic Conference on Composite Materials (NCCM-21). Beijing, 2020: 359-364(in Chinese).
[18] 朱怡臻, 王瑛, 陈鸣亮, 等. 先进树脂基复合材料RTM成型工艺研究及应用进展[J]. 塑料工业, 2020, 48(5):18-22, 128. DOI: 10.3969/j.issn.1005-5770.2020.05.002 ZHU Yizhen, WANG Ying, CHEN Mingliang, et al. Research progress and application of RTM for advanced resin matrix composites[J]. China Plastics Industry,2020,48(5):18-22, 128(in Chinese). DOI: 10.3969/j.issn.1005-5770.2020.05.002
[19] 赵安安, 杨文凯, 于飞, 等. 大型高性能复合材料构件RTM工艺进展[J]. 南京航空航天大学学报, 2020, 52(1):39-47. DOI: 10.16356/j.1005-2615.2020.01.004 ZHAO An'an, YANG Wenkai, YU Fei, et al. RTM progress for large-scale high-performance composite components[J]. Journal of Nanjing University of Aeronautics and Astronautics,2020,52(1):39-47(in Chinese). DOI: 10.16356/j.1005-2615.2020.01.004
[20] 祝庆君. 环氧树脂增韧改性及RTM工艺制备碳纤维复合材料性能研究[D]. 哈尔滨: 黑龙江省科学院石油化学研究院, 2018. ZHU Qingjun. Study on the toughening modification of epoxy resin and the properties of carbon fiber composites prepared by RTM process[D]. Harbin: Institute of Petrochemistry, Heilongjiang Academy of Sciences, 2018(in Chinese).
[21] GAJJAR T, SHAH D, JOSHI S, et al. Investigation on dimensional accuracy for CFRP antenna reflectors using autoclave and VARTM processes[C]//Advances in Computational Methods in Manufacturing. Lecture Notes on Multidisciplinary Industrial Engineering. Singapore: Springer, 2019: 693-702.
[22] AMIRKHOSRAVI M, PISHVAR M, CENGIZ ALTAN M. Fabricating high-quality VARTM laminates by magnetic consolidation: Experiments and process model[J]. Composites Part A: Applied Science and Manufacturing,2018,114:398-406. DOI: 10.1016/j.compositesa.2018.09.003
[23] 轩立新, 王茂源, 束长朋, 等. RTM成型工艺用改性硅炔杂化树脂性能研究[J]. 中国胶粘剂, 2020, 29(6):1-6. DOI: 10.13416/j.ca.2020.06.002 XUAN Lixin, WANG Maoyuan, SHU Changpeng, et al. Study on properties of modified silicone alkyne hybrid resin for RTM molding process[J]. China Adhesives,2020,29(6):1-6(in Chinese). DOI: 10.13416/j.ca.2020.06.002
[24] 中国国家标准化管理委员会. 纤维增强塑料弯曲性能试验方法: GB/T 1449−2005[S]. 北京: 中国标准出版社, 2005. Standardization Administration of the People's Republic of China. Fiber-reinforced plastic composites: Determination of flexural properties: GB/T 1449−2005[S]. Beijing: China Standards Press, 2005(in Chinese).
[25] 中国国家标准化管理委员会. 纤维增强塑料层间剪切强度试验方法: GB/T 1450.1−2005[S]. 北京: 中国标准出版社, 2005. Standardization Administration of the People's Republic of China. Fiber-reinforced plastic composites: Determination of interlaminar shear strength: GB/T 1450.1−2005[S]. Beijing: China Standards Press, 2005(in Chinese).
-
期刊类型引用(2)
1. 陈向明,李新祥,柴亚南,陈普会,孙侠生. 复合材料壁板后屈曲设计与分析技术研究进展. 复合材料学报. 2024(09): 4673-4700 . 本站查看
2. 吴菁,章向明,朱子旭,胡明勇,王安稳. 复合材料船体帽型加筋板面内刚度. 海军工程大学学报. 2023(04): 1-6+28 . 百度学术
其他类型引用(1)
-
目的
RTM成型工艺是一种液体闭模成型工艺,成品表面光滑、安全环保、成本较低,极具发展潜力。硅炔杂化树脂是一种耐高温性能优异的有机无机杂化树脂,加工性能良好、介电性能突出,在航空航天、核电等军用民用领域均有广泛应用,但其与纤维的结合性能较差,复合材料力学性能不佳,成为制约其发展应用的主要因素。本文从硅炔杂化树脂(PSA)出发,采用苯并噁嗪PBZ和氨基稀释剂PAD共混改性硅炔杂化树脂PSA,制备了适用于RTM成型工艺的PPA树脂,以提高复合材料的力学性能。
方法以苯并噁嗪树脂(PBZ)和氨基稀释剂(PAD)共混改性PSA树脂,制备三元树脂体系PSA/PBZ/PAD(PPA树脂),并采用RTM工艺制备复合材料。通过DSC、FTIR分析改性硅炔杂化树脂的结构变化,研究固化过程与反应机理;利用旋转流变仪、旋转黏度计测试改性硅炔杂化树脂黏度及流变特性,分析PPA树脂的RTM成型工艺窗口;采用TGA测定改性硅炔杂化树脂耐热性能变化情况;通过力学试验测试对比PPA树脂复合材料的室温和400 ℃热老化后的力学性能变化;利用DMA及热老化测试研究复合材料的热性能。
结果(1)PPA树脂在常温下是具有一定黏度的液体状树脂,40 ℃黏度降至800 mPa·s以下,且在40 ℃能够长时间保持黏度稳定。(2)PPA树脂在250 ℃附近存在单一的的固化峰,随着PBZ的含量增加,固化峰值温度逐渐降低。固化过程中PSA树脂、PBZ树脂和PAD之间存在相互反应,包括噁嗪环开环、炔基与氨基的氢胺化反应、Si−H与炔基的硅氢加成反应以及炔基之间的Diels-Alder反应和环三聚反应。(3)PPA树脂固化物质量损失5%的温度(T)和质量保留率均低于PSA树脂固化物,并且T和质量保留率均随着PBZ树脂的增加而逐渐下降,其中PPA-1(PSA:PBZ:PAD质量比为5:1:2)树脂在氮气和空气中的T分别为585.4 ℃ 和568.3 ℃,1000 ℃质量保留率分别为88.3% 和26.0%。(4)石英纤维增强PPA树脂复合材料(QF/PPA)力学性能随PBZ含量的增加而逐渐升高,其中QF/PPA-1常温下弯曲强度为346.2 MPa,层间剪切强度为21.4 MPa,较QF/PSA复合材料分别提高了120.6%和72.6%。400 ℃热老化2 h后,QF/PPA-1复合材料的弯曲强度和层间剪切强度分别为256.5 MPa和17.1 MPa,400 ℃热老化2 h后的力学性能保留率超过70%。(6)QF/PPA复合材料的储能模量随着温度升高明显提高,玻璃化转变温度均高于500 ℃。(7)QF/PPA复合材料在400 ℃高温下长时间热老化后,整体颜色发白,结构完好,层间结合牢固,其中在热老化10 h后,QF/PPA复合材料质量损失率小于3%,低于QF/PSA复合材料。
结论本文制备了改性硅炔杂化树脂PPA,PPA树脂加工窗口较宽,40 ℃的黏度满足RTM成型工艺的要求,但固化温度高于PSA树脂。PBZ/PAD引入改善了PSA树脂与纤维的界面结合性能,复合材料的力学性能均有大幅提升,室温弯曲强度较QF/PSA复合材料提高了120%以上。PPA树脂和QF/PPA复合材料的耐热性能均非常优异,400 ℃热老化10 h后QF/PPA复合材料质量损失率低于3%,玻璃化转变温度高于500 ℃。
-
硅炔杂化树脂是一种兼具优异加工性、介电性能、耐热性能、低密度和低吸水率的有机无机杂化树脂,广泛应用于航空航天及核电等军用民用领域。但硅炔杂化树脂固化之后交联密度过高,脆性大,与纤维界面结合性差,导致复合材料力学性能较差,限制了其在各领域的应用。
本文从硅炔杂化树脂PSA出发,采用含乙炔基苯并噁嗪PBZ和氨基稀释剂PAD共混改性硅炔杂化树脂PSA,得到改性硅炔杂化树脂PPA,并制备石英纤维增强PPA树脂复合材料(QF/PPA)。引入PBZ和PAD后,PPA树脂不仅黏度降低至600 mPa·s且能够保持足够长的时间,满足RTM成型工艺的要求,同时提高了树脂与纤维之间的结合性,改善了复合材料的界面性能,并且PPA树脂中三种组分相互之间可以交联固化,耐热性能优异。因此,PPA树脂可以在40 ℃注射树脂,浸润纤维,注射温度较低,安全环保。所制备的QF/PPA复合材料力学性能明显提高,其中QF/PPA-1常温下弯曲强度为346.2 MPa,层间剪切强度为21.4 MPa,较QF/PSA复合材料分别提高了120.6%和72.6%,400 ℃热处理2 h后的弯曲强度和层间剪切强度分别为256.5 MPa和17.1 MPa,高温力学强度保留率超过70%。
QF/PPA复合材料的弯曲强度和层间剪切强度