• Composites Research
• /
• v.36 no.1
• /
• pp.42-47
• /
• 2023

In this paper, we present a fatigue life prediction methodology using multiple S-N curves according to the different stress states of laminated composites. The stress states of the plies of the laminated composites are classified into five modes: longitudinal tension or compression and transverse tension or compression, and shear according to the maximum stress criterion and Puck's criterion with a scaling factor K. This methodology has advantages in computational cost, and it can also consider microstructural characteristics of the composites by applying different S-N curves. The S-N curves for the fatigue analysis are obtained by experimental fatigue test. The proposed methodol is implemented into commercial software, ABAQUS user material subroutine and therefore, the fatigue analysis is conducted using the structural analysis results. The finite element (FE) simulation results are presented for unidirectional composites with and without open-hole. The FE simulation results show that the stress condition is different depending on the fiber orientation of the unidirectional composite, so the fatigue life is calculated with different S-N curves.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.48 no.9
• /
• pp.663-670
• /
• 2020

A study on mechanical property characterization and modeling technique was carried out to approximate the behaviour of structures with 2.5D C/SiC material. Several tensile tests were performed to analyze the behaviour characteristics of the 2.5D C/SiC material and elastic property was characterized by applying a mathematical homogenization and a modified rule of mixture. SiC matrix representing the elasto-plastic behavior approximates as a bilinear function. Then the equivalent yield strength and equivalent plastic stiffness were calculated by minimizing errors in experiment and approximation. RVE(Representative Volume Element)was defined from the fiber and matrix configuration of 2.5D C/SiC and a process of calculating the effective stiffness matrix by applying the modified rule of mixture to RVE was implemented in the ABAQUS User-defined subroutine. Finite element analysis was performed by applying the mechanical properties of fiber and matrix calculated based on the proposed process, and the results were in good agreement with the experimental results.

• International Journal of Highway Engineering
• /
• v.11 no.1
• /
• pp.165-175
• /
• 2009

With the current move towards adopting mechanistic-empirical concepts in the design of pavement structures, state-of-the-art mechanistic analysis methodologies are needed to determine accurate pavement responses, such as stress, strain, and deformation. Previous laboratory studies of pavement foundation geomaterials, i.e., unbound granular materials used in base/subbase layers and fine-grained soils of a prepared subgrade, have shown that the resilient responses followed by nonlinear, stress-dependent behavior under repeated wheel loading. This nonlinear behavior is commonly characterized by stress-dependent resilient modulus material models that need to be incorporated into finite element (FE) based mechanistic pavement analysis methods to predict more realistically predict pavement responses for a mechanistic pavement analysis. Developed user material subroutine using aforementioned resilient model with nonlinear solution technique and convergence scheme with proven performance were successfully employed in general-purpose FE program, ABAQUS. This numerical analysis was investigated in predicted critical responses and domain selection with specific mesh generation was implemented to evaluate better prediction of pavement responses. Results obtained from both axisymmetric and three-dimensional (3D) nonlinear FE analyses were compared and remarkable findings were described for nonlinear FE analysis. The UMAT subroutine performance was also validated with the instrumented full scale pavement test section study results from the Federal Aviation Administration's National Airport Pavement Test Facility (FAA's NAPTF).

• Journal of Mechanical Science and Technology
• /
• v.20 no.9
• /
• pp.1399-1409
• /
• 2006

An improved friction model was proposed with consideration of the effect of the sliding speed, the contact pressure and the temperature, and it was implemented into a user subroutine of a commercial FEM code, ABAQUS/Explicit. Then, a smooth tire was simulated for free rolling, driving, braking and cornering situations using the improved friction model and the Coulomb friction model, and the effect of the friction models on the slip, the frictional energy distribution and the cornering force and moment was analyzed. For the free rolling, the driving and the braking situations, the improved friction model and the Coulomb friction model resulted in similar profiles of the slip and the frictional energy distributions although the magnitudes were different. The slips obtained from the simulations were in a good correlation with experimental data. For the cornering situation, the Coulomb friction model with the coefficient of friction of 1 or 2 resulted in lower or higher cornering forces and moments than experimental data. In addition, in contrast to experimental data it did not result in a maximum cornering force and a decrease of the cornering moment for the increase of the speed. However, the improved friction model resulted in similar cornering forces and moments to experimental data, and it resulted in a maximum cornering force and a decrease of the cornering moment for the increase of the speed, showing a good correlation with experimental data.

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

• 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 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.22 no.6
• /
• pp.519-525
• /
• 2009

In this study a cohesive zone model was used to simulate the delamination phenomena which occurs by a successive crack initiation and propagation in composite laminates. The cohesive zone model was incorporated to the classical finite element method via cohesive element formulation and then implemented into the user-subroutine UEL of a commercial finite element program Abaqus. To validate the formulation and implementation of the cohesive element the finite element results were compared with the experimental data of double cantilever beam and end notched flexure tests. The numerical results well agree with the experimental load-displacement curves. Also the effect of the elastic stiffness and the size of the cohesive element on the global load-displacement curves were studied numerically. To minimize the mesh-dependency of the crack propagation path and eliminate the zig-zag patterns in the load-displacement curve, cohesive elements should be refined at the crack-tip.