L.A. de Lacerda¹, J. M. da Silva¹

CMES-Computer Modeling in Engineering & Sciences, Vol.14, No.3, pp. 153-160, 2006, DOI:10.3970/cmes.2006.014.153

Abstract A two-dimensional boundary element formulation is presented and coupled to a genetic algorithm to identify polarization curves of buried slender structures. The dual boundary element method is implemented to model the cathodic protection of the metallic body and the genetic algorithm is employed to deal with the inverse problem of determining the non-linear polarization curve, which describes the relation between current density and electrochemical potential at the soil metal interface. In this work, this non-linear relation resulting from anodic and cathodic reactions is represented by a classical seven parameters expression. Stratified soil resistivity is modeled More >

Z. D. Han¹, H. T. Liu¹, A. M. Rajendran², S. N. Atluri³

CMES-Computer Modeling in Engineering & Sciences, Vol.14, No.2, pp. 119-128, 2006, DOI:10.3970/cmes.2006.014.119

Abstract This paper presents the implementation of a three-dimensional dynamic code, for contact, impact, and penetration mechanics, based on the Meshless Local Petrov-Galerkin (MLPG) approach. In the current implementation, both velocities and velocity-gradients are interpolated independently, and their compatibility is enforced only at nodal points. As a result, the time consuming differentiations of the shape functions at all integration points is avoided, and therefore, the numerical process becomes more stable and efficient. The ability of the MLPG code for solving high-speed contact, impact and penetration problems with large deformations and rotations is demonstrated through several computational More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.14, No.2, pp. 77-90, 2006, DOI:10.3970/cmes.2006.014.077

Abstract In this paper we derive the first-order and second-order one-step GPS applied to the estimation of thermophysical properties. Solving the resultant algebraic equations, which usually converges within ten iterations, it is not difficult to estimate the unknown temperature-dependent thermal conductivity and heat capacity simultaneously, if some supplemented data of measured temperature at a time T is provided. When the measured temperature in the conducting slab is contaminated by noise, our estimated results are also good. The new method does not require any prior information on the functional forms of thermal conductivity and heat capacity. Numerical examples More >

F. Han¹, E. Pan¹, A.K. Roy², Z.Q. Yue³

CMES-Computer Modeling in Engineering & Sciences, Vol.14, No.1, pp. 15-30, 2006, DOI:10.3970/cmes.2006.014.015

Abstract In this paper, an analytical solution is presented to study the response of piezoelectric, transversely isotropic, functionally graded, and multilayered half spaces to uniform circular surface loadings (pressure or negative electric charge). The inhomogeneous material is exponentially graded in the vertical direction and can have multiple discrete layers. The propagator matrix method and cylindrical system of vector functions are used to first derive the solution in the transformed domain. In order to find the responses in the physical-domain, which are expressed in one-dimensional infinite integrals of the Bessel function products, we introduced and adopted an… 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 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 >

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 More >

L. Gao¹, K. Liu^1,2, Y. Liu³

CMES-Computer Modeling in Engineering & Sciences, Vol.12, No.3, pp. 181-196, 2006, DOI:10.3970/cmes.2006.012.181

Abstract A new numerical algorithm based on the Meshless Local Petrov-Galerkin approach is presented for analyzing the dynamic fracture problems in elastic media. To simplify the treatment of essential boundary condition, a novel modified Moving Least Square (MLS) procedure is proposed by introducing Lagrange multiplier into MLS procedure, which can perform both MLS approximation and interpolation in one approximation domain. The compact spline function is used as the test function in the local form of elasto-dynamic equations. For the feature of stress wave propagation, the coupled second-order ODEs respect to the time are solved by the… More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.12, No.3, pp. 161-180, 2006, DOI:10.3970/cmes.2006.012.161

Abstract Materials' internal spacetime may bear certain similarities with the external spacetime of special relativity theory. Previously, it is shown that material hardening and anisotropy may cause the internal spacetime curved. In this paper we announce the third mechanism of mixed-control to cause the curvedness of internal spacetime. To tackle the mixed-control problem for a Prager kinematic hardening material, we demonstrate two new formulations. By using two-integrating factors idea we can derive two Lie type systems in the product space of M^m+1⊗Mⁿ⁺¹. The Lie algebra is a direct sum of so(m,1)⊕so(n,1), and correspondingly the symmetry group is… More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.12, No.2, pp. 83-108, 2006, DOI:10.3970/cmes.2006.012.083

Abstract The system we consider consists of two parts: a purely algebraic system describing the manifold of constraints and a differential part describing the dynamics on this manifold. For the constrained dynamical problem in its engineering application, it is utmost important to developing numerical methods that can preserve the constraints. We embed the nonlinear dynamical system with dimensions n and with k constraints into a mathematically equivalent n + k-dimensional nonlinear system, which including k integrating factors. Each subsystem of the k independent sets constitutes a Lie type system of X^˙_i = A_iX_i with A_i ∈ so(n_i,1) and n₁ +···+n_k = n.… More >