复合材料织物层板压溃分析（二）

波纹板压溃过程中的分层效果图

4

单元类型

层压板面内与continuum shell截面属性匹配的的单元类型也为连续壳单元，如SC8R。

为了模拟压溃过程中的分层损伤，在原始层板单元两层之间插入0厚度cohesive单元，Abaqus2016版本及以后版本可以直接在CAE中插入0厚度cohesive 单元了，不需要写程序咯，具体参见之前的文章“Abaqus 中创建零厚度cohesive单元的几种方法”，插入的零厚度cohesive单元默认是COH3D8P类型，是做水力压裂分析用的，切记要把它切换到COH3D8哈，此处不再赘述。

5

装配

创建完层压板，同时载创建一个离散刚体作为底部承载板和顶部加载板，每个刚体上均设置参考点，并赋予质量及惯性矩等属性，然后将零件装配到一起。

6

分析步类型

因为压溃过程中存在严重的刚度退化，隐式分析是几乎不收敛的，所以采用显式动力学求解器进行准静态分析，分析步类型选择“Dynamic/Explicit”，。

至于准静态分析中计算效率和计算准确度的矛盾，在之前的文章中也介绍过了，各位自行考量吧，鱼和熊掌不可兼得，有高性能计算机就尽量不用笔记本台式机。

7

接触设置

整个压溃过程中，除了上下加载板与波纹板之间有接触外，当某些层间cohesive单元或者层板单元删除后，裸露出来的新的单元面之间以及新的单元面与加载版之间都会产生新的接触，而这些内面接触默认是不考虑的，因此在设置通用接触的时候要特别注意创建内部接触。

8

边界载荷

底部支撑板固支，顶部加载板参考点上施加向下的位移或者速度，如果是位移加载的话，记得指定幅值曲线amplitude，准静态分析推荐smooth step类型哈。

9

创建任务

创建任务，提交计算就可以啦。出来的效果图就是这个样子的，条件允许的把加载速率放慢，把网格画细一点，效果会更好哈。

需要注意的是，即使采用了显式分析，如果单元删除参数不合理的话，中途也会出现因单元畸变而跳出的问题，这个自己去悟吧。

整个加载过程中的载荷位移曲线如下图所示，后期会趋于稳定，震荡有些剧烈哈，加载速度放缓慢一些会好很多哈，感兴趣的自己去尝试咯。

题外话：今年的“上纬杯”比赛题目是真的难，要求载荷不能高也不能低，对设计水平、计算能力都有很高的要求，看着这个题目自己就觉得脑袋疼，要么拿大量实验去测，要么把计算搞准确，然鹅，想算得准又是何其的难，本文及实例仅供参考，欢迎批评指正，欢迎探讨。

10

如何获取模型

以往的很多模型直接就把链接公布了，这次呢，需要大家小小付出一下。

需要上述原始CAE模型或inp的朋友，需要搜集一种公开发布的T700或T800或T1000或其他单向板的力学性能（T300除外哈，已经有很多了），要求至少有E1、E2 、NU12 、G12四个参数哈，搜集好之后，将原始pdf发管理员微信“jm19961996”，并备注材料名，然后就可以拿到下载链接了。本活动有效期截止到2019.8.10。