An integral equation approach to the inclusion-crack interactions in three-dimensional infinite elastic domain

 In this paper, an integral equation method to the inclusion-crack interaction problem in three-dimensional elastic medium is presented. The method is implemented following the idea that displacement integral equation is used at the source points situated in the inclusions, whereas stress integral equation is applied to source points along crack surfaces. The displacement and stress integral equations only contain unknowns in displacement (in inclusions) and displacement discontinuity (along cracks). The hypersingular integrals appearing in stress integral equation are analytically transferred to line integrals (for plane cracks) which are at most weakly singular. Finite elements are adopted to discretize the inclusions into isoparametric quadratic 10-node tetrahedral or 20-node hexahedral elements and the crack surfaces are decomposed into discontinuous quadratic quadrilateral elements. Special crack tip elements are used to simulate the variation of displacements near the crack front. The stress intensity factors along the crack front are calculated. Numerical results are compared with other available methods. Received: 28 January 2002 / Accepted: 4 June 2002 The work described in this paper was partially supported by a grant from the Research Grant Council of the Hong Kong Special Administration Region, China (Project No.: HKU 7101/99 E).     

Interaction between cracks and rigid-line inclusions by an integral equation approach

 An integral equation approach is presented to investigate the interaction between cracks and rigid-line inclusions embedded in an infinite isotropic elastic matrix subject to remote loading. The relevant fundamental solutions in the integral formulation are presented. Special tip elements are used to simulate the variation of the discontinuous displacements over the crack surfaces, and the axial and shear forces along the rigid-line inclusions. The stress intensity factors at the crack tips and at the ends of the rigid-line inclusions are computed and compared with available solutions. Received: 6 August 2002 / Accepted: 3 February 2003 The work described in this paper was partially supported by a grant from the Research Grant Council of the Hong Kong Special Administration Region, China (Project No.: HKU 7011/01E). The authors would like to thank two reviewers for their constructive comments and suggestions to the paper. The comments of Professor H.P. Hong at the Department of Civil and Environmental Engineering of the University of Western Ontario of Canada are also appreciated.     

Boundary integral equation method for two dissimilar elastic inclusions in an infinite plate

This paper provides a numerical solution for an infinite plate containing two dissimilar elastic inclusions, which is based on complex variable boundary integral equation (CVBIE). The original problem is decomposed into two problems. One is an interior boundary value problem (BVP) for two elastic inclusions, while other is an exterior BVP for the matrix with notches. After performing discretization for the coupled boundary integral equations (BIEs), a system of algebraic equations is formulated. The inverse matrix technique is suggested to solve the relevant algebraic equations, which can avoid using the assembling of some matrices. Several numerical examples are carried out to prove the efficiency of suggested method and the hoop stress along the interface boundary is evaluated.     

Numerical analysis of doubly periodic array of cracks/rigid-line inclusions in an infinite isotropic medium using the boundary integral equation method

In this paper, the boundary integral equation approaches are used to study the doubly periodic array of cracks/rigid-line inclusions in an infinite isotropic plane medium. For the doubly periodic rigid-line inclusion problems, the special integral equation containing the axial and shear forces within the rigid-line inclusion is used. The doubly periodic crack problems are dealt with using the displacement discontinuous integral equation approach. Stress intensity factors, effective elastic properties for doubly periodic array of cracks/rigid-line inclusions are calculated and compared with the available numerical solutions.     

Numerical analysis of the inclusion-crack interactions using an integral equation

 A general-purpose integral formulation is proposed for the analysis of the interaction between inclusions and cracks embedded in an elastic isotropic homogeneous infinite medium subjected to a remote loading. This formulation is tailored for the inclusions of arbitrary shapes with the presence of cracks. The discretization is limited to the inclusions (with continuous quadratic triangular and quadrilateral elements) and the cracks (using discontinuous quadratic elements). For the calculation of the stress intensity factors at the crack tips, special crack tip elements are used to model the variation of the displacements near the crack tips. Maximum circumferential stress criterion is adopted to determine the crack propagating direction. Numerical results of benchmark examples are compared with other available methods. Received: 8 January 2002 / Accepted: 24 September 2002     

Diffraction of elastic waves by four rigid strips embedded in an infinite orthotropic medium

In the present paper, the elastodynamic response of four coplanar rigid strips embedded in an infinite orthotropic medium due to elastic waves incident normally on the strips is analyzed. The resulting mixed boundary-value problem is solved by an integral-equation method. The normal stress and the vertical displacement are derived in closed analytic form. Numerical values of stress-intensity factors at the edges of the strips and vertical displacements at point in the plane of the strips for several orthotropic materials are calculated and plotted graphically to show the effect of material orthotropy.     

Coupling of the BEM with the MFS for the numerical simulation of frequency domain 2-D elastic wave propagation in the presence of elastic inclusions and cracks

This paper proposes a coupling formulation between the boundary element method (BEM displacement and TBEM traction formulations) and the method of fundamental solutions (MFS) for the transient analysis of elastic wave propagation in the presence of multiple elastic inclusions to overcome the specific limitations of each of these methods. The full domain of the original problem is divided into sub-domains, which are handled separately by the BEM or the MFS. The coupling is enforced by imposing the required boundary conditions.The accuracy, efficiency and stability of the proposed algorithms, using different combinations of BEM and MFS, are verified by comparing the solutions against reference solutions. The computational efficiency of the proposed coupling formulation is illustrated by computing the CPU time and the error at high frequencies.The potential of the proposed procedures is illustrated by simulating the propagation of elastic waves in the vicinity of an empty crack, with null thickness placed close to an elastic inclusion.     

The elastic problem for an infinite solid containing a circular hole with a pair of radial edge cracks of different lengths

By using the principle of superposition and the results of an earlier paper the problem is reduced to that of solving a pair of singular integral equations. A number of special cases are considered and numerical results given.     

Doubly periodic arrays of cracks and thin inhomogeneities in an infinite magnetoelectroelastic medium

A two-dimensional boundary element method for the analysis of a magnetoelectroelastic medium containing doubly periodic sets of cracks or thin inclusions is developed in this paper. The integral equations and closed-form expressions for corresponding kernels are obtained. Based on the quasi-periodicity of extended displacement and stress function, the integral representations for average stress, strain, electric displacement, magnetic induction etc. are developed. The algorithm of effective properties determination is given. The numerical examples prove the efficiency and high accuracy of the proposed approach in determination of stress, electric displacement and magnetic induction intensity factors and effective properties of the material containing doubly periodic arrays of cracks or thin inclusions.     

The DRM‐MD integral equation method: an efficient approach for the numerical solution of domain dominant problems

This work presents a multi‐domain decomposition integral equation method for the numerical solution of domain dominant problems, for which it is known that the standard Boundary Element Method (BEM) is in disadvantage in comparison with classical domain schemes, such as Finite Difference (FDM) and Finite Element (FEM) methods. As in the recently developed Green Element Method (GEM), in the present approach the original domain is divided into several subdomains. In each of them the corresponding Green's integral representational formula is applied, and on the interfaces of the adjacent subregions the full matching conditions are imposed. In contrast with the GEM, where in each subregion the domain integrals are computed by the use of cell integration, here those integrals are transformed into surface integrals at the contour of each subregion via the Dual Reciprocity Method (DRM), using some of the most efficient radial basis functions known in the literature on mathematical interpolation. In the numerical examples presented in the paper, the contour elements are defined in terms of isoparametric linear elements, for which the analytical integrations of the kernels of the integral representation formula are known. As in the FEM and GEM the obtained global matrix system possesses a banded structure. However in contrast with these two methods (GEM and non‐Hermitian FEM), here one is able to solve the system for the complete internal nodal variables, i.e. the field variables and their derivatives, without any additional interpolation. Finally, some examples showing the accuracy, the efficiency, and the flexibility of the method for the solution of the linear and non‐linear convection–diffusion equation are presented. Copyright © 1999 John Wiley & Sons, Ltd.     

A monolithic approach for interaction of incompressible viscous fluid and an elastic body based on fluid pressure Poisson equation

This paper describes a new monolithic approach based on the fluid pressure Poisson equation (PPE) to solve an interaction problem of incompressible viscous fluid and an elastic body. The PPE is derived so as to be consistent with the coupled equation system for the fluid‐structure interaction (FSI). Based on this approach, we develop two kinds of efficient monolithic methods. In both methods, the fluid pressure is derived implicitly so as to satisfy the incompressibility constraint, and all other unknown variables are derived fully explicitly or partially explicitly. The coefficient matrix of the PPE for the FSI becomes symmetric and positive definite and its condition is insensitive to inhomogeneity of material properties. The arbitrary Lagrangian–Eulerian (ALE) method is employed for the fluid part in order to take into account the deformable fluid‐structure interface. To demonstrate fundamental performances of the proposed approach, the developed two monolithic methods are applied to evaluate the added mass and the added damping of a circular cylinder as well as to simulate the vibration of a rectangular cylinder induced by vortex shedding. Copyright © 2005 John Wiley & Sons, Ltd.     

Constitutive equation for a class of isotropic, perfectly elastic solids using a new measure of finite strain and corresponding stress

The definition of a measure of strain, referred to as the bi-configuration strain tensor, centres on the difference between the left Cauchy-Green deformation tensor and its inverse. A new measure of stress, coined the bi-configuration stress tensor, has been defined. This measure of stress refers the traction in the current configuration jointly to the referential and spatial configurations, that is, to an effective element of area identified as an element of bi-configuration area. The stress and strain tensors are assumed to be constitutively related by a finite strain form of a generalised Hooke's law. The predictions obtained from the proposed constitutive equation are compared with the observed mechanical behaviour of various test materials. Comparison with experiment centres on biaxial stress measurements in various simple modes of deformation identified by way of a generalised stress-strain relation. The predictions from the proposed constitutive theory are in good accord with the results of experiment.     

The impact of human consideration,schedule types and product mix on scheduling objectives for unpaced mixed-model assembly lines

There is an increasing awareness in scheduling research that human behaviour needs to be considered explicitly in scheduling models. Although most scheduling literature ignores human behaviour, especially sequence-dependent processing times form a good basis for explicit consideration. Hence, a processing time function is derived that considers the effects of learning, forgetting, fatigue and recovery. The necessity for explicit human consideration can be regarded as most urgent for unpaced highly-manual mixed-model assembly lines. Based on real data a simulation study is conducted to determine the effect of explicit human consideration while also taking into account the effects of different idealised schedule types and the product mix. The results strongly indicate that the product mix has a consistently high impact on scheduling objectives, the schedule type affects lower-level objectives like starving and blocking times to a greater extent than higher-level objectives like makespan and flow time, and that for certain objectives the height of the objective values and the relative favourability of schedule types depends on human consideration.