0 引 言

复合材料应用从飞行器次承力结构向主承力结构发展，这对其设计技术提出了更高的挑战，要求设计部门能够全面掌握其在各种服役环境下满足性能要求的设计方法。在飞行器多次使用的任务要求下，飞行器复合材料结构设计应满足损伤容限与耐久性要求。

目前金属结构已经有定量的、可以作为设计依据的损伤容限分析方法，常用的断裂判据如：应力强度因子K判据、能量释放率G判据、J积分判据[1]等，成熟的疲劳裂纹扩展速率方程如 Pairs模型、Walker模型、Forman模型[1]等。而复合材料结构损伤容限设计在目前还没有普遍适用的计算方法，设计人员通常在细节设计及设计概念上进行定性设计，通过一系列的试验确定设计许用值，以验证性试验证明满足损伤容限耐久性要求[2]。

高昂的试验代价将严重制约复合材料结构的高效设计，因此，国内外不断完善对复合材料结构的损伤容限分析方法和技术，在有效指导损伤容限设计的同时最大程度减少试验成本。

复合材料结构损伤容限分析主要针对低速冲击损伤分析。层合板受冲击损伤后可以被等效为含圆（椭圆）孔或软化夹杂的含缺陷层合板。许希武等[3]提出了弹性核模型，由子层的屈曲应力和压缩破坏应力比值确定弹性核的材料性能降，从而导致载荷的重新分配，重新分配的载荷在损伤区域边缘产生应力集中。庄茁等[4]提出复合材料层合板在低速冲击损伤下主要存在纤维断裂、基体开裂、基体挤裂、分层4种损伤模式，并给出每种损伤模式的损伤准则以及相应的刚度折减方案，将这个损伤分析模型写成Abaqus子程序并实现应用。Chang等[5]较早提出了含孔复合材料层合板的二维逐渐损伤模型，并且分析了拉伸载荷作用下含孔层合板的破坏过程，失效分析中采用了 Hashin准则及Yamada-Sun准则。

本文将层合板的冲击损伤区域假设为圆孔，圆孔的半径为最大损伤尺寸，对冲击后含应力集中的复合材料层合板进行剩余强度分析方法研究，分别采用平均应力准则和有限元方法对带孔复合材料平板的破坏应力进行了分析，并与试验结果进行了比对。

1 带孔复合材料平板结构分析模型

本文研究对象尺寸如图1所示，平板为T300/5208层压板，单层性能如表 1所示，平板铺层顺序为[0/±45/90]2s。

图1 研究对象尺寸 Fig.1 Research Object Size

表1 带孔复合材料平板单层性能[7] Tab.1 Lamina Properties of the Composite Laminates with a Hole Open

变量 参数 变量 参数E1/GPa 181.00 G13/GPa 7.17 E2/GPa 10.30 G23/GPa 3.78 E3/GPa 10.30 Xt/MPa 1500.00 υ12 0.28 Xc/MPa 1500.00 υ13 0.28 Yt/MPa 40.00 υ23 0.30 Yc/MPa 246.00 G12/GPa 7.17 S/MPa 68.00

承受均匀拉应力带孔复合材料平板有限元模型采用参数化对称模型计算，如图2所示。

“国……国亡了！我……我也……老了！你们还年青，你们去救国吧！我的老骨头再……再也不中用了！我是个老亡国奴，我不会眼见你们把日本旗撕碎，等着我埋在坟里……也要把中国旗子插在坟顶，我是中国人！我要中国旗子。我不当亡国奴，生是中国人，死是中国鬼……不……不是亡……亡国奴……“

图2 带孔复合材料平板有限元模型 Fig.2 Finite Element Models of the Composite Laminates with a Hole Open

2 含应力集中复合材料结构损伤容限分析

2.1 平均应力准则

Whitney-Nuismer平均应力准则[6]认为，孔边界某一特征距离处的平均应力达到无缺口强度 σ0时，断裂发生。平均应力准则如图3所示，平均应力准则的特征长度a0的大小与材料和铺层情况有关，实践中a0的值通过试验根据强度降低数据确定。a为椭圆孔的长轴半径，b为短轴半径，当a=b=R时，平均应力准则表示含圆孔平板的失效准则。

图3 平均应力准则 Fig.3 Whitney-Nuismer Average Stress Criterion

在无限远处施加平行于 y轴的均匀应力，在孔边沿x轴的正应力yσ近似为

式中 TK为正交各向异性应力集中系数。

式中 xE，yE，xyυ和xyG为层合板的各项等效模量。

工程用复合材料层合板一般为对称均衡铺层，Ex，Ey，υxy和Gxy是板的面内特性，根据正交各向异性板的面内刚度矩阵A计算：

未采用刚度降级的线性分析，取各层的最大失效因子，计算得到层合板破坏应力。

表4和图5以试验数据为基准，对这3种计算方法的精度进行比较。平均应力准则计算结果与试验结果相比，最大误差为8.7%，几乎精确估计了带孔复合材料平板的破坏应力。未采用刚度降级的分析方法没有考虑单层初始破坏造成损伤后结构的承载能力，计算结果很保守，计算得到的破坏应力为试验结果的50%～59%。采用刚度降级的非线性分析结果比试验结果高10%～19%，该分析方法基于经典的层合板理论考虑了损伤发生后理想的结构刚度降级，没有考虑到层间剪切应力等因素对刚度降级产生的影响，由此导致计算结果偏高。

北京市人大常委会相关负责人表示,条例实施以来,北京市人大常委会每年都确定一个重点难点问题,持续开展监督。聚焦居家养老健康服务和医养结合工作,是今年北京市人大常委会确定的重点。

其中，

2.2 有限元方法

较之传统节庆表演的固定形式化、民俗活动较单一的内容选择，社区教育能够整合调度多种教育资源，吸收各级政府、公共文化机构、学校、商业组织等多方的聚合力量，充实非物质文化遗产传承与创新活动的实施内容、形式和手段，更大范围地扩大社会凝聚力和影响力。

文章编号：1004-7182(2018)02-0105-05DOI：10.7654/j.issn.1004-7182.20180221

CGGA 全转录组测序的 mRNA 抽提和质检步骤同上，测序平台为 Illumina HiSeq platform（美国 Illumina 公司）。

式中 X为材料主方向1方向许用应力；Y为2方向许用应力；Z为 3方向许用应力；S12为平面内剪切方向许用应力；S13为13方向剪切许用应力；S23为23方向剪切许用应力；F为失效因子系数。各许用应力值参见表1，本文中失效因子系数F取1。

摊铺过程的质量控制应首先从摊铺设备抓起，如果条件许可，应尽量选用性能稳定、技术先进的摊铺机，全套设备尽量使用相同品牌。面层施工必须全幅摊铺，采用浮动基准梁控制摊铺厚度和平整度。当路面滞水时，应暂停施工；气温低于10℃时，应停止摊铺；气温低于15℃时，应控制沥青混凝土的最大运距，以保证碾压温度在规定的范围之内；未经压实即遭雨淋的沥青混凝土应全部清除，更换新料。

当积分点单层应力状态满足失效准则后，该层材料模量降低一定百分比，本文中模量剩余百分比为1%。具体的刚度降级方法与失效准则相对应，Hill失效准则只包含一个失效因子，所有模量都被降低，如表2所示， dQ为降阶后的刚度，Q为降阶前的刚度。

复合材料层合板的受力特性与各单层密切相关。但一层甚至几层单层板的破坏，虽然会引起层合板刚度的显著变化，但是层合板仍可能由余下的各个单层板来承受更大的载荷，直到全部单层破坏导致层合板的总体破坏[7]。将层合板刚度降级引入有限元方法，具体过程为，每一个载荷步更新刚度矩阵，判断积分点应力状态是否满足失效准则，如果满足则根据相应的刚度降准则改变刚度矩阵，直到计算程序终止。

表2 刚度降级准则 Tab.2 Stiffness Reduction Criterion

失效模式 刚度降级准则1方向应力达到许用值 ( )Q QQEEGυ d = =0.01 , , ,11 22 12 12 2方向应力达到许用值 ( )Q QQEEGυ d = =0.01 , , ,11 22 12 12面内剪切应力达到许用值 ( )Q QQEEGυ d = =0.01 , , ,11 22 12 12

3 分析方法比较

本文对 3种不同的分析方法进行比较，分别为：Whitney-Nuismer平均应力准则法、未采用刚度降级的有限元方法、采用刚度降级的有限元方法。

本文选用 Hill失效准则[8]，该准则计算失效因子FI，强度比 1=SR ，失效准则具体为

在Whitney-Nuismer平均应力准则法中，从文献[9]可得到层合板的平均应力准则特征长度a0=3.8 mm，由式（1）～（5）得到各带孔复合材料平板的破坏应力。

式中 h为层合板厚度。

通过采用刚度降级的非线性分析得到各带孔复合材料平板的力-位移曲线，如图4所示，分析曲线得到层合板的破坏应力。

由表2可以看出,复合混沌序列有均衡的0-1比,并且满足E<0.02的要求,满足Golomb伪随机假设的第一条件。

图4 各带孔复合材料平板的力-位移曲线 Fig.4 Force-deformation Curve of Each Composite Laminates

续图4

3种方法破坏应力计算结果如表3所示，试验结果借用文献[9]中试验数据，对各分析方法所得的数据进行比较。

表3 各方法破坏应力计算结果 Tab.3 Ultimate Stress Analysed by Different Methods

孔半径mm试验结果MPa平均应力准则MPa未采用刚度降级MPa采用刚度降级MPa 1.25 355.68 384.21 168.70 395.83 3.25 306.28 305.65 128.59 364.58 3.75 296.4 293.82 123.41 354.17 7.5 266.76 243.33 109.05 312.50 12.5 227.24 216.06 95.04 250.00

即当时，产生破坏，得到破坏应力：

表4 各计算方法结果与试验结果的比较误差 Tab.4 Comparation of Errors between Analysis Results and Test Results

孔半径/mm平均应力准则 未采用刚度降级 采用刚度降级1.25 8.02% -52.57% 11.29%3.25 -0.21% -58.02% 19.03%3.75 -0.87% -58.36% 19.49%7.5 -8.78% -59.12% 17.15%12.5 -4.92% -58.18% 10.02%

图5 分析结果与试验结果的对比 Fig.5 Comparation of Analysis Results and Test Results

4 工程应用

对航天器一块主承力的复合材料壁板拟开孔造成的强度下降进行预估，该壁板及拟开孔情况如图 6所示，拟开孔半径2.5 mm。

图6 复合材料壁板及拟开孔情况 Fig.6 Composite Plate to be Opened a Hole

该壁板为有曲率壁板，曲率大小为 10-3，同时拟开孔尺寸与壁板的最短尺寸相比不足3%，该壁板主要受载方式为纵向受拉。因此可以将该壁板拟开孔导致强度降低问题简化为带孔复合材料平板的强度评估，应用 Whitney-Nuismer平均应力准则和刚度降有限元分析方法分别进行计算。拟开孔前后的破坏应力如表5所示，开孔前破坏应力由材料手册中相同铺层层合板性能得到。

两组患者术前SAI评分差异无统计学意义（P＞0.05）；镇静组患者术中SAI评分程度低，但差异无统计学意义（P=0.05）；镇静组患者对麻醉的满意度较高，差异有统计学意义（P＜0.05，表2）。

表5 拟开孔前后的破坏应力 Tab.5 Ultimate Stress before and after Opening a Hole

开孔后破坏应力/MPa开孔前破坏应力/MPa平均应力准则 刚度降有限元分析483 321 393

根据第3节各分析方法与试验结果的误差范围，预测开孔后的实际强度极限范围为330～352 MPa。此结果为进一步分析开孔的可行性和结构设计提供了依据，应用实践表明分析结果正确有效。

5 结 论

本文对复合材料结构的损伤容限分析方法进行了研究并应用于工程实例。同时得到了以下结论，可以用于飞行器结构损伤容限设计。

a）从本文分析结果得到，平均应力准则可以精确估计孔所导致的强度下降，可以用于航空航天复合材料结构损伤容限分析；

b）未考虑刚度降级过程的有限元分析方法过于保守，不能用于结构减重设计；

推荐理由:本书由陈志武、韦森、高全喜、赖建诚诚挚推荐。与伟大的经济学家同行，洞悉经济学的本质与规律。本书既是一部经济史，又是一部思想史。

c）为得到更高精度的有限元分析结果，需要综合层间剪切应力、残余应力、工艺偏差等影响对刚度降级过程进行修正。

参 考 文 献

[1] 殷之平. 结构疲劳与断裂[M]. 西安: 西北工业大学出版社, 2012.

Yin Zhiping. Structural fatigue and fracture[M]. Xi`an: Northwestern Polytechnical University Press, 2012.

[2] 陈业标, 等. 飞机复合材料结构强度分析[M]. 上海: 上海交通大学出版社, 2011.

Chen Yebiao, et al. Structural strength analysis of aircraft composite materials[M]. Shanghai: Shanghai Jiao Tong University Press, 2011.

[3] 毛春见, 许希武, 林智育. 含冲击损伤复合材料层板及加筋壁板剩余强度研究[J]. 应用力学学报, 2010, 27(2): 286-292.

Mao Chunjian, Xu Xiwu, Lin Zhiyu. A research for residual strength of stiffened and non-stiffened composite panels with impact damage[J].Chinese Journal of Applied Mechanics, 2010, 27(2): 286-292.

[4] 庄茁, 等. 复合材料低速冲击损伤研究及等效模型的应用[J]. 工程力学,2012, 29(增Ⅱ): 15-22.

Zhuang Zhuo, et al. Damage of composite under low-velocity impact and application of mechanical equivalent model[J]. Engineering Mechanics,2012, 29(S2): 15-22.

[5] Chang F, Chang K. A progressive damage modal for laminated composites containing stress concentrations[J]. J. Comps Mater, 1987(21): 834-855.

[6] Nuismer R J, Whitney J M. Uniaxial failure of composite laminates containing stress concentrations[R]. ASTM STP 593, 1975.

[7] 沈观林, 等. 复合材料力学[M]. 北京: 清华大学版社, 2006.

Shen Guanlin, et al. Mechanics of composite materials[M]. Beijing:Tsinghua University Press, 2006.

[8] 陈火红, 等. MSC Patran/Marc培训教程和实例[M]. 北京: 科学出版社,2004.

Chen Huohong, et al. MSC Patran/Marc training courses and examples[M].Beijing: Science Press, 2004.

[9] Nuismer R J, Lzbor J D. Application of the average stress failure criterion,Part 1: tension[J]. Journal of Composite Materials, 1978(12): 238.