• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.12 no.2
• /
• pp.43-54
• /
• 1992

For shallow arches with large dynamic loading, linear analysis is no longer considered as practical and accurate. In this study, a method is presented for the dynamic analysis of shallow arches in which geometric nonlinearity must be considered. A program is developed for the analysis of the nonlinear dynamic behavior and for evaluation of critical buckling loads of shallow arches. Geometric nonlinearity is modeled using Lagrangian description of the motion. The finite element analysis procedure is used to solve the dynamic equation of motion and Newmark method is adopted in the approximation of time integration. A shallow arch subject to radial step loads is analyzed. The results are compared with those from other researches to verify the developed program. The behavior of arches is analyzed using the non-dimensional time, load, and shape parameters. It is shown that geometric nonlinearity should be considered in the analysis of shallow arches and probability of buckling failure is getting higher as arches are getting shallower. It is confirmed that arches with the same shape parameter have the same deflection ratio at the same time parameter when arches are loaded with the same parametric load. In addition, it is proved that buckling of arches with the same shape parameter occurs at the same load parameter. Circular arches, which are under a single or uniform normal load, are analyzed for comparison. A parabolic arch with radial step load is also analyzed. It is verified that the developed program is applicable for those problems.

• Journal of Navigation and Port Research
• /
• v.34 no.3
• /
• pp.195-203
• /
• 2010

This paper concerns about total dynamic analysis of integrated mining system. This system consists of vertical steel pipe, intermediate buffer station, flexible pipe and self-propelled miner. The self-propelled miner and buffer are assumed as rigid-body of 6-dof. Discrete models of vertical steel pipe and flexible pipe are adopted, which are obtained by means of lumped-parameter method. The motion of mining vessel is not considered. Instead, the motion of mining vessel is taken into account in form of various boundary conditions (e.g. forced excitation in slow motion and/or fast oscillation and so on). A terramechanics model of extremely cohesive soft soil is applied to the self-propelled miner. Hinged and ball constraints are used to define the connections between sub-systems (vertical steel pipe, buffer, flexible pipe, self-propelled miner). Equations of motion of the coupled model are derived with respect to the each local coordinates system. Four Euler parameters are used to express the orientations of the sub-systems. To solve the equations of motion of the total dynamic model, an incremental-iterative formulation is employed. Newmark-${\beta}$ method is used for time-domain integration. The total dynamic responses of integrated mining system are investigated.

• Journal of Korean Society of Steel Construction
• /
• v.23 no.2
• /
• pp.199-210
• /
• 2011

The cables in cable-stayed bridges are under high stress and are very sensitive to vibration due to their small section areas compared with other members. Therefore, it is reasonable to evaluate the cable impact factor by taking into account the dynamic effect due to moving-vehicle motion. In this study, the cable impact factors were evaluated via moving-vehicle-load analysis, considering the design parameters, i.e., vehicle weight, cable model, road surface roughness, vehicle speed, longitudinal distance between vehicles. For this purpose, two steel composite cable-stayed bridges with 230- and 540-m main spans were selected. The results of the analysis were then compared with those of the influence line method that is currently being used in design practice. The road surface roughness was randomly generated based on ISO 8608, and the convergence of impact factors according to the number of generated road surfaces was evaluated to improve the reliability of the results. A9-d.o.f. tractor-trailer vehicle was used, and the vehicle motion was derived from Lagrange's equation. 3D finite element models for the selected cable-stayed bridges were constructed with truss elements having equivalent moduli for the cables, and with beam elements for the girders and the pylons. The direct integration method was used for the analysis of the bridge-vehicle interaction, and the analysis was conducted iteratively until the displacement error rate of the bridge was within the specified tolerance. It was acknowledged that the influence line method, which cannot consider the dynamic effect due to moving-vehicle motion, could underestimate the impact factors of the end-cables at the side spans, unlike moving-vehicle-load analysis.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.50 no.12
• /
• pp.889-898
• /
• 2022

In this paper, parallel computing method on the three-dimensional electromagnetic field is proposed. The present electromagnetic scattering analysis is conducted based on the time-harmonic vector wave equation and the finite element method. The edge-based element and 2^nd -order absorbing boundary condition are used. Parallelization of the elemental numerical integration and the matrix assemblage is accomplished by allocating the partitioned finite element subdomain for each processor. The graph partitioning library, METIS, is employed for the subdomain generation. The large sparse matrix computation is conducted by MUMPS, which is the parallel computing library based on the multi-frontal method. The accuracy of the present program is validated by the comparison against the Mie-series analytical solution and the results by ANSYS HFSS. In addition, the scalability is verified by measuring the speed-up in terms of the number of processors used. The present electromagnetic scattering analysis is performed for a perfect electric conductor sphere, isotropic/anisotropic dielectric sphere, and the missile configuration. The algorithm of the present program will be applied to the finite element and tearing method, aiming for the further extended parallel computing performance.

• Journal of Korea Water Resources Association
• /
• v.42 no.12
• /
• pp.1053-1067
• /
• 2009

The objective of this study is to develop an efficient and accurate quadratic finite element model based on Streamline Upwind/Petrov Galerkin (SU/PG) scheme for analyzing and predicting two dimensional flow features in complex natural rivers. For a development of model, quadratic tin, quadrilateral and mixed elements as well as linear tin, quadrilateral and mixed elements were used in the model. Also, this model was developed through reinforcement of Gauss Quadrature which was necessary to integral of governing equation. Several tests for bottom-rising channel and U-type channel were performed for the purpose of validation and verification of the developed model. Such results showed that solutions of second order elements are better accurate and improved than those of linear elements. Results obtained by the developed model and RMA-2 model are compared, and the results for the developed model were better accurate than those of RMA-2 model. In the future if the developed model is applied in natural rivers, it can provide better accurate results than those of existing model.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.9 no.4
• /
• pp.159-168
• /
• 2005

This study carries out a finite difference nonlinear analysis of anisotropic advanced composite plate structures with various layer sequences. In the numerical analysis of various mechanical problems involving complex partial differential equations, the finite difference method (FDM) developed in this study has an advantage over the finite element method in its ability to avoid mesh generation and numerical integration. Many studies in FDM have been made on clamped or simple boundary conditions using merely an energy approach. These approaches cannot be satisfied, however, with pivotal points along the free boundary. Therefore, this study addresses the nonlinear problem of anisotropic plates by adopting a finite difference modeling elimination of pivotal difference points in the case of a free boundary condition. Complex nonlinear behaviors of composite plate structures for various parameters, especially for layer sequences, are analyzed using the proposed approach.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.2
• /
• pp.317-324
• /
• 1994

The circulation of Lake Paldang is analysed numerically as an ultimate goal to develop a vehicle predicting the dispersion and concentration of pollutants and sediment flowed into the lake. In finite difference formulation of 2-D depth averaged governing equations. Abbott's 3-time level scheme is employed and for nonlinear terms time centering iteration technique in time and space is used. Model parameters for shear stresses and eddy diffusivities are determined through measured data in and near the lake. Predicted velocities for steady flow are shown to be close to the measured velocities and further improved by taking into account of wind effect. This indicates that the wind effect is needed for proper circulation analysis and it calls for the inclusion of the wind effect. Simulated results of unsteady flow caused by flood inflows and release through Paldang dam are found to characterize the flow features quite well as expected. This implies that the developed model can be used as a tool to analysing the circulation in the lake.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.2 no.2
• /
• pp.47-55
• /
• 1998

Steady-state and transient performance for the medium scale civil aircraft turbofan engine was analyzed. Steady-state performance was analyzed on maximum take-off condition, maximum climb condition, and cruise condition. At 90％RPM of the low pressure compressor, the partload performance was economized. The transient performance was analyzed with cases of the step increase, the ramp increase, the ramp decrease, and the step increase and ramp decrease for the input fuel flow. For the transient performance analysis, work matching between compressor and turbine was needed. Modified Euler method was used the integration of residual torque in work matching equation. At all flight condition, the overshoot of the high pressure turbine inlet temperature was appeared in the step and ramp increase case, and the surge of high pressure compressor was appeared in the step increase case and the ramp increase case within 5.5 seconds of maximum climb condition.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.16 no.2
• /
• pp.207-216
• /
• 2003

In this paper, the finite element method for the elastodynamic axisymmetric fracture analysis is presented in matrix form through the application of the Galerkin method to the time integral equations of motion with no inertia forces. Isoparametric quadratic quadrilateral element and triangular crack tip singular elements with one-quarter node are used in the mesh division of the finite element model. To show the validity and accuracy of the proposed method, the infinite elastic medium with the penny shaped crack is solved first and compared with the analytical solution and the numerical results by the finite difference method and the boundary element method existing in the published literatures, and then the dynamic stress intensity factors of solid and hollow cylinders of finite dimensions haying penny-shaped cracks and internal and external circumferential tracks are computed in detail.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.29 no.3
• /
• pp.167-170
• /
• 2018

The method of moments(MoM) is one of the most popular integral-equation-based full-wave simulation methods, and the multi-level fast multipole method(MLFMM) algorithm can be used for its efficient calculation. When calculating the surface current on the large scatterer in the MoM or MLFMM, iterative methods for the final matrix inversion are used. Among them, BiCGstab(l) has been widely adopted due to its good convergence rate. The number of iterations can be reduced when l becomes larger, but the number of operations per iteration is increased. Herein, we analyze the computational complexity of BiCGstab(l) in the MLFMM method and propose an optimum choice of l.