• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.33 no.5
• /
• pp.447-454
• /
• 2009

Rate-independent crystal plasticity model suffers from the non-uniqueness of activated slip systems and the determination of the shear slip rates on the active slip systems. In this paper, a time-integration algorithm which circumvents the problem of the multiplicity of the slip systems was developed and implemented into the user subroutine VUMAT of a commercial finite element program ABAQUS. The magnitude of the slip shears on the active slip systems in f.c.c Cu single crystal aligned with the specific crystallographic orientation was investigated to validate our solution procedure. Also, texture developments under various deformation modes such as simple compression, simple tension and plane strain compression were compared with the results of the rate-dependent model by using the rate-independent crystal plasticity model. The computation time employing the rate-independent model is much more reduced than the those of the rate-dependent model.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2007.10a
• /
• pp.67-69
• /
• 2007

Unusual mechanical constitutive behavior of magnesium alloy sheets has been implemented into the finite element program ABAQUS via user material subroutine. For the verification purpose, the springback of AZ31B magnesium alloy sheet was measured using the unconstrained cylindrical bending test of Numisheet'2002. In addition to the developed constitutive models, the other two models based on isotropic constitutive equations with tensile and compressive properties were also considered. Preliminary comparisons have been made between simulated results by the finite element analysis and corresponding experiments and the newly proposed model showed enhanced prediction capability in springback prediction.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.29 no.2 s.233
• /
• pp.253-261
• /
• 2005

In decades, a substantial body of work on a unified viscoplastic model which considers the mechanism of plastic deformation and creep deformation has developed. The systematic scheme for numerical analysis of unified model is necessary because the dominant failure mechanism is the defect growth and coalescence in materials. In the present study, the unified viscoplastic model for materials with defects suggested by Suquet and Michel was employed for numerical analysis. The constitutive equations are integrated based on the generalized mid-point rule and implemented into a finite element program (ABAQUS) by means of user-defined subroutine (UMAT). To evaluate the validity of the developed UMAT code and the assessment of the adopted viscoplastic model, the results obtained from the UMAT code was compared with the numerical reference solution and experimental data. The unit cell analysis also has been investigated to study the effect of strain rate, temperature, stress triaxiality and initial defect volume fraction on the growth and coalescence of the defect.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2006.05a
• /
• pp.151-155
• /
• 2006

Elastic-viscoplastic structural analysis is performed for regenerative cooling chamber of liquid rocket thrust chamber using Bodner-Partom visco-plastic model. Strain rate test was also conducted for a copper alloy at various temperatures in order to get material constants of visco-plastic model used in the structural analysis. Material constants of visco-plastic model were obtained from strain rate test results and visco-plsstic model was incorporated into finite element program, Marc, by means of user subroutine. The structural analysis results indicate that the deformation of cooling channel is mostly caused by thermal loading rather than pressure loading and confirmed structural stability of the cooling channel under operating condition.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.10 no.4
• /
• pp.69-76
• /
• 2006

Elastic-viscoplastic structural analysis has been performed for regenerative cooling chamber of liquid rocket thrust chamber using Bodner-Partom visco-plastic model. Strain rate test was conducted for a copper alloy at various temperatures in order to get material constants of visco-plastic model used in the structural analysis. Material constants of visco-plastic model were obtained from strain rate test results and visco-plastic model was incorporated into finite element program, Marc, by means of a user subroutine. The structural analysis results indicated that the deformation of cooling channel is mostly caused by thermal loading rather than pressure loading and confirmed structural stability of the cooling channel under the operating condition.

• Journal of Ocean Engineering and Technology
• /
• v.36 no.6
• /
• pp.390-402
• /
• 2022

Offshore wind energy has become a major energy source, and various studies are underway to increase the economic feasibility of floating offshore wind turbines (FOWT). In this study, the characteristics of wave-induced motion of a combined wind-wave energy platform were analyzed to reduce the variability of energy extraction. A user subroutine was developed, and numerical analysis was performed in connection with the ANSYS-AQWA hydrodynamic program in the time domain. A platform combining the TLP-type FOWT and the Wavestar-type wave energy converter (WEC) was proposed. Each motion response of the platform on the second-order wave load, the effect of WEC attachment and Power take-off (PTO) force were analyzed. The mooring line tension according to the installation location was also analyzed. The vertical motion of a single FOWT was increased approximately three times due to the second-order sum-frequency wave load. The PTO force of the WEC played as a vertical motion damper for the combined platform. The tension of the mooring lines in front of the incident wave direction was dominantly affected by the pitch of the platform, and the mooring lines located at the side of the platform were mainly affected by the heave of the platform.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.3 no.4
• /
• pp.95-107
• /
• 1983

The general purpose programs are in their fixed algorithm and theory of mechanics which can not be altered without painful program modifications. Users are usually guided by user's manual for data input. The several symbolic manipulation programs for structural analysis are introduced recently. These programs allow users to include a wide class of solution algorithm and to specify, by means of some symbolic manipulation, a combination of analytical steps to suit a particular problem. As they can solve a single domain problem, a large computer is usually needed. The scope of this study is to develop an efficient symbolic manipulation program with space beam element, plate bending element and eigen value routines. The incorporated Substructure capability and generation capability of finite element characteristic arrays (e.g., stiffness matrix, mass matrix) enables users to analyse multidomain problem with small computer. The program consists of modulized independent processors, each having its own specific function and is easily modified, eliminated and added. The processors are efficiently handling data by the Data base approach which is the concept of integrated program network(IPN).

• Journal of the Korean Geotechnical Society
• /
• v.20 no.4
• /
• pp.89-102
• /
• 2004

The generated blasting pressure wave initiated under decoupled-charge condition is a function of peak blasting pressure, rise time, and wave-shape function. The peak blasting pressure and the rise time are also the function of explosive and rock properties. The probabilistic distributions of explosive and rock properties are derived from the results of their property tests. Since the probabilistic distributions of explosive and rock properties displayed a normal distribution, the peak blasting pressure and the rise time can also be regarded as a normal distribution. Parameter analysis and uncertainty analysis were performed to identify the most influential parameter that affects the peak blasting pressure and the rise time. Even though the explosive properties were found to be the most influential parameters on the peak blasting pressure and the rise time from the parameter analyses, the result of uncertainty analysis showed that rock properties constituted major uncertainties in estimating the peak blasting pressure and the rise time rather than explosive properties. Damage and overbreak of the remaining rock around the excavation line induced by blasting were evaluated by dynamic numerical analysis. A user-subroutine to estimate the rock damage was coded based on the continuum damage mechanics. This subroutine was linked to a commercial program called 'ABAQUS/Explicit'. The results of dynamic numerical analysis showed that the rock damages generated by the initiation of stopping hole were larger than those from the initiation of contour hole. Several methods to minimize those damages were proposed such as relocation of stopping hole, detailed subdivision of rock classification, and so on. It was found that fracture probability criteria and fractured zones could be distinctively identified by applying fuzzy-random probability.

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