裂纹问题的非局部弹性力学分析

求解并给出非局部弹性力学平面问题的单位集中不连续位移基本解,基于这些基本解和经典弹性力学中的不连续位移边界积分方程＿边界元方法,提出了一种非局部弹性力学平面问题的一般解法·利用该解法,研究分析了Ｇｒｉｆｉｔｈ裂纹、边缘裂纹等断裂力学中基本的但又很重要的问题·结果表明,裂纹前沿的应力集中系数与裂纹长度有关,给出了裂纹长度对断裂韧性ＫⅠｃ的影响·所得结果与已有实验结果一致·     

任意厚度具有自由边叠层板的精确解析解

自由边问题一直是三维弹性力学中的难题,通常很难满足自由边上一个正应力和两个剪应力都等于0．基于三维弹性力学基本方程和状态空间方法,引入自由边界位移函数并考虑全部弹性常数,建立了正交异性具有自由边单层和叠层板的状态方程．对状态方程中的变量以级数形式展开,通过边界条件的满足精确求解任意厚度具有自由边叠层板的位移和应力,此解满足层间应力和位移的连续条件．算例计算表明,采用引入的位移函数形式,简化了计算过程并且采用较少的级数项可以获得收敛解．与有限元方法计算结果进行了对比,可以得到较高精度的数值结果．其解可以作为其它数值方法和半解析方法的参考解．     

平面弹性裂纹分析的一种有效边界元方法

提出了一种简单而有效的平面弹性裂纹应力强度因子的边界元计算方法．该方法由Crouch与Starfield建立的常位移不连续单元和闫相桥最近提出的裂尖位移不连续单元构成A·D2在该边界元方法的实施过程中,左、右裂尖位移不连续单元分别置于裂纹的左、右裂尖处,而常位移不连续单元则分布于除了裂尖位移不连续单元占据的位置之外的整个裂纹面及其它边界．算例(如单向拉伸无限大板中心裂纹、单向拉伸无限大板中圆孔与裂纹的作用)说明平面弹性裂纹应力强度因子的边界元计算方法是非常有效的．此外,还对双轴载荷作用下有限大板中方孔分支裂纹进行了分析．这一数值结果说明平面弹性裂纹应力强度因子的边界元计算方法对有限体中复杂裂纹的有效性,可以揭示双轴载荷及裂纹体几何对应力强度因子的影响．     

两种材料组成空间的弹性力学基本解

本文研究了由两种不同材料的半空间所组成弹性体的弹性力学基本解。应用三维弹性理论中的Papkovich-Neuber通解以及Kelvin特解,求解出了在空间内部作用有集中力时空间的弹性力学位移场。该位移场在两个半空间内部分别满足各自的位移平衡方程,在其交接面上满足位移及面力的连续条件。作为本文结果的几种特殊情况,半空间的Lorentz问题与Mindlin问题的解,以及Stokes流中类似问题的解均可从该解答中导出。     

非局部非对称弹性固体理论

本文基于非局部连续统场论和非线性连续体力学理论,建立了非局部非对称弹性固体的非线性理论.它完善和发展了Eringen等人所建立的非局部弹性场论.将文献[1]中所建立的非局部非对称弹性力学的线性理论推广到有限变形.证明了在非局部弹性固体中存在着非局部体力矩.非局部体力矩引起应力的非对称性,而非局部体力矩则由原子晶格相互作用形成的共价键所产生的.并应用本文建立的理论合理地解释了平面横波和纵波色散系关的不相似性.     

弹性力学中一类新的边界积分方程

本文首先从弹性力学中的守恒积分出发,提出了弹性力学中某点的应力用边界上的位移和应力表达出来的一个积分公式。接着,依据这个公式提出了一类新型的边界积分方程。在已知位移的边界上它是第二类积分方程,在已知外力的边界上它是第一类积分方程。     

两相材料空间问题基本解的显式张量表示*

本文应用张量运算将文献中的三维两相无限体的集中力基本解表示为张量形式,从而使其能够直接用于边界积分方程和边界元方法,以分析两相材料空间弹性力学问题.本文结果包括了Mindlin问题、Lorentz问题和均质体空间问题的基本解.     

求解平片裂纹问题的有限部积分与边界元法

本文利用位移的Somigliana公式和有限部积分的概念,导出了求解三维弹性力学中的任意形状平片裂纹问题的超奇异积分方程组,进而联合使用有限部积分法与边界元法对所得方程建立了数值法.为验证本文的方法,计算了若干数值例子的裂纹面的位移间断及裂纹前沿的应力强度因子,它们与理论值相比符合很好.     

新的三维力学GELD正演和反演算法

在本文中 ,我们提出了新的整体积分和局部微分GILD的力学正演和反演方法 .我们建立了弹性和塑性力学的体积分微分方程 .我们证明了这个体积分方程和伽辽金虚功原理等价 .新的GILD方法是基于这个体积分微分方程 .GL方法是进一步的发展 ,GL是一种整体场和局部场相互作用的全新方法 .在这个方法中 ,仅仅需要解 3× 3或者 6 × 6的局部小矩阵 .特别是 ,用GL方法求解无限域的偏微分方程时 ,不需要任何人工边界 ,不需要任何吸收边界条件和不需要任何边界积分方程 .新的三维力学GILD正演和反演算法已被应用研究奈米材料的力学性质的模拟计算 .我们获得非常好的奈米材料的力学变形的超拉力的力学性质 .我们提出了新的奈米地球物理新概念和发现了GILD数值量子     

非齐次弹性力学问题双互易边界元方法研究北大核心CSCD

基于弹性力学边界元方法理论,将边界元法与双互易法结合,采用指数型基函数对非齐次项进行插值得到双互易边界积分方程.将边界积分方程离散为代数方程组,利用已知边界条件和方程特解求解方程组,得出域内位移和边界面力.指数型基函数的形状参数是由插值点最近距离的最小值决定,采用这种形状参数变化方案,分析径向基函数(RBF)插值精度以及插值稳定性.再次将指数型基函数应用到双互易边界元法中,分析双互易边界元方法下计算精度及稳定性,验证了指数型插值函数作为双互易边界元方法的径向基函数解决弹性力学域内体力项问题的有效性.     

双I—型裂纹断裂动力学问题的非局部理论解

研究了非局部理论双中I－型裂纹弹性波散射的力学问题,并利用富里叶变换使本问题的求解转换为三重积分方程的求解,进而采用新方法和利用一维非局部积分核代替二维非局部积分核来确定裂纹尖端的应力状态,这种方法就是Schmidt方法,所得结是比艾林根研究断裂静力学问题的结果准确和更加合理,克服了艾林根研究断裂静力学问题时遇到的数学困难,与经典弹性解相比,裂纹尖端不再出现物理意义下不合理的应力奇异性,并能够解释宏观裂纹与微观裂纹的力学问题。     

采用新方法研究非局部理论中Ⅰ-型裂纹的断裂问题

采用新的方法研究非局部理论中Ⅰ_型裂纹的断裂问题,进而确定裂纹尖端的应力状态,这种方法就是Schmidt方法·　所得结果比艾林根研究同样问题的结果准确和更加合理,克服了艾林根研究同样问题时遇到的数学困难·　与经典弹性解相比,裂纹尖端不再出现物理意义上不合理的应力奇异性,并能够解释宏观裂纹与微观裂纹的力学问题·     

Analysis of arbitrarily shaped planar cracks in two-dimensional hexagonal quasicrystals with thermal effects. Part II: Numerical solutions

The extended displacement discontinuity (EDD) boundary element method is developed to analyze an arbitrarily shaped planar crack in two-dimensional (2D) hexagonal quasicrystals (QCs) with thermal effects. The EDDs include the phonon and phason displacement discontinuities and the temperature discontinuity on the crack face. Green's functions for uniformly distributed EDDs over triangular and rectangular elements for 2D hexagonal QCs are derived. Employing the proposed EDD boundary element method, a rectangular crack is analyzed to verify the Green's functions by discretizing the crack with rectangular and triangular elements. Furthermore, the elliptical crack problem for 2D hexagonal QCs is investigated. Normal, tangential, and thermal loads are applied on the crack face, and the numerical results are presented graphically.     

A boundary element analysis for stress intensity factors of multiple circular arc cracks in a plane elasticity plate

In this paper, a numerical approach for analyzing interacting multiple cracks in infinite linear elastic media is presented. By extending Bueckner’s principle suited for a crack to a general system containing multiple interacting cracks, the original problem is divided into a homogeneous problem (the one without cracks) subjected to remote loads and a multiple crack problem in an unloaded body with applied tractions on the crack surfaces. Thus, the results in terms of the stress intensity factors (SIFs) can be obtained by considering the latter problem, which is analyzed easily by means of the displacement discontinuity method with crack-tip elements proposed recently by the author. Test examples are given to illustrate that the numerical approach is very accurate for analyzing interacting multiple cracks in an infinite linear elastic media under remote uniform stresses. In addition, the displacement discontinuity method with crack-tip elements is used to analyze a multiple crack problem in a finite plate. It is found that the boundary element method is also very accurate for investigating interacting multiple cracks in a finite plate. Specially, a generalization of Bueckner’s principle and the displacement discontinuity method with crack-tip elements are used to analyze multiple circular arc crack problems in infinite plate in tension (including: Two Collinear Circular Arc Cracks, Three Collinear Circular Arc Cracks, Two Parallel Circular Arc Cracks, Three Parallel Circular Arc Cracks and Two Circular Arc Cracks) in a plane elasticity plate. Many results are given.     

Analysis of arbitrarily shaped planar cracks in two-dimensional hexagonal quasicrystals with thermal effects. Part I: Theoretical solutions

An extended displacement discontinuity (EDD) boundary integral equation method is proposed for analysis of arbitrarily shaped planar cracks in two-dimensional (2D) hexagonal quasicrystals (QCs) with thermal effects. The EDDs include the phonon and phason displacement discontinuities and the temperature discontinuity on the crack surface. Green's functions for unit point EDDs in an infinite three-dimensional medium of 2D hexagonal QC are derived using the Hankel transform method. Based on the Green's functions and the superposition theorem, the EDD boundary integral equations for an arbitrarily shaped planar crack in an infinite 2D hexagonal QC body are established. Using the EDD boundary integral equation method, the asymptotic behavior along the crack front is studied and the classical singular index of 1/2 is obtained at the crack edge. The extended stress intensity factors are expressed in terms of the EDDs across crack surfaces. Finally, the energy release rate is obtained using the definitions of the stress intensity factors.     

Construction of discontinuous solutions of the equations of plane elasticity theory by the method of generalized functions

A method of constructing integral representations of discontinuous solutions of the equations of plane elasticity theory based on the use of the apparatus of the theory of generalized functions is described. The representations obtained for the discontinuous displacement and stress field components are utilized to formulate sufficient conditions ensuring continuous continuation of these quantities at almost all the points of the line of discontinuity.     

Analysis solution method for 3D planar crack problems of two-dimensional hexagonal quasicrystals with thermal effects

An analysis solution method (ASM) is proposed for analyzing arbitrarily shaped planar cracks in two-dimensional (2D) hexagonal quasicrystal (QC) media. The extended displacement discontinuity (EDD) boundary integral equations governing three-dimensional (3D) crack problems are transferred to simplified integral-differential forms by introducing some complex quantities. The proposed ASM is based on the analogy between these EDD boundary equations for 3D planar cracks problems of 2D hexagonal QCs and those in isotropic thermoelastic materials. Mixed model crack problems under combined normal, tangential and thermal loadings are considered in 2D hexagonal QC media. By virtue of ASM, the solutions to 3D planar crack problems under various types of loadings for 2D hexagonal QCs are formulated through comparison to the corresponding solutions of isotropic thermoelastic materials which have been studied intensively and extensively. As an application, analytical solutions of a penny-shaped crack subjected uniform distributed combined loadings are obtained. Especially, the analytical solutions to a penny-shaped crack subjected to the anti-symmetric uniform thermal loading are first derived for 2D hexagonal QCs. Numerical solutions obtained by EDD boundary element method provide a way to verify the validity of the presented formulation. The influences of phonon-phason coupling effect on fracture parameters of 2D hexagonal QCs are assessed.     

A method of solving three-dimensional problems of elasticity for laminated composite plates

An efficient method of solving 3D elasticity problems for thick and thin laminated composite plates is presented. It is based on a new concept of reference surfaces inside the plate. According to this concept, into each nth layer, In arbitrary reference surfaces parallel to the midsurface are introduced, and the displacement vectors of these surfaces are chosen as unknown functions. Such a choice allows one to represent the governing equations of the high-order theory of plates proposed in a very compact form and to derive strain–displacement relationships correctly describing all rigid-body motions of laminated plates.     

Numerical Simulation of Crack Growth

In recent years the X-FEM based on the partition of unity method and the strong discontinuity approach (SDA) have shown to be powerful tools to model crack growth. Both methods model the crack surface by introducing additional d.o.f.. In the X-FEM the nodes in the mesh around a crack are globally enhanced with new d.o.f. while in in the SDA the new d.o.f. are commonly introduced as internal ones. Thus the jump displacement fields are constant across elements. Therefore the d.o.f. can be condensed on element level which results in jumps in the displacement field at element edges. In this contribution the strong discontinuity approach is used approximating the displacement jump linearly across the crack length similar as e.g. in [3]. New additional nodes of the cracked elements that lie on the element edges are introduced but are not considered as internal nodes but remain global. Thus crack path continuity is automatically given. These global d.o.f. approximate the discontinuous part of the displacement field. The sum of the aforementioned part and the continuous displacement field represent the total displacement field including a possible jump. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)