• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.3
• /
• pp.479-486
• /
• 2002

In this study, the progressive inelastic deformation, so called, thermal ratchet phenomenon which can occur in high temperature structures of liquid metal reactor was simulated with thermal ratchet structural test facility and 316L stainless steel test cylinder. The thermal ratchet deformation at the reactor baffle cylinder of the liquid metal reactor can occur due to the moving temperature distribution along the axial direction as the sodium free surface moves up and down under the cyclic heat-up and cool-down transients. The ratchet deformation was measured with the laser displacement sensor and LVDTs after cooling the structural specimen which is heated up to 55$0^{\circ}C$ with steep temperature gradients along the axial direction. The temperature distribution of the test cylinder along the axial direction was measured with 28 channels of thermocouples and was used for the ratchet analysis. The thermal ratchet deformation was analyzed with the constitutive equation of nonlinear combined hardening model which was implemented as ABAQUS user subroutine and the analysis results were compared with those of the test. Thermal ratchet load was applied 9 times and the residual displacement after 9 cycles of thermal load was measured to be 1.79mm. The ratcheting deformation shapes obtained by the analysis with the combined hardening model were in reasonable agreement with those of the structural tests.

• Geomechanics and Engineering
• /
• v.5 no.5
• /
• pp.463-483
• /
• 2013

Considering the influences of fluid penetration, casing, excavation processes of wellbore and perforation tunnels, the seepage-deformation finite element model of oil and gas well coupled with perforating technique is established using the tensile strength failure criterion, in which the user-defined subroutine is developed to investigate the dynamic evolvement of the reservoir porosity and permeability. The results show that the increases of perforation angle and decreases of perforation density lead to a higher fracture initiation pressure, while the changes of the perforation diameter and length have no evident influences on the fracture initiation pressure. As for initiation location for the fracture in wellbore, it is on the wellbore face while considering the presence of the casing. By contrast, the fractures firstly initiate on the root of the tunnels without considering casing. Besides, the initial fracture position is also related with the perforation angle. The fracture initiation position is located in the point far away from the wellbore face, when the perforation angle is around $30^{\circ}$; however, when the perforation angle is increased to $45^{\circ}$, a plane fracture is initiated from the wellbore face in the maximum horizontal stress direction; no fractures was found around perforation tunnels, when the angel is close to $90^{\circ}$. The results have been successfully applied in an oilfield, with the error of only 1.1% comparing the fracture initiation pressure simulated with the one from on-site experiment.

• Structural Engineering and Mechanics
• /
• v.67 no.6
• /
• pp.643-658
• /
• 2018

In this paper, we introduce a new framework for running the finite element (FE) packages inside an online Loop together with MATLAB. Contrary to the Hardware-in-the-Loop techniques (HiL), in the proposed Software-in-the-Loop framework (SiL), the FE package represents a simulation platform replicating the real system which can be out of access due to several strategic reasons, e.g., costs and accessibility. Practically, SiL for sophisticated structural design and multi-physical simulations provides a platform for preliminary tests before prototyping and mass production. This feature may reduce the new product's costs significantly and may add several flexibilities in implementing different instruments with the goal of shortlisting the most cost-effective ones before moving to real-time experiments for the civil and mechanical systems. The proposed SiL interconnection is not limited to ABAQUS as long as the host FE package is capable of executing user-defined commands in FORTRAN language. The focal point of this research is on using the compiled FORTRAN subroutine as a messenger between ABAQUS/CAE kernel and MATLAB Engine. In order to show the generality of the proposed scheme, the limitations of the available SiL schemes in the literature are addressed in this paper. Additionally, all technical details for establishing the connection between FEM and MATLAB are provided for the interested reader. Finally, two numerical sub-problems are defined for offline and online post-processing, i.e., offline optimization and closed-loop system performance analysis in control theory.

• 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.

• Steel and Composite Structures
• /
• v.35 no.6
• /
• pp.739-752
• /
• 2020

This paper aims to investigate cyclic plasticity of a new type of high-performance austenitic stainless steel with both high strength and high ductility. The new stainless steel termed as QN1803 has high nitrogen and low nickel, which leads to reduction of cost ranging from 15% to 20%. Another virtue of the new material is its high initial yield strength and tensile strength. Its initial yield strength can be 40% to 50% higher than conventional stainless steel S30408. Elongation of QN1803 can also achieve approximately 50%, which is equivalent to the conventional one. QN1803 also has a corrosion resistance as good as that of S30408. In this paper, both experimental and numerical studies on the new material were conducted. Full-range true stress-true strain relationships under both monotonic and cyclic loading were obtained. A cyclic plasticity model based on the Chaboche model was developed, where a memory surface was newly added and the isotropic hardening rule was modified. A user-defined material subroutine was written, and the proposed cyclic plasticity model can well evaluate full-range hysteretic properties of the material under various loading histories.

• Composites Research
• /
• v.25 no.5
• /
• pp.136-140
• /
• 2012

This paper describes the bone healing process of fractured long bones such as a tibia applied by composite IM rods using finite element analysis. To simulated tissue differentiation process mechano-regulation theory with a deviatoric strain was implemented and a user's subroutine programmed by a Python code for an iterative calculation was used. To broadly find the appropriate rod modulus for healing bone fractures, composite IM rods were analyzed considering the stacking sequence. To compare mechanical stimulation at fracture gap, two kinds of initial loading conditions were applied. As a result, it was found that the initial loading condition was the most sensitive factor for the healing performance. In case a composite IM rod made of a plain weave carbon fiber/epoxy (WSN3k) had a stacking sequence of $[{\pm}45]_{nT}$, the healing efficiency was the most effective under a initial load of 10%BW.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.29 no.5
• /
• pp.437-445
• /
• 2016

Recently, polyurethane foam has been used in various industry fields to preserve temperature environment of structures, and a wide range of loads from the static to the dynamic are imposed on the material during a life period. The biggest characteristic of polyurethane foam is porosity as being polymeric material, and it is generally known that insulation performance of the material strongly depends on internal void size. In addition, polyurethane foam's mechanical behavior has high dependence on strain rate and temperature as well as being highly non-linear ductile for compression. In the non-linear compressive behavior, volume fraction of voids and elastic modulus decrease as strain increases. Therefore, in this study, temperature-dependent viscoplastic-damage constitutive model was developed to describe the non-linear compressive behavior with the aforementioned features of polyurethane foam.

• Composites Research
• /
• v.26 no.2
• /
• pp.91-98
• /
• 2013

A three-dimensional finite analysis method was proposed to predict the failure of composite double-lap bolted joints, which is based on the stiffness degradation method using damage variables and Hashin's three-dimensional failure criteria. Ladeveze's theory using damage variables to consider the matrix/shear damage was combined with stiffness degradation in fiber direction. Four different failure modes were considered including matrix compression/shear, matrix tension/shear, fiber compression, and tension failures. The friction between bolt and composite and the clamping force were considered using a commercial finite element software ABAQUS. The damage model was incorporated using the user-defined subroutine of the software. The predicted result was verified with the existing test result for bearing tension double shear and showed the deviation ranging 7~16% from test results.

• Journal of Ocean Engineering and Technology
• /
• v.30 no.6
• /
• pp.474-483
• /
• 2016

Accidental events such as collisions, groundings, and hydrocarbon explosions in marine structures can cause catastrophic damage. Thus, it is extremely important to predict the extent of such damage, which determines the total amount of oil spills and the residual hull girder strength. Punching fracture tests were conducted by Choung (2009b), where various sizes of indenters and circular unstiffened steel plates with different thicknesses were used to quasi-statically realize damage extents. A three-dimensional fracture strain surface was developed based on a reference (Choung et al., 2015b), where the average stress triaxiality and average normalized Lode angle were used as the parameters governing the fracture of ductile steels. In this study, new numerical analyses were performed using very fine axisymmetric elements in combination with an Abaqus user-subroutine to implement the three-dimensional fracture strain surface. Conventional numerical analyses were also conducted for the tests to identify the best fit fracture strain values by changing the fracture strains. Based on the phenomenon of the average normalized Lode angle starting out positive and then becoming slightly negative, it was inferred that the shear stress primarily dominates in determining the fractures locations, with a partial contribution from the compressive stress. It should be stated that the three-dimensional fracture surface effectively predicted at least the shear stress-dominant fracture behavior of a mild steel.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.28 no.5
• /
• pp.543-551
• /
• 2015

Polyurethane Foam(PUF), a outstanding thermal insulation material, is used for various structures as being composed with other materials. These days, PUF composed with glass fiber, Reinforced PUF(R-PUF), is used for a insulation system of LNG Carrier and performs function of not only the thermal insulation but also a structural member for compressive loads like a sloshing load. As PUF is a porous material made by mixing and foaming, mechanical properties depend on volume fraction of voids which is a dominant parameter on density. Thus, In this study, density is considered as the effect parameter on mechanical properties of Polyurethane Foam, and mechanical behavior for compression of the material is described by using modified Gurson damage model.