• Journal of the Korean Geotechnical Society
• /
• v.27 no.3
• /
• pp.75-83
• /
• 2011

Emerging issues related with fully coupled Thermo-Hydro-Mechanical (THM) behavior of unsaturated soil demand the development of a numerical tool in diverse geo-mechanical and geo-environmental areas. This paper presents general governing equations for coupled THM processes in unsaturated porous media. Coupled partial differential equations are derived from three mass balances equations (solid, water, and air), energy balance equation, and force equilibrium equation. With Galerkin formulation and time integration of these governing equations, finite element code is developed to find nonlinear solution of four main variables (displacement-u, gas pressure-$P_g$), liquid pressure-$P_1$), and temperature-T) using Newton's iterative scheme. Three cases of numerical simulations are conducted and discussed: one-dimensional drainage experiments (u-$P_g-P_1$), thermal consolidation (u-$P_1$-T), and effect of pile on surrounding soil due to surface temperature variation (u-$P_1$-T).

• Structural Engineering and Mechanics
• /
• v.66 no.1
• /
• pp.15-25
• /
• 2018

Post-construction subgrade settlement especially differential settlement, has become a key issue in construction and operation of non-ballasted track on high-speed railway soil subgrade, which may also affect the dynamic performance of passing trains. To estimate the effect of differential subgrade settlement on the mechanical behaviors of the vehicle-slab track system, a detailed model considering nonlinear subgrade support and initial track state due to track self-weight is developed. Accordingly, analysis aiming at a typical high-speed vehicle coupled with a deteriorated slab track owing to differential subgrade settlement is carried out, in terms of two aspects: (i) determination of an initial mapping relationship between subgrade settlement and track deflections as well as contact state between track and subgrade based on a semi-analytical method; (ii) simulation of dynamic performance of the coupled system by employing a time integration approach. The investigation indicates that subgrade settlement results in additional track irregularity, and locally, the contact between the concrete track and the soil subgrade is prone to failure. Moreover, wheel-rail interaction is significantly exacerbated by the track degradation and abnormal responses occur as a result of the unsupported areas. Distributions of interlaminar contact forces in track system vary dramatically due to the combined effect of track deterioration and dynamic load. These may not only intensify the dynamic responses of the coupled system, but also have impacts on the long-term behavior of the track components.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2000.04a
• /
• pp.126-130
• /
• 2000

In this study, flutter characteristics of a composite wing have been studied for the variation of laminate angles in the subsonic, transonic and supersonic flow regime. The laminate angles are selected by the aspect of engineering practice such as 0, $\pm$45 and 90 degrees. To calculate the unsteady aerodynamics for flutter analysis, the Doublet Lattice Method(DLM) in subsonic flow and the Doublet Point Method(DPM) in supersonic flow are applied in the frequency domain. In transonic flow, transonic small disturbance(TSD) code is used to calculate the nonlinear unsteady aerodynamics in the time domain. Aeroelastic governing equation has been solved by v-g method in the frequency domain and also by Coupled Time-Integration Method(CTIM) in the time domain. from the results of present study, characteristics of free vibration responses and aeroelastic instabilities of a composite wing are presented for the set of various lamination angles in the all flow range.

• Journal of Ocean Engineering and Technology
• /
• v.24 no.5
• /
• pp.16-21
• /
• 2010

During the loading/offloading operation of a liquefied natural gas carrier (LNGC) that is moored in a side-by-side configuration with an offshore plant, sloshing that occurs due to the partially filled LNG tank and the interactive effect between the two floating bodies are important factors that affect safety and operability. Therefore, a time-domain software program, called CHARM3D, was developed to consider the interactions between sloshing and the motion of a floating body, as well as the interactions between multiple bodies using the potential-viscous hybrid method. For the simulation of a floating body in the time domain, hydrodynamic coefficients and wave forces were calculated in the frequency domain using the 3D radiation/diffraction panel program based on potential theory. The calculated values were used for the simulation of a floating body in the time domain by convolution integrals. The liquid sloshing in the inner tanks is solved by the 3D-FDM Navier-Stokes solver that includes the consideration of free-surface non-linearity through the SURF scheme. The computed sloshing forces and moments were fed into the time integration of the ship's motion, and the updated motion was, in turn, used as the excitation force for liquid sloshing, which is repeated for the ensuing time steps. For comparison, a sloshing motion coupled analysis program based on linear potential theory in the frequency domain was developed. The computer programs that were developed were applied to the side-by-side offloading operation between the offshore plant and the LNGC. The frequency-domain results reproduced the coupling effects qualitatively, but, in general, the peaks were over-predicted compared to experimental and time-domain results. The interactive effects between the sloshing liquid and the motion of the vessel can be intensified further in the case of multiple floating bodies.

• Journal of computational fluids engineering
• /
• v.13 no.4
• /
• pp.33-38
• /
• 2008

For study on the unsteady blockage effect, flows around a rotationally oscillating circular cylinder with relatively low forcing frequency in closed test-section wind tunnels have been numerically investigated by solving compressible Navier-Stokes equations. The numerical scheme is based on a node-based finite-volume method with the Roe's flux-difference splitting and an implicit time-integration method coupled with dual time-step sub-iteration. The computed results of the oscillating cylinder in the test section showed that the fluctuations of lift and drag are augmented by the blockage effects. The drag further increases because of low base pressure. The pressure on the test section wall shows the harmonics having the oscillating and the shedding frequencies contained in the blockage effect.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2005.04a
• /
• pp.295-300
• /
• 2005

The freeplay, one of the concentrated structural nonlinearities, is inevitable for control surfaces of a real air vehicle due to normal wear of components and manufacturing mismatches. Also aerodynamic nonlinearities caused by a shock wave occur in transonic region. In practice, these nonlinearities induce the limit cycle oscillation (LCO) and decrease the transonic flutter speed. In this study, the fictitious mass method is used to apply a modal approach to nonlinear structural models due to freeplay. The transonic small-disturbance (TSD) equation is used to calculate unsteady aerodynamic forces in transonic region. Nonlinear aeroelastic time responses are predicted by the coupled time integration method (CTIM). This method was also applied to a 3D all-movable control wing to investigate its nonlinear aeroelastic responses. The angle of attack effect on the LCO characteristics has been found to be closely related with the initial pitching moment.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.05a
• /
• pp.66-72
• /
• 2000

Aerodynamic sensitivity analysis is performed for the Navier-Stokes equations coupled with two-equation turbulence models using a discrete adjoint method and a direct differentiation method respectively. Like the mean flow equations, the turbulence model equations are also hand-differentiated to accurately calculate the sensitivity derivatives of flow quantities with respect to design variables in turbulent viscous flows. Both the direct differentiation code and the adjoint variable code adopt the same time integration scheme with the flow solver to efficiently solve the differentiated equations. The sensitivity codes are then compared with the flow solver in terms of solution accuracy, computing time and computer memory requirements. The sensitivity derivatives obtained from the sensitivity codes with different turbulence models are compared with each other. Using two-equation turbulence models, it is observed that a usual assumption of constant turbulent eddy viscosity in adjoint methods may lead to seriously inaccurate results in highly turbulent flows.

• Advances in concrete construction
• /
• v.6 no.3
• /
• pp.269-283
• /
• 2018

This study presents the effects of geometrical dimensions of concrete gravity dams on the seismic response considering different base width/dam height (L/H) ratios. In the study, a concrete gravity dam with the height of 200 m is selected and finite element models of the dam are constituted including five different L/H ratios such as 0.25, 0.5, 0.75, 1.00, 1.25. All dams are modeled in ANSYS software considering dam-reservoir-foundation interaction. 1989 Loma Prieta earthquake records are applied to models in upstream-downstream direction and linear time history analyses are performed. Dynamic equilibrium equations of motions obtained from the finite element models of the coupled systems are solved by using Newmark time integration algorithm. The seismic response of the models is evaluated from analyses presenting natural frequencies, mode shapes, displacements and principal stresses. The results show that the L/H ratios considerably affect the seismic response of gravity dams. Also, the model where L/H ratio is 1.00 has more desirable results and most appropriate representation of the seismic response of gravity dams.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.12
• /
• pp.2511-2518
• /
• 2002

Dynamic behaviors and stability of an optical disk drive coupled with an automatic ball balancer (ABB) are analyzed by a theoretical approach. The feeding system is modeled a rigid body with six degree-of-freedom. Using Lagrange's equation, we derive the nonlinear equations of motion for a non -autonomous system with respect to the rectangular coordinate. To investigate the dynamic stability of the system in the neighborhood of the equilibrium positions, the monodromy matrix technique is applied to the perturbed equations. On the other hand, time responses are computed by the Runge -Kutta method. We also investigate the effects of the damping coefficient and the position of ABB on the dynamic behaviors of the system.

• 제어로봇시스템학회:학술대회논문집
• /
• 1992.10a
• /
• pp.165-170
• /
• 1992

In this paper, the Hardware-In-The-Loop Simulation(HILS) of missile systems are studied. The HILS is an effective performance evaluation technique that bridges the simulation fidelity gap between analytic all-digital simulations and actual flight tests of missile systems. The HILS may be required to perform system integration tests, performance evaluation at system or subsystem level. Major elements of this HILS facility will include the flight table, simulation computers, I/O computer and peripheral equipments. HILS of missile systems typically involve computer modeling of flight dynamics coupled with a hardware guidance and control(G/C) systems. This paper describes a real time performance evaluation technique of a G/C system, Development of a HILS for a Autopilot of SAM G/C will be used as an example.