S. Okamoto¹, Y. Omura¹

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.2, pp. 287-300, 2003, DOI:10.3970/cmes.2003.004.287

Abstract The purpose of this study is to develop a procedure for performing a dynamic analysis in the case that a structure undergoes large translational and rotational displacements when moving along a nonlinear trajectory at variable velocity. Finite-element equations of motion that include the inertial force of the structure's motion have been derived. The equations also account for the geometric nonlinearity that has to be considered in a problem of finite translational and rotational displacements. A finite rotational matrix was used to transfer vectors or matrices measured in a certain coordinate frame to those measured in another coordinate frame. The computational… More >

L. Briseghella¹, C. Majorana¹, P. Pavan¹

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.2, pp. 273-286, 2003, DOI:10.3970/cmes.2003.004.273

Abstract Some aspects of the application of a conservative time integration scheme to the non-linear dynamics of elasto-damaged thin shells are presented. The main characteristic of the scheme is to be conservative, in the sense that it allows the time-discrete system to preserve the basic laws of continuum, namely the balance of the linear and angular momenta as well as the fulfilment of the second law of thermodynamic. Here the method is applied to thin shells under large displacements and rotations. The constitutive model adopted is built coupling the linear elastic model of De Saint Venant-Kirchhoff with a scalar damage function… More >

Wen Yi Lin¹, Kuo Mo Hsiao²

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

Abstract The buckling and postbuckling behavior of spatial rods under different types of end torque and compressive axial force is investigated using finite element method. All coupling among bending, twisting, and stretching deformations for beam element is considered by consistent second-order linearization of the fully geometrically nonlinear beam theory. However, the third order term of the twist rate is also considered. An incremental-iterative method based on the Newton-Raphson method combined with constant arc length of incremental displacement vector is employed for the solution of nonlinear equilibrium equations. The zero value of the tangent stiffness matrix determinant of the structure is used… More >

M. Iura¹, Y. Suetake², S. N. Atluri³

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.2, pp. 249-258, 2003, DOI:10.3970/cmes.2003.004.249

Abstract An accuracy of finite element solutions for 3-D Timoshenko's beams, obtained using a co-rotational formulation, is discussed. The co-rotational formulation has often been used with an assumption that the relative deformations are small. A fundamental question, therefore, has been raised as to whether or not the numerical solutions obtained approach the solutions of the exact theory. In this paper, from theoretical point of view, we investigate the accuracy of the co-rotational formulation for 3-D Timoshenko's beam undergoing finite strains and finite rotations. It is shown that the use of the conventional secant coordinates fails to give satisfactory numerical solutions. We… More >

K. Ijima¹, H. Obiya¹, S. Iguchi², S. Goto²

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.2, pp. 239-248, 2003, DOI:10.3970/cmes.2003.004.239

Abstract Defining element coordinates in space frame, element end deformations become statically clear from the energy principle. Therefore, the deformations can be expressed by nodal displacement without any approximation. The paper indicates that the exact expressions of the deformations and the geometrical stiffness strictly based on the equations makes large displacement analysis of space frame possible with robustness on the computation. More >

P.B. Béda¹

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.2, pp. 231-238, 2003, DOI:10.3970/cmes.2003.004.231

Abstract The paper studies the behavior of a spatial Euler-Bernoulli beam loaded by a terminal thrusting force and a couple. The classical Clebsch-Kirchhoff equilibrium equations are written by using appropriate angular coordinates describing the finite rotations of the local frames attached to each cross-sections of the beam with respect to a fixed system. When we have geometric boundary conditions at one end and dynamic boundary conditions (a force and a couple) at the other the set of equilibrium equations form and initial value probem which can easily be solved with standard Runge-Kutta method. More >

Y. Başar, O. Kintzel¹

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.2, pp. 217-230, 2003, DOI:10.3970/cmes.2003.004.217

Abstract The objective of this contribution is the development of a finite element model for finite rotation and large strain analysis of thin walled shells involving geometry intersections. The shell configuration is described by a linear polynomial in the thickness coordinate. The director of the shell is multiplicatively decomposed into a stretching parameter and an inextensible unit vector whose rotation is accomplished by an updated-rotation formulation. A rotation vector with three independent components is used throughout the shell which permits advantageously to consider smooth shells and compound shells by a unified procedure. This formulation is introduced into an isoparametric four-node element.… More >

M. Iura, S. N. Atluri

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.2, pp. 213-216, 2003, DOI:10.3970/cmes.2003.004.213

Abstract This article has no abstract. More >

J. P. Agnantiaris¹, D. Polyzos^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.1, pp. 197-212, 2003, DOI:10.3970/cmes.2003.004.197

Abstract A Boundary Element Method (BEM), for the three-dimensional solution of both non-axisymmetric and axisymmetric coupled acoustic-elastic problems in the frequency domain, is presented. The present BEM makes use of the Burton and Miller integral equation for infinite acoustic spaces, while elastic structures are dealt with the standard boundary integral equation of elastodynamics. The axisymmetric formulation involves the use of the fast Fourier transform algorithm. Highly accurate numerical algorithms are used for the evaluation of singular integrals, while nearly singular integrals are treated, also with high accuracy, through the use of practical numerical techniques, for both the axisymmetric and non-axisymmetric cases.… More >

Z.Tang, S. Shen, S.N. Atluri¹

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.1, pp. 177-196, 2003, DOI:10.3970/cmes.2003.004.177

Abstract A meshless numerical implementation is reported of the 2-D Fleck-Hutchinson phenomenological strain-gradient theory, which fits within the framework of the Toupin-Mindlin theories and deals with first-order strain gradients and the associated work-conjugate higher-order stresses. From a mathematical point of view, the two-dimensional Toupin-Mindlin strain gradient theory is a generalization of the Poisson-Kirchhoff plate theories, involving, in addition to the fourth-order derivatives of the displacements, also a second-order derivative. In the conventional displacement-based approaches in FEM, the interpolation of displacement requires C$^{1}$ --continuity (in order to ensure convergence of the finite element procedure for 4$^{th}$ order theories), which inevitably involves very… More >