发明内容
本发明建立了一种预测平面斜交编织复合材料双向拉伸模量与强度的新方法,该方法具有计算简便且精度高等优点,其技术方案如下:
步骤一、根据纤维束的对称平面斜交编织方式(如周期性和重复性等),选择最小的周期性单元作为代表性体积元,由此确定其胞体单元,定义胞体单元内斜交编织几何形状与尺寸。
根据图1所示的平面对称斜交编织布的编织形状,考虑编织的周期性和重复性,选择代表性体积单元。图2为所选的代表性体积单元(胞体单元),单胞内包含了4条对称斜交的纤维束,定义由该四条纤维束组成的菱形的对角线方向分别为横向和纵向。编织角α指两个方向纤维束的夹角之半(如图2所示)。由于纤维束编织的中心对称性,可以仅取其中一条纤维束进行内力分析。因为纤维束自身的轴向横截面尺寸相比与其长度来说很小,将其视为受力情况是面内拉伸,面外弯曲的弯曲梁。由于纤维束编织导致的相互挤压作用,将其波动状态用正弦波进行描述。可以得到纤维束的中心线Z坐标表达式:
纤维束截面形状考虑纤维束间的相互挤压作用后认为如图3所示,截面两端为圆弧过渡,纤维束的截面最大宽度为w,厚度为h。可由此计算纤维束截面的截面面积A和截面惯性矩Iy、Iz、Ip的表达式,得到
Ip=Iy+Iz(5)
式中w和h分别为纤维束截面的宽度和厚度。
步骤二、对步骤一中的胞体单元施加双向面内拉伸外载,对胞体单元内织布的纤维束进行受力分析,从而,建立平面斜交编织复合材料胞体单元纤维织布的细观力学模型,确定胞体单元内纤维织布的总应变余能,并利用最小余能原理,求解胞体单元纤维织布的内力分布。
在单胞上作用横向及纵向拉伸载荷P1,P2。由单胞的受力分析(如图4所示)可以看出这是典型的对称结构在对称载荷作用下的静不定问题。由受力分析可知(如图5所示)梁的内力主要有轴力N(x),面内剪力FS1,面外剪力FS2(x),面内弯矩M1(x),面外弯矩M2(x),扭矩T(x)。
在对称载荷作用下,对称结构对称轴处横截面上的反对称内力为零,可知梁端面上的轴力与剪力为
由此可得到弯曲梁的所有内力表达式,即
可知这个静不定结构中有四个未知内力M,M1,Q,T。采用最小余能原理求解多余未知力,首先得到单胞的总余能表达式,Ji代表计算中间变量,可参见附录A。
由最小余能原理,可知,
计算未知内力,即
步骤三、根据卡式定理,求解胞体单元纤维织布的变形,再根据应力-应变本构方程,得到平面斜交编织复合材料织布的双向拉伸模量。在此基础上,根据混合定理,可以得到平面斜交编织复合材料的双向拉伸模量。
将未知内力代入单胞总余能表达式,认为材料满足线弹性假设,根据能量法定义可计算得到双向拉伸外载作用下的单胞双向位移。
根据应力和应变定义,可得到单胞的应力和应变表达式
根据应力和应变关系表达式,可以得到纤维织布双向拉伸模量公式:
式中E1和E2分别表示纤维织布的横向及纵向拉伸模量,σ1和σ2分别表示横向及纵向拉伸应力,ε1和ε2分别表示横向及纵向拉伸应变。
取单胞厚度为t,纤维编织布占单胞的体积分数可表示为
根据混合定理,可以得到平面斜交编织复合材料的双向弹性模量为
步骤四、根据组分材料(纤维束和基体)的破坏强度,采用最大拉应力准则和畸变能准则分别作为纤维束和基体破坏判据,并利用变形协调条件,分析平面斜交编织复合材料的双向拉伸强度,得到其解析解。
在已知单胞的内力分布后,对危险点进行破坏强度计算,并根据纤维束与基体的强度极限计算平面斜交编织复合材料的拉伸强度。采用最大拉应力准则和畸变能准则分别作为纤维束和基体破坏判据,即
当平面斜交编织复合材料横向和纵向上有拉伸外载T1和T2作用时,纤维编织布和基体各自承担一部分外载Tf,Tm,即
由于单胞内纤维束与基体的变形协调,在已知两者模量的基础上,可得
联立方程(19)和(20),可得
令T2为定值,研究T1满足强度准则的最大值,即横向极限载荷。根据强度准则,基体和纤维编织布存在极限载荷[Ff1]和[Fm1],即
对于平面斜交编织复合材料的横向极限载荷[F1],横向拉伸强度σx的解析解变为
根据前面单胞的内力分析以及强度准则,可分别求解[Ff1]和[Fm1]。
(1)求解[Ff1]
由于纤维编织布可看作由弯曲梁组成的编织结构,其梁内正应力由弯矩与轴力共同作用引起,因此,
因为编织交叉处(x=0或x=2L)变形最大,强度最值得考虑,则
根据第一强度理论,令σb1=σmax,则
根据上式,可得到纤维编织布横向及纵向的极限载荷[Ff1]:
(2)求解[Fm1]
基体上的双向拉伸应力可由名义应力法得到,即
取基体内部内的任一点的应力状态为双向拉应力状态,根据第四强度理论,可得
求解方程(31),可得[Fm1]。
将[Ff1]和[Fm1]的计算结果代入式(23)和式(24),可得到平面斜交编织复合材料的横向拉伸强度。同理,按照此方法可以得到纵向拉伸强度。
本发明是一种预测平面斜交编织复合材料双向拉伸模量与强度的新方法,其特点是根据编织几何形状和尺寸以及纤维束和基体的性能参数,便可方便快捷地预测平面斜交编织复合材料的双向拉伸模量与强度。
具体实施方式
步骤一、根据纤维束的对称平面斜交编织方式(如周期性和重复性等),选择最小的周期性单元作为代表性体积元,由此确定其胞体单元,定义胞体单元内斜交编织几何形状与尺寸。
根据图1所示的平面对称斜交编织布的编织形状,考虑编织的周期性和重复性,选择代表性体积单元。图2为所选的代表性体积单元(胞体单元),单胞内包含了4条对称斜交的纤维束,定义由该四条纤维束组成的菱形的对角线方向分别为横向和纵向。编织角α指两个方向纤维束的夹角之半(如图2所示)。由于纤维束编织的中心对称性,可以仅取其中一条纤维束进行内力分析。因为纤维束自身的轴向横截面尺寸相比与其长度来说很小,将其视为受力情况是面内拉伸,面外弯曲的弯曲梁。由于纤维束编织导致的相互挤压作用,将其波动状态用正弦波进行描述。可以得到纤维束的中心线Z坐标表达式:
纤维束截面形状考虑纤维束间的相互挤压作用后认为如图3所示,截面两端为圆弧过渡,纤维束的截面最大宽度为w,厚度为h。可由此计算纤维束截面的截面面积A和截面惯性矩Iy、Iz、Ip的表达式,得到
Ip=Iy+Iz(5)式中w和h分别为纤维束截面的宽度和厚度。
步骤二、对步骤一中的胞体单元施加双向面内拉伸外载,对胞体单元内织布的纤维束进行受力分析,从而,建立平面斜交编织复合材料胞体单元纤维织布的细观力学模型,确定胞体单元内纤维织布的总应变余能,并利用最小余能原理,求解胞体单元纤维织布的内力分布。
在单胞上作用横向及纵向拉伸载荷P1,P2。由单胞的受力分析(如图4所示)可以看出这是典型的对称结构在对称载荷作用下的静不定问题。由受力分析可知(如图5所示)梁的内力主要有轴力N(x),面内剪力FS1,面外剪力FS2(x),面内弯矩M1(x),面外弯矩M2(x),扭矩T(x)。
在对称载荷作用下,对称结构对称轴处横截面上的反对称内力为零,可知梁端面上的轴力与剪力为
由此可得到弯曲梁的所有内力表达式,即
可知这个静不定结构中有四个未知内力M,M1,Q,T。采用最小余能原理求解多余未知力,首先得到单胞的总余能表达式,Ji代表计算中间变量,可参见附录A。
由最小余能原理,可知,
计算未知内力,即
步骤三、根据卡式定理,求解胞体单元纤维织布的变形,再根据应力-应变本构方程,得到平面斜交编织复合材料织布的双向拉伸模量。在此基础上,根据混合定理,可以得到平面斜交编织复合材料的双向拉伸模量。
将未知内力代入单胞总余能表达式,认为材料满足线弹性假设,根据能量法定义可计算得到双向拉伸外载作用下的单胞双向位移。
根据应力和应变定义,可得到单胞的应力和应变表达式
根据应力和应变关系表达式,可以得到纤维织布双向拉伸模量公式:
式中E1和E2分别表示纤维织布的横向及纵向拉伸模量,σ1和σ2分别表示横向及纵向拉伸应力,ε1和ε2分别表示横向及纵向拉伸应变。
取单胞厚度为t,纤维编织布占单胞的体积分数可表示为
根据混合定理,可以得到平面斜交编织复合材料的双向弹性模量为
步骤四、根据组分材料(纤维束和基体)的破坏强度,采用最大拉应力准则和畸变能准则分别作为纤维束和基体破坏判据,并利用变形协调条件,分析平面斜交编织复合材料的双向拉伸强度,得到其解析解。
在已知单胞的内力分布后,对危险点进行破坏强度计算,并根据纤维束与基体的强度极限计算平面斜交编织复合材料的拉伸强度。采用最大拉应力准则和畸变能准则分别作为纤维束和基体破坏判据,即
当平面斜交编织复合材料横向和纵向上有拉伸外载T1和T2作用时,纤维编织布和基体各自承担一部分外载Tf,Tm,即
由于单胞内纤维束与基体的变形协调,在已知两者模量的基础上,可得
联立方程(19)和(20),可得
令T2为定值,研究T1满足强度准则的最大值,即横向极限载荷。根据强度准则,基体和纤维编织布存在极限载荷[Ff1]和[Fm1],即
对于平面斜交编织复合材料的横向极限载荷[F1],横向拉伸强度σx的解析解变为
根据前面单胞的内力分析以及强度准则,可分别求解[Ff1]和[Fm1]。
(1)求解[Ff1]
由于纤维编织布可看作由弯曲梁组成的编织结构,其梁内正应力由弯矩与轴力共同作用引起,因此,
因为编织交叉处(x=0或x=2L)变形最大,强度最值得考虑,则
根据第一强度理论,令σb1=σmax,则
根据上式,可得到纤维编织布横向及纵向的极限载荷[Ff1]:
(2)求解[Fm1]
基体上的双向拉伸应力可由名义应力法得到,即
取基体内部内的任一点的应力状态为双向拉应力状态,根据第四强度理论,可得
求解方程(31),可得[Fm1]。
将[Ff1]和[Fm1]的计算结果代入式(23)和式(24),可得到平面斜交编织复合材料的横向拉伸强度。同理,按照此方法可以得到纵向拉伸强度。
附录A
J34=J1+J15+J31 2(J2+J16)+J31(J3+J19)+J30 2J14+J30J17+J30J31J18
J35=J4+J7+J10+J20+J22+J26+J32 2(J5+J23)+J32(J6+J12+J25+J28)+J33 2(J8+J21)+J33(J9+J11+J24+J27)+J32J33(J13+J29)。