• Composites Research
• /
• v.14 no.6
• /
• pp.16-23
• /
• 2001

Analysis of mechanical behavior for a pultruded-wound hollow rod is presented. For this purpose, the pultruded-wound hollow rod is manufactured by the new winder attached to the conventional pultrusion system. And the conventional pultrusion process is newly altered to manufacture pultruded-wound specimens. A computer program, POST II, is modified to perform this study, In the nonlinear finite element formulation, the updated Lagrangian description method based on the second Piolar-Kirchhoff stress tensor and the Green strain tensor are used. For the finite element modeling of the composite hollow rod, the eight-node degenerated shell element is utilized. In order to estimate the failure, the maximum stress criterion is adopted to the averaged stress in the each layer of the finite elements. As numerical examples, the behavior of glass/up composite hollow rod is investigated from the initial loading to the final collapse. Present finite element results considering stiffness degradation and stress unload due to failure shows excellent agreement with experiments in the ultimate load, failure and deformations.

• Atmosphere
• /
• v.23 no.1
• /
• pp.1-11
• /
• 2013

The angular momentum transport of an idealized axisymmetric vortex in the developing stage was investigated using the Weather Research and Forecast (WRF) model. The balanced axisymmetric vortex was constructed based on an empirical function for tangential wind, and the temperature, geopotential, and surface pressure were obtained from the balanced equation. The numerical simulation was carried out for 6 days on the f-plane with the Sea Surface Temperature (SST) set as constant. The weak vortex at initial time was intensified with time, and reached the strength of tropical cyclone in a couple of days. The Absolute Angular Momentum (AAM) was transported along with the secondary circulation of the vortex. Total AAM integrated over a cylinder of radius of 2000 km decreased with simulation time, but total kinetic energy increased rapidly. From the budget analysis, it was found that the surface friction is mainly responsible for the decrease of total AAM. Also, contribution of the surface friction to the AAM loss was about 90% while that of horizontal advection was as small as 8%. The trajectory of neutral numerical tracers following the secondary circulation was presented for the Lagrangian viewpoint of the transports of absolute angular momentum. From the analysis using the trajectory of tracers it was found that the air parcel was under the influence of the surface friction continuously until it leaves the boundary layer near the core. Then the air parcel with reduced amount of angular momentum compared to its original amount was transported from boundary layer to upper level of the vortex and contributed to form the anti-cyclone. These results suggest that the tropical cyclone loses angular momentum as it develops, which is due to the dissipation of angular momentum by the surface friction.

• Journal of the Korea Concrete Institute
• /
• v.20 no.3
• /
• pp.317-323
• /
• 2008

In the current study, A model for the flow analysis of fresh and highly flowable concrete has been developed using a particle method, the moving particle semi-implicit (MPS) method. The phenomena on the flow of concrete has been considered as a visco-plastic flow problem, and the basic governing equation of concrete particle dynamics has been based on the Navier-Stokes equation in Lagrangian form and the conservation of mass. In order to formulate a visco-plastic flow constitutive law of fresh concrete, concrete is modeled as a highly viscous material in the state of non-flow and as a visco-plastic material in the state of flow after reaching the yield stress of fresh concrete. A flow test of fresh concrete in the L-box was simulated and the predicted flow was well matched with the experimental result. The developed method was well showed the flow motion of concrete particles because it was formulated to be based on the motion of visco-plastic fluid dynamics.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.34 no.9
• /
• pp.817-824
• /
• 2010

An improved spray model of a pressure-swirl atomizer was developed and the grid dependency of the model was investigated. Since the Lagrangian-Eulerian approach was adopted for tracking droplets, very small grids could not be used. However, in order to detect swirl flow accurately, small grids were needed because of the consideration of swirl injection. In order to overcome these limitations, numerical studies were performed by using various grids with cell sizes ranging from 10.0 $\times$ 10 mm to 0.625 $\times$ 0.625 mm. From these calculated results, it was observed that the most efficient grid cell size was 1.25 $\times$ 1.25 mm.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.42 no.5
• /
• pp.151-159
• /
• 2000

Dynamic loading of structures often causes excursions of stresses will into the inelastic range and the influence of geometry changes on the response is also significant in may cases. In general , the shell structures designed according to quasi-Static analysis may collapse under condition of dynamic loading. Therefore, for a more realistic prediction on the lad carrying capacity of these shell. both material and geometric nonlinear effects should be considered. In this study , the material nonlinearity effect on the dynamic response is formulated by the elasto-viscoplastic model highly corresponding to the real behavior of the material. Also, the geometrically nonlinear behavior is taken into account using a Total Lagrangian formulation. the reinforcing bars are modeled by the equivalent steel layer at the location of reinforcements, and Von Mises yield criteria is adopted for the steel layer behavior. Also, Drucker-Prager yield criteria is applied for the behavior of concrete. the shape imperfection of dome is assumed as 'dimple type' which can be expressed Wd1=Wd0(1-(r-a)m)n while the shape imperfection of wall is assumed as sinusoidal curve which is Wwi =Wwo sin(n $\pi$y/l). In numerical test, three cases of shape imperfection of 0.0 -5.0cm(opposite direction to loading ; inner shape imperfection)and 5cm (direction to loading : outward shape imperfection) and thickness of steel layer determined by steel ratio of 0,3, and 5% were analyzed. The effect of shape imperfection and steel ratio and behavior characteristics of perfect shape shell and imperfect shape shell are identified through analysis of above mentioned numerical test. Dynamic behaviors of dome and wall according toe combination of shape imperfection and steel ratio are also discussed in this paper.

• Journal of the Korean GEO-environmental Society
• /
• v.7 no.6
• /
• pp.133-143
• /
• 2006

A fully coupled effective stress dynamic analysis procedure for modeling seismic liquefaction on slope is presented. An elasto-plastic formulation is used for the constitutive model UBCSAND in which the yield loci are radial lines of constant stress ratio and the flow rule is non-associated. This is incorporated into the 2D version of Fast Lagrangian Analysis of Continua (FLAC) by modifying the existing Mohr-Coulomb model. This numerical procedure is used to simulate centrifuge test data from the Rensselaer Polytechnic Institute (RPI). UBCSAND is first calibrated to cyclic direct simple shear tests performed on Nevada sand. Both pre- and post-liquefaction behaviour is captured. The centrifuge test is then modeled and the predicted accelerations, excess porewater pressures, and displacements are compared with the measurements. The results are shown to be in general agreement. The procedure is currently being used in the design of liquefaction remediation measures for a number of dam, bridge, tunnel, and pipeline projects in Western Canada.

• Structural Engineering and Mechanics
• /
• v.24 no.6
• /
• pp.647-664
• /
• 2006

In this study, an improved finite element formulation with a scheme of solution for the large deflection analysis of inextensible prismatic and nonprismatic slender beams is developed. For this purpose, a three-noded Lagrangian beam-element with two dependent degrees of freedom per node (i.e., the vertical displacement, y, and the actual slope, $dy/ds=sin{\theta}$, where s is the curved coordinate along the deflected beam) is used to derive the element stiffness matrix. The element stiffness matrix in the global xy-coordinate system is achieved by means of coordinate transformation of a highly nonlinear ($6{\times}6$) element matrix in the local sy-coordinate. Because of bending with large curvature, highly nonlinear expressions are developed within the global stiffness matrix. To achieve the solution after specifying the proper loading and boundary conditions, an iterative quasi-linearization technique with successive corrections are employed considering these nonlinear expressions to remain constant during all iterations of the solution. In order to verify the validity and the accuracy of this study, the vertical and the horizontal displacements of prismatic and nonprismatic beams subjected to various cases of loading and boundary conditions are evaluated and compared with analytic solutions and numerical results by available references and the results by ADINA, and excellent agreements were achieved. The main advantage of the present technique is that the solution is directly obtained, i.e., non-incremental approach, using few iterations (3 to 6 iterations) and without the need to split the stiffness matrix into elastic and geometric matrices.

• Structural Engineering and Mechanics
• /
• v.15 no.1
• /
• pp.83-114
• /
• 2003

A new simple 3-node triangular flat-shell element with standard nodal DOF (6 DOF per node) is proposed for the linear and geometrically nonlinear analysis of very thin to thick plate and shell structures. The formulation of element GT9 (Long and Xu 1994), a generalized conforming membrane element with rigid rotational freedoms, is employed as the membrane component of the new shell element. Both one-point reduced integration scheme and a corresponding stabilization matrix are adopted for avoiding membrane locking and hourglass phenomenon. The bending component of the new element comes from a new generalized conforming Kirchhoff-Mindlin plate element TSL-T9, which is derived in this paper based on semiLoof constrains and rational shear interpolation. Thus the convergence can be guaranteed and no shear locking will happen. Furthermore, a simple hybrid procedure is suggested to improve the stress solutions, and the Updated Lagrangian formulae are also established for the geometrically nonlinear problems. Numerical results with solutions, which are solved by some other recent element models and the models in the commercial finite element software ABAQUS, are presented. They show that the proposed element, denoted as GMST18, exhibits excellent and better performance for the analysis of thin-think plates and shells in both linear and geometrically nonlinear problems.

• Proceedings of the SAREK Conference
• /
• 2008.11a
• /
• pp.546-550
• /
• 2008

In this study, 3-dimensional Computational Fluid Dynamics (CFD) analysis was induced to simulate air flow and particle motion in the axial flow cyclone separator. The commercialized CFD code FLUENT was used to visualize pressure drop and particle collection efficiency inside the cyclone. We simulated 4 cyclone models with different shape of vane, such as turning angle or shape of cross section. For the air flow simulation, we calculated the flow field using standard ${\kappa}-{\varepsilon}$ turbulence viscous model. Each model was simulated with different inlet or outlet boundary conditions. Our major concern for the flow filed simulation was pressure drop across the cyclone. For the particle trajectory simulation, we adopted Euler-Lagrangian approach to track particle motion from inlet to outlet of the cyclone. Particle collection efficiencies of various conditions are calculated by number based collection efficiency. The result showed that the rotation angle of the vane plays major roll to the pressure drop. But the smaller rotation angle of vane causes particle collection efficiency difference with different inlet position.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2012.05a
• /
• pp.71-71
• /
• 2012

현재 전 세계적으로 극한강우의 발생빈도가 점차 높아지고 있으며 홍수량 또한 강도가 커지고 있는 것이 현실이다. 하지만 과거의 홍수발생 빈도에 따라 설계된 홍수방어시설들이 점차 한계를 보이고 있으므로 이를 대비하기위한 구조적 대책뿐만 아니라 홍수피해 발생 가능지역에 사전 예경보를 시행하는 비구조적 대책마련 또한 필요하다. 기존의 홍수예측은 확정적인 하나의 유량예측값만을 제공함으로써 신속하고 편리하였지만 이에 대한 불확실성이 큰 경우 예상치 못한 큰 인적 물적 피해를 가져올 수 있다. 이처럼 확률론적 홍수예측의 필요성이 대두되어 지면서 유럽이나 미국등 선진국에서는 EFFS(European Flood Forecasting System)과 NWSRFS(National Water Service River Forecast System)같이 이미 확률론적 홍수예측에 대한 연구 및 기술개발이 활발하게 진행되어지고 있다. 하지만 홍수예측의 확률론적 접근에 있어서는 많은 불확실성들이 내포되어 있으므로 예측시스템에서 생성된 앙상블 유량예측 결과의 신뢰도 분석과 올바른 불확실성 정보의 제공이 필요하다. 본 연구는 확률론적 홍수예측 방법을 국내에 적용시켜서 기상청의 예측시스템 KLAPS(Korea Local Analysis and Prediction System), MAPLE(McGill Algorithm for Precipitation Nowcasting by Lagrangian Extrapolation), UM(Unified Model) 그리고 MOGREPS(Met Office Global Regional Ensemble Prediction System)으로부터 생성된 기상앙상블을 현재 국토해양부 홍수통제소에서 사용하고 있는 강우-유출모형인 저류함수모형(Storage Function Method)의 입력 자료로 사용한다. 확률론적 홍수예측에서 오는 불확실성을 분석하기 위해서 첫 번째로 제공되는 기상예측 시스템의 시 공간적 스케일 및 대상유역의 공간특성에 따라 어떠한 형태로 전파되어지는지를 분석하였다. 두 번째는 각각의 예측시스템들이 선행기간(Lead time)에 따라 불확실성의 특성이 어떻게 나타나게 되는지를 확인하였다. 이러한 불확실성의 특성을 정확하게 파악하게 된다면 예측에 있어서 현재 갖고 있는 문제점들로부터 개선해 나가야 할 방향을 제시해주어 향후연구에 유용하게 활용될 수 있을 것이다.