S. Wong, Y. Shie

CMES-Computer Modeling in Engineering & Sciences, Vol.36, No.2, pp. 97-118, 2008, DOI:10.3970/cmes.2008.036.097

Abstract In this paper, we propose a Galerkin-based smoothed particle hydrodynamics (SPH) formulation with moving least-squares meshless approximation, applied to solid mechanics and large deformation. Our method is truly meshless and based on Lagrangian kernel formulation and stabilized nodal integration. The performance of the methodology proposed is tested through various simulations, demonstrating the attractive ability of particle methods to handle severe distortions and complex phenomena. More >

H. Nguyen-Van¹, N.Mai-Duy¹, T. Tran-Cong¹

CMES-Computer Modeling in Engineering & Sciences, Vol.36, No.1, pp. 65-96, 2008, DOI:10.3970/cmes.2008.036.065

Abstract A novel node-based smoothing element for triangular and quadrilateral meshes is presented for static analysis of planar piezoelectric structures. In contrast to the smoothed finite element formulation that was based on sub-cells within an original quadrilateral element, this new method transforms a general original finite element mesh into a mesh of new smoothing cells individually associated with a single node which is termed as node-based elements. The displacement fields of the element are approximated by the linear interpolation functions of the original mesh while the approximations of mechanical strains and electric potential fields are normalized using the stabilized conforming nodal… More >

Giorgio Pini¹, Annamaria Mazzia¹, Flavio Sartoretto²,

CMES-Computer Modeling in Engineering & Sciences, Vol.36, No.1, pp. 43-64, 2008, DOI:10.3970/cmes.2008.036.043

Abstract Meshless methods have been explored in many 2D problems and they have been shown to be as accurate as Finite Element Methods (FEM). Compared to the extensive literature on 2D applications, papers on solving 3D problems by meshless methods are surprisingly few. Indeed, a main drawback of these methods is the requirement for accurate cubature rules. This paper focuses on the so called Meshless Local Petrov Galerkin (MLPG) methods. We show that accurate solutions of 3D potential problems can be attained, provided suitable cubature rules are identified, sparse data structures are efficiently stored, and strategies are devised in order to… More >

A. Alaimo¹, A. Milazzo², C. Orlando³

CMES-Computer Modeling in Engineering & Sciences, Vol.36, No.1, pp. 23-42, 2008, DOI:10.3970/cmes.2008.036.023

Abstract In this paper a Hierarchical approach for the analysis of advanced aerospace structures is presented. The proposed Global/Local model uses two kind of numerical methods. The first step of the Hierarchical procedure is performed by the Finite Element code Patran/Nastran\texttrademark , using a coarse mesh to study the global structure, then the local region is analyzed by using a Boundary Element code based on the multidomain anisotropic technique. This code accurately predicts stress concentrations at crack tips with a reduction of the modeling efforts and of the computational time. The Global/Local interface code implemented allows an intuitive extraction of the… More >

Cheng-Hung Huang¹, Chi-An Chen¹

CMES-Computer Modeling in Engineering & Sciences, Vol.36, No.1, pp. 1-22, 2008, DOI:10.3970/cmes.2008.036.001

Abstract A three-dimensional shape identification problem (or inverse geometry problem) in estimating the unknown irregular shape of internal cavity by using the steepest descent method (SDM) and a general purpose commercial code CFD-RC is examined in this study based on the simulated measured temperature distributions on the outer surface by infrared thermography. The advantage of calling CFD-RC as a subroutine in the present shape identification problem lies in its characteristics of easily-handling the problem considered here since the auto mesh function of CFD-RC enables the handling of this moving boundary problem. \newline Three test cases are performed to test the validity… More >

H.X. Wang, S.X. Wang

CMES-Computer Modeling in Engineering & Sciences, Vol.35, No.3, pp. 253-274, 2008, DOI:10.3970/cmes.2008.035.253

Abstract In the meshfree cohesive crack method, the discrete crack is modeled by a set of cohesive crack segments which can be arbitrarily oriented. Propagation of the crack is achieved by activation of crack surfaces at individual nodes, so no representation of the crack surface is needed. The crack is modeled by a local enrichment of the test and trial functions with sign function, so that discontinuities are along the direction of the crack. A set of cracking rules is developed to avoid spurious cracking. More >

De-hao Yu¹ ,Hong-ying Huang²

CMES-Computer Modeling in Engineering & Sciences, Vol.35, No.3, pp. 227-252, 2008, DOI:10.3970/cmes.2008.035.227

Abstract In this paper, we apply the artificial boundary method to solve a three-dimensional nonlinear interface problem on an unbounded domain. A spherical or ellipsoidal surface as the artificial boundary is introduced. The exact artificial boundary conditions are derived explicitly in terms of an infinite series and then the well-posedness of the coupled weak formulation in a bounded domain, which is equivalent to the original problem in the unbounded domain, is obtained. The error estimate depends on the mesh size, the term after truncating the infinite series and the location of the artificial boundary. Some numerical examples are presented to demonstrate… More >

Yongfeng Zhang¹, Hanchen Huang^1,2, Satya N. Atluri³

CMES-Computer Modeling in Engineering & Sciences, Vol.35, No.3, pp. 215-226, 2008, DOI:10.3970/cmes.2008.035.215

Abstract Molecular dynamics simulations reveal the asymmetrical yield strength of twinned copper nanowires under tension and compression. The simulation results show that the strength of nanowires depends on loading conditions, morphologies, and twin spacing. Under tensile loading condition the Schmidt factor of the leading partial is larger than that under compression. Effectively, the yield strength under tension is smaller than that under compression. When the cross-section is circular in morphology, dislocation nucleation requires larger stress, and the asymmetry of yield strength depends on the nucleation stress. When the cross section is square in morphology, dislocation nucleation requires smaller stress, and the… More >

Chih-Ping Wu^1,2, Kuan-Hao Chiu², Yung-Ming Wang²

CMES-Computer Modeling in Engineering & Sciences, Vol.35, No.3, pp. 181-214, 2008, DOI:10.3970/cmes.2008.035.181

Abstract A mesh-free collocation method based on differential reproducing kernel (DRK) approximations is developed for the three-dimensional (3D) analysis of simply-supported, doubly curved functionally graded (FG) magneto-electro-elastic shells under the mechanical load, electric displacement and magnetic flux. The material properties of FG shells are firstly regarded as heterogeneous through the thickness coordinate and then specified to obey an identical power-law distribution of the volume fractions of the constituents. The novelty of the present DRK-based collocation method is that the shape functions of derivatives of reproducing kernel (RK) approximants are determined using a set of differential reproducing conditions without directly taking the… More >

Chein-Shan Liu ¹

CMES-Computer Modeling in Engineering & Sciences, Vol.35, No.2, pp. 157-180, 2008, DOI:10.3970/cmes.2008.035.157

Abstract For an inverse vibration problem of nonlinear mechanical system to estimate displacement- and velocity-dependent restoring force, we transform the equation of motion into a parabolic type partial differential equation (PDE). Then by a semi-discretization of the PDE, the inverse vibration problem is formulated as a multi-dimensional two-point boundary value problem with unknown sources, allowing a closed-form estimation through a Lie-group shooting method to construct the restoring force surface over phase plane. Only one set of displacements measured at sampling time points is used in the estimation. The new method does not require to assume a priori the functional form of… More >