• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.37 no.1
• /
• pp.9-16
• /
• 2024

Using the proper orthogonal decomposition (POD) based intrusive reduced order model (ROM), the total degrees of freedom of the structural system can be significantly reduced and the critical time step satisfying the conditional stability increases in the explicit time integrations. In this study, therefore, the changes in the critical time step in the explicit time integrations are investigated using both the POD-ROM and Voronoi-cell lattice model (VCLM). The snapshot matrix is composed of the data from the structural response under the arbitrary dynamic loads such as seismic excitation, from which the POD-ROM is constructed and the predictive capability is validated. The simulated results show that the significant reduction in the computational time can be achieved using the POD-ROM with sufficiently ensuring the numerical accuracy in the seismic analyses. In addition, the validations show that the POD based intrusive ROM is compatible with the Voronoi-cell lattice based explicit dynamic analyses. In the future study, the research results will be utilized as an elemental technology for the developments of the real-time predictive models or monitoring system involving the high-fidelity simulations of structural dynamics.

• Structural Engineering and Mechanics
• /
• v.37 no.2
• /
• pp.197-213
• /
• 2011

Concrete is a heterogeneous material exhibiting quasi-brittle behaviour. While homogenization of concrete is commonly accepted in general engineering applications, a detailed description of the material heterogeneity using a mesoscale model becomes desirable and even necessary for problems where drastic spatial and time variation of the stress and strain is involved, for example in the analysis of local damages under impact, shock or blast load. A mesoscale model can also assist in an investigation into the underlying mechanisms affecting the bulk material behaviour under various stress conditions. Extending from existing mesoscale model studies, where use is often made of specialized codes with limited capability in the material description and numerical solutions, this paper presents a mesoscale computational model developed under a general-purpose finite element environment. The aim is to facilitate the utilization of sophisticated material descriptions (e.g., pressure and rate dependency) and advanced numerical solvers to suit a broad range of applications, including high impulsive dynamic analysis. The whole procedure encompasses a module for the generation of concrete mesoscale structure; a process for the generation of the FE mesh, considering two alternative schemes for the interface transition zone (ITZ); and the nonlinear analysis of the mesoscale FE model with an explicit time integration approach. The development of the model and various associated computational considerations are discussed in this paper (Part 1). Further numerical studies using the mesoscale model for both quasi-static and dynamic loadings will be presented in the companion paper (Part 2).

• Earthquakes and Structures
• /
• v.7 no.2
• /
• pp.159-178
• /
• 2014

Although the family methods with unconditional stability and numerical dissipation have been developed for structural dynamics they all are implicit methods and thus an iterative procedure is generally involved for each time step. In this work, a new family method is proposed. It involves no nonlinear iterations in addition to unconditional stability and favorable numerical dissipation, which can be continuously controlled. In particular, it can have a zero damping ratio. The most important improvement of this family method is that it involves no nonlinear iterations for each time step and thus it can save many computationally efforts when compared to the currently available dissipative implicit integration methods.

• Structural Engineering and Mechanics
• /
• v.65 no.1
• /
• pp.91-98
• /
• 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.

• Computational Structural Engineering
• /
• v.8 no.3
• /
• pp.5-10
• /
• 1995

공학분야에서 흔히 접하게 되는 접촉충격 문제들의 해석을 위해 접촉충격의 지배방정식과 효율적인 해석방법인 pinball contact-impact 알고리즘에 대해 기술하였다. 예제로 보여준 충돌하는 두 구조계의 접촉충격 해석과 자동차의 crashworthiness 해석결과들은 전술한 알고리즘이 매우 효율적임을 입증하고 있다. 자동차의 crashworthiness 해석과 같이 해석대상 모델이 크고 복잡한 경우에는 짧은 경과시간 동안의 거동변화를 파악하고자 하더라도 엄청난 계산시간이 소요된다. 따라서 효율적인 crashworthiness 해석을 위해서는 전술한 접촉충격 알고리즘 외에도, 논외의 대상이라 본문에서 소개하지는 않았으나, 강력한 요소(one-point integration element), 효과적인 비선형 동적해석 알고리즘, 해의 정확도를 높일 수 있는 adaptive method, 효율적인 시간적분을 위한 mixed time integration method(subcycling)와 explicit method 등을 복합적으로 혼합하여 사용하여야 한다.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2000.10a
• /
• pp.136-140
• /
• 2000

It is desirable to predict springback quantitatively and accurately for the tool and process design in sheet stamping operations, however, it is blown very difficult. The result of springback analysis by the finite element method is sensitively influenced by numerical factors such as blank element size, number of integration point, punch velocity, contact algorithm etc. In the present work, a parametric study by Taguchi method is performed in order to evaluate the influence of numerical factors on springback Quantitatively and to obtain the combination of numerical factors which yields the best approximation to experimental data. Since springback is determined by the residual stress after forming process, it is important to evaluate stress distribution accurately. The oscillation in the time history curve of stress obtained by explicit FEM says that the stress solution at termination time is in very unstable state. Therefore, a variability study is also carried out in this study in order to assess the stability of implicit springback analysis starting from the stress solution by explicit forming simulation. The 2D draw bending process, one of the NUMISHEET '93 benchmark problems, is adopted as an application model.

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

A stepped specimen which is subjected to step loading is modeled to study the initiation and growth of adiabatic shear band using explicit time integration finite element code. The material model for specimen includes effects of thermal softening, strain hardening and strain rate hardening. Various mesh sizes are tested to check whether they are small enough to model highly localized discontinuous phenomena reasonably well. It is shown that the number of adiabatic shear band depends on impact velocity and it is also shown that the initiation and growth of adiabatic shear band inversely depends on prescribed velocity at the top of specimen.

• The Journal of the Korea Contents Association
• /
• v.11 no.9
• /
• pp.1-8
• /
• 2011

A fuzzy inference technique for real-time textile animation without integration at textile model based Mass-Spring model is introduced. Until now many techniques have used the Mass-Spring model to describe elastically deformable objects like textile. A textile object is able to represent as a deformable surface composed of spring and masses, the movement of textile surface which is analysed through the numerical integration by the fundamental law of dynamics such as Hooke's law. However, the integration methods have 'instability problems' if the explicit Euler's method is applied or 'large amounts of calculation' if the implicit Euler's method is applied. A simple and fast animation technique for Mass-Spring model of a textile with fuzzy inference is proposed. The stabilized simulation result is obtained the state of each mass-point in real-time for the n of mass-points by a relatively simple calculation.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.5
• /
• pp.1551-1562
• /
• 1991

최근까지 발표된 유한 요소법을 이용한 압연 해석 관련 주요 논문들을 정리해 보면 다음과 같다. Li와 Kobayashil는 강소성 유한 요소법(rigidplastic finite element method)을 여러가지 마찰조건에 대하여 해석하였다. 이때 압연롤은 강체 (rigid body)로 시편은 가공경화(workhardening)를 동반한 강소성체로 모델링하였다. Hwang과 Kobayashi는 강소성 유한 요소법을 이용한 평면 변형 압연에서 재료 손실을 최소화하는 예비 성형체(preform)의 설계에 대한 연구를 수행하였다. 이 경우에도 역시 압연롤은 강체로 시편은 가공 경화를 동반한 강소성체와 완전 소성체로 모델링 되었으나, 고착(sticking) 마찰 조건에 대해서만 해석을 수행하였다. Mori와 Osak- ada 그리고 Oda는 약간 압축성이 있는 재료의 평면 변형 압연에 대하여 연구하였다. 이때 압연롤은 강체로 시편은 가공 경화를 동반한 강소성체로 모델링 되었으며 경계 면에서는 Coulomb 마찰을 고려하였다. 이밖에도 오일러(Eulerian) 수식화를 이용한 Dawson과 Thompson, Berman의 해석 결과가 있으며, 또 폭 방향의 변형까지를 고려한 Li와 Kobayashi, Mori와 Osakada의 3차원 해석 결과가 있다.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2002.02a
• /
• pp.274-284
• /
• 2002

Deep drawing process for rectangular drawn section is different with that for axisymmetric circular one. Therefore deep drawing process for rectangular drawn section requires several intermediate steps to generate the final configuration without any significant defect. In this study, finite element analysis for multi-stage deep drawing process for high precision rectangular cases is carried out especially for an extreme aspect ratio. The analysis is performed using rigid-plastic finite element method with an explicit time integration scheme of the commercial program, LS-DYNA3D. The sheet blank is modeled using eight-node continuum brick elements. The results of analysis show that the irregular contact condition between blank and die affects the occurrence of failure, and the difference of aspect ratio in the drawn section leads to non-uniform metal flow, which may cause failure. A series of experiments for multi-stage deep drawing process for the rectangular cases are conducted, and the deformation configuration and the thickness distribution of the drawn rectangular cases are investigated by comparing with the results of the numerical analysis. The numerical analysis with an explicit time integration scheme shows good agreement with the experimental observation.