断裂力学论文2.doc

断裂力学是近半个世纪才发展起来的一门新兴学科，他是对经典连续介质的一个重要贡献。断裂力学主要研究带裂纹固定的强度和传播的规律。它的主要任务是研究裂纹尖端在载荷作用下的扩展规律，了解带裂纹体的承载能力，从而提出抗裂纹的设计方法，以保证构件的安全工作。 断裂力学产生于人们对各种工程断裂事故的思考，为了避免断裂事故，人们与之进行了长期的，艰苦的和卓有成效的斗争，起初凭经验，后来发展成为理论。在断裂力学出现以前，传统的控制构件不发生断裂而能够安全工作的理论，称为强度条件或安全设计，其基本思想是保证构件的工作应力不超过材料的许用应力。即： 安全设计对确保构件安全工作起了重大作用，至今仍然是必不可少的，但人们在长期的生产实践中逐步认识到在某种情况下，安全设计出的构件并不安全，断裂事故仍不断发生，特别是对于高强度材料构件，焊接构件，处在低温。或腐蚀环境中的结构等，断裂事故就更加频繁，例如，比利时阿尔伯运河上几座大桥的断裂，1949年东俄亥俄煤气公司的圆柱形液体天然气罐爆炸使周围街市变为废墟。这些接连不断的工程断裂事故引起了人们的高度警觉，这些事故发生在工作应力低于材料屈服极限的条件下，用传统的安全设计观点是无法解释的。从大量的工程事故中分析，断裂接起源于构件有缺陷，传统的设计思想的一个严重是把材料视为无缺陷的均匀连续体，而实际上构件总是存在着形式不同的缺陷，因而实际材料的强度大大低于理论模型的强度。断裂力学正好弥补了传统设计思想的不足。 断裂力学新进展： 目前常规的有限元模拟裂纹扩展问题有网格重构法，打开节点法，消去单元法，断裂过程区法（CZM）等，而这些方法受制于网格的划分，有限元模拟不连续（如裂纹）的另一个思路就是在单元常规的连续位移上迭加上有应变间断的位移项，在一个单元范围内构建应变不连续面，因此不连续面可以不受制于网格的限制。 扩展有限元(XFEM）是目前分析不连续力学问题最有吸引力的数值方法，该方法通过在单元常规的连续位移场上迭加上不连续的间断附加项和列间的应力奇异性，但目前的XFEM要引入附加自由度，系统总的自由度随裂纹扩展而不断增加，对程序编写中的自由度编程造成一定的困难，而更为关键的是构造裂纹尖端附件的弹塑性场非常困难，因此目前该方法只限于弹性分析。现在有的商业软件虽已加入了XFEM模块，但还只能针对各向同性线弹性材料，裂纹扩展用的的是类似CZM（裂纹过程区法）准则。还不能很好的表现裂纹尖端奇异性，不能使用应力强度因子准则。此外，裂纹每次扩展都必须穿过一个完整的单元。有网敏感性，尚不能模拟裂纹大偏角扩展和多裂纹交汇的情况。 问题提出与选模分析： 1，单元子划分和子结构 本文方法中，将结构中的单元分为三种，不含裂纹单元即普通单元；包含裂尖的单元；被裂纹贯穿的单元需要进行自划分，如图，斜阴影单元为包含裂尖的单元，竖阴影单元为裂尖 贯穿单元，其他为普通单元单元节点上的空心圆圈为初始网格的结点，单元内部和边界上的是实心点为子划分得到的内部结点，虚线为子划分单元的边。 为了反映裂尖的奇异性包含裂尖的子划分退化为三角形的六结点奇异单元，对于界面裂纹，裂尖在两个单元的界面上，此时需要对两个单元同时进行子划分（6和7）对于被裂纹贯穿的单元，有三种可能的子划分形式：裂纹贯穿单元的对边，此时原单元自然子划分为两个子单元（2）;裂纹穿过单元的相邻边，原单元子可划分为三个子单元（3）当裂纹遇到界面时，裂纹可能沿单元边偏转，此时原单元也可以子划分为两个子单元（4），对于贯穿的裂纹单元（5），也可以作为单元2子划分的特例，可令其中退化为一条线的子划分单元刚度矩阵的元素退化为零。 覆盖一条裂纹的所有子划分单（上图中的浅、深灰色单元）可以组成一个子结构，子结构内部结点和子结构边界上的子划分结点（上图中单元4边上粗线条正方形表示的点）的自由度，可以凝聚到子结构的表面结点（上图中粗线条圆圈）的出口自由度。子结构的平衡方程为： 为子结构的刚度矩阵、总结点位移和总的右端项，展开写为： 其中 分别为出口节点位移，内部节点位移，边界节点位移， 为刚度矩阵中与上述各位移相对应的项， 为交叉相， 分别为出口结点荷载矢量、内部节点荷载矢量和边界结点荷载矢量，当裂纹表面不受载荷时有: 子结构边界节点（上图粗线条正方向表示的点）同时也是子结构外的边界上的结点，因此，其位移可以由所加他的两个出口结点的位移线性插值得到： 将(3)带入（2）的： 将子结构的内部结点和边界结点的自由度凝聚到出口结点自由度之后，子结构的刚度矩阵为：