• Structural Engineering and Mechanics
• /
• v.87 no.4
• /
• pp.363-373
• /
• 2023

Since the first family of structure-dependent methods can simultaneously integrate unconditional stability and explicit formulation in addition to second order accuracy, it is very computationally efficient for solving inertial problems except for adopting auto time-stepping techniques due to no nonlinear iterations. However, an unusual stability property is first found herein since its unconditional stability interval is drastically different for zero and nonzero damping. In fact, instability might occur for solving a damped stiffness hardening system while an accurate result can be obtained for the corresponding undamped stiffness hardening system. A technique of using a stability factor is applied to overcome this difficulty. It can be applied to magnify an unconditional stability interval. After introducing this stability factor, the formulation of this family of structure-dependent methods is changed accordingly and thus its numerical properties must be re-evaluated. In summary, a large stability factor can result in a large unconditional stability interval but also lead to a large relative period error. As a consequence, a stability factor must be appropriately chosen to have a desired unconditional stability interval in addition to an acceptable period distortion.

• Structural Engineering and Mechanics
• /
• v.91 no.3
• /
• pp.325-334
• /
• 2024

This paper deals with the dynamic buckling behavior of truncated conical shells composed of carbon nanotube composites, an important area of study in view of their very wide engineering applications in aerospace industries. In this regard, the effective material properties of the nanocomposite have been computed using the Mori-Tanaka model, which has already been established for such analyses. The motion equations ruling the structure's behavior are derived using first order shear deformation theory, Hamilton's principle, and energy method. This will provide adequate background information on its dynamic response. In an effort to probe the dynamic instability region of the structure, differential quadrature method combined with Bolotin's method will be adopted to tackle the resulting motion equations, which enables efficient and accurate analysis. This work considers the effect of various parameters in the geometrical parameters and the volume fraction of CNTs on the structure's DIR. Specifically, it became clear that increasing the volume fraction of CNTs shifted the frequency range of the DIR to higher values, indicating the significant role of nanocomposite composition regarding structure stability.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.10 no.5
• /
• pp.545-565
• /
• 2018

A numerical code is developed based on potential flow theory to investigate nonlinear sloshing in rectangular Liquefied Natural Gas (LNG) tanks under forced excitation. Using this code, internal free-surface elevation and sloshing loads on liquid tanks can be obtained both in time domain and frequency domain. In the mathematical model, acceleration potential is solved in the calculation of pressure on tanks and the artificial damping model is adopted to account for energy dissipation during sloshing. The Boundary Element Method (BEM) is used to solve boundary value problems of both velocity potential and acceleration potential. Numerical calculation results are compared with published results to determine the efficiency and accuracy of the numerical code. Sloshing properties in partially filled rectangular and membrane tank under translational and rotational excitations are investigated. It is found that sloshing under horizontal and rotational excitations share similar properties. The first resonant mode and excitation frequency are the dominant response frequencies. Resonant sloshing will be excited when vertical excitation lies in the instability region. For liquid tank under rotational excitation, sloshing responses including amplitude and phase are sensitive to the location of the center of rotation. Moreover, experimental tests were conducted to analyze viscous effects on sloshing and to validate the feasibility of artificial damping models. The results show that the artificial damping model with modifying wall boundary conditions has better applicability in simulating sloshing under different fill levels and excitations.

• Journal of Ocean Engineering and Technology
• /
• v.22 no.5
• /
• pp.52-60
• /
• 2008

This study investigates the applicability of CADMAS-SURF (Super Roller Flume for Computer Aided Design of a MArtime Structure) code basal on the Navier-Stokes solver to predict bore phenomena. The time variation of ware levels and velocities due to the bore propagation were computed for the different bore strength conditions. In order to verify the numerical results by CADMAS-SURF, laboratory experiments were also performed, using the DPIV and LDV measuring system. The numerical results were compared to the experimental data and the analytical predictions by the NSC method basal on fully nonlinear shallow-water theory by the method of characteristics. It appears that the CADMAS-SURF slightly overestimated the water-surface level measured by the laboratory experiments and its discrepancy becomes prominent as the bore strength increases. The predicted propagation speed for a bore was also slaver than that by the experiment and NSC method. However, the temporal variations in velocities were in relatively good agreement with the experimental data for all cases, except for overshooting and undershooting in the front face of a bore, which may be derived from the numerical instability. Further, CADMAS-SURF successfully simulated the decrease in the water level and velocity caused by the effects of negative waves reflected from the upstream end wall.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.18 no.1
• /
• pp.1-7
• /
• 2010

The paper proposes a multi-body dynamic simulation to numerically evaluate the generated axial force(G.A.F) and plunging resistant force(P.R.F) practically related to the shudder and idling vibration of an automobile. A numerical analysis of two plunging types of CV joints, tripod joint(TJ) and very low axial tripod joint(VTJ), is conducted using the commercial program DAFUL. User-defined subroutines of a friction model illustrating the contacted parts of the outboard and inboard joint are subsequently developed to overcome the numerical instability and improve the solution performance. The Coulomb friction effect is applied to describe the contact models of the lubricated parts in the rolling and sliding mechanisms. The numerical results, in accordance with the joint articulation angle variation, are validated with experimentation. The offset between spider and tulip housing is demonstrated to be the critical role in producing the 3rd order component of the axial force that potentially causes the noise and vibration in vehicle. The VTJ shows an excellent behavior for the shudder when compared with TJ. In addition, a flexible nonlinear contact analysis coupled with rigid multi-body dynamics is also performed to show the dynamic strength characteristics of the rollers, housing, and spider.

• Geomechanics and Engineering
• /
• v.22 no.6
• /
• pp.475-488
• /
• 2020

The recent increase in the use of Expanded Polystyrene (EPS) geofoam in construction and geotechnical projects has driven researchers to investigate its behavior, more deeply. In this paper, a series of experimental tests to investigate the stress-strain behavior and the mechanical properties of EPS blocks, under monotonic axial loading are presented. Four different densities of cylindrically shaped EPS with different dimensions are used to investigate the effects of loading rate, height and diameter, as well as the influence of the density of EPS on the stress-strain response. The results show that increasing the height of the EPS samples leads to instability of the sample and consequent lower resistance to the applied pressure. Large EPS samples show higher Young's modulus and compressive resistance due to some boundary effects. An increase in the rate of loading can increase the elastic moduli and compressive resistance of the EPS geofoam samples, which also varies depending on the density of the samples. It was also determined that the elastic modulus of EPS increases with increasing EPS density. By implementing an efficient numerical procedure, the stress-strain response of EPS geofoam samples can be reproduced with great accuracy. The numerical analysis based on the proposed method can used to evaluate the effect of different factors on the behavior of EPS geofoam.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.1
• /
• pp.1-11
• /
• 2010

The definition of the speed of sound is reexamined since it is crucial in the numerical analysis of compressible real gas flows. The thermodynamic speed of sound (TSS), $a_{th}$, and the characteristic speed of sound (CSS), $a_{ch}$, are derived using generalized equation of state (EOS). It is found that the real gas EOS, for which pressure is not linearly dependent on density and temperature, results in slightly different TSS and CSS. in this formalism, Roe's approximate Riemann solver was derived again with corrections for real gases. The results show a little difference when the speeds of sound are applied to the Roe's scheme and Advection Upstream Splitting Method (AUSM) scheme, but a numerical instability is observed for a special case using AUSM scheme. It is considered reasonable to use of CSS for the mathematical consistency of the numerical schemes. The approach is applicable to multi-dimensional problems consistently.

• Journal of the Korean GEO-environmental Society
• /
• v.10 no.7
• /
• pp.81-91
• /
• 2009

Recently, twin-tunnel is often designed in the aspects of disaster prevention and economical reasons. However, the design cases and the studies are relatively insufficient. By the twin-tunnel excavation, deviate stresses of pillar between tunnels are increased and the increased stresses induce the instability of the twin-tunnel. In this study, numerical analyses about the twin-tunnel behaviour were conducted with varying ground strength, width of pillar and depth of earth cover and a series of regression analyses were carried out by using the results of numerical analyses for the twin-tunnel. Based on the numerical analyses, an estimation method of derived stresses is suggested through the regression analyses. Also, based on the results of regression analyses, an quantitative estimation method considering the reinforcement effects is also suggested. Then various parametric studies were conducted to be considered the reinforcement type and various design parameters. Finally, the efficiency of the suggested method based on the Hoek-Brown Failure Criterion is verified through the results of parametric studies.

• Journal of the Korean Society of Safety
• /
• v.21 no.2 s.74
• /
• pp.143-149
• /
• 2006

This paper deals with the space-time finite element analysis of two-dimensional vibration problem with a single variable. The method of space-time finite elements enables the simpler solution than the usual finite element analysis with discretization in space only. We present a discretization technique in which finite element approximations are used in time and space simultaneously for a relatively large time period. The weighted residual process is used to formulate a finite element method for a space-time domain. A stability problem is described and some investigations for chosen type of rectangular space-time finite elements are carried out. Instability is caused by a too large time step of successive time steps in the traditional time-dependent problems. It has been shown that the numerical stability of time-stepping on the larger time steps is quite good. The unstructured space-time finite element not only overcomes the shortcomings of the stability in the traditional numerical methods, but it is also endowed with the features of an effective computational technique. Some numerical examples have been presented to illustrate the efficiency of the described method.

• Journal of the Korean Geosynthetics Society
• /
• v.13 no.4
• /
• pp.45-56
• /
• 2014

Frozen ground in cold region consists of an upper active layer and lower permafrost which is permanently frozen land. During the summer season, the air temperature is high enough to make the frozen ground melt, which causes the reduction of soil strength and thaw settlement. These phenomena result in structural instability, so it is necessary to apply frozen ground stability techniques. Thermosyphon is a closed natural two-phase convection device to maintain the ground temperature below $0^{\circ}C$ by extracting heat from the ground and discharges it into the atmosphere. Experimental and numerical investigation has been performed to estimate the effect of the refrigerant filling ratio in thermosyphon using R-134a refrigerant and the thermal conductance of the thermosyphon.