孔洞裂纹问题的一种数值方法

研究孔洞与裂纹的相互作用问题,通过把适于单一裂纹的Bueckner原理扩充到含有多孔洞多裂纹的一般体系,将原问题分解为承受远处载荷不含裂纹不含孔洞的均匀问题,和在远处不承受载荷但在裂纹面上和孔洞表面上承受面力的多孔洞多裂纹问题.于是,以应力强度因子作为参量的问题可以通过考虑后者来解决,而利用笔者提出的杂交位移不连续法,这种多孔洞多裂纹问题是容易数值求解的.算例说明该数值方法对分析平面弹性介质中孔洞与裂纹的相互作用既简单又有效.     

主裂纹与微裂纹相互作用问题的有效数值计算方法

提出了平面弹性介质中主裂纹与微裂纹相互作用问题的有效数值计算 方法. 通过把适于单一裂纹的Bueckner原理扩充到含有多裂纹的一般体系，将原问题分解 为承受远处载荷不含裂纹的均匀问题，和在远处不承受载荷但在裂纹面上承受面力的多裂纹 问题. 于是，以应力强度因子作为参量的问题可以通过考虑后者(多裂纹问题)来解决，而 利用提出的杂交位移不连续法，这种多裂纹问题是容易数值求解的. 列举 Cai和 Faber为评价主裂纹与微裂纹相互作用问题的近似方法而列举的算例，说明 该数值方法对分析平面弹性介质中主裂纹与微裂纹相互作用问题既简单又非常有效.     

缺陷孔对平面裂纹应力强度因子的影响

带有裂纹和缺陷孔洞的板的问题是一个多连域的边值问题，这类问题适合用边界元法所具有的高精度特性来求解．采用子域边界元法，在平面应变的条件下对存在中心裂纹的平面板受远处拉伸和剪切裁荷的作用进行了数值分析．研究了圆形孔洞对Ⅰ型和Ⅱ型应力强度因子的影响，与有限元法进行了对比，求解结果更加精确．计算了椭圆形孔对Ⅰ型应力强度因子的影响，得到了一些有意义的结果，并对移动接触弹性体作用下的带裂纹板进行了钻孔研究．     

含单排周期微裂纹平板结构中的弹性波传播

受损伤固体中含有的微裂纹或微孔洞往往具有周期性，对含周期性缺陷结构中的弹性波分析是力学研究中的重要课题，它直接关系到结构的强度和使用寿命。目前对损伤固体中弹性波散射与透射研究结果主要是弹性动力学平面问题。1995年。Scarpetta和Sumbatyan采用解析法研究了平面波在双周期裂纹弹性介质中的传播问题。并推出显式分析结果。本文基于弹性动力学理论，分析研究了含有单排横向周期裂纹的平板中弯曲波的反射与透射问题。给出了含单排裂纹时反射波与透射波系数的数值结果。对于多排裂纹情况，可采用具有退化核第一类Fredholm积分方程方法分析求解，在求解中给出相应的无量纲数，例如无量纲波数、裂纹尺寸比等。本文分析结果可望能在工程振动控制中应用。     

界面裂纹问题中的通用权函数法

热载荷和机械载荷共同作用下复合材料中的裂纹扩展往往发生在界面处.传统求解热冲击及机械载荷共同作用下界面裂纹尖端的应力强度因子的数值方法(如有限元、边界元法等),计算工作量大、效率低.通用权函数与时间无关,运用通用权函数法可以免除对每个时刻的应力分析,计算效率可得到很大提高.本文将通用权函数法推广到求解热载荷和机械载荷共同作用下界面裂纹尖端的应力强度因子过渡过程的问题中,推导出求解平面双材料界面裂纹问题应力强度因子的通用权函数法计算格式.基于此格式,计算热载荷和机械载荷共同作用下界面裂纹尖端的应力强度因子.通过实例计算比较,表明此方法得到的结果可以达到与相互作用积分法相当的工程应用精度.最后,应用此方法研究了热障涂层受热冲击及表面力共同作用时裂纹长度以及涂层厚度对应力强度因子的影响.结果表明:在一定边界条件下,当热障涂层中存在边缘裂纹时,随着涂层厚度的增加,更容易导致裂纹的扩展和涂层的剥落.     

缝洞型多孔介质中多相流的有限体积法数值模拟

本文研究的碳酸盐岩油藏储集体属于缝洞型多孔介质.这类缝洞型多孔介质由裂缝、溶蚀孔洞和低孔隙度低渗透率的基岩组成.裂缝是空隙流体流动的主要通道;溶蚀孔洞大小从几厘米到数米不等,渗透率和孔隙度都很高,是流体主要的储集空间.由于缝洞型多孔介质空隙空间的复杂性和强非均质性,数值计算中基本控制方程的空间离散应采用非结构化网格的计算模型.本文采用有限体积法模拟缝洞型多孔介质中多相流体的流动,并给出了相应的单元中心格式有限体积法的计算公式.裂缝介质和溶洞介质中单元间多相流体的流动考虑为高速非达西流,其质量通量采用Forchheimer定律计算.非线性方程的离散选取全隐式格式,并采用Newton-Raphson迭代进行求解.通过两个二维模型注水驱油的数值模拟,验证了本文方法的有效性.     

拉伸载荷作用下单边缺陷-边裂纹应力强度因子研究

利用杂交位移不连续法研究拉伸载荷作用下矩形板中单边缺陷-边裂纹(半圆孔裂纹和半方孔裂纹)问题,给出了这三种平面弹性裂纹问题的应力强度因子的详细数值解。通过半圆孔裂纹问题和半方孔裂纹问题与单边裂纹问题的应力强度因子的比较,发现半圆孔和半方孔对单边裂纹有屏蔽影响。此外,本文的研究结果表明,杂交位移不连续法用于分析平面弹性有限体中复杂裂纹问题的应力强度因子简单且又准确。     

含中心裂纹压电体的二维精确解

解析地研究了含中心裂纹的压电体，它在无穷远处承受机电载荷，并在裂面上满足由Parton和Kudryavtsev以及Hao和Shen提出的绝对电边界条件。在平面应变假设下，给出其二维精确解，并提供了机械应变能释放率和裂尖能量释放率等数值结果。考虑工业应用范围之内常用的远场载荷时，由绝对电边界条件得出的能量释放率表现出明显的非线性特征及载荷相关性，而不是完全与电场无关，这一结论与Xu和Rajapakse在较小载荷下得到的规律不同。     

磁电弹性体中纳米孔边任意位置贯穿裂纹的解析解

本文研究了反平面机械载荷、面内电载荷和面内磁载荷作用下磁电弹材料中含有纳米尺度孔边任意位置贯穿裂纹的Ⅲ型断裂力学性能.基于Gurtin-Murdoch表面弹性理论考虑纳米缺陷(孔洞和裂纹)的表面效应,利用磁电弹理论和复变弹性理论获得了纳米缺陷表面为磁电不可通条件下磁电弹场的精确解,给出了贯穿裂纹两端裂尖的磁电弹场强因子的解析表达.所得结果与已有研究比较说明了本文方法的有效性.讨论了裂纹位置、裂纹相互作用与施加多物理场载荷对无量纲磁电弹场强因子的影响.结果表明：贯穿裂纹裂尖的无量纲磁电弹场强因子尺寸效应显著;缺陷表面效应对裂纹耦合尖端场的影响受裂纹位置的制约;无量纲磁电弹场强因子受贯穿裂纹两端的裂纹长度比与施加力电磁载荷的显著影响.     

一维六方压电准晶三角形孔边裂纹反平面问题

通过引入合适的数值保角映射,利用Stroh型公式研究一维六方压电准晶中正三角形孔边裂纹的反平面问题,给出在电非渗透边界条件下三角形孔边裂纹尖端的场强度因子和能量释放率。通过数值算例,讨论场强度因子和能量释放率随缺陷几何尺寸和力电荷载的变化规律。结果表明:随孔边裂纹长度的增加,场强度因子先急剧增加后减小,并趋于定值1,正三角形孔洞的尺寸对其影响可忽略不计;声子场和相位子场机械载荷总是促进裂纹扩展,而电位移对裂纹的扩展极大地依赖于声子场和相位子场载荷的大小。     

Automated numerical simulation of the propagation of multiple cracks in a finite plane using the distributed dislocation method

In this paper, an automated numerical simulation of the propagation of multiple cracks in a finite elastic plane by the distributed dislocation method is developed. Firstly, a solution to the problem of a two-dimensional finite elastic plane containing multiple straight cracks and kinked cracks is presented. A serial of distributed dislocations in an infinite plane are used to model all the cracks and the boundary of the finite plane. The mixed-mode stress intensity factors of all the cracks can be calculated by solving a system of singular integral equations with the Gauss–Chebyshev quadrature method. Based on the solution, the propagation of multiple cracks is modeled according to the maximum circumferential stress criterion and Paris' law. Several numerical examples are presented to show the accuracy and efficiency of this method for the simulation of multiple cracks in a 2D finite plane.     

Interacting multiple cracks in piezoelectric materials

In this paper, a method is presented to calculate the plane electro–elastic fields in piezoelectric materials with multiple cracks. The cracks may be distributed randomly in locations, orientations and sizes. In the method, each crack is treated as a continuous distributed dislocations with the density function to be determined according to the conditions of external loads and crack surfaces. Some numerical examples are given to show the interacting effect among multiple cracks.     

The stress distribution in an infinite anisotropic plate with co-linear cracks

A general solution of the plane problem of a finite number of co-linear cracks in an anisotropic material is presented. The solution is obtained by reducing the problem to four very simple Riemann-Hilbert problems. From the solution it is concluded that if the loads acting on the cracks have the resultant zero for each of the cracks, then the normal and shear stresses created on the line of the cracks are independent of the elastic constants. Expressions for the stress intensity factors are derived, and some examples are presented.     

Opening-Function Simulation of the Three-Dimensional Nonstationary Interaction of Cracks in an Elastic Body

Integral relations between three-dimensional dynamic displacements (stresses) in an infinite elastic body with arbitrarily located plane cracks and discontinuities in the displacements of the opposite crack faces are presented. The influence of opening cracks on each other is considered in the problem on crack faces loaded by pulse forces. This problem is reduced to a system of boundary integral equations of the wave-potential type in a time domain. The dynamic mode I stress intensity factors are determined for two coplanar elliptic cracks under forces in the form of the Heaviside function     

Modeling of voids/cracks and their interactions

Numerical modeling of a two-dimensional elastic body containing multiple voids/cracks to study the interaction between these defects can be significantly simplified by developing special finite elements, each containing an internal circular/elliptic hole or a slit crack. These finite elements are developed using complex potentials and the conformal mapping technique. The elements developed can be divided into two categories, namely, the semi-analytic-type and hybrid-type elements. The latter element type is an improved version of the former due to the implementation of displacement continuity along the inter-element boundary. All the proposed elements can be easily combined with the conventional displacement elements, such as isoparametric elements, to analyze the above-mentioned problems without using complicated finite element meshes. Numerical examples have been employed to illustrate the modeling of voids/cracks and their interactions. The results obtained using the semi-analytic-type elements are in good agreement with the theoretical results, and the corresponding results obtained using the hybrid-type elements show an improvement of the agreement with the theoretical results. However, the former element type is much easier to construct.     

Integrodifferential equations for plane filtration in the presence of cracks

A system of singular integr Differential equations is derived for the plane problem of steady-state filtration in a plate cut by a system of cracks. We consider an arbitrary set of cracks, and also monoperiodic and biperiodic systems of cracks, in an infinite plane. In the case of a system of infinite parallel rectilinear cracks, the general solution is obtained in explicit form-in quadratures. As an example, we find the complex potential and the formula for the output from a borehole for a linear system of tiered, flooded plates, cut by a system of rectilinear parallel cracks.     

Fissuration multiple par corrosion sous contrainte : Modélisation du comportement mécanique par une méthode d'homogénéisation

Stress Corrosion Cracking (SCC) phenomenon is characterised by initiation and propagation of surface cracks frequently multiple. In order to model the mechanical behaviour of such materials we propose, as for composites, the use of homogenization techniques. Two materials are considered, first one corresponding to the cracked external volume and the second to the internal safe material. The cracked volume is considered as a two phase material, i.e., elastic matrix containing elliptical voids. The overall behaviour of the equivalent material is obtained applying the usual homogenization rules. Comparison between simulations and experimental results is done.     

The seismic pulse in a semi-infinite medium

Summary In a foregoing paper the present author developed methods for studying the transient field from a vertical electric antenna placed in the vicinity of the plane boundary of two semi-infinite dielectric media.As the theory involved is applicable to the comparable elastodynamic pulse problem the present paper deals with the field from a buried transient longitudinal source in an elastic half space.The method appears to be relatively simple and is also applicable to the more general problem in which two elastic semi-infinite solids are separated by a plane boundary.