S. S. Chen¹, Y. H. Liu^1,2, Z. Z. Cen¹

CMES-Computer Modeling in Engineering & Sciences, Vol.28, No.1, pp. 39-56, 2008, DOI:10.3970/cmes.2008.028.039

Abstract In most engineering applications, solutions derived from the lower-bound theorem of plastic limit analysis are particularly valuable because they provide a safe estimate of the load that will cause plastic collapse. A solution procedure based on the meshless local Petrov-Galerkin (MLPG) method is proposed for lower-bound limit analysis. This is the first work for lower-bound limit analysis by this meshless local weak form method. In the construction of trial functions, the natural neighbour interpolation (NNI) is employed to simplify the treatment of the essential boundary conditions. The discretized limit analysis problem is solved numerically with the reduced-basis technique. The self-equilibrium… More >

B Radhika¹, S S P,a¹, C S Manohar^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.27, No.1&2, pp. 79-110, 2008, DOI:10.3970/cmes.2008.027.079

Abstract Reliability of nonlinear vibrating systems under stochastic excitations is investigated using a two-stage Monte Carlo simulation strategy. For systems with white noise excitation, the governing equations of motion are interpreted as a set of Ito stochastic differential equations. It is assumed that the probability distribution of the maximum in the steady state response belongs to the basin of attraction of one of the classical asymptotic extreme value distributions. The first stage of the solution strategy consists of selection of the form of the extreme value distribution based on hypothesis tests, and the next stage involves the estimation of parameters of… More >

Jyoti Ranjan Majhi, Ranjan Ganguli¹

CMES-Computer Modeling in Engineering & Sciences, Vol.27, No.1&2, pp. 25-36, 2008, DOI:10.3970/cmes.2008.027.025

Abstract The flapping equation for a rotating rigid helicopter blade is typically derived by considering 1) small flap angle, 2) small induced angle of attack and 3) linear aerodynamics. However, the use of nonlinear aerodynamics can make the assumptions of small angles suspect. A general equation describing helicopter blade flap dynamics for large flap angle and large induced inflow angle of attack is derived in this paper with nonlinear aerodynamics . Numerical simulations are performed by solving the nonlinear flapping ordinary differential equation for steady state conditions and the validity of the small angle approximations are examined. It is shown that… More >

Zdenko Tonković¹, Jurica Sorić^1,2, Ivica Skozrit¹

CMES-Computer Modeling in Engineering & Sciences, Vol.26, No.2, pp. 75-90, 2008, DOI:10.3970/cmes.2008.026.075

Abstract An efficient numerical algorithm for modeling of cyclic elastoplastic deformation of shell structures is derived. The constitutive model includes highly nonlinear multi-component forms of kinematic and isotropic hardening functions in conjunction with von Mises yield criterion. Therein, the closest point projection algorithm employing the Reissner-Mindlin type kinematic model, completely formulated in tensor notation, is applied. A consistent elastoplastic tangent modulus ensures high convergence rates in the global iteration approach. The integration algorithm has been implemented into a layered assumed strain isoparametric finite shell element, which is capable of geometrical nonlinearities including finite rotations. Numerical examples, considering the symmetric and nonsymmetric… More >

Q.W. Ma¹

CMES-Computer Modeling in Engineering & Sciences, Vol.23, No.2, pp. 75-90, 2008, DOI:10.3970/cmes.2008.023.075

Abstract In the MLPG_R (Meshless Local Petrove-Galerkin based on Rankine source solution) method, one needs a meshless interpolation scheme for an unknown function to discretise the governing equation. The MLS (moving least square) method has been used for this purpose so far. The MLS method requires inverse of matrix or solution of a linear algebraic system and so is quite time-consuming. In this paper, a new scheme, called simplified finite difference interpolation (SFDI), is devised. This scheme is generally as accurate as the MLS method but does not need matrix inverse and consume less CPU time to evaluate. Although this scheme… More >

P. H. Wen¹, Y. C. Hon²

CMES-Computer Modeling in Engineering & Sciences, Vol.21, No.3, pp. 177-192, 2007, DOI:10.3970/cmes.2007.021.177

Abstract In this paper, we perform a geometrically nonlinear analysis of Reissner-Mindlin plate by using a meshless collocation method. The use of the smooth radial basis functions (RBFs) gives an advantage to evaluate higher order derivatives of the solution at no cost on extra-interpolation. The computational cost is low and requires neither the connectivity of mesh in the domain/boundary nor integrations of fundamental/particular solutions. The coupled nonlinear terms in the equilibrium equations for both the plane stress and plate bending problems are treated as body forces. Two load increment schemes are developed to solve the nonlinear differential equations. Numerical verifications are… More >

Amit Shaw¹, D Roy²

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.1, pp. 69-98, 2007, DOI:10.3970/cmes.2007.019.069

Abstract A novel discretization strategy and derivative reproduction based on reproducing kernel (RK) particle approximations of functions are proposed. The proposed scheme is in the form of an RK interpolation that offers significant numerical advantages over a recent version of the strategy by Chen et al. (2003), wherein the authors added a set of primitive functions to the reproducing kernel (enrichment) functions. It was also required that the support size of the primitive function be less than the smallest distance between two successive grid points. Since the primitive function was required to vary from 0 to 1 within half of this… More >

Q.W. Ma¹

CMES-Computer Modeling in Engineering & Sciences, Vol.18, No.3, pp. 223-234, 2007, DOI:10.3970/cmes.2007.018.223

Abstract Two methods have been recently developed by the author and his group: one called MLPG_R (Meshless Local Petrov-Galerkin method based on Rankine source solution) and the other called QALE-FEM (Quasi Arbitrary Lagrangian-Eulerian Finite Element Method). The former is a meshless method developed from a general MLPG (Meshless Local Petrov-Galerkin) method and is more computationally efficient than the general one when applied to modelling nonlinear water waves. The later is a mesh-based method similar to a conventional finite element method (FEM) when discretizing the governing equations but different from the conventional one in managing the mesh. In this paper, they are… More >

T. Nishioka¹, Y. Kobayashi¹, T. Fujimoto¹

CMES-Computer Modeling in Engineering & Sciences, Vol.17, No.3, pp. 231-238, 2007, DOI:10.3970/cmes.2007.017.231

Abstract First, for growing cracks in elastic-plastic materials, an incremental variational principle is developed to satisfy the boundary conditions near newly created crack surfaces. Then using this variational principle, a moving finite element method is formulated and developed, based on the Delaunay automatic triangulation. Furthermore, theoretical backgrounds on numerical prediction for fracture path of curving crack using T* integral are explained. Using the automatic moving finite element method, fracture-path prediction simulations are successfully carried out. More >

Ivan Christov¹, C.I. Christov², P.M. Jordan³

CMES-Computer Modeling in Engineering & Sciences, Vol.17, No.1, pp. 47-54, 2007, DOI:10.3970/cmes.2007.017.047

Abstract Two widely-used weakly-nonlinear models of acoustic wave propagation --- the inviscid Kuznetsov equation (IKE) and the Lighthill--Westervelt equation (LWE) --- are investigated numerically using a Godunov-type finite-difference scheme. A reformulation of the models as conservation laws is proposed, making it possible to use the numerical tools developed for the Euler equations to study the IKE and LWE, even after the time of shock-formation. It is shown that while the IKE is, without qualification, in very good agreement with the Euler equations, even near the time of shock formation, the same cannot generally be said for the LWE. More >