复合材料结构强度⑤ 应力分析及强度校核
本章是全书的篇幅重心，也是复合材料结构强度分析的方法主干。它从能量原理出发搭起有限元基础，再用经典层板理论把层压板写成一个可计算的整体，最后落到强度校核：当量强度、蔡-吴二阶张量理论与"层压板效应"。读完应能讲清一条链路——节点位移如何变成单层应力，单层应力又如何变成一张层板的破坏结论。 一、能量原理：有限元方程的理论根
一切从功与能的关系开始。在保守系统中，外力所做的功全部转为应变能，而内力功为负，所以应变能恒为正值。写成矩阵形式，应变能就是位移的二次型： 由此得到两条主线。其一，卡氏第一定理：节点力等于应变能对相应位移的偏导数。 上式正是位移法（刚度法）的基本方程，也是后续有限元列式的起点。 由求导次序可交换立刻推出刚度互等定理，总刚的对称性不是假设，而是能量函数光滑性的必然结果。 势能一侧同样清楚：总势能等于内力势能与外力势能之和，一维弹簧的算例最能说明问题——势能取驻值（通常取极小）时结构处于平衡，这就是最小势能原理。原书进一步指出，势能极值原理、虚功原理与平衡方程三者可以互相导出，对简单问题近乎同义反复，但面对复杂结构，不同提法的解题效率差别很大。 二、有限元基础：从单元刚度到总刚
冯康先生的概括很精炼：有限元法的数学基础是变分原理加分割近似。方砖可砌圆井，直锯可锯弯板，把连续曲线分段取直得到的折线，分割越细逼真度越高——这就是分片插值。整体能量泛函因此被简化为单元能量的累加，无限自由度问题离散为有限自由度的多元二次函数极值问题，最终等价于线性代数方程组求解。 在单元自身坐标系中写出位移函数，由单元节点位移反求待定系数，得到形状函数N与几何矩阵B=C[x̄i]⁻¹； 由 {ε}=B{δe}、{σ}=[D]{ε} 得到应力应变表达式，进而写出单元势能，用最小势能原理建立单元平衡方程； 由 [Ke]=∫Bᵀ[D]BdΩ 得到单元刚度，杆单元则直接化为 [K′e]=(EA/L)[[1,−1],[−1,1]]； 在自身坐标系得到 [K′e]后，用变换矩阵统一到总体坐标：[Ke]=[T]ᵀ[K′e][T]，节点力同步变为 {F}=[T]ᵀ[K′e][T]{u}； 各单元刚度按节点自由度编号"对号入座"扩阶、重排后叠加为总刚 [K]，这种早期做法被称为直接刚度法。 注意一个热应力的细节：含热膨胀的总应变为弹性应变与热应变之和，本构关系须写成 {σ} = [G]·({ε} − ΔT·{α})，温度载荷按等效节点力加入，而不是当作一种新的约束。 若干工程问题值得记住。插值函数分为Lagrange型与Hermite型，必须满足协调性，例如任意四边形不能选用双线性插值；总刚具有三带特性，大量零元素应去零存放；总刚各行（列）之和为零，因此奇异，必须删去约束点对应的行与列。气动载荷则按面积元或广义插值面样条函数折算到节点，网格粗细直接影响求解精度，计算结果必须做可信性确认，不能只看云图颜色。 三、单元选型与网格畸变控制
模型先要选对单元。杆/梁元用于承轴力与弯曲的加筋，板元承担面内与弯曲载荷，体元用于厚截面与连接区，弹簧元用来替代连接的柔性，刚体元（RBAR/RBE2/RBE3/RSPLINE/RSSCON）用于不协调部位与网格过渡连接。 板元中QUAD4与QUAD8的选用有明确的曲率代价：圆筒曲面用QUAD4时一般占5°-10° 圆心角，用QUAD8时可放大到10°-25°。 网格畸变会直接污染刚度矩阵，必须守住形状限制：长宽比一般 < 4:1；后掠不超过20°-30°；梯度小于20°-30°；翘曲h/a控制在5%之内。高应力区或应力梯度大的区域要加密。 节点协调靠过渡单元或端部补点实现：板-板、板-梁、梁-梁各有对应做法，梁与板刚度不匹配时用加强板过渡，加筋板按筋条与蒙皮刚度叠加求等效刚度，疏密网格交界可用RSPLINE元协调。 四、层压板应力分析：从ABD矩阵到单层应力
复合材料结构应力分析走进5.2节，主体是经典层板理论加薄板理论，把层压板当作均质各向异性板单元处理。单层板是正交异性材料，正应力不产生剪应变，因此必须先建立单层的应力-应变关系，再逐层讨论层压板的本构关系。 层压板的变形用中面应变 {ε0} 与曲率 {κ} 描述，任意点的应变由二者线性组合得到；层压板内力则由各层应力沿厚度积分求得，且面内力互等Nxy=Nyx。把逐层积分结果写成矩阵，就得到整套方法的骨架： A、B、D的物理意义决定了铺层设计的自由度：B矩阵非零意味着拉剪、拉弯耦合。对称均衡铺层的B为零，可由ABD矩阵反求层压板工程常数；非均衡铺层则必须正视耦合效应，热膨胀系数也按各层Q̄、α 与 [A]⁻¹ 求和得到。 值得说明的是，B的存在是"层压板效应"的代数根源——它把面内载荷与弯曲变形绑在一起，也是逐层破坏顺序之所以复杂的力学原因。 五、强度校核：当量强度、蔡-吴准则与二次破坏
强度指结构的最大承载能力，工程上更要用设计许用值代替材料强度。分析分两级：总体强度给出内力与安全余量MS，局部强度补齐总体模型无法细化的细节——钉载、胶接、圆角、自由边。建议局部强度计算数据库化、自动化，并自动生成计算报告。 当量强度的思路是把复杂受力折算为便于表征的应力：连接件破坏时的挤压应力Fbr=P/dh称为当量强度；胶接强度则由单搭接试件的破坏拉力算出"当量剪切强度"。它对应力状态的简化是显式的，使用时要与控制该值的试验件构型一致。剩余强度则定义为损伤后强度与完好强度之比，是损伤容限评价的核心指标。 层压板破坏用"逐层破坏"理论：对每一层单独判断失效，失效判据是蔡-吴二阶张量多项式理论。它由五个基本强度参数定出强度张量系数——纵向拉伸Xt、纵向压缩Xc、横向拉伸Yt、横向压缩Yc与面内剪切S，F1、F2、F11、F22、F66全部由这五个试验确定。破坏标注取f<1为未破坏、f≥1为破坏发生。 工程上更关心裕度，于是引入强度/应力比 R：把应力按比例放大R倍后代入判据，得到一个关于R的一元二次方程，取正根即为该层的强度比，R同时在破坏包线上给出破坏模式。应变空间中的破坏包线只是把应力准则代入应力-应变关系后的重新描述，本质仍是平面应力准则，不是平面应变准则。 层压板与金属最大的差别在于破坏是渐进的，因此层压板效应必须显式建模。层板强度试验的载荷-位移曲线多为线性且突然破坏，据此提出二次破坏分析：逐层卸载已破坏单层的刚度，重新分配载荷，循环计算直至层板丧失承载能力。单步实施顺序是：逐层计算应力并代入蔡-吴判据；找出最先破坏的层，按其失效模式退化刚度（现场单层强度还要按层厚与相邻层夹角修正，以体现层间约束）；重新求解内力分配，跳过已破坏层继续判断余下各层；统计破坏层数占总层数的比例 Ψ，与特征数 Ψc比较。 判据的核心不是"某一层坏了"，而是比例 Ψ。当 Ψ 超过特征数 Ψc时判定层板失效。下表对比三种判据的差别： 特征数 Ψc是层压板形式与外载（Nx、Ny、Nxy）的函数，应按铺层比与载荷工况分别近似确定，取值应小于0.3。把它当作常数套用会系统性高估强度。 T300/双马层板的验证算例给出了完整对比：12种铺层（另一轮9种）的计算值与试验值中除个别铺层外吻合较好；本方法所得强度高于二阶张量逐层破坏法，原因是计及了层间约束。程序LAMSTRENG.EXE可输出层板的纵向与横向拉压强度、剪切强度以及联合载荷下的承载能力，例如某联合加载工况下x向应力达到255.6MPa即判破坏，并给出明确的破坏或不破坏结论——这正是强度校核报告最需要的一句话。 这期的分享就到这里啦！如果内容对你有帮助，欢迎点赞、转发、收藏。你们的每一次鼓励，都是我持续输出的动力。感谢大家的支持，我们下期见！ #复合材料#有限元分析#应力分析#强度校核#蔡-吴准则#层压板 【免责声明】本文系对《飞机复合材料结构强度分析》（陈业标、汪海、陈秀华编著，上海交通大学出版社，大飞机出版工程，2011年12月第1版，ISBN 978-7-313-07882-7）相关章节的解读与重新组织，非原书转载，仅供技术交流与学习参考；文中观点、数据与公式以原书为准。版权归原作者与出版社所有，所述观点为笔者对资料的个人理解与归纳，不代表原作立场；若认为引用不当或涉版权问题，请及时告知，我们将在核实后第一时间处理。感谢理解与支持！ 著作权归作者所有，欢迎分享，未经许可，不得转载
首次发布时间：2026-10-15
最近编辑：53分钟前
作者推荐
