• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1996.04a
• /
• pp.95-102
• /
• 1996

The stress intensity factor for an interface crack in dissimilar elastic and viscoelastic materials is derived and the time-domain boundary element analysis is performed. Numerical results show that the proposed method is very useful for the analysis of the interface crack in elastic and viscoelastic bimaterials.

• Structural Engineering and Mechanics
• /
• v.65 no.1
• /
• pp.121-128
• /
• 2018

Three structure-dependent integration methods with no numerical dissipation have been successfully developed for time integration. Although these three integration methods generally have the same numerical properties, such as unconditional stability, second-order accuracy, explicit formulation, no overshoot and no numerical damping, there still exist some different numerical properties. It is found that TLM can only have unconditional stability for linear elastic and stiffness softening systems for zero viscous damping while for nonzero viscous damping it only has unconditional stability for linear elastic systems. Whereas, both CEM and CRM can have unconditional stability for linear elastic and stiffness softening systems for both zero and nonzero viscous damping. However, the most significantly different property among the three integration methods is a weak instability. In fact, both CRM and TLM have a weak instability, which will lead to an adverse overshoot or even a numerical instability in the high frequency responses to nonzero initial conditions. Whereas, CEM possesses no such an adverse weak instability. As a result, the performance of CEM is much better than for CRM and TLM. Notice that a weak instability property of CRM and TLM might severely limit its practical applications.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2009.10a
• /
• pp.281-284
• /
• 2009

Flexible forming process for sheet metal using reconfigurable die is introduced based on numerical simulation. Numerical simulation of sheet metal forming process is carried out by using flexible dies model instead of conventional matched die set. Elastic cushion which has high resilience behavior from excessive deformation are inserted between forming punches and blank material for smoothing the forming surface which has discrete due to characteristics of the flexile die. As an elastic cushion, urethane pads are utilized using hyperelastic material model in the simulation. Formability in view of surface defect such as onset of dimple is compared with regard to various punch sizes. Consequently, it is confirmed that the flexible forming process for sheet material has appropriate capability and feasibility for manufacturing of smoothly curved surface instead of conventional die forming process.

• The Journal of the Acoustical Society of Korea
• /
• v.17 no.3E
• /
• pp.35-42
• /
• 1998

Transient acoustic and elastic wave propagation in inhomogeneous media are studied by using the Fourier method. It is known that the fourier method has advantages in memory requirements and computing speed over conventional methods such as FDM and FEM, because the Fourier method needs less grid points for achieving the same accuracy. To verify the proposed numerical scheme, several examples having analytic solutions are considered, where two different semi-infinite media are in contact along a plane boundary. The comparisons of numerical results by the Fourier method and analytic solutions show good agreements. In addition, the fourier method is applied to a layered half-plane, in which an elastic semi-infinite medium is covered by an elastic layer of finite thickness. It is showed how to derive the analytic solutions by using the Cagniard-de Hoop method. The numerical solutions are in excellent agreements with analytic results.

• Computational Structural Engineering
• /
• v.9 no.4
• /
• pp.155-163
• /
• 1996

In the analysis of thin elastic structures such as beam-, arch-, plate- and shell-like bodies using standard finite element schemes, there may occur deterioration of approximation quality owing to shear and membrane lockings. Moreover, a recognition of this phenomenon in the computed numerical results is not easy without comparing with other available reference numerical data. This paper analyses briefly this phenomenon and introduces one inexpensive but reliable a posteriori locking detection method. Numerical examples are given supporting the theoretical results.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 1998.10a
• /
• pp.399-406
• /
• 1998

Transient acoustic and elastic wave propagation in inhomogeneous media are studied by using the Fourier method. To verify the proposed numerical scheme, several examples having analytic solutions are considered, where two different semi-infinite media are in contact along a plane boundary. The comparisons of numerical results by the Fourier method and analytic solutions show good agreements. In addition, the Fourier method is applied to a layered half-plane, in which an elastic semi-infinite medium is covered by an elastic layer of finite thickness. It is showed how to derive the analytic solutions by using the Cagniard-de Hoop method. The numerical solutions are in excellent agreements with analytic results.

• Structural Engineering and Mechanics
• /
• v.22 no.2
• /
• pp.197-222
• /
• 2006

In this paper, the non-linear time-dependent closed-form, discrete and combined solutions for the post-elastic response of a geometrically and physically non-linear suspended cable to a uniformly distributed load considering the creep effects, are presented. The time-dependent closed-form method for the particularly straightforward determination of a vertical uniformly distributed load applied over the entire span of a cable and the accompanying deflection at time t corresponding to the elastic limit and/or to the elastic region, post-elastic and failure range of a suspended cable is described. The actual stress-strain properties of steel cables as well as creep of cables and their rheological characteristics are considered. In this solution, applying the Irvine's theory, the direct use of experimental data, such as the actual stress-strain and strain-time properties of high-strength steel cables, is implemented. The results obtained by the closed-form solution, i.e., a load corresponding to the elastic limit, post-elastic and failure range at time t, enable the direct use in the discrete non-linear time-dependent post-elastic analysis of a suspended cable. This initial value of load is necessary for the non-linear time-dependent elastic and post-elastic discrete analysis, concerning incremental and iterative solution strategies with tangent modulus concept. At each time step, the suspended cable is analyzed under the applied load and imposed deformations originated due to creep. This combined time-dependent approach, based on the closed-form solution and on the FEM, allows a prediction of the required load that occurs in the post-elastic region. The application of the described methods and derived equations is illustrated by numerical examples.

• Structural Engineering and Mechanics
• /
• v.26 no.4
• /
• pp.441-457
• /
• 2007

A new numerical model based on the spline finite strip method is presented here for the analysis of buckling of built-up columns with and without end stay plates. The channels are modelled with spline finite strips while the connecting elements are represented by a 3D beam finite element, for which the stiffness matrix is modified in order to ensure complete compatibility with the strips. This numerical model has the advantage to give all possible failure modes of built-up columns for different boundary conditions. The end stay plates are also taken into account in this method. To validate the model a comparative study was carried out. First, a general procedure was chosen and adopted. For each numerical analysis, the lowest buckling loads and modes were calculated. The basic or "pure" buckling modes were identified and their critical loads were compared with solutions obtained using analytical methods and/or other numerical methods. The results showed that the proposed numerical model can be used in practice to study the elastic buckling of built-up columns. This model is considered accurate and efficient for the local buckling of short columns and global buckling for slender columns.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1998.10a
• /
• pp.397-404
• /
• 1998

A plate on the elastic half-space may generally be analyzed by the finite element method. However, there are some difficulties to obtain the flexibility matrix of the foundation based on the Boussinesq's theory. In this study, an efficient numerical procedure which uses the analysis results of the vertical displacements due to the uniformly distributed loading in a circular area is presented. Some numerical examples represent better results than those of numerical integration technique or subsection method especially in the case of irregular mesh pattern.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.34 no.4
• /
• pp.501-509
• /
• 2010

A floating crane is a crane-mounted ship and is used to assemble or to transport heavy blocks in shipyards. In this paper, the static and dynamic response of a floating crane and a heavy block that are connected using elastic booms and wire ropes are described. The static and dynamic equations of surge, pitch, and heave for the system are derived on the basis of flexible multibody system dynamics. The equations of motion are fully coupled and highly nonlinear since they involve nonlinear mass matrices, elastic stiffness matrices, quadratic velocity vectors, and generalized external forces. A floating frame of reference and nodal coordinates are employed to model the boom as a flexible body. The nonlinear hydrostatic force, linear hydrodynamic force, wire-rope force, and mooring force are considered as the external forces. For numerical analysis, the Hilber-Hughes-Taylor method for implicit integration is used. The dynamic responses of the cargo are analyzed with respect to the results obtained by static and numerical analyses.