• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.18 no.4 s.70
• /
• pp.385-393
• /
• 2005

This paper established the dynamic model of a flexible Timoshenko beam capable of geometrical nonlinearities subject to large overall motions by using the finite element method. Equations of motion are derived by using Hamilton principle and are formulated in terms of finite elements in which the nonlinear constraint equations are adjoined to the system using Lagrange multipliers. The Newmark direct integration method and the Newton-Raphson iteration are employed here for the numerical study which is to demonstrate the efficiency of the proposed formulation.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.4
• /
• pp.642-651
• /
• 1988

The present paper develops finite element procedures to calculate displacements, strains and stresses in initially stressed elastic solids subjected to static or time-dependent loading conditions. As a point of departure, we employ Hamilton's principle to obtain nonlinear equations of motion characterizing the displacement in a solid. The equations of motion reduce to linear equations of motion if incremental stresses are assumed to be infinitesimal. In the case of linear problem, finite element solutions are obtained by Newmark's direct integration method and by modal analysis. An analytic solution is referred to compare with the linear finite element solution. In the case of nonlinear problem, finite element solutions are obtained by Newton-Raphson iteration method and compared with the linear solution. Finally, the effect of the order of Gauss-Legendre numerical integration on the nonlinear finite element solution, has been investigated.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2006.05a
• /
• pp.388-391
• /
• 2006

In the development of sheet-handling machinery, it is important to predict the static and dynamic behavior of the sheets with a high degree of reliability because the sheets are fed and stacked at such a high speed. Flexible media behaves geometric nonlinearity of large displacement and small strain. In this paper, static and dynamic analyses of flexible media are performed by FEM considering geometric nonlinearity. Linear stiffness matrix and geometric nonlinear stiffness matrix based on the Co-rotational(CR) approach are derived and numerical simulations are performed by Updated Newton-Raphson(UNR) method and Newmark integration scheme.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2006.11a
• /
• pp.721-724
• /
• 2006

In the development of sheet-handling machinery, it is important to predict the static and dynamic behavior of the sheets with a high degree of reliability. Flexible media is very thin, very light and very flexible so it behaves geometric nonlinearity of large displacement and large rotation but small strain. In this paper, static and dynamic analyses of flexible media are performed by dynamic FEM considering geometric nonlinearity. Mass and tangent stiffness matrices based on the Co-rotational(CR) approach are derived and numerical simulations are performed by full Newton-Raphson(FNR) method and Newmark integration scheme.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.7
• /
• pp.44-50
• /
• 2002

Low-velocity impact on composite sandwich panel has been investigated. Contact force is computed from a proposed modified Hertzian contact law. The Hertzian contact law is constructed by adjusting numerical value of the exponent and reducing the through-the- thickness elastic constant of honeycomb core. The equivalent transverse elastic constant is calculated from the rule of mixture. Nonlinear equation to calculate the contact force is solved by the Newton-Raphson method and time integration is done by the Newmark-beta method. A finite element program for the low-velocity impact analysis is coded by implementing these techniques and an 18-node assumed strain solid element. Behaviors of composite sandwich panels subjected to low-velocity impact are analyzed for various cases with different geometry and lay-ups. It has been found that the present code with the proposed contact law can predict measured contact forces and contact times for most cases within reasonable error bounds.

• Structural Engineering and Mechanics
• /
• v.2 no.4
• /
• pp.395-402
• /
• 1994

An improved finite element-transfer matrix method is applied to the transient analysis of plates with large displacement under various excitations. In the present method, the transfer of state vectors from left to right in a combined finite element-transfer matrix method is changed into the transfer of generally incremental stiffness equations of every section from left to right. Furthermore, in this method, the propagation of round-off errors occurring in recursive multiplications of transfer and point matrices is avoided. The Newmark-${\beta}$ method is employed for time integration and the modified Newton-Raphson method for equilibrium iteration in each time step. An ITNONDL-W program based on this method using the IBM-PC/AT microcomputer is developed. Finally numerical examples are presented to demonstrate the accuracy as well as the potential of the proposed method for dynamic large deflection analysis of plates with random boundaries under various excitations.

• Smart Structures and Systems
• /
• v.19 no.4
• /
• pp.403-413
• /
• 2017

The present article encompasses a nonlinear finite element (FE) and genetic algorithm (GA) based optimal vibration energy harvesting from nonprismatic piezo-laminated cantilever beams. Three cases of cross section profiles (such as linear, parabolic and cubic) are modelled to analyse the geometric nonlinear effects on the output responses such as displacement, voltage, and power. The simultaneous effects of taper ratios (such as breadth and height taper) on the output power are also studied. The FE based nonlinear dynamic equation of motion has been solved by an implicit integration method (i.e., Newmark method in conjunction with the Newton-Raphson method). Besides this, a real coded GA based constrained optimization scheme has also been proposed to determine the best set of design variables for optimal harvesting of power within the safe limits of beam stress and PZT breakdown voltage.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2000.04a
• /
• pp.38-42
• /
• 2000

In this study, the low velocity impact behavior of the composite laminates has been described by using 3 dimensional nonlinear finite elements. To describe the geometric nonlinearity due to large deformation, the dynamic contact problem is formulated using the exterior penalty finite element method on the base of Total Lagrangian formulation. The incremental decomposition is introduced, and the converged solution is attained by Newton-Raphson Method. The Newmark's constant-acceleration time integration algorithm is used. To make verification of the finite element program developed in this study, the solution of the nonlinear static problem with occurrence of large deformation is compared with ABAQUS, and the solution of the static contact problem with indentation is compared with the Hertz solution. And, the solution of low velocity impact problem for isotropic material is verificated by comparison with that of LS-DYNA3D. Finally the contact force of impact response from the nonlinear analysis are compared with those from the linear analysis.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.2 no.2
• /
• pp.87-94
• /
• 1998

Behavioral characteristics of shallow circular arches with dynamic loading and different end conditions are analysed. Geometric nonlinearity is modelled using Lagrangian description of the motion. The finite element analysis procedure is used to solve the dynamic equation of motion, and the Newmark method is adopted in the approximation of time integration. The behavior of arches is analysed using the buckling criterion and non-dimensional time, load and shape parameters which Humphreys suggested. But a new deflection-ratio formula including the effect of horizontal displacement plus vertical displacement is presented to apply for the non-symmetric buckling problems. Through the model analysis, it's confirmed that fix-ended arches have higher buckling stability than hinge-ended arches, and arches with the same shape parameter have the same deflection ratio at the same time parameter when loaded with the same parametric load.

• Coupled systems mechanics
• /
• v.1 no.4
• /
• pp.361-379
• /
• 2012

This paper presents a fully coupled three-dimensional solver for the analysis of interaction between pulsatile flow and large deformation structure. A partitioned time marching algorithm is employed for the solution of the time dependent coupled discretised problem, enabling the use of highly developed, robust and well-tested solvers for each field. Conservative transfer of information at the fluid-structure interface is combined with an effective multi-predict-correct iterative scheme to enable implicit coupling of the interacting fields at each time increment. The three-dimensional unsteady incompressible fluid is solved using a powerful implicit time stepping technique and an ALE formulation for moving boundaries with second-order time accurate is used. A full spectrum of total variational diminishing (TVD) schemes in unstructured grids is allowed implementation for the advection terms and finite element shape functions are used to evaluate the solution and its variation within mesh elements. A finite element dynamic analysis of the highly deformable structure is carried out with a numerical strategy combining the implicit Newmark time integration algorithm with a Newton-Raphson second-order optimisation method. The proposed model is used to predict the wave flow fields of a particular flow-induced vibrational phenomenon, and comparison of the numerical results with available experimental data validates the methodology and assesses its accuracy. Another test case about three-dimensional biomedical model with pulsatile inflow is presented to benchmark the algorithm and to demonstrate the potential applications of this method.