S. Beyer¹, Ch. Zhang², S. Hirose³, J. Sladek, V. Sladek⁴

Structural Durability & Health Monitoring, Vol.3, No.3, pp. 177-190, 2007, DOI:10.3970/sdhm.2007.003.177

Abstract This paper presents a hypersingular time-domain boundary element method (BEM) for transient dynamic crack analysis in two-dimensional (2-D), homogeneous, anisotropic and linear elastic solids. A finite crack in an infinite or a finite solid subjected to impact loading conditions is investigated. A combination of the classical displacement boundary integral equations (BIEs) on the external boundary and the hypersingular traction BIEs on the crack-faces is applied. The present BEM uses the time-domain dynamic fundamental solutions for anisotropic solids derived by Wang and Achenbach (1994). An explicit time-stepping scheme based on collocation method is developed. Numerical examples More >

P.D. Shah¹, C.L. Tan¹, X. Wang¹

The International Conference on Computational & Experimental Engineering and Sciences, Vol.1, No.2, pp. 75-80, 2007, DOI:10.3970/icces.2007.001.075

Abstract The path-independent mutual- or M-integral for the computation of the T-stress for interface cracks lying between dissimilar anisotropic, linear elastic solids is developed in this paper. For the numerical stress analysis, the Boundary Element Method (BEM) is employed and special crack-tip elements with the proper oscillatory traction singularity are used. The successful implementation of the scheme for evaluating the T-stress of an interface crack between anisotropic bi-materials with BEM is demonstrated by numerical examples. More >

Dalin Tang^*, Chun Yang^†, Tal Geva^‡,§, Pedro J. del Nido^¶

Molecular & Cellular Biomechanics, Vol.4, No.3, pp. 159-176, 2007, DOI:10.3970/mcb.2007.004.159

Abstract A single-layer isotropic patient-specific right/left ventricle and patch (RV/LV/Patch) combination model with fluid-structure interactions (FSI) was introduced in our previous papers to evaluate and optimize human pulmonary valve replacement/insertion (PVR) surgical procedure and patch design. In this paper, an active anisotropic model with two-layer structure for ventricle wall and tissue fiber orientation was introduced to improve previous isotropic model for more accurate assessment of RV function and potential application in PVR surgery and patch design. A material-stiffening approach was used to model active heart contraction. The computational models were used to conduct ``virtual (computational)'' surgeries More >

J. S. Ren^*

Molecular & Cellular Biomechanics, Vol.4, No.1, pp. 55-66, 2007, DOI:10.3970/mcb.2007.004.055

Abstract Formation and rupture of aneurysms due to the inflation of an artery with collagen fibers distributed in two preferred directions, subjected to internal pressure and axial stretch are examined within the framework of nonlinear elasticity. A two layer tube model with a fiber-reinforced composite based incompressible anisotropic hyperelastic constitutive material is employed to model the stress-strain behavior of the artery wall with distributed collagen fibers. The artery wall takes up a uniform inflation deformation, and there are no aneurysms in the artery under the normal condition. But an aneurysm may be formed in arteries when… More >

J.N. Johnson¹, J.M. Owen²

CMES-Computer Modeling in Engineering & Sciences, Vol.22, No.3, pp. 165-188, 2007, DOI:10.3970/cmes.2007.022.165

Abstract In this paper, we propose a Meshless Local Petrov-Galerkin method for studying the diffusion of a magnetic field within a non-magnetic (μ = μ₀) conducting medium with non-homogeneous and anisotropic electrical resistivity. We derive a local weak form for the magnetic diffusion equation and discuss the effects of different trial/test functions and nodal spacings on its solution. We then demonstrate that the method produces convergent results for several relevant one-dimensional test problems for which solutions are known. This method has the potential to be combined with other mesh-free methods such as Smoothed Particle Hydrodynamics (SPH) to More >

Ping’en Li¹, Xianyue Su^1,2, Youquan Yin¹

CMC-Computers, Materials & Continua, Vol.6, No.2, pp. 93-102, 2007, DOI:10.3970/cmc.2007.006.093

Abstract Based on the perturbation method, for flexural wave in cased hole in anisotropic formation, the alteration in the phase velocity caused by the differences in elastic constants between anisotropic formation of interest and a reference, or unperturbed isotropic formation is obtained. Assuming the cased hole is well bonded, the Thomson-Haskell transfer matrix method is applied to calculate the dispersion relation of flexural wave in cased hole in unperturbed isotropic formation. Both the cases of a fast and slow formation are considered where the symmetry axis of a transversely isotropic (TI) formation makes an angle with… More >

P.D. Shah¹, Ch. Song², C.L. Tan¹, X. Wang¹

Structural Durability & Health Monitoring, Vol.2, No.4, pp. 225-238, 2006, DOI:10.3970/sdhm.2006.002.225

Abstract Numerical T-stress solutions in two dimensional anisotropic cracked bodies are very scarce in the literature. Schemes to evaluate this fracture parameter in anisotropy have been reported only fairly recently. Among them are those developed in conjunction with two different computational techniques, namely, the Boundary Element Method (BEM) and the Scaled Boundary Finite-Element Method (SBFEM). This paper provides a review of the respective schemes using these techniques and demonstrates their efficacy with three examples. These examples, which are of engineering importance, involve cracks lying in a homogeneous medium as well as at the interface between dissimilar media. More >

Y.C. Shiah¹, Y.C. Lin¹, C. L. Tan²

CMES-Computer Modeling in Engineering & Sciences, Vol.16, No.1, pp. 15-26, 2006, DOI:10.3970/cmes.2006.016.015

Abstract In this paper, the order of singularity of the integrals appearing in the boundary integral equation for two-dimensional BEM analysis in anisotropic elasticity is reduced using integration by parts. The integral containing the traction fundamental solution is then analytically integrated to give an exact formulation for a general element of n-order interpolation of the variables. This allows the integrals to be very accurately evaluated even for very thin, slender bodies without the need for excessively refined meshes as in conventional BEM analysis. Three example problems involving thin, layered materials are presented to demonstrate the veracity and More >

B. Yang², V. K. Tewary³

CMES-Computer Modeling in Engineering & Sciences, Vol.15, No.3, pp. 165-178, 2006, DOI:10.3970/cmes.2006.015.165

Abstract Green's function (GF) modeling of defects may take effect only if the GF as well as its various integrals over a line, a surface and/or a volume can be efficiently evaluated. The GF is needed in modeling a point defect, while integrals are needed in modeling line, surface and volumetric defects. In a matrix of multilayered, generally anisotropic and linearly elastic and piezoelectric materials, the GF has been derived by applying 2D Fourier transforms and the Stroh formalism. Its use involves another two dimensions of integration in the Fourier inverse transform. A semi-analytical scheme has… More >

P.D. Shah¹, C.L. Tan^1,2, X. Wang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.13, No.3, pp. 185-198, 2006, DOI:10.3970/cmes.2006.013.185

Abstract In this paper, the path independent mutual or M-integral for the computation of the T-stress for interface cracks between dissimilar anisotropic, linear elastic solids, is developed. The required auxiliary field solution is derived from the solution of the problem of an anisotropic composite wedge subjected to a point force at its apex. The Boundary Element Method (BEM) is employed for the numerical stress analysis in which special crack-tip elements with the proper oscillatory traction singularity are used. The successful implementation of the procedure for evaluating the T-stress in a bi-material interface crack and its application are demonstrated More >