G. M. Kulikov¹, S. V. Plotnikova¹

CMES-Computer Modeling in Engineering & Sciences, Vol.28, No.1, pp. 15-38, 2008, DOI:10.3970/cmes.2008.028.015

Abstract This paper presents a robust non-linear geometrically exact four-node solid-shell element based on the first-order seven-parameter equivalent single-layer theory, which permits us to utilize the 3D constitutive equations. The term "geometrically exact" reflects the fact that geometry of the reference surface is described by analytically given functions and displacement vectors are resolved in the reference surface frame. As fundamental shell unknowns six displacements of the outer surfaces and a transverse displacement of the midsurface are chosen. Such choice of displacements gives the possibility to derive strain-displacement relationships, which are invariant under arbitrarily large rigid-body shell More >

Adan Vega¹, Sherif Rashed², Yoshihiko Tango³, Morinobu Ishiyama³, Hidekazu Murakawa²

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

Abstract Experimental observations have shown that the inherent deformation produced by multi-heating lines is not a simple addition of the inherent deformation produced by single heating lines. Therefore, to accurately predict inherent deformation, the method of superposing inherent deformation of single heating lines is not appropriate. To overcome this difficulty, the authors investigate the influence of multi-heating lines on line heating inherent deformation. First, the influence of previous heating lines on inherent deformation of overlapping, parallel and crossing heating lines is clarified. The influence of the proximity to plate side edge on inherent deformation is also More >

F. Cosmi¹

CMES-Computer Modeling in Engineering & Sciences, Vol.27, No.3, pp. 187-196, 2008, DOI:10.3970/cmes.2008.027.187

Abstract The Cell Method is a recent numerical method that can be applied in several fields of physics and engineering. In this paper, the elastodynamics formulation is extended to include system internal damping, highlighting some interesting characteristics of the method. The developed formulation leads to an explicit solving system. The mass matrix is diagonal (without lumping) and in the most general case a time-dependent damping coefficient can be defined for each node. \newline Accuracy and convergence rate have been tested with reference to the classical problem of a particle free vibration with viscous damping.
An application More >

Chih-Ping Wu¹, Kuan-Hao Chiu, Yun-Ming Wang

CMES-Computer Modeling in Engineering & Sciences, Vol.27, No.3, pp. 163-186, 2008, DOI:10.3970/cmes.2008.027.163

Abstract A differential reproducing kernel particle (DRKP) method is proposed and developed for the analysis of simply supported, multilayered elastic and piezoelectric plates by following up the consistent concepts of reproducing kernel particle (RKP) method. Unlike the RKP method in which the shape functions for derivatives of the reproducing kernel (RK) approximants are obtained by directly taking the differentiation with respect to the shape functions of the RK approximants, we construct a set of differential reproducing conditions to determine the shape functions for the derivatives of RK approximants. On the basis of the extended Hellinger-Reissner principle, More >

Chia-Cheng Tsai¹

CMES-Computer Modeling in Engineering & Sciences, Vol.27, No.3, pp. 151-162, 2008, DOI:10.3970/cmes.2008.027.151

Abstract In this paper we develop analytical particular solutions for the polyharmonic and the products of Helmholtz-type partial differential operators with Chebyshev polynomials at right-hand side. Our solutions can be written explicitly in terms of either monomial or Chebyshev bases. By using these formulas, we can obtain the approximate particular solution when the right-hand side has been represented by a truncated series of Chebyshev polynomials. These formulas are further implemented to solve inhomogeneous partial differential equations (PDEs) in which the homogeneous solutions are complementarily solved by the method of fundamental solutions (MFS). Numerical experiments, which include More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.27, No.3, pp. 137-150, 2008, DOI:10.3970/cmes.2008.027.137

Abstract For the inverse vibration problem, a Lie-group shooting method is proposed to simultaneously estimate the time-dependent damping and stiffness functions by using two sets of displacement as inputs. First, we transform these two ODEs into two parabolic type PDEs. Second, we formulate the inverse vibration problem as a multi-dimensional two-point boundary value problem with unknown coefficients, allowing us to develop the Lie-group shooting method. For the semi-discretizations of PDEs we thus obtain two coupled sets of linear algebraic equations, from which the estimation of damping and stiffness coefficients can be written out explicitly. The present More >

Mira Mitra¹, S. Gopalakrishnan²

CMES-Computer Modeling in Engineering & Sciences, Vol.27, No.1&2, pp. 125-136, 2008, DOI:10.3970/cmes.2008.027.125

Abstract In this paper, the wave characteristics, namely, the spectrum and dispersion relations of multi-wall carbon nanotubes (MWNTs) are studied. The MWNTs are modeled as multiple thin shells coupled through van der Waals force. Each wall of the MWNT has three displacements, i.e, axial, circumferential and radial with variation along the axial and circumferential directions. The wave characteristics are obtained by transforming the governing differential wave equations to frequency domain via Fourier transform. This transformation is first done in time using fast Fourier transform (FFT) and then in one spatial dimension using Fourier series. These transformed equations More >

K. Jayabal, A. Arockiarajan, S.M. Sivakumar¹

CMES-Computer Modeling in Engineering & Sciences, Vol.27, No.1&2, pp. 111-124, 2008, DOI:10.3970/cmes.2008.027.111

Abstract A three dimensional micromechanically motivated model is proposed here based on firm thermodynamics principles to capture the nonlinear dissipative effects in the polycrystal ferroelectrics. The constraint imposed by the surrounding grains on a subgrain at its boundary during domain switching is modeled by a suitable modification of the switching threshold in a subgrain. The effect of this modification in the dissipation threshold is studied in the polycrystal behavior after due correlation of the subgrain behavior with the single crystal experimental results found in literature. Taking into consideration, all the domain switching possibilities, the volume fractions 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 More >

Sanjay Mittal¹

CMES-Computer Modeling in Engineering & Sciences, Vol.27, No.1&2, pp. 63-78, 2008, DOI:10.3970/cmes.2008.027.063

Abstract Flow past a circular cylinder looses stability at a Reynolds number,Re~47. It has been shown, in the past, that the linear stability analysis (LSA) of the steady state solution can predict not only the critical Re, but also the non-dimensional frequency, St, of the associated instability. For larger Re the non-linear effects become important and the LSA of the steady-state flow does not predict the correct St. It is shown that, in general, the LSA applied to the time-averaged flow can result in useful information regarding its stability. This idea is applied to the Re = 100 flow past More >