• 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.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.15 no.4
• /
• pp.567-578
• /
• 2002

The main purpose of this study is to investigate the dynamic behaviour of containment building in nuclear power plant excited by aircraft impact loading using a lower order 8-node solid element. The yield and failure surfaces for concrete material model is formulated on the basis of Drucker-Prager yield criteria and are assumed to be varied by taking account of the visco-plastic energy dissipation. The standard 8-node solid element has prone to exhibit the element deficiencies and the so-called B bar method proposed by Hughes is therefore adopted in this study. The implicit Newmark method is adopted to ensure the numerical stability during the analysis. Finally, the effect of different levels of cracking strain and several types of aircraft loading are examined on the dynamic behaviour of containment building and the results are quantitatively summarized as a future benchmark.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.9 no.4
• /
• pp.15-25
• /
• 1989

This paper presents the results of the numerical analysis on the behavior of cable-stayed bridge considering geometric nonlinearity of cables. Finite element method is used and geometric nonlinearities are considered on the analysis of cable-stayed bridge. The governing equilibrium equations are derived by the principle of virtual work, and modified Newton-Raphson method and Newmark-${\beta}$ method are employed in response calculations. The validity of this study is demonstrated by comparing the examples with analytical results by other method and testing results.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.8 no.6 s.40
• /
• pp.1-11
• /
• 2004

In the previous study(1), the finite element method was used for the vertical vibration analysis of suspension bridge considering the effects of the shear deformation and the rotary inertia under moving load. This study firstly performs the eigenvalue analysis for the free vertical vibration of suspension bridge using FEM analysis. Next the equations of motion considering interaction between suspension bridge and vehicles/train are derived using mode superposition method. And dynamic analysis was performed using the Newmark $\beta$ Method. Finally through the numerical examples, the dynamic responses of bridges by this study are investigated.

• Earthquakes and Structures
• /
• v.16 no.3
• /
• pp.359-368
• /
• 2019

This paper focuses on the study of dynamic analysis of thick plates resting on Winkler foundation. The governing equation is derived from Mindlin's theory. This study is a parametric analysis of the reflections of the thickness / span ratio, the aspect ratio and the boundary conditions on the earthquake excitations are studied. In the analysis, finite element method is used for spatial integration and the Newmark-${\beta}$ method is used for the time integration. While using finite element method, a new element is used. This element is 17-noded and it's formulation is derived from using higher order displacement shape functions. C++ program is used for the analyses. Graphs are presented to help engineers in the design of thick plates subjected to earthquake excitations. It is concluded that the 17-noded finite element is used in the earthquake analysis of thick plates. It is shown that the changes in the aspect ratio are more effective than the changes in the aspect ratio. The center displacements of the reinforced concrete thick clamped plates for b/a=1, and t/a=0.2, and for b/a=2, and t/a=0.2, reached their absolute maximum values of 0.00244 mm at 3.48 s, and of 0.00444 mm at 3.48 s, respectively.

• Advances in nano research
• /
• v.12 no.6
• /
• pp.567-581
• /
• 2022

The free and forced vibration in addition to electric energy harvesting of a piezoelectric disk resting on two-parameter foundation modeled by modified couple stress as well as Kirchhoff plate theory is probed. The governing equations and boundary conditions are obtained using Hamilton's principle. Then, the free and forced vibration are solved using numerical solutions, generalized differential quadrature method (GDQM) and Newmark-beta method. The forced vibration is resulted from a base excitation load. Also, the possible voltage which can be harvested from this system is obtained using generalized integral quadrature method. The validity of the formulation and solution procedure is confirmed using a compassion study. The impact of parameters such as length effect, inner to outer radius ratio, and foundations parameters on the free and forced vibration as well as energy harvesting is investigated in detail. This paper can be a basis for future studies in the area of piezoelectric harvesters in small scales.

• Journal of the Korean Geotechnical Society
• /
• v.35 no.4
• /
• pp.15-26
• /
• 2019

The stability of the gravity quay wall against earthquakes is evaluated on the basis of the allowable displacement of the wall. To estimate the displacement caused by external forces, empirical equations based on the Newmark sliding block method or numerical analysis are widely used. In numerical analysis, it is possible to analyze precisely a complicated site and structure, but difficult to set the appropriate parameters and environments; there are limitations in obtaining reliable results, depending on one's level of expertise. The Newmark method, with only seismic motions, is widely used because it is simpler than numerical simulations when estimating permanent displacement. However, the empirical equations do not have any parameters for the response characteristics and sliding block of the structure, and sliding blocks being assumed as rigid bodies does not consider the nonlinear behavior of the soil and interaction with the structure. Therefore, in order to evaluate the seismic stability of the gravity quay wall, a newly-developed empirical equation is needed to overcome the above-mentioned limitations. In this study, numerical simulations are performed to analyze the response characteristics of the backfill of the structure, and to propose an optimal method of calculating the active area. For this purpose, finite element analyses were performed to analyze the response characteristics, and stress-strain relationships for various seismic motions. As a result, the response characteristics, sliding block, and failure surface of the backfill vary depending on the input seismic motions.

• 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.

• Journal of the Korean Geotechnical Society
• /
• v.33 no.7
• /
• pp.29-41
• /
• 2017

There are uncertainties about the seismic load caused by seismic waves, which cannot be predicted due to the characteristics of the earthquake occurrence. Therefore, it is necessary to consider these uncertainties by probabilistic analysis. In this paper, procedures to develop a fragility curve that is a representative method to evaluate the safety of a structure by stochastic analysis were proposed for cut slopes. Fragility curve that considers uncertainties of soil shear strength parameters was prepared by Monte Carlo Simulation using pseudo static analysis. The fragility curve considering the uncertainty of the input ground motion was developed by performing time-history seismic analysis using selected 30 real ground input motions and the Newmark type displacement evaluation analysis. Fragility curves are represented as the cumulative probability distribution function with lognormal distribution by using the maximum likelihood estimation method.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.9 no.1
• /
• pp.11-23
• /
• 1989

This paper presents the results of the numerical analysis on the behaviour of mooring system of offshore guyed tower. Finite element method is used and geometric nonlinearities are considered in the analysis of mooring line. The governing equilibrium equations are derived by the principle of virtual work, and modified Newton-Raphson method and Newmark-${\beta}$ method are employed in response calculations. The drag and inertia effects of fluid are included using a Morrison type equation. The influences of changing typical parameters like initial inclination and tension of line at the guy attachment point, the length of clump weight, its unit weight and the anchor line length are examined. The effects of idealising the clump weight as a point load(lumped clump weight) on the behaviour of mooring lines are also discussed. Numerical examples demonstrate the validity and capability of the mathematical formulation.