9-第九章-层间应力解析

第9章层间应力经典层合板的理论假定了所有单层处于平面应力状态，不考虑层间应力（面外应力）分量，层合板的应力分析比较简单。但无论是机械加载还是湿热加载，层合板中都会产生层间应力，尤其是在板边缘附近层间应力分布复杂，变化梯度大。层间应力往往引起层合板边缘脱黏，形成层间裂纹，造成整个层合板的刚度和强度下降，使结构过早失效。经典层合板的理论不能完全确定引起复合材料破坏的应力，它无法解决层间应力这类三维各向异性弹性力学问题。本章主要介绍层合板产生层间应力的原因和基于弹性力学的一些层间应力的分析方法。9.1层间应力的定性分析层合板一般由不同铺设方向的单层组成，各单层的弹性性能不同，受力下的变形也不同。但是层合板中的各单层相互黏结成一体，层和层之间变形相互制约和协调，于是在层间产生相应的正应力和剪应力，即层间应力。以下通过对层合板和[0/90]层合板的拉剪耦合变形协调和泊松耦合变形协调分析以及力矩平衡原理，解释层间应力产生的原因。一、拉剪耦合变形协调引起的层间剪应力一块材料和厚度相同的斜交铺设层合板，承受拉伸应力x，若将单层分别考虑，各层除了有线应变x和y外，由于拉剪耦合效应的存在，还会出现剪切应变xy，如图9.1所示。由式（3.26）可知，+和–层的12Q相差一个负号，因此–层xy1与+层的xy2大小相等，方向相反。黏合在一起的层合板的A16等于零，无拉剪耦合，层合板无剪切变形，因此层的剪切变形必须相互协调为零。这一效果是通过各单层相互对对方施加剪应力来实现的，又因为在层合板的自由边缘无剪力，因此只能靠层间剪切应力zx来提供。相等，因此两个相反方向铺层的正应变x，y相同，+和–层的16Q图9.1拉剪耦合变形协调引起层间剪应力示意图同理，一块材料和厚度相同的斜交铺设层合板，只承受剪切应力xy，由于剪拉耦合效应的存在，+和–层中会出现耦合线应变x和y，黏合在一起的层合板无轴向变形，为协调轴向变形，在层间会产生层间剪切应力zx和zy。层间应力也会在正交铺设层合板中出现。设有一块[0/90]层合板承受有平均轴向拉应力,如图9.2所示。假如将0°单层和90°单层分别考虑，各自只受有x方向的正应力，且沿x方向的变形相同。由于0°单层和90°单层在y方向变形的泊松耦合效应不同，由式（3.12）可知0°单层沿y方向收缩较多，90°单层沿y方向收缩较少。为了保证两单层黏层在一起后y方向变形协调一致，就需要通过0°单层和90°单层相互施加y方向的力，将0°单层往外拉，90°单层往里压，得到相同的y方向变形，这样就会在两板中产生沿y方向的正应力y1（拉应力）和y2（压应力）。因为层合板两侧是自由边界，不能提供沿y方向的作用力，所以两层中沿y方向的内力就只能由层间相互作用来提供，这就形成了层间剪应力zy。正交铺设层合板没有拉剪耦合相应，0°单层和90°单层都不会出现面内剪切变形，因而没有剪切变形需要协调，各层之间也不会出现协调剪切变形的层间剪应力zx。二、泊松耦合变形协调引起的层间剪应力x图9.2泊松耦合变形协调引起层间剪应力示意图三、力矩平衡引起的层间正应力[0/90]层合板中的y方向的正应力y1和y2与层间剪应力zy没有作用在同一平面内，从而形成一个附加力矩，为了平衡该力矩必须产生层间正应力z。已有的研究表明，该层间正应力z在层合板靠近自由边缘处可能达到无穷大，z的分布特征是形成的力矩正好与y方向的正应力形成的力矩平衡。yzzyyzz9.2单向拉伸下对称层合板的弹性力学基本方程考虑一个有限宽度的对称层合板，选取xOy坐标面为对称面，而z轴为材料主轴，如图9.3所示。把层合板的每一个单层视为宏观匀质的各向异性体，各个单层之间存在一个理想的物理非连续界面。当在层合板的两端沿x方向承受均匀拉伸时，界面上产生层间应力。一般情况下，界面上有层间正应力和切应力三个分量：z，zx，zy。因此对各个单层要从三维应力状态出发进行弹性力学分析。图9.3有限宽度的层合板在以xOy坐标面为对称的层合板中，任一个单层中可视为以z轴为弹性主方向的单对称材料，其应力—应变关系可按角度铺设单对称材料表达为xyxzyzzyxxyxzyzzyxCCCCCCCCCCCCCCCCCCCC66362616554545443633231326232212161312110000000000000000（9.1）应变与位移的几何关系为yuxvzuxwzvywzwyvxuxyxzyzzyx（9.2）由于在层合板中各个单层的几何、弹性特性和受力形式沿x方向都是均匀分布的，即不随x而变化，因此在层合板两端承受均匀轴向拉伸作用力时，所有应力与x无关，应变也与x无关。假设层合板各单层的任意点沿x方向的应变为常数，即0x并设层合板中的其他应变分量与x无关。因此，各个铺层的位移场可表示为zyWwzyVvzyUxu,,,0（9.3）利用应变与位移的几何关系可得应变场为yUzUzVyWzWyVxyxzyzzyx0（9.4）再根据应力—应变关系式（9.1）可得应力分量为yUCzWCyVCCyUCzWCyVCCyUCzWCyVCCzyx363323013262322012161312011（9.5a）yUCzWCyVCCzUCzVyWCzUCzVyWCxyxzyz66362601655454544（9.5b）不考虑体积力，并注意到所有应力分量与x无关，可把静力平衡微分方程简化为000zyzyzyzyzyzyyzxy（9.6）把式（9.5）代入式（9.6）就得到用位移表示任意一个单层的平衡微分方程为000223322442442324536244232244222222452226245362245222622552266zWCyWCzyVCCzyUCCzyWCCzVCyVCzUCyUCzyWCCzVCyVCzUCyUC（9.7）再考虑层合板的应力边界条件。在层合板的最外层，上、下两表面上边界条件为000zzyzx（9.8）在层合板的两侧自由边缘上边界条件为000yyzyx（9.9）另外，还要考虑相邻两个单层在其界面上满足静力和位移连续条件，即111kzkzkzykzykzxkzx（9.10）111kkkkkkwwvvuu（9.11）从单层的位移控制微分方程式（9.7）可知，它仅包含对y和z两个坐标的偏微分，所以此模型属于准三维力学模型。由于在一般情况下联立偏微分方程（9.7）得不到封闭解析解，因此就提出了求解层间应力问题的一些方法，主要有：①直接解法。它是利用解析法和数值法直接求解微分方程的边值问题，包括复变函数解法、级数解法、摄动法和有限差分法等。②变分解法。它是建立在变分原理基础上的近似解法，它包括有限元法、瑞利—李兹法、伽辽金法等。③混合解法。它是将直接解法与变分法结合起来求解问题，例如边界层法与瑞利—李兹法的组合等。这些近似解法比较实用，但也存在一些问题，主要是不能精确满足所有给定的边界条件和界面连续条件，对于一边界值问题，若采用不同近似解法有时将得到不同的结果。9.3斜交对称层合板层间应力的近似解法斜交对称层合板只沿x方向承受均匀拉伸时（见图9.3），各单层在空间应力状态中只有x,xy,xz，和三个应力分量起主要作用，而y,z,yz，和可以忽略。因此对这种层合板进行层间应力分析时间可假设000yzzy（9.12）这种情况下平衡微分方程式（9.6）中第2，3式自然满足，有效方程只剩第1式。由应力—应变关系式（9.5）中的第2，3，4式，可得0002454423633230132623222012zWCzVyWCyUCzWCyVCCyUCzWCyVCC（9.13）由此可解得zUCCzVyWyUCCCCCCCCzWyUCCCCCCCCyV444536222632022132312332636320331223131111（9.14）式中32233322CCCC（9.15）将式（9.15）代入应变表达式（9.4）中，可把应变分量表示如下：yUzUzUCCyUyUxyxzyzzyx44454032010（9.16）式中3622262342213231233326362323312231311111CCCCCCCCCCCCCCCC（9.17）利用以上公式可将应力分量式（9.4）表示为000445455556636426203632611616134122013212111yzzyxzxyxzUCCCCyUCCCCCCyUCCCCCC（9.18）将以上结果代入平衡微分方程式（9.6）或式（9.7）中的第1式，则得一个位移函数U(y，z)的偏微分方程为02222zUyU（9.19）式中663642624454455544CCCCCCCC（9.20）由此可见，只要从方程式（9.19）求解出适当的位移函数U，就可求得层合板的位移和应力。这就是说，把解决斜交对称层合板沿x方向承受载荷时的空间问题归结为在边界条件和连续条件式（9.8）~式（9.11）的约束下求解方程式（9.19）。为了研究层间应力的变化规律，考虑四层斜交对称层合板S45单层厚度为t1，板厚h=4t1，板宽度为2b，且取b=8t1。采用高模量石墨/环氧复合材料，弹性常数为E11=138GPa，E22=E33=14.5GPa9.5231312GGGGPa21.0231312xxz/按近似弹性理论解和有限差分解法进行数值计算，可以得到各单层的应力。图9.4表示了在单层+45°和–45°之间的层间应力随y/b的变化情况。可以看出，自由边缘效应的边界层宽度相当于层合板的厚度。by/图9.445°层间应力的变化9.4斜交对称和正交层合板层间应力的完全解对于斜交对称和正交层合板，可以用解析法得到层间应力场的完全解。按三维应力状态分析层合板，建立整体力和力矩的平衡关系，假定层间应力的表达形式，应用最小余能原理进行数值计算，确定有关参数。考虑一种对称的斜交或正交的层合板，其形状和尺寸标记如图9.5所示，两侧面为自由边界，在x方向承受均匀拉伸。1为了简化分析三维应力状态下的平衡问题，提出以下假设：①每个单层作为均质宏观各向异性体，其弹性性质可用三维有效弹性模量表征；②远离自由边缘的中间区域，用经典