• Journal of the Korean Society for Advanced Composite Structures
• /
• v.3 no.3
• /
• pp.22-30
• /
• 2012

This study deals with an enhanced assumed strain (EAS) three-dimensional element for free vibration analysis of laminated composite and sandwich structures. The three-dimensional finite element (FE) formulation based on the EAS method for composite structures shows excellence from the standpoints of computational efficiency, especially for distorted element shapes. Using the EAS FE formulation developed for this study, the effects of side-to-thickness ratios, aspect ratios and ply orientations on the natural frequency are studied and compared with the available elasticity solutions and other plate theories. The numerical results obtained are in good agreement with those reported by other investigators. The new approach works well for the numerical experiments tested, especially for complex structures such as sandwich plates with laminated composite faces.

• Smart Structures and Systems
• /
• v.18 no.1
• /
• pp.155-181
• /
• 2016

In the present paper we discuss the stability and the post-buckling behaviour of thin piezoelastic plates. The first part of the paper is concerned with the modelling of such plates. We discuss the constitutive modelling, starting with the three-dimensional constitutive relations within Voigt's linearized theory of piezoelasticity. Assuming a plane state of stress and a linear distribution of the strains with respect to the thickness of the thin plate, two-dimensional constitutive relations are obtained. The specific form of the linear thickness distribution of the strain is first derived within a fully geometrically nonlinear formulation, for which a Finite Element implementation is introduced. Then, a simplified theory based on the von Karman and Tsien kinematic assumption and the Berger approximation is introduced for simply supported plates with polygonal planform. The governing equations of this theory are solved using a Galerkin procedure and cast into a non-dimensional formulation. In the second part of the paper we discuss the stability and the post-buckling behaviour for single term and multi term solutions of the non-dimensional equations. Finally, numerical results are presented using the Finite Element implementation for the fully geometrically nonlinear theory. The results from the simplified von Karman and Tsien theory are then verified by a comparison with the numerical solutions.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.8
• /
• pp.1429-1437
• /
• 1992

A new finite element formulation based on the relaxation type hereditary integral is presented for a time-domain analysis of isotropic, linear viscoelastic problems. The semi-discrete variational approximation and elastic-viscoelastic correspondence principle are used in the theoretical development of the proposed method. In a time-stepping procedure of final, linear algebraic system equations, only a small additional computation for past history is required since the equivalent stiffness matrix is constant. The viscoelasticity matrices are derived and the stress computation algorithm is given in matrix form. The effect of time increment and Gauss point numbers on the numerical accuracy is examined. Two dimensional numerical examples of plane strain and plane stress are solved and compared with the analytical solutions to demonstrate the versatility and accuracy of the present method.

• Journal of Mechanical Science and Technology
• /
• v.19 no.12
• /
• pp.2187-2196
• /
• 2005

In this paper, the stiffness and the mass matrices for the in-plane motion of a thin circular beam element are derived respectively from the strain energy and the kinetic energy by using the natural shape functions of the exact in-plane displacements which are obtained from an integration of the differential equations of a thin circular beam element in static equilibrium. The matrices are formulated in the local polar coordinate system and in the global Cartesian coordinate system with the effects of shear deformation and rotary inertia. Some numerical examples are performed to verify the element formulation and its analysis capability. The comparison of the FEM results with the theoretical ones shows that the element can describe quite efficiently and accurately the in-plane motion of thin circular beams. The stiffness and the mass matrices with respect to the coefficient vector of shape functions are presented in appendix to be utilized directly in applications without any numerical integration for their formulation.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.27 no.4
• /
• pp.313-320
• /
• 1991

The paper describes the formulation for the analysis of the flexural free vibration of rectangular plate structure by the transfer influence coefficient method, which was developed on the base of the concept of the successive transmission of dynamic influence coefficients. For the analysis of rectangular plate which two opposite sides are simply-supported edge condition, the results of simple numerical examples demonstrate the validity of the present method, that is, the numerical high accuracy, the high speed and the flexibility for programming, compared with results of the transfer matrix method and exact solution or Leissa's method.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.23 no.3
• /
• pp.111-116
• /
• 1987

The accurate prediction of the ship wave resistance is very important to design ships which operate satisfactorily in a wave environment. Thus, work should continue on development and validation of methods to compute ship wave patterns and wave resistance. Research efforts to improve the prediction of ship waves and wavemaking resistance are categorized in two major areas. First is the development of higher-order theories to take account of the nonlinear effect of the free surface condition and improved analytical treatment of the body boundary condition. Second is the development of direct numerical methods aimed at solving body and free-surface boundary conditions as accurately as possible. A new formulation of the slender body theory for a ship with constant speed is developed by Maruo. It is quite different from the existing slender ship theory by Vossers, Maruo and Tuck. It may be regarded as a substitute for the Neumann-Kelvin approximation. In present work, the method of asymptotic expansion of the Kelvin source is applied to obtain a new wave resistance formulation in fluid of finite depth. It takes a simple form than existing theory.

• Tunnel and Underground Space
• /
• v.10 no.2
• /
• pp.173-179
• /
• 2000

The influence of the plastic flow rules on the elasto-plastic behaviour of a discrete joint element was investigated by performing the numerical direct shear tests under both constant normal displacement and normal displacement conditions. The finite interface elements obeying Plesha’s joint constitutive law was used to allow the relative motion of the rock blocks on the joint surface. Realistic results were obtained in the tests adopting the non-associated flow rule, while the associated flow rule overestimated the joint dilation. To overcome the computational drawbacks coming from the non-symmetric element stiffness matrix in the conventional non-associated plasticity, the symmetric formulation of the tangential stiffness matrix for a non-associated joint element was proposed. The symmetric elasto-plastic matrix it derived by assuming an imaginary equivalent joint with associated flow rule which shows the same plastic response as that of original Joint with non-associated flow rule. The validity of the formulation was confirmed through the numerical direct shear tests under constant normal stress condition.

• The Journal of the Acoustical Society of Korea
• /
• v.10 no.1
• /
• pp.12-19
• /
• 1991

A numerical analysis is carried out for the nonlinear phenomena of the bubble oscillator. The model is based on the Keller's formulation for the bubble dynamics. Interpretation of the bubble interior is based on the formulation by Prosperetti. His formulation adopts the energy equation for the analysis of the bubble interior. The numerical simulation Shows typical nonlinear phenomena in its frequency response. Among such nonlinear aspects are the jump phenomenon, the shift of natural frequency of the system, and the appearance of superharmonic resonances. It is deduced that the nonlinear frequency response is dependent upon the initial condition of the bubble oscillator and some multi-valued frequency region can appear in the response curve. Nonlinear phenomena appeared in the bubble oscillator is compared with those of the Duffing equation and it may be said that the bubble dynamic equation has similar nonlinear aspects to the Duffing equation.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.5A
• /
• pp.457-465
• /
• 2009

This study presents a moving least squares (MLS) finite difference method for solving elastic crack problems with stress singularity at the crack tip. Near-tip functions are intrinsically employed in the MLS approximation to model near-tip field inducing singularity in stress field. employment of the functions does not lose the merit of the MLS Taylor polynomial approximation which approximates the derivatives of a function without actual differentiating process. In the formulation of crack problem, computational efficiency is considerably improved by taking the strong formulation instead of weak formulation involving time consuming numerical quadrature Difference equations are constructed on the nodes distributed in computational domain. Numerical experiments for crack problems show that the intrinsically enriched MLS finite difference method can sharply capture the singular behavior of near-tip stress and accurately evaluate stress intensity factors.

• Coupled systems mechanics
• /
• v.4 no.1
• /
• pp.67-98
• /
• 2015

Numerical prediction of dynamic behavior of fully coupled saturated porous media is of great importance in many engineering problems. Specifically, static and dynamic response of soils - porous media with pores filled with fluid, such as air, water, etc. - can only be modeled properly using fully coupled approaches. Modeling and simulation of static and dynamic behavior of soils require significant Verification and Validation (V&V) procedures in order to build credibility and increase confidence in numerical results. By definition, Verification is essentially a mathematics issue and it provides evidence that the model is solved correctly, while Validation, being a physics issue, provides evidence that the right model is solved. This paper focuses on Verification procedure for fully coupled modeling and simulation of porous media. Therefore, a complete Solution Verification suite has been developed consisting of analytical solutions for both static and dynamic problems of porous media, in time domain. Verification for fully coupled modeling and simulation of porous media has been performed through comparison of the numerical solutions with the analytical ones. Modeling and simulation is based on the so called, u-p-U formulation. Of particular interest are numerical dispersion effects which determine the level of numerical accuracy. These effects are investigated in detail, in an effort to suggest a compromise between numerical error and computational cost.