断裂力学数值计算方法研究.docxVIP

断裂力学数值计算方法研究

0断裂力学数值计算方法的发展

断裂力学是固体力学最新发展的一个分支。为了控制和防止结构对象的断裂破坏，研究了在设计结构中确定断裂边缘的应力场和适应性场，分析了断裂扩展的条件和规律，并在结构、机械、岩石、振动等工程领域得到广泛应用。在断裂力学中，裂缝边缘的应力场、适应场和表示断裂边缘应力场强度的应力场是一个重要的研究内容。然而，只有少数简单、特殊的故障分析和解决方案，大多数情况都可以通过数值分析方法来解决。事实上，数值计算是科学研究的三大支柱。数值计算对各种问题适应性强，应用广泛。详细描述了不同的值结果，并可以通过图像显示直观地描述力学过程。数值模拟可以是比实验更有效的。由于断裂边缘附近的应力场非常不同，因此传统的数值方法不能直接用于分析断裂力学问题的效果。因此，有必要结合断裂力学的特点，制定更有效的数值计算方法。随着断裂力学研究的深入，需要解决的问题变得越来越复杂。因此，如何建立有效的数值计算方法，已成为科学家研究的重点。随着计算机、数学和力学等学科的发展，数值计算方法不断发展，从最初的有限差分法、有限元素法到现在的无网络法、数值流形法、小波值法、非连续变形分析等。

本文将在详细分析近年来应用于断裂力学的各类数值计算方法特点的基础上,着重概述它们解决断裂力学问题的研究成果与进展,并讨论断裂力学数值计算方法今后的研究发展方向.

1断裂力学主要数值计算方法

1.1自适应有限元法

普遍认为,20世纪60年代初出现的有限元法(FiniteElementMethod)是计算力学诞生的标志.有限元法是建立在传统的Ritz法的基础上,利用变分原理导出代数方程组求解.它将连续介质离散成有限个单元来进行数值计算,通过对连续体的离散化,在每个单元上建立插值函数,从而建立整个求解区域上的函数,然后利用节点位移求出应力分量.有限元法实现了统一的计算模型、离散方法、数值求解和程序设计方法,从而能广泛地适应求解复杂结构的力学问题.所以,该方法自问世至今已得到了迅猛发展,从最初的用于结构和固体力学的计算分析不断向其他领域扩展,也成为分析断裂力学问题的首选数值方法.

由于裂纹尖端的奇异性,用普通单元求解应力强度因子往往需要过多的自由度,而自由度数要受计算机的限制,同时解的收敛性也没有保证.为此,学者们提出在裂纹尖端采用特殊的有限单元,使其位移模式反映奇异性,这样就不需要细分网格,解的精度也得到了保证.这种用于裂纹奇异性模拟的有限单元称为奇异单元,其中典型的两种是奇应变单元和等参数单元.

对于一个已取定的特定网格,由于其准确解未知,无法对其解做出可靠性判断.为此,人们提出自适应有限元法,在已求得有限元解的情况下,通过应力恢复、误差估计和新网格自动划分,进而再进行有限元求解,重复这一过程直至得到满意的有限元解.这种基于自适应分析的数值方法,在有效误差估计及初始定义可接受的误差水平的网格基础上,可自动形成高精度的数值分析模型,已经成为数值分析的热点.杨庆生和杨卫在研究断裂过程的有限元模拟时,基于自适应有限元法的一般原理,提出了一种简化的高精度和高效率有限元网格的动态重新划分策略,并已在各种裂纹问题中得以实现与应用.文献采用自适应有限元法来确定裂纹尖端塑性区.文献研究了自适应有限元模拟裂纹扩展的网格生成技术,通过修改裂纹周围单元的形状及单元间的邻接关系,实现网格动态划分对裂纹扩展进行跟踪.

随机分析是断裂力学发展的一个重要方向,也是结构可靠性评估的基础.随机有限元法在有限元法的基础上,采用随机参数来描述工程实际问题,主要研究内容包括随机变分原理、随机有限元控制方程的建立及其求解.尽管关于随机有限元法的研究起步较晚,但是它在工程实际中却有着广泛的应用,并取得了一定的成果.文献用幂律非线性随机有限元所提供的位移随机场、应力、应变随机场,进行弹、塑性断裂参数及其变化率的计算.

另外,用于分析断裂力学问题的有限元法研究成果还有:用于动态断裂问题的空—时有限元法,解析法与有限元法相结合的半解析有限元法等.弹塑性有限元法和刚塑性有限元法也有一些应用.

有限元法的不足之处在于随着计算精度要求的提高,有限元网格的划分十分困难,计算工作量十分庞大,在输入数据的准备上很费事.所以,进一步提高计算效率,降低计算机存储的要求仍然是计算力学算法今后研究的焦点.在断裂问题的数值分析中,有限元法求解裂纹扩展的局限是要么预先知道裂纹的扩展路径,要么必须重新划分网格,这给求解造成了很大的困难.

1.2基于边界元法的数值计算方法

边界元法(BoundaryElementMethod)是继有限元法之后发展起来的一种求解力学问题的数值方法.其构成包含如下三个主要部分:(1)基本解的特性及其选用.大多数基本解是前人已经得出的带奇异性的某个特殊问题的解.例如,在无限