Weight functions for interface cracks in dissimilar anisotropic materials

Bueckner‘s work conjugate integral customarily adopted for linear elastic materials is established for an interface crack in dissimilar anisotropic materials. The difficulties in separating Stroh‘s six complex arguments involved in the integral for the dissimilar materials are overcome and then the explicit function representations of the integral are given and studied in detail. It is found that the pseudo-orthogonal properties of the eigenfunction expansion form (EEF) for a crack presented previously in isotropic elastic cases, in isotopic bimaterial cases, and in orthotropic cases are also valid in the present dissimilar arbitrary anisotropic cases. The relation between Bueckner‘s work conjugate integral and the J-integral in these cases is obtained by introducing a complementary stressdisplacement state. Finally, some useful path-independent integrals and weight functions are proposed for calculating the crack tip parameters such as the stress intensity factors.     

A New Method for Establishing Pseudo Orthogonal Properties of Eigenfunction Expansion Form in Fracture Mechanics

A new and simple method is developed to establish the pseudo orthogonal properties (POP) of the eigenfunction expansion form (EEF) of crack-tip stress complex potential functions for cracked anisotropic and piezoelectric materials, respectively. Di?erent from previous research, the complex argument separation technique is not required so that cumbersome manipulations are avoided. Moreover, it is shown, di?erent from the previous research too, that the orthogonal …     

Closed-form weight function for edge crack problems

An accurate crack surface opening displacement formula for edge cracks in finite bodies has been derived, and was subsequently used to develop wide-range closed-form weight functions for edge cracks. The analytical weight function is general and can be easily used to determine accurate stress intensity factors and crack surface opening displacements for arbitrarily loaded edge cracks with high efficiency. Examples have been given for problems of edge crack(s) with a number of typical load conditions for illustration.     

Variational formula for a body with a crack and its applications (review)

A review is made of studies on the variational formula for a body with a crack, asymptotic expansions, weight functions, and displacements of crack faces __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 8, pp. 30–46, August 2007.     

Hypersingular integral equation method for a three-dimensional crack in anisotropic electro-magneto-elastic bimaterials

Using Green’s functions, the extended general displacement solutions of a three-dimensional crack problem in anisotropic electro-magneto-elastic (EME) bimaterials under extended loads are analyzed by the boundary element method. Then, the crack problem is reduced to solving a set of hypersingular integral equations (HIE) coupled with boundary integral equations. The singularity of the extended displacement discontinuities around the crack front terminating at the interface is analyzed by the main-part analysis method of HIE, and the exact analytical solutions of the extended singular stresses and extended stress intensity factors (SIFs) near the crack front in anisotropic EME bimaterials are given. Also, the numerical method of the HIE for a rectangular crack subjected to extended loads is put forward with the extended crack opening dislocation approximated by the product of basic density functions and polynomials. At last, numerical solutions of the extended SIFs of some examples are obtained.     

On waves due to rolling of a ship in water of finite depth

Summary The problem of the generation of waves due to small rolling oscillations of a thin vertical plate partially immersed in uniform finite-depth water is investigated here by utilizing two mathematical methods assuming the linearised theory of water waves. In the first method, the use of eigenfunction expansion of the velocity potentials on the two sides of the plate produces the amplitude of wave motion at infinity in terms of an integral involving the unknown horizontal velocity across the gap, and also in terms of another integral involving the unknown difference of the potential across the plate. These unknown functions satisfy two integral equations. Any one of these, when solved numerically, can be used to compute the amplitude of the wave motion set up at either infinity on the two sides of the plate for various values of the wave number.In the second method, the problem is formulated in terms of a hypersingular integral equation involving the difference of the potential function across the plate. The hypersingular integral equation is solved numerically, and its numerical solution is used to compute the wave amplitude at infinity. The two methods produce almost the same numerical results. The results are illustrated graphically, and a comparison is made with the deep-water result. It is observed that the deep-water result effectively holds good if the plate is partially immersed to the order of one-tenth of the bottom depth.This work was initiated when the first Author was visiting Mathematics Department, Indian Institute of Science, Bangalore. It was partially supported by DST, and by CSIR. The authors take this opportunity to thank the Managing Editor for his suggestions to improve the paper in the present form.     

Green's functions for anisotropic piezoelectric bimaterials with a slipping interface

An explicit full-field expression of the Green's functions for anisotropic piezoelectric bimaterials with a slipping interface is derived. When the electro-elastic singularity reduces to a pure dislocation in displacement and electric potential, interaction energy between the dislocation and the bimaterials is obtained explicitly while the generalized force on the dislocation is given in a real form which is also valid for degenerate materials. The investigation demonstrates that the boundary conditions at the slipping interface between two piezoelectric materials will exert a prominent influence on the mobility of the dislocation. Project supported by the National Natural Science Foundation of China (No. 59635140).     

一种新型三角形裂尖单元及其在结构裂纹分析中的应用

构造了一种适合边界元分析裂纹问题的三角形单元,该单元中的形函数包含两部分,主要部分用于捕捉裂纹尖端上位移分布的陡峭特性(性质),另一部分为常规的拟合函数,体现裂纹尖端位置附近的物理量在其他方向上的连续分布。形函数主要部分的构造充分利用了已有理论研究获得的结论,在裂纹表面,随着距离远离尖端,位移分布与■函数保持同阶变化。在传统形函数的基础上,通过先乘以一项同阶于■的变量项,再在系数中将其在形函数所在点上的值除去,便得到新型的用于拟合裂纹尖端附近位移和面力分布的形函数。新的形函数能够满足形函数的delta性质,但归一性不再满足,因此,新的形函数只用于物理量的拟合,而几何量的拟合依然采用传统方案。通过对偶边界元方法计算裂纹尖端的张开位移后,利用一种位移外插方法计算获得应力强度因子。数值算例关注了一种无限域内的圆盘裂纹,应用新构造的三角形单元于对偶边界元中计算结构在受到斜拉力时裂纹尖端的三种应力强度因子。通过与参考解进行对比,验证了该插值方案用于对偶边界元分析裂纹问题时的正确性和高精度。     

轴对称热冲击SIF过渡过程分析的热权函数法

热权函数法利用温度与热权函数的乘积的积分来直接计算热冲击过程中裂纹尖端的应力强度因子过渡过程。热权函数与时间τ无关,由于免除了对每一时刻τ所需作的有限元或边界元分析,计算过程大大简化,计算效率得到极大提高,本文将热权函数与有限元法直接耦合,给出了基于刚度阵导数法的轴对称问题的热权函数计算格式,实例计算表明,本文给出的热权函数计算格式具有满意的计算精度。     

Contact zone assessment for a fast growing interface crack in an anisotropic bimaterial

Summary A plane strain problem for a crack with a frictionless contact zone at the leading crack tip expanding stationary along the interface of two anisotropic half-spaces with a subsonic speed under the action of various loadings is considered. The cases of finite and infinite-length interface cracks under the action of a moving concentrated loading at its faces are considered. A finite-length crack for a uniform mixed-mode loading at infinity is considered as well. The associated combined Dirichlet-Riemann boundary value problems are formulated and solved exactly for all above-mentioned cases. The expressions for stresses and the derivatives of the displacement jumps at the interface are presented in a closed analytical form for an arbitrary contact zone length. Transcendental equations are obtained for the determination of the real contact zone length, and the associated closed form asymptotic formulas are found for small values of this parameter. It is found that independently of the types of the crack and loading, an increase of the crack tip speed leads to an increase of the real contact zone length and the correspondent stress intensity factor. The latter increase significantly for an interface crack tip speed approaching the Ragleigh wave speed.     

焊接节点半椭圆表面裂纹应力强度因子的权函数计算

本文采用一种改进的权函数法来计算焊接节点半椭圆表面裂纹应力强度因子ＫＩ值，并给出了相应的数值处理方法，就Ｔ型板节点进行了数值验算。     

Green’s functions for bending problem of a thin anisotropic plate with an elliptic hole

The Green’s functions have not been studied in open literatures for the bending problem of an anisotropic plate with an elliptic hole subjected to a normal concentrated force and a concentrated moment. In this paper, the problem is investigated and the Green’s functions are first obtained by using the complex potential approach. The techniques of conformal mapping transformation and analytic continuation are used to derive the closed-form complex stress functions. The Green’s functions obtained have some potential applications in the analysis of composite structures such as the modification of the displacement compatibility model for notched stiffened composite panels and the formulation of a new method for interlaminar stress analysis around holes of laminates.     

An internal penalty discontinuous Galerkin method for simulating conjugate heat transfer in a closed cavity

Using the discontinuous Galerkin (DG) method for conjugate heat transfer problems can provide improved accuracy close to the fluid‐solid interface, localizing the data exchange process, which may further enhance the convergence and stability of the entire computation. This paper presents a framework for the simulation of conjugate heat transfer problems using DG methods on unstructured grids. Based on an existing DG solver for the incompressible Navier‐Stokes equation, the fluid advection‐diffusion equation, Boussinesq term, and solid heat equation are introduced using an explicit DG formulation. A Dirichlet‐Neumann partitioning strategy has been implemented to achieve the data exchange process via the numerical flux of interface quadrature points in the fluid‐solid interface. Formal h and p convergence studies employing the method of manufactured solutions demonstrate that the expected order of accuracy is achieved. The algorithm is then further validated against 3 existing benchmark cases, including a thermally driven cavity, conjugate thermally driven cavity, and a thermally driven cavity with conducting solid, at Rayleigh numbers from 1000 to 100 000. The computational effort is documented in detail demonstrating clearly that, for all cases, the highest‐order accurate algorithm has several magnitudes lower error than first‐ or second‐order schemes for a given computational effort.     

Improved preconditioned conjugate gradient method and its application in F. E. A. for engineering

In this paper two theorems with theoretical and practical significance are given in respect to the preconditioned conjugate gradient method(PCCG).The theorems discuss respectively the qualitative property of the iterative solution and the construction principle of the iterative matrix.The authors put forward a new incompletely LU factorizing technique for non-M-matrix and the method of constructing the iterative matrix.This improved PCCG is used to calculate the ill-conditioned problems and large-scale three-dimensional finite element problems,and simultaneously contrasted with other methods.The abnormal phenomenon is analyzed when PCCG is used to solve the system of ill-conditioned equations,It is shown that the method proposed in this paper is quite effective in solving the system of large-scale finite element equations and the system of ill-conditioned equations.     

Symmetric and skew-symmetric weight functions in 2D perturbation models for semi-infinite interfacial cracks

In this paper we address the vector problem of a 2D half-plane interfacial crack loaded by a general asymmetric distribution of forces acting on its faces. It is shown that the general integral formula for the evaluation of stress intensity factors, as well as high-order terms, requires both symmetric and skew-symmetric weight function matrices. The symmetric weight function matrix is obtained via the solution of a Wiener–Hopf functional equation, whereas the derivation of the skew-symmetric weight function matrix requires the construction of the corresponding full field singular solution.The weight function matrices are then used in the perturbation analysis of a crack advancing quasi-statically along the interface between two dissimilar media. A general and rigorous asymptotic procedure is developed to compute the perturbations of stress intensity factors as well as high-order terms.     

Three-dimensional flow driven pore-crack networks in porous composites: Boltzmann Lattice method and hybrid hypersingular integrals

This study introduces a hybrid hypersingular integral equation-Lattice Boltzmann method (HHIE-LBM) for analyzing extended 3D flow driven pore-crack networks problem in various porosity composites. First, the extended hybrid electronic-ionic, thermal, magnetic, electric and force coupled fields’ pressure and velocity boundary conditions for HHIE-LBM model are established, and the closed form solutions of extended distribution functions are given. Second, an extended 3D flow driven pore-crack networks problem in various porosity composites is translated into a coupled of HHIE-LBM equations, and the pore-crack networks propagation parameters are analyzed. Third, the extended dynamic stress intensity factors (EDSIFs) are calculated by using the parallel numerical technology and the visualization results are presented. Last, the relationship between the EDSIFs and the differential porosity is discussed, and several rules have been found, which can be utilized to understand the extended fluid flow mechanism in various porosity composites and analyze the extended fluid flow varying mechanism on coseismal slip.     

Explicit conditions for the existence of exceptional body waves and subsonic surface waves in anisotropic elastic solids

It is known that a subsonic surface (Rayleigh) wave exists in an anisotropic elastic half-space x₂  0 if the first transonic state is not of Type 1. If the first transonic state is of Type 1 but the limiting wave is not exceptional, a subsonic surface wave exists. If the first transonic state is of Type 1 and the limiting wave is exceptional, a subsonic surface wave exists when . It is shown that an exceptional body wave is necessarily an exceptional transonic wave, and could be an exceptional limiting wave. Only two wave speeds are possible for an exceptional body wave. We present explicit conditions in terms of the reduced elastic compliances for the existence of an exceptional body wave. If an exceptional body wave exists, conditions are given for identifying whether the transonic state is of Type 1. Hence, through the existence of an exceptional body wave we provide explicit conditions for the existence of a subsonic surface wave with the exception when needs to be computed.     

Analysis of an interfacial crack in a piezoelectric bi-material via the extended Green’s functions and displacement discontinuity method

Based on the extended Stroh formalism, we first derive the extended Green’s functions for an extended dislocation and displacement discontinuity located at the interface of a piezoelectric bi-material. These include Green’s functions of the extended dislocation, displacement discontinuities within a finite interval and the concentrated displacement discontinuities, all on the interface. The Green’s functions are then applied to obtain the integro-differential equation governing the interfacial crack. To eliminate the oscillating singularities associated with the delta function in the Green’s functions, we represent the delta function in terms of the Gaussian distribution function. In so doing, the integro-differential equation is reduced to a standard integral equation for the interfacial crack problem in piezoelectric bi-material with the extended displacement discontinuities being the unknowns. A simple numerical approach is also proposed to solve the integral equation for the displacement discontinuities, along with the asymptotic expressions of the extended intensity factors and J-integral in terms of the discontinuities near the crack tip. In numerical examples, the effect of the Gaussian parameter on the numerical results is discussed, and the influence of different extended loadings on the interfacial crack behaviors is further investigated.     

Passages to the limit in the dynamic problem for a crack at the interface between elastic media

The paper analyzes numerically the passages to the limit in the dynamic problem for a penny-shaped crack at the interface between dissimilar linear elastic, homogeneous, isotropic materials as either the frequency of harmonic load or the difference between the properties of the materials decreases. It is shown that as the frequency decreases, the solution of the dynamic problem tends to that of the static problem, and as the physical and mechanical properties of the materials become less different, the original problem goes into the dynamic problem for a crack in a homogeneous body __________ Translated from Prikladnaya Mekhanika, Vol. 44, No. 7, pp. 26–34, July 2008.     

Lubrication theory for micropolar fluids and its application to a journal bearing with finite length

In this paper, the field equation of micropolar fluid with general lubrication theory assumptions is simplified into two systems of coupled ordinary differential equation. The analytical solutions of velocity and microrotation velocity are obtained. Micropolar fluid lubrication Reynolds equation is deduced. By means of numerical method, the characteristics of a finitely long journal bearing under various dynamic parameters, geometrical parameters and micropolar parameters are shown in curve form. These characteristics are pressure distribution, load capacity, coefficient of flow flux and coefficient of friction. Practical value of micropolar effects is shown, so micropolar fluid theory further closes to engineering application.