• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.9 no.3
• /
• pp.65-73
• /
• 1989

Continuum formulations for the expressions of dynamic energy release rates and computational methods for dynamic stress intensity factors are developed for the analysis of dynamic fracture problems subjected to stress wave loading. Explicit volume integral expressions for instantaneous dynamic energy release rates are derived by modeling virtual crack extensions with the dynamic Eulerian-Lagrangian kinematic description. In the finite element applications a finite region around a crack-tip is modeled by using quarter-point singular isoparametric elements, and the volume integrals are evaluated for each crack-tip element during virtual crack extensions while the singularity is maintained. It is shown that the use of the present method is more reliable and accurate for the dynamic fracture analysis than that of other path-independent integral methods when the effects of stress waves are significant.

• Journal of the Korean Society for Precision Engineering
• /
• v.29 no.6
• /
• pp.632-639
• /
• 2012

A fully coupled thermo-mechanical model is adopted to study the temperature distribution and the material deformation in friction stir welding(FSW) process. Rotational speed is most important parameters in this research. Three dimension results under different process parameters were presented. Result indicate that the maximum temperature is lower than the melting point of the welding material. The higher temperature gradient occurs in the leading side of the workpiece. The maximum temperature can be increased with increasing the tool angular velocity, rpm in the current numerical modeling. In this research ABAQUS Ver.6.7 is to analyze a fully coupled thermo-mechanical model. ALE(Arbitrary Lagrangian-Eulerian) finite element formulation is used for the large deformation in FSW process and using the Mass scaling for the analysis time efficiency.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.43 no.3
• /
• pp.85-93
• /
• 2001

In this study, a computer program considering shape imperfections of arch under dynamic loading was developed. The shape imperfection of arch was assumed as higher degree polynomial expressed as $\omega$$_{i}$ = $\omega$$_{o}$ (1-(2$\chi$/L)$^{m}$ )$^n$and sinusoidal curve such as $\omega$$_{i}$ = $\omega$$_{o}$ sin(η$\pi$$\chi$/L). In finite element formulation, the material nonlinear behavior was assumed the elasto-viscoplastic model highly corresponding to the real behavior of the material and the geometrically nonlinear behavior was modeled using Lagrangian description of motion. Also, the behavior of steel was modeled by applying yield criteria of Von Mises. The developed program was applied to the analysis of the dynamic behavior for the clamped beam subjected to the concentrated load at midspan and the results were compared with those from other research to investigate accuracy of the presented finite element program. In numerical examples, the shape imperfections of L/500, L/1,000 and L/2,000 were considered and the modes of shape imperfections of the symmetric and antisymmetric were adopted. The effects of the shape imperfections on the dynamic behavior of arch were conspicuous and results of analysis indicate that the reasonable values of arch rise to arch span ratio ranged between 0.1 and 0.3.

• Particle and aerosol research
• /
• v.5 no.3
• /
• pp.123-131
• /
• 2009

This paper presents simulation results of particle transport in low-pressure, low-temperature plasma environment. The size dependent transport of particles in the plasma is investigated with a two-dimensional simulation tool developed in-house for plasma chamber analysis and design. The plasma model consists of the first two and three moments of the Boltzmann equation for ion and electron fluids respectively, coupled to Poisson's equation for the self-consistent electric field. The particle transport model takes into account all important factors, such as gravitational, electrostatic, ion drag, neutral drag and Brownian forces, affecting the motion of particles in the plasma environment. The particle transport model coupled with both neutral fluid and plasma models is simulated through a Lagrangian approach tracking the individual trajectory of each particle by taking a force balance on the particle. The size dependant trap locations of particles ranging from a few nm to a few ${\mu}m$ are identified in both electropositive and electronegative plasmas. The simulation results show that particles are trapped at locations where the forces acting on them balance. While fine particles tend to be trapped in the bulk, large particles accumulate near bottom sheath boundaries and around material interfaces, such as wafer and electrode edges where a sudden change in electric field occurs. Overall, small particles form a "dome" shape around the center of the plasma reactor and are also trapped in a "ring" near the radial sheath boundaries, while larger particles accumulate only in the "ring". These simulation results are qualitatively in good agreement with experimental observation.

• Computational Structural Engineering
• /
• v.11 no.4
• /
• pp.155-168
• /
• 1998

보강된 판 및 쉘구조의 안정성 및 후좌굴을 포함하는 기하학적 비선형 해석을 수행하기 위하여, total Lagrangian formulation에 근거한 연속체의 증분평형방정식으로부터 변형된 쉘요소인 유한요소이론을 제시하였다. 쉘구조의 곡률이 불연속적으로 변하거나 쉘부재들이 유한한 각도로 만나는 보강된 판 및 쉘구조의 비선형 해석이 가능하도록 주부재와 보강재 간의 연결점에 대한 일반적인 변환관계를 제시하였으며 좌굴해석 및 기하학적 비선형해석의 경우에 해의 정확성 및 수렴성을 개선시키기 위하여 접선강도행렬 산정시 회전각의 2차항을 포함시켰다. 또한, shear locking 현상을 극복하기 위하여 감차적분을 적용하였고 쉘구조의 좌굴해석에서는 power method를 적용하여 해석의 효율을 높였으며, 후좌굴해석에서는 변위 및 하중증분법을 적절히 결합시켜 보강된 쉘구조의 후좌굴 거동추적이 용이하였다. 또한, 입력자료를 손쉽게 준비하고 좌굴모드 및 후좌굴거동을 효율적으로 분석하기 위하여 전, 후 처리 프로그램을 개발하였고 다양한 해석예제를 통하여 다른 문헌의 해석결과를 비교함으로써 본 연구에서 개발된 유한요소 해석프로그램의 타당성 및 정확성을 입증하였다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.16 no.1
• /
• pp.77-86
• /
• 2003

Dynamic response of baffled fuel-storage tank in turnaround motion is simulated using the ALE finite element method. Fuel-storage tank undergoes abrupt impact load caused by inertia force of internal fuel in turnaround motion. Also, large dynamic force and moment caused by this load influence structural stability and control system. In this paper, ring-type baffles are adopted to suppress the dynamic influence. Through the parametric analysis with respect to the baffle number and location, the effects of baffle on the dynamic response of baffled fuel-storage tank is analyzed. The ALE finite element method is adopted for the accurate and effective simulation of the hydrodynamic interaction between fluid and structure.

• Proceedings of the KSR Conference
• /
• 2003.05a
• /
• pp.376-381
• /
• 2003

This study is focused on the dynamic response of curved bridge when the rubber tired AGT vehicles is running with alternative articulations. For the analytic approach, there is necessary for the three dimensional vehicle model with 11 degree of freedom and the three dimensional curved bridge model by means of finite element method. It can be described by conventional Lagrangian formula with respect to the dynamic interactions between vehicles and its met bridge. The formula is implemented by Fortran language on the simulation program designated BADIA II(Bridge-AGT Dynamic Interaction Analysis II). The solutions of the formula are derived by Newmark- ${\beta}$ method. The BADIA II is for the dynamic interactions between vehicle and curved bridge in terms of the roughness of running surface and guide rail. The applicability of the BADIA II is verified in terms of displacement and modal frequency. This study is described that the dynamic interactive behaviors between the rubber tired AGT vehicle and curved bridge in terms of the radius of curvatures of curved bridge, vehicle articulations, vehicle speeds, vehicle weights, flatness of running surface and roughness of guide rail using BADIA II.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.04a
• /
• pp.319-326
• /
• 2002

A p-version finite element model based on degenerate shell element is proposed for the analysis of orthotropic laminated plates. In the nonlinear formulation of the model, the total Lagrangian formulation is adopted with large deflection and moderate rotation being accounted for in the sense of von Karman hypothesis. The material model Is based on the Huber-Mises yield criterion and Prandtl-Reuss flow rule in accordance with the theory of strain hardening yield function, which is generalized for anisotropic materials by introducing the parameters of anisotropy. The model is also based on extension of equivalent-single layer laminate theory(ESL theory) with shear deformation, leading to continuous shear strain at the interface of two layers. The Integrals of Legendre Polynomials we used for shape functions with p-level varying from 1 to 10. Gauss-Lobatto numerical quadrature is used to calculate the stresses at the nodal points instead of Gauss points. The validity of the proposed p-version finite element model is demonstrated through several comparative points of view in terms of ultimate load, convergence characteristics, nonlinear effect, and shape of plastic zone

• Structural Engineering and Mechanics
• /
• v.4 no.1
• /
• pp.25-35
• /
• 1996

This paper describes a continuum variational formulation for design optimization of nonlinear structures in the elastic-plastic domain, where unloading and reloading of the structures are allowed to occur. The Total Lagrangian procedure is used for the description of the structural deformation. The direct differentiation approach is used to derive the sensitivities of the various structural response measures with respect to the design parameters. Since the material goes into the inelastic range and unloading and reloading of the structure are allowed to occur, the structural response is path dependent and an additional step is needed to integrate the constitutive equations. It can be shown, consequently, that design sensitivity analysis is also path-dependent. The theory has been discretized by the finite element technique and implemented in a structural analysis code. Mathematical programming approach is used for the optimization process. Numerical applications on trusses are performed, where cross-sectional areas and nodal point coordinates are treated as design variables. Optimal designs have been obtained and compared by using two different strategies: a two level strategy where the levels are defined accordingly the type of design variables, cross sectional areas or node coordinates, and optimizing simultaneously with respect to both types of design variables.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.30 no.12 s.255
• /
• pp.1518-1525
• /
• 2006

Finite element analyses of structures made of the fiber reinforced composites require an adequate method to characterize the high anisotropic behavior induced by one or several layers of fiber cords with different spatial orientation embedded in a rubber matrix. This paper newly proposes a continuum based rebar element considering change of the orientation of the fiber during deformation of the composite. The mechanical behavior of the embedded fiber is modeled using two-node bar elements in order to consider the relative deformation and spatial orientation of the embedded fiber. For improvement of the analysis accuracy, the load-displacement curve of fiber is applied to the stiffness matrix of fiber. A finite element program is constructed based on the total Lagrangian formulation considering both geometric and material nonlinearity. Finite element analyses of the tensile test are carried out in order to evaluate the validity of the proposed method. Analysis results obtained with the proposed method provides realistic representation of the fiber reinforced rubber composite compared to results of other two models by the Halpin-Tsai equation and a rebar element in ABAQUS/Standard.