• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.4
• /
• pp.191-197
• /
• 2018

This paper details the development of an isotropic elastic-damage constitutive model for polymeric foam based on irreversible thermodynamics to consider the growth and coalescence of voids. The constitutive equations describe the material behavior sustaining unilateral damage. To facilitate finite element analysis, the material properties for specific types of polymeric foams are applied to the developed model; the model is then implemented in ABAQUS as a user-defined material subroutine. To validate the developed damage model, the simulated results are compared to the results of a series of tensile tests on various polymeric foams. The proposed damage model can be utilized to further research on continuum damage mechanics and finite element analysis of polymeric foams in computational engineering.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.30 no.3A
• /
• pp.309-316
• /
• 2010

This paper proposed a nonlocal anisotropic damage model to simulate the behavior of plain and reinforced concrete structures that are predominantly tensile and compressive load. This model based on continuum damage mechanics, used a symmetric second-order tensor as the damage variable. For quasi-brittle materials, such as concrete, the damage patterns were different in tension and in compression. These two damage states were modeled by damage evolution laws ensuring a damage tensor rate proportional to the total strain tensor in terms of principal components. To investigate the effectiveness of proposed model, the double edge notched specimen experimented by nooru-mohamed and reinforced concrete bending beam were analyzed using the implementation of the proposed model. As the results for the simulation, the nonlocal anisotropic damage model with an adequate control of rupture correctly represented the crack propagation for mixed mode fracture. In the structural failure of reinforced concrete bending beam, the proposed model can be showed up to a very high damage level and yielding of the reinforcements.

• Interaction and multiscale mechanics
• /
• v.5 no.3
• /
• pp.229-265
• /
• 2012

Hygroexpansion of wood is a known and undesired characteristic in civil engineering. When wood is exposed to changing environmental humidity, it adsorbs or desorbs moisture and warps. The resulting distortions or - at restrained conditions - cracks are a major concern in timber engineering. We herein present a multiscale model for prediction of the macroscopic hygroexpansion behavior of individual pieces of softwood from their microstructure, demonstrated for spruce. By applying poromicromechanics, we establish a link between the swelling pressure, driving the hygroexpansion of wood at the nanoscale, and the resulting macroscopic dimensional changes. The model comprises six homogenization steps, which are performed by means of continuum micromechanics, the unit cell method and laminate theory, all formulated in a poromechanical framework. Model predictions for elastic properties of wood as functions of the moisture content closely approach corresponding experimental data. As for the hygroexpansion behavior, the swelling pressure has to be back-calculated from macroscopic hygroexpansion data. The good reproduction of the anisotropy of wood hygroexpansion, based on only a single scalar calibration parameter, underlines the suitability of the model. The multiscale model constitutes a valuable tool for studying the effect of microstructural features on the macroscopic behavior and for assessing the hygroexpansion behavior at smaller length scales, which are inaccessible to experiments. The model predictions deliver input parameters for the analysis of timber at the structural scale, therewith enabling to optimize the use of timber and to prevent moisture-induced damage or failure.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2006.05a
• /
• pp.349-352
• /
• 2006

The plastic deformation of polycrystalline materials is induced by changes of the microstructure when the loading is beyond the critical state of stress. Constitutive models for the crystal plasticity have the common objective which relates microscopic single crystals in the crystallographic texture to the macroscopic continuum point. In this paper, a new consistent tangent stiffness for the anisotropic elasto-viscoplastic analysis of polycrystalline deformation is developed, which can be used in the finite element analysis for the slip-dominated large deformation of polycrystalline materials. In order to calculate the consistent tangent stiffness, the state function is defined based on the consistency condition between the elastic and plastic stress. The rate of shearing increment($\Delta{\gamma}^{\alpha}$) is calculated with satisfying the consistency condition. The consistency condition becomes zero when the trial resolved shear stress($\tau^{{\alpha}^*}$) becomes resolved shear stress($\tau^{\alpha}$) at every step. Iterative method is utilized to calculate the rate of shearing increment based on the implicit backward Euler method. The consistent tangent stiffness can be formulated by differentiating the rate of shearing increment with total strain increment after the instant rate of shearing increment converges. The proposed tangent stiffness is applied to the ABAQUS/Standard by implementing in the ABAQUS/UMAT.

• International Journal of Precision Engineering and Manufacturing
• /
• v.5 no.4
• /
• pp.50-60
• /
• 2004

A very flexible beam can be used to model various types of continuous mechanical parts such as cables and wires. In this paper, the dynamic properties of a very flexible beam, included in a multibody system, are analyzed using absolute nodal coordinates formulation, which is based on finite element procedures, and the general continuum mechanics theory to represent the elastic forces. In order to consider the dynamic interaction between a continuous large deformable beam and a rigid multibody system, a combined system equations of motion is derived by adopting absolute nodal coordinates and rigid body coordinates. Using the derived system equation, a computation method for the dynamic stress during flexible multibody simulation is presented based on Euler-Bernoulli beam theory, and its reliability is verified by a commercial program NASTRAN. This method is significant in that the structural and multibody dynamics models can be unified into one numerical system. In addition, to analyze a multibody system including a very flexible beam, formulations for the sliding joint between a very deformable beam and a rigid body are derived using a non-generalized coordinate, which has no inertia or forces associated with it. In particular, a very flexible catenary cable on which a multibody system moves along its length is presented as a numerical example.

• Structural Engineering and Mechanics
• /
• v.45 no.6
• /
• pp.759-777
• /
• 2013

In this paper, a topology optimization method based on the element independent nodal density (EIND) is developed for continuum solids with multiple load cases and multiple constraints. The optimization problem is formulated ad minimizing the volume subject to displacement constraints. Nodal densities of the finite element mesh are used a the design variable. The nodal densities are interpolated into any point in the design domain by the Shepard interpolation scheme and the Heaviside function. Without using additional constraints (such ad the filtering technique), mesh-independent, checkerboard-free, distinct optimal topology can be obtained. Adopting the rational approximation for material properties (RAMP), the topology optimization procedure is implemented using a solid isotropic material with penalization (SIMP) method and a dual programming optimization algorithm. The computational efficiency is greatly improved by multithread parallel computing with OpenMP to run parallel programs for the shared-memory model of parallel computation. Finally, several examples are presented to demonstrate the effectiveness of the developed techniques.

• Tunnel and Underground Space
• /
• v.31 no.4
• /
• pp.270-288
• /
• 2021

We proposed a numerical method to simulate the hydro-mechanical behavior of rock fracture using a grain-based distinct element model (GBDEM) in the paper. As a part of DECOVALEX-2023 Task G, we verified the method via benchmarks with analytical solutions. DECOVALEX-2023 Task G aims to develop a numerical method to estimate the coupled thermo-hydro-mechanical processes within the crystalline rock fracture network. We represented the rock sample as a group of tetrahedral grains and calculated the interaction of the grains and their interfaces using 3DEC. The micro-parameters of the grains and interfaces were determined by a new methodology based on an equivalent continuum approach. In benchmark modeling, a single fracture embedded in the rock was examined for the effects of fracture inclination and roughness, the boundary stress condition and the applied pressure. The simulation results showed that the developed numerical model reasonably reproduced the fracture slip induced by boundary stress condition, the fracture opening induced by fluid injection, the stress distribution variation with fracture inclination, and the fracture roughness effect. In addition, the fracture displacements associated with the opening and slip showed good agreement with the analytical solutions. We expect the numerical model to be enhanced by continuing collaboration and interaction with other research teams of DECOVALEX-2023 Task G and validated in further study experiments.

• Tunnel and Underground Space
• /
• v.31 no.6
• /
• pp.593-609
• /
• 2021

We proposed a numerical method for the thermo-mechanical behavior of rock fracture using a grain-based distinct element model (GBDEM) and simulated thermally induced fracture slip. The present study is the benchmark simulation performed as part of DECOVALEX-2023 Task G, which aims to develop a numerical method to estimate the coupled thermo-hydro-mechanical processes within the crystalline rock fracture network. We represented the rock sample as an assembly of Voronoi grains and calculated the interaction of the grains (blocks) and their interfaces (contacts) using a distinct element code, 3DEC. Based on an equivalent continuum approach, the micro-parameters of grains and contacts were determined to reproduce rock as an elastic material. Then, the behavior of the fracture embedded in the rock was characterized by the contacts with Coulomb shear strength and tensile strength. In the benchmark simulation, we quantitatively examined the effects of the boundary stress and thermal stress due to heat conduction on fracture behavior, focusing on the mechanism of thermally induced fracture slip. The simulation results showed that the developed numerical model reasonably reproduced the thermal expansion and thermal stress increment, the fracture stress and displacement and the effect of boundary condition. We expect the numerical model to be enhanced by continuing collaboration and interaction with other research teams of DECOVALEX-2023 Task G and validated in further study experiments.

• Journal of the Korean Magnetics Society
• /
• v.20 no.3
• /
• pp.100-105
• /
• 2010

For predicting magnetic signals due to the remanent magnetization distributed on a ferromagnetic ship hull, this paper presents an efficient methodology for solving inverse problems, where the material sensitivity analysis based on the continuum mechanics is combined with the equivalent magnetic models. To achieve this, the 3D magnetic charge model and the magnetic dipole moment model are introduced and material sensitivity formulae applicable to each equivalent model are derived. The formulae offer the first-order gradient information of an objective function with respect to the variation of the magnetic charge or magnetic dipole and so an optimal solution can be easily obtained regardless of the number of design variables. To validate the proposed method, the numerical results are comparison with the real measurements of a mock-up model.

• Interaction and multiscale mechanics
• /
• v.3 no.1
• /
• pp.95-110
• /
• 2010

In recent years, study of the static response of pavements to moving vehicle and aircraft loads has received significant attention because of its relevance to the design of pavements and airport runways. The static response of beams resting on an elastic foundation and subjected to moving loads was studied by several researchers in the past. However, most of these studies were limited to steady-state analytical solutions for infinitely long beams resting on Winkler-type elastic foundations. Although the modelling of subgrade as a continuum is more accurate, such an approach can hardly be incorporated in analysis due to its complexity. In contrast, the two-parameter foundation model provides a better way for simulating the underlying soil medium and is conceptually more appealing than the one-parameter (Winkler) foundation model. The finite element method is one of the most suitable mathematical tools for analysing rigid pavements under moving loads. This paper presents an improved solution algorithm based on the finite element method for the static analysis of rigid pavements under moving vehicular or aircraft loads. The concrete pavement is discretized by finite and infinite beam elements, with the latter for modelling the infinity boundary conditions. The underlying soil medium is modelled by the Pasternak model allowing the shear interaction to exist between the spring elements. This can be accomplished by connecting the spring elements to a layer of incompressible vertical elements that can deform in transverse shear only. The deformations and forces maintaining equilibrium in the shear layer are considered by assuming the shear layer to be isotropic. A parametric study is conducted to investigate the effect of the position of moving loads on the response of pavement.