• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.1
• /
• pp.145-153
• /
• 1991

A numerical simulation of the nonsteady state rolling process in the plane strain condition is presented in the basis of the rigid-plastic finite element method by considering large deformation. In order to apply the large deformation theory to the numerical method for sheet rolling problems, constitutive equation relating 2nd-Piola Kirchhoff stress and Lagrangian strain which reflect geometrical nonlinearity is used. To confirm the validity of the developed algorithm, the analysis of the neutral flow region, roll separating force, torque, pressure and stress/strain distributions on the workpiece is conducted from the bite of the material until the steady state is reached. The computed results of the roll force and torque in the present finite element analysis are lower than those corresponding to small strain theory. The pressure distribution at the work piece-roll interface is found to show the typical 'friction hill' type only. The peak value in near the neutral region, however, is good agrements with the existing results. the neutral region, however, is good agrements with the existing results. The frictional force at the roll interface provide detailed information about the neutral point where the shear forces change direction. In addition, the analysis also includes the effect and influence of material condition, strip thickness, work roll diameter, as well as roll speed and lubricant on each deformation process.

• Korea-Australia Rheology Journal
• /
• v.19 no.4
• /
• pp.233-242
• /
• 2007

We performed a numerical investigation to find out the optimal choice of the spatial discretization in the distributed-Lagrangian-multiplier/fictitious-domain (DLM/FD) method for the solid/fluid interaction problem. The elastic solid bar attached on the bottom in a pressure-driven channel flow of a Newtonian fluid was selected as a model problem. Our formulation is based on the scheme of Yu (2005) for the interaction between flexible bodies and fluid. A fixed regular rectangular discretization was applied for the description of solid and fluid domain by using the fictitious domain concept. The hydrodynamic interaction between solid and fluid was treated implicitly by the distributed Lagrangian multiplier method. Considering a simplified problem of the Stokes flow and the linearized elasticity, two numerical factors were investigated to clarify their effects and to find the optimum condition: the distribution of Lagrangian multipliers and the solid/fluid interfacial condition. The robustness of this method was verified through the mesh convergence and a pseudo-time step test. We found that the fluid stress in a fictitious solid domain can be neglected and that the Lagrangian multipliers are better to be applied on the entire solid domain. These results will be used to extend our study to systems of elastic particle in the Stokes flow, and of particles in the viscoelastic fluid.

• Tunnel and Underground Space
• /
• v.9 no.4
• /
• pp.296-305
• /
• 1999

In ordinary back analysis it if hardly possible to obtain the mechanical properties of tunnel lining by using commonly measured displacements of tunnel lining, because only a few displacements could be measured at the site. Therefore, it is necessary to develop a new method which can evaluate the state of stresses of tunnel by using measured data. In this study, in order to assess tunnel lining stability by estimating its stresses with a few measured displacements, a formulation of back analysis method was proposed. The accuracy of results were investigated through the parametric study for several types of measurement model of two dimensional elastic lining. This new back analysis method to assess tunnel lining stresses and strains with a few numbers of measured displacements showed high accuracy and good applicability when compared to the results of numerical experiments by FEM. The method has been tested on subway tunnel and its applicability has been confirmed by comparing field and analytical data. It is verified that the stress on the tunnel lining can be obtained by only more than 3 point of input displacements without any condition of external loads.

• Coupled systems mechanics
• /
• v.1 no.4
• /
• pp.361-379
• /
• 2012

This paper presents a fully coupled three-dimensional solver for the analysis of interaction between pulsatile flow and large deformation structure. A partitioned time marching algorithm is employed for the solution of the time dependent coupled discretised problem, enabling the use of highly developed, robust and well-tested solvers for each field. Conservative transfer of information at the fluid-structure interface is combined with an effective multi-predict-correct iterative scheme to enable implicit coupling of the interacting fields at each time increment. The three-dimensional unsteady incompressible fluid is solved using a powerful implicit time stepping technique and an ALE formulation for moving boundaries with second-order time accurate is used. A full spectrum of total variational diminishing (TVD) schemes in unstructured grids is allowed implementation for the advection terms and finite element shape functions are used to evaluate the solution and its variation within mesh elements. A finite element dynamic analysis of the highly deformable structure is carried out with a numerical strategy combining the implicit Newmark time integration algorithm with a Newton-Raphson second-order optimisation method. The proposed model is used to predict the wave flow fields of a particular flow-induced vibrational phenomenon, and comparison of the numerical results with available experimental data validates the methodology and assesses its accuracy. Another test case about three-dimensional biomedical model with pulsatile inflow is presented to benchmark the algorithm and to demonstrate the potential applications of this method.

• Proceedings of the Korean Powder Metallurgy Institute Conference
• /
• 1997.10a
• /
• pp.7-7
• /
• 1997

During sintering of very porous green bodies, as obtained by compaction of hard powders - such as tungsten carbide or ceramics - or by injection moulding, important shrinkage occurs. Due to heterogeneous green density field, gravity effects, friction on the support, thermal gradients, etc., this shrinkage is often non-uniform, which' may induce significant shape changes. As the ratio of compact dimension to powder size is very high, the mechanics of continuum is relevant to model such phenomena. Thus numerical techniques, such as the finite element method can be used to simulate the sintering process and predict the final shape of the sintered part. Such type of simulation has much been developed in the last decade firstly for hot isostatic pressing and next for die compaction. Finite element modelling has been recently applied to free sintering. The simulation of sintering should be based on constitutive equations describing the thermo-mechanical behaviour of the material under any state of stress and any temperature which may arise within the sintering body. These equations can be drawn either from experimental data or from micromechanical models. The experiments usually consist in free sintering and sinter-forging tests. Indeed applying more complex loading conditions at high temperature under controlled atmosphere is delicate. Micromechanical models describe the constitutive behaviour of aggregates of spheres from the deformation of two-sphere contact either by viscous flow or grain boundary diffusion. Such models are not able to describe complex microstructure and mechanisms as observed in real materials but they can give some basic information on the formulation of constitutive equations. Practically both experimental and theoretical approaches can be coupled to identify the constitutive equations. Such procedure has been performed for modelling the sintering of compacts obtained by die pressing of a mixture of tungsten carbide and cobalt powders. The constitutive behaviour of this material during sintering has been described by a linear viscous constitutive model, whose functions have been fitted from results of free sintering and sinter-forging experiments. This model has next been introduced in ABAQUS finite element code to simulate the sintering of heterogeneous green compacts of various geometries at constant temperature. Examples of simulations are shown and compared with experiments.

• Structural Engineering and Mechanics
• /
• v.68 no.3
• /
• pp.299-312
• /
• 2018

In this paper, an extended finite element method is proposed to analyze both geometric and material non-linear behavior of general Functionally Graded Material (FGM) plate-shell type structures. A user defined subroutine (UMAT) is developed and implemented in Abaqus/Standard to study the elastoplastic behavior of the ceramic particle-reinforced metal-matrix FGM plates-shells. The standard quadrilateral 4-nodes shell element with three rotational and three translational degrees of freedom per node, S4, is extended in the present study, to deal with elasto-plastic analysis of geometrically non-linear FGM plate-shell structures. The elastoplastic material properties are assumed to vary smoothly through the thickness of the plate-shell type structures. The nonlinear approach is based on Mori-Tanaka model to underline micromechanics and locally determine the effective FGM properties and self-consistent method of Suquet for the homogenization of the stress-field. The elasto-plastic behavior of the ceramic/metal FGM is assumed to follow Ludwik hardening law. An incremental formulation of the elasto-plastic constitutive relation is developed to predict the tangent operator. In order to to highlight the effectiveness and the accuracy of the present finite element procedure, numerical examples of geometrically non-linear elastoplastic functionally graded plates and shells are presented. The effects of the geometrical parameters and the volume fraction index on nonlinear responses are performed.

• Journal of the Society of Naval Architects of Korea
• /
• v.33 no.2
• /
• pp.1-12
• /
• 1996

Hydrodynamic characteristics of shipping water on deck are investigated by using a simplified two-dimensional model. Formulation of the shipping water on deck leads to a nonlinear hyperbolic system of equations based on the shallow-water wave theory. Time-domain solution of these equations are obtained numerically using a finite difference method which adopts predictor-corrector method for time-marching and 2nd upwind differencing method for convection term calculation. To confirm the validity of the present numerical method, we calculated some shallow-water wave problems accompanying a bore and compared the obtained results with the analytic solutions. We found good agreements between them. Though the calculation results of shipping water on deck, we show that the shipping water flows into the deck as a rarefying wave arid grows into a bore after colliding with a deck structure. Also we examined the effects of acceleration and slope of deck and found that they have significant influences on the flow of shipping water.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.18 no.4 s.70
• /
• pp.395-404
• /
• 2005

This paper presents a reliable numerical procedure for nonlinear time-history analysis of space steel frames subjected to dynamic loads. Geometric nonlinearities of member (P-$\delta$) and frame (P-$\Delta$) are taken into account by the use of stability functions in framed stiffness matrix formulation. The gradual yielding along the member length and over the cross section is included by using a tangent modulus concept and a softening plastic hinge model based on the New-Orbison yield surface. A computer program utilizing the average acceleration method for the integration scheme is developed to numerically solve the equation of motion of framed structure formulated in an incremental form. The results of several numerical examples are compared with those derived from using beam element model of ABAQUS program to illustrate the accuracy and the computational efficiency of the proposed procedure.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2003.11a
• /
• pp.244-251
• /
• 2003

A two-dimensional BEM code, $FRACOD^{2D}$, was applied to simulate fracture initiation and propagation processes in a rock pillar during an in situ heater test of a rock pillar planned at the $\"{A}sp\"{o}$ Underground Rock laboratory of SKB, in Southern Sweden. To take the advantage of conventional BEM for simulating fracturing processes, but without efforts for domain integral transformation, a hybrid approach is developed to simulate the fracturing processes in rock pillar under coupled thermo-mechanical loading. The code FRACOD was used for simulating the fracture initiation and propagation processes with its boundary tractions reflecting the effects of the initial and redistributed thermomechanical stresses in the domain of interest at multiple excavation and heating steps were produced by a special algorithm of stress inversion, based on resultant thermo-mechanical stress fields at each excavation and heat loading step by a FEM code without considering fracturing processes. This hybrid approach can take the advantages of both types of numerical methods and avoids their shortcomings for fracturing process simulation and domain effects, respectively. In this paper, we present the hybrid approach for the stress, displacements, and fracturing processes at sequential excavation and heating steps of the in situ heater test as a predictive modelling, the formulation of the fracturing models and the predictive results. Two sections of borehole depth, 0.5 m and 1.5 m below the tunnel floor are considered. The pillar area is modelled with the FRACOD and the stress field produced by excavation and heating is transferred with corresponding boundary stresses. From the modelling results, the degree of fracturing and damage are evaluated for 120 days of heating. Dominated shear fracturing in the vicinity of the central pillar was observed from the models at both sections, but spalled area appears to be limited. Based on the modelling results, a sensitivity study for the effect of pre-existing fractures in the vicinity of the holes is also conducted, and the initiation and evolution of EDZ around the deposition holes are investigated using this particular numerical technique.

• Journal of computational fluids engineering
• /
• v.18 no.4
• /
• pp.25-34
• /
• 2013

A numerical study of a laminar natural convection of the CuO-water nanofluid in a square cavity using the Buongiorno's nonhomogeneous model is presented. All the governing equations including the volume fraction equation are discretized on a cell-centered, non-uniform grid employing the finite-volume method with a primitive variable formulation. Calculations are performed over a range of Rayleigh numbers and volume fractions of the nanopartile. From the computed results, it is shown that both the homogeneous and nonhomogeneous models predict the deterioration of the natural convection heat transfer well with an increase of the volume fraction of nanoparticle at the same Rayleigh number, which was observed in the previous experimental studies. It is also shown that the differences in the computed results of the average Nusselt number at the wall between the homogeneous and nonhomogeneous models are very small, and this indicates that the slip mechanism of the Brown diffusion and thermophoresis effects are negligible in the laminar natural convection of the nanofluid. The degradation of the heat transfer with an increase of the volume fraction of the nanoparticle in the natural convection of nanofluid is due to the increase of the viscosity and the decrease of the thermal expansion coefficient and the specific heat. It is clarified in the present study that the previous controversies between the numerical and experimental studies are owing to the different definitions of the Nusselt number.