论文部分内容阅读
新型复合材料在工程实际应用中正发挥着日渐重要的作用,特别是在土木工程、机械制造、航空航天和海上平台等领域得到了广泛地应用。复合材料通常由两种或两种以上具有不同性能的材料设计加工而成,具有优良的物理力学性能,可以完成许多传统材料无法承担的任务。但像层合板、编织物等复合材料由于都具有非均质性和各向异性的材料属性,因此在复杂的加载条件下将会导致其出现裂纹任意合并和分叉,从而影响结构的有效力学性能和可靠性。材料失效的高精度数值仿真一直以来是工程科学领域的一个巨大挑战。近几年来,得到快速发展和改进的增强有限元法已被证明可以精确和高效地模拟准静态非线性断裂问题。然而,当前的增强有限元法主要应用于线弹性材料和裂纹稳定扩展问题,尚未针对弹塑性材料和裂纹不稳定扩展问题开展详细的研究。为此,本文提出了一种处理弹塑性固体中任意裂纹萌生和扩展的新型弹塑性增强有限元法。同时,针对准静态模拟中出现局部或者全局数值失稳这一问题,本文提出了一种新型的基于惯性力的数值稳定性方法。基于弹塑性理论,推导得到了可以处理塑性变形和非线性断裂的一维和二维新型弹塑性增强有限元法,借助通用商业有限元软件ABAQUS,结合用户自定义单元子程序功能,开发了一维和二维弹塑性增强有限元计算程序。新型弹塑性单元引入线性等向强化模型和von Mises屈服函数描述开裂之前的弹塑性效应,结合分段线性内聚力关系描述裂纹萌生和扩展。在发生内聚力裂纹断裂之后,借助单元内部节点描述材料从连续状态到不连续状态的演化过程,并通过内聚力裂纹张开位移自洽算法把内部节点自由度在每个单元中进行凝聚。通过一系列单元测试和数值算例证明,该方法与其它方法相比,具有网格不敏感性、计算高效性、模拟精确性和数值稳定性等优势。针对裂纹快速和动态传播的断裂问题,基于惯性力效应推导得到了可以克服准静态模拟当中发生计算不收敛的数值稳定性法。该方法只需要通过逐渐缩短时间步长,并且不需要借助任何数值阻尼或其它增强求解参数,就可以无条件地保证计算结果的收敛性。通过对几个局部或者全局数值不稳定问题的数值模拟,很好地证明了该方法可以高效和无条件地通过数值不稳定点,同时在之后顺利和快速地回到裂纹稳定扩展状态。在所有的数值测试中,这个新方法的求解精度和数值稳定性均优于像数值阻尼法、弧长法、隐式动态模拟法等已有的增强求解技术。