断裂力学讲义ch10-断裂力学数值模拟_67409446导论.ppt

扩展有限元计算原理 * 应力和位移哪种匹配好？为什么？ 计算J积分的一种方法（区域积分，domain-integral） 可能是最有前景的裂纹模拟算法——扩展有限元 （Belytschko and Moes，1999） 允许裂纹穿过单元，在规则网格上可计算复杂形状裂纹。 通过适当的形函数捕捉裂纹尖端奇异场，可在粗网格上获得精确解答。 裂纹扩展时无需重新划分网格。 * 第十章 断裂力学数值模拟 为什么要进行数值模拟？ 解析解太少，不能满足工程需要（应力强度因子手册） 没有解析解！ 对于稍微复杂的断裂过程，理论解就只是近似，而实验又难以观测  问题二：如何准确高效地模拟裂纹扩展的过程。 解决方案一：在裂纹扩展路径上预先加密网格 优点：实现起来简单。 缺点：不能处理未知路径情形；在裂纹未到达处和裂纹尾岸的过密网格造成较大的计算量 问题二：如何准确高效地模拟裂纹扩展的过程。 解决方案二：将自适应网格划分技术用于裂纹扩展——网格在裂尖处自动加密和在裂纹尾岸自动弥合。 优点：随裂纹扩展合理分配自由度 缺点：网格改变的前处理仍需一些计算量；目前针对二维自适应网格划分技术有一定发展，三维还没有较成熟的方案。 问题二：如何准确高效地模拟裂纹扩展的过程。 最大的瓶颈是自适应网格划分技术，该处理过程本身有时就相当繁琐。无网格法和物质点法被一部分学者认为可以解决网格划分技术问题，但是这些新兴的计算方法还距成熟可靠地应用于工程计算和被人们接受，还有一定差距。  问题一：如何根据断裂准则准确判断裂纹何时起裂 问题二：如何准确高效地模拟裂纹扩展的过程。 能否同时解决这两个问题？虚拟内键模型（VIB=Virtual Internal Bond）提出一个新的思路。将物质的断裂准则通过原子间的相互作用自然引入连续介质的本构。Gao and Kliein Deformation Gradient Green Strain Bond Length Bond Potential Energy Density Stress VIB应用示例  虚拟内键模型VIB的不足： 目前主要还是模拟弹性或脆性断裂 泊桑比固定，各向同性材料 单元尺度如何定？（稍后讲） ABAQUS等商业有限元软件中有一种模拟裂纹扩展的方法： 当判断单元中应力超过一临界应力水平时，将该单元设为死单元，即将刚度设为零，来模拟裂纹扩展。 同学们认为该方法怎么样？  Bazant等人提出crack band approach方法最先应用于混凝土的有限元计算中 争取做到计算结果与单元尺寸无关。 1 Needleman等人的界面元模拟裂纹扩展 缺陷和不足： 两倍的结点，计算量大 计算结果对单元尺寸敏感 断裂力学的数值模拟要解决两个问题： 问题一：如何根据断裂准则准确判断裂纹何时起裂 问题二：如何准确高效地模拟裂纹扩展的过程。  问题一：如何根据断裂准则准确判断裂纹何时起裂 方法一：应力和位移匹配法来确定应力强度因子（最直接的想法） 应力匹配：， 位移匹配：， 优点：想法直接简单 缺点：需要在在裂尖过分加密网格来保证计算精度，当计算裂纹扩展时，计算量很大 除了渐近展开的主项，再多包括几项可进一步改善精度  问题一：如何根据断裂准则准确判断裂纹何时起裂 方法三：环路积分来计算J积分 利用J积分的路径无关性，考虑到在远离裂尖处的应力和位移的计算有较好精度，在那里采用环路积分来计算J积分。 优点：裂尖网格无需象应力位移匹配法那么密；可以表征线性和非线性断裂 缺点：路径上的积分精度还是不够高。 之后的改进方法会考察一个面积上的积分来提高精度(而虚拟裂纹扩展应变能法相当于计算两个面积积分之差) 问题一：如何根据断裂准则准确判断裂纹何时起裂 方法四：虚拟裂纹扩展刚度导数法来计算J积分或能量释放率 上式的推导中考虑到有限元的解满足，且恒定外载条件，所以其中的第一和第三项为零。  问题一：如何根据断裂准则准确判断裂纹何时起裂 方法四(续)：虚拟裂纹扩展刚度导数法来计算J积分或能量释放率 == 在采用刚度导数法时，可以不必对所有单元刚度求导数，只需计算一个环形域。 优点：将J积分回路的积分转化为一个环形面域的积分，将这种思想连续化并进一步推广，可以得到一个准确性更好的面域的积分，对网格在裂尖的精细化要求进一步降低。 缺点：当前所示的方法仍需数值微分，但在改进推广版中可以克服。 问题一：如何根据断裂准则准