F. Cosmi¹

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.3, pp. 365-372, 2001, DOI:10.3970/cmes.2001.002.365

Abstract The aim of this paper is to present a methodology for solving the plane elasticity problem using the Cell Method. It is shown that with the use of a parabolic interpolation in a vectorial problem, a convergence rate of 3.5 is obtained. Such a convergence rate compares with, or is even better than, the one obtained with FEM with the same interpolation – depending on the integration technique used by the FEM application. The accuracy of the solution is also comparable or better. More >

KwangsooHo¹

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.3, pp. 351-364, 2001, DOI:10.3970/cmes.2001.002.351

Abstract Negative, zero or positive rate sensitivity of the flow stress can be observed in metals and alloys over a certain range of strain, strain rate and temperature. It is believed that negative rate sensitivity is an essential feature of dynamic strain aging, of which the Portevin-Le Chatelier effect is one other manifestation. The viscoplasticity theory based on overstress (VBO), one of the unified state variable theories, is generalized to model zero (rate independence) and negative as well as positive rate sensitivity in a consistent way. The present model does not have the stress rate term… More >

Xiaolin Chen, Yijun Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.3, pp. 337-350, 2001, DOI:10.3970/cmes.2001.002.337

Abstract An advanced boundary element method (BEM) is developed in this paper for analyzing thin layered structures, such as thin films and coatings, under the thermal loading. The boundary integral equation (BIE) formulation for steady-state thermoelasticity is reviewed and a special case, that is, the BIE for a uniform distribution of the temperature change, is presented. The new nearly-singular integrals arising from the applications of the BIE/BEM to thin layered structures under thermal loading are treated in the same way as developed earlier for thin structures under the mechanical loading. Three 2-D test problems involving layered More >

R. Vignjevic, T. De Vuyst, M. Gourma¹

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.3, pp. 319-336, 2001, DOI:10.3970/cmes.2001.002.319

Abstract The work presented provides an overview of different types of kernel interpolation used in the SPH method: conventional SPH, normalised SPH (NSPH), corrected kernel SPH (CSPH) and normalised corrected kernel SPH (NCSPH). These four methods are considered in a fully mesh-free form (using no background mesh). To illustrate the effect of using different interpolation methods two problems were simulated: a 1D symmetric elastic impact problem, and a shock-tube. An overview of the simulation results for the two problems is given. Shortcomings for the interpolation schemes tested were identified and discussed. It is concluded that NCSPH More >

Jianming Zhang, Zhenhan Yao¹

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.3, pp. 307-318, 2001, DOI:10.3970/cmes.2001.002.307

Abstract The Meshless Local Boundary Integral Equation (MLBIE) method is a truly meshless one as it does not need a 'finite element or boundary element mesh', either for variable interpolation or for 'energy' integration. The Boundary Node Method (BNM) further reduces the dimensionality of the problem by one, i.e. it only requires nodes constructed on the surface. However, the BNM is not truly meshless, as a background mesh is needed for boundary integration; and the MLBIE does not have the advantage of reduced dimensionality as the BNM. A new Regular Hybrid Boundary Node method based on… More >

Yu.A. Melnikov, M.Yu. Melnikov¹

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.2, pp. 291-306, 2001, DOI:10.3970/cmes.2001.002.291

Abstract The use of potential (integral) representations is studied when computing Green's functions for boundary value problems stated for Laplace and biharmonic equations over regions of complex configuration in two dimensions. The emphasis is on the non-traditional potentials, whose observation and source points occupy different sets. Such potentials reduce the original boundary value problems to functional (integral) equations with smooth kernels. Special integral representations are studied, the ones whose kernels are built not of the fundamental solutions of governing differential equations but of the Green's functions for simply shaped regions, which are associated with boundary value More >

H.-K. Ching, R. C. Batra¹

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.2, pp. 273-290, 2001, DOI:10.3970/cmes.2001.002.273

Abstract It is shown that the MLPG method augmented with the enriched basis functions and either the visibility criterion or the diffraction criterion successfully predicts the singular stress fields near a crack tip. Results are presented for a single edge-cracked plate and a double edge-cracked plate with plate edges parallel to the crack axis loaded in tension, the single edge-cracked plate with one plate edge parallel to the crack axis clamped and the other loaded by tangential tractions, and for a double edge-notched plate with the side between the notches loaded in compression. For the first More >

R. Sedaghati, A. Suleman¹, S. Dost, B. Tabarrok²

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.2, pp. 259-272, 2001, DOI:10.3970/cmes.2001.002.259

Abstract A structural analysis and optimization method is developed to find the optimal topology of adaptive determinate truss structures under various impact loading conditions. The objective function is based on the maximization of the structural strength subject to geometric constraints. The dynamic structural analysis is based on the integrated finite element force method and the optimization procedure is based on the Sequential Quadratic Programming (SQP) method. The equilibrium matrix is generated automatically through the finite element analysis and the compatibility matrix is obtained directly using the displacement-deformation relations and the Single Value Decomposition (SVD) technique. By More >

Enzo TONTI¹

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.2, pp. 237-258, 2001, DOI:10.3970/cmes.2001.002.237

Abstract We present a new numerical method for the solution of field equations. The essence of the method is to directly provide a discrete formulation of field laws, without using and requiring a differential formulation. It is proved that, for linear interpolation, the stiffness matrix so obtained coincides with the one of the Finite Element Method. For quadratic interpolation, however, the present stiffness matrix differs from that of FEM; moreover it is unsymmetric. It is shown that by using a parabolic interpolation, a convergence of the fourth order is obtained. This is greater than the one More >

S.-P. Zhu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.2, pp. 227-236, 2001, DOI:10.3970/cmes.2001.002.227

Abstract Problems defined on infinite domains must be treated on a finite computational domain. The treatment of the artificially placed boundaries (usually referred to as open boundaries) of such domain truncations can be quite subtle; an over truncation would normally result in large, undesirable reflection of signals back to the computational domain whereas an under truncation would imply an injudicious use of computational resources. In particular, problems occur when strongly nonlinear free surface waves generated in a numerical wave tank are passing through such an open boundary.
In this paper, some recent numerical test results of… More >