Xianwu Ling¹, S.N. Atluri¹

CMES-Computer Modeling in Engineering & Sciences, Vol.13, No.3, pp. 219-228, 2006, DOI:10.3970/cmes.2006.013.219

Abstract In this paper, two matrix algebraic tools are provided for studying the solution-stabilities of inverse heat conduction problems. The propagations of the computed temperature errors, as caused by a noise in temperature measurement, are presented. The spectral norm analysis reflects the effect of the computational time steps, the sensor locations and the number of future temperatures on the computed error levels. The Frobenius norm analysis manifests the dynamic propagations of the computed errors. As an application of the norm analysis, we propose a method for the best positioning of the thermocouples. More >

H. K. Ching^1,2, J. K. Chen²

CMES-Computer Modeling in Engineering & Sciences, Vol.13, No.3, pp. 199-218, 2006, DOI:10.3970/cmes.2006.013.199

Abstract The Meshless Local Petrov-Galerkin (MLPG) method is a novel numerical approach similar to finite element methods, but it allows the construction of the shape function and domain discretization without defining elements. In this study, the MLPG analysis for transient thermomechanical response of a functionally graded composite heated by Gaussian laser beams is presented. The composite is modeled as a 2-D strip which consists of metal and ceramic phases with the volume fraction varying over the thickness. Two sets of the micromechanical models are employed for evaluating the effective material properties, respectively. Numerical results are presented for the thermomechanical responses in… 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 by numerical examples. More >

V. Vavourakis¹, E. J. Sellountos², D. Polyzos³

CMES-Computer Modeling in Engineering & Sciences, Vol.13, No.3, pp. 171-184, 2006, DOI:10.3970/cmes.2006.013.171

Abstract Comparison studies on the accuracy provided by five different elastostatic Meshless Local Petrov-Galerkin (MLPG) type formulations, based on Local Boundary Integral Equation (LBIE) considerations, are made. The main differences of these MLPG(LBIE) formulations, as they compared to each other, are concentrated on the treatment of tractions on the local and global boundaries and the way of imposing the boundary conditions of the elastostatic problem. Both the Moving Least Square (MLS) approximation scheme and the Radial Basis Point Interpolation Functions (RBPIF) are exploited for the interpolation of the interior and boundary variables. Two representative elastostatic problems are solved and the relative… More >

A. N. Guz, V. A. Dekret¹

CMES-Computer Modeling in Engineering & Sciences, Vol.13, No.3, pp. 165-170, 2006, DOI:10.3970/cmes.2006.013.165

Abstract Stability problem of composite material reinforced by two parallel short fibers is solved. The problem is formulated with application of equations of linearized three-dimensional theory of stability. The composite is modeled as piecewise-homogeneous medium. The influence of geometrical and mechanical parameters of the material on critical strain is investigated. More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.13, No.2, pp. 149-164, 2006, DOI:10.3970/cmes.2006.013.149

Abstract The present paper provides a Lie-group shooting method for the numerical solutions of second order nonlinear boundary value problems exhibiting multiple solutions. It aims to find all solutions as easy as possible. The boundary conditions considered are classified into four types, namely the Dirichlet, the first Robin, the second Robin and the Neumann. The two Robin type problems are transformed into a canonical one by using the technique of symmetric extension of the governing equations. The Lie-group shooting method is very effective to search unknown initial condition through a weighting factor r ∈ (0,1) Furthermore, the closed-form solutions are derived… More >

S.Y. Wang^1,2, M.Y. Wang³

CMES-Computer Modeling in Engineering & Sciences, Vol.13, No.2, pp. 119-148, 2006, DOI:10.3970/cmes.2006.013.119

Abstract In this paper, an implicit free boundary parametrization method is presented as an effective approach for simultaneous shape and topology optimization of structures. The moving free boundary of a structure is embedded as a zero level set of a higher dimensional implicit level set function. The radial basis functions (RBFs) are introduced to parametrize the implicit function with a high level of accuracy and smoothness. The motion of the free boundary is thus governed by a mathematically more convenient ordinary differential equation (ODE). Eigenvalue stability can be guaranteed due to the use of inverse multiquadric RBF splines. To perform both… More >

J. Sladek¹, V. Sladek¹, P. H. Wen², M.H. Aliabadi³

CMES-Computer Modeling in Engineering & Sciences, Vol.13, No.2, pp. 103-118, 2006, DOI:10.3970/cmes.2006.013.103

Abstract A meshless local Petrov-Galerkin (MLPG) method is applied to solve bending problems of shear deformable shallow shells described by the Reissner theory. Both static and dynamic loads are considered. For transient elastodynamic case the Laplace-transform is used to eliminate the time dependence of the field variables. A weak formulation with a unit test function transforms the set of governing equations into local integral equations on local subdomains in the mean surface of the shell. Nodal points are randomly spread on that surface and each node is surrounded by a circular subdomain to which local integral equations are applied. The meshless… More >

Radek Pecher¹, Steve Elston, Peter Raynes

CMES-Computer Modeling in Engineering & Sciences, Vol.13, No.2, pp. 91-102, 2006, DOI:10.3970/cmes.2006.013.091

Abstract Meshfree techniques for solving partial differential equations in physics and engineering are a powerful new alternative to the traditional mesh-based techniques, such as the finite difference method or the finite element method. The elimination of the domain mesh enables, among other benefits, more efficient solutions of nonlinear and multi-scale problems. One particular example of these kinds of problems is a Q-tensor based model of nematic liquid crystals involving topological defects.
This paper presents the first application of the meshless local Petrov-Galerkin method to solving the Q-tensor equations of nematostatics. The theoretical part introduces the Landau -- de Gennes free-energy… More >

Zai You Yan¹

CMES-Computer Modeling in Engineering & Sciences, Vol.13, No.2, pp. 81-90, 2006, DOI:10.3970/cmes.2006.013.081

Abstract Boundary element method in acoustics for Neumann boundary condition problems including sharp edges & corners is investigated. In previous acoustic boundary element method, acoustic pressure and normal velocity are the two variables at sharp edges & corners. However, the normal velocity at sharp edges & corners is discontinuous due to the indefinite normal vector. To avoid the indefinite normal vector and the discontinuous normal velocity at sharp edges & corners, normal vector of elemental node is defined and applied in the numerical implementation. Then the normal velocity is transformed to velocity which is unique even at sharp edges & corners.… More >