• 한국전산구조공학회논문집
• /
• 제23권6호
• /
• pp.635-640
• /
• 2010

When using finite element pogram ABAQUS to compute material characteristics, one builds a user material subroutine if unique constitutive feature needs to be included. In ABAQUS/Standard, UMAT subroutine should be built, and in ABAQUS/Explicit, VUMAT should be used. Although two subroutines carry out the same type of task, two different programs should be made depending on the working environment, and it is not easy to program the subroutines following the manual without enough understanding of solid mechanics. In this paper, difference between UMAT and VUMAT subroutines is epitomized, and a conversion scheme from UMAT to VUMAT is discussed. An example shows that the two programs give the same stress computation result.

• Steel and Composite Structures
• /
• 제18권5호
• /
• pp.1291-1303
• /
• 2015

Damage caused by low velocity impact is so dangerous in composites because although in most cases it is not visible to the eye, it can greatly reduce the strength of the composite material. In this paper, damage development in U-section glass/polyester pultruded beams subjected to low velocity impact was considered. Different failure criteria such as Maximum stress, Maximum strain, Hou, Hashin and the combination of Maximum strain criteria for fiber failure and Hou criteria for matrix failure were programmed and implemented in ABAQUS software via a user subroutine VUMAT. A suitable degradation model was also considered for reducing material constants due to damage. Experimental tests, which performed to validate numerical results, showed that Hashin and Hou failure criteria have better accuracy in predicting force-time history than the other three criteria. However, maximum stress and Hashin failure criteria had the best prediction for damage area, in comparison with the other three criteria. Finally in order to compare numerical model with the experimental results in terms of extent of damage, bending test was performed after impact and the behavior of the beam was considered.

• 한국지반환경공학회 논문집
• /
• 제18권2호
• /
• pp.53-58
• /
• 2017

본 연구에서는 ABAQUS(2014)를 이용한 다공질 매체의 혼합유한요소해석에서 요소 간의 분리를 모사할 수 있는 방법을 제안한다. ABAQUS에서는 변위과 간극수압(u-p모델)의 자유도를 갖는 혼합유한요소의 분리를 standard(implicit) 버전 상에서 cohesive element와 함께 해석을 제안하지만, 요소 간의 이탈, 강체운동, 접촉 등과 같은 분리현상에 대해서는 경계조건 문제로 수치 해석상 한계가 있다. ABAQUS-explicit 해석에서는 경계조건 문제에 대해 자유롭지만 지금까지의 혼합요소 간의 분리를 제공하고 있지 않다. 그러므로, 본 연구에서는 ABAQUS-explicit 상에서 u-p 모델에 대한 분리를 모사할 수 있는 새로운 접근방법이 제안된다. VUMAT 서브루틴을 통하여 구성모델이 적용되고, 간극수압 변화에 따른 요소의 분리 조건을 판단한다. 그리고 VDISP 서브루틴을 통하여 요소의 분리를 발생시킨다. 이렇게 제안된 알고리즘은 간단한 2차원 다공질 매체 예제를 통하여 구현된다.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2008년도 추계학술대회A
• /
• pp.668-673
• /
• 2008

In this paper, we developed a physically-based micromechanical model for inelastic deformation of ceramic powders. The aggregate response of ceramic particles was modeled using the two-surface yield function which considered the shear-induced dilatancy caused by friction, rolling resistance and cohesion between powder particles and consolidation caused by plastic deformation of powder themselves under high compression. The constitutive equations were implemented into the user-subroutine VUMAT of finite element program ABAQUS/Explicit. The material parameters in the constitutive model were identified by calibrating the model to reproduce data from triaxial compression tests and simple compression tests. The density distribution obtained by using the proposed model was in good quantitative agreement with the experimental results of the triaxial compression and cold isostaic compression as well.

• 대한기계학회논문집A
• /
• 제33권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.

• Steel and Composite Structures
• /
• 제17권2호
• /
• pp.159-172
• /
• 2014

A finite element model with the consideration of damage initiation and evolution has been developed for the analysis of the dynamic response of a composite sandwich panel subject to low velocity impact. Typical damage modes including fiber breakage, matrix crushing and cracking, delamination and core crushing are considered in this model. Strain-based Hashin failure criteria with stiffness degradation mechanism are used in predicting the initiation and evolution of intra-laminar damage modes by self-developed VUMAT subroutine. Zero-thickness cohesive elements are adopted along the interface regions between the facesheets and the foam core to simulate the initiation and propagation of delamination. A crushable foam core model with volumetric hardening rule is used to simulate the mechanical behavior of foam core material at the plastic state. The time history curves of contact force and the core collapse area are obtained. They all show a good correlation with the experimental data.

• 소성∙가공
• /
• 제15권3호
• /
• pp.220-225
• /
• 2006

This paper investigates the characteristics of a hydro-mechanical punching process. The hydro-mechanical punching process is divided into two stages: the first stage is the mechanical half piercing in which an upper punch goes down before the initial crack is occurred; the second stage is the hydro punching in which a lower punch goes up until the final fracture is occurred. Ductile fracture criteria such as the Cockcroft, Brozzo and Oyane are adopted to predict the fracture of sheet material. The index values of ductile fracture criteria are calculated with a user material subroutine, VUMAT in the ABAQUS Explicit. The hydrostatic pressure retards the initiation of a crack in the upper region of the blank and induces another crack in the lower region of the blank during the punching process. The final fracture zone is placed at the middle surface of the blank to the thickness direction. The result demonstrates that the hydro-mechanical punching process makes a finer shearing surface than the conventional one as hydrostatic pressure increases.

• 소성∙가공
• /
• 제21권5호
• /
• pp.298-304
• /
• 2012

The shielded slot plates for a molten carbonate fuel cell(MCFC) have a sheared corrugated trapezoidal pattern. In the FEM simulations for the production of the shielded slot plate, the user material subroutine VUMAT in the commercial FEM software ABAQUS was used to implement a ductile fracture criterion. The critical damage value for the ductile fracture criterion was determined by comparing the experimental results of the shearing process with the simulation results. Using the ductile fracture criterion, the FEM simulation of the three-dimensional forming process of the shielded slot plate was conducted. The effects of the shearing process on the forming process were examined through FEM simulation and experiments. The forming simulation of nine unit cells was conducted. Using the simulation results of the forming process, the deformed shape after springback was calculated. The experimental result shows good agreement with the simulation.

• Structural Engineering and Mechanics
• /
• 제81권1호
• /
• pp.39-50
• /
• 2022

Layer separation (delamination) is an essential threat to fiber-reinforced polymer (FRP) plates under dynamic, static, and fatigue loads. Under compressive load, the growth of delamination will lead to structural instability. The aim of this paper is to present a method using shape memory alloy (SMA) stitches to suppress the delamination growth in a FRP plate and to improve the buckling behavior of the plate with a rectangular hole. The present paper is divided into two parts. Firstly, a closed-form (CF) formulation for evaluating the buckling load of the FRP plate is presented. Secondly, the finite element method (FEM) will be employed to calculate the buckling loads of the plates which serves to validate the results obtained from the closed-form method. The novelty of this work is the development of the closed-form solution using the p-Ritz energy approach regarding the stress-dependent phase transformation of SMA to trace the equilibrium path. For the FEM, the Lagoudas constitutive model of the SMA material is implemented in FORTRAN programming language using a user material subroutines (VUMAT). The model is simulated in ABAQUS/Explicit solver due to the nature of the loading type. The cohesive zone model (CZM) is applied to simulate the delamination growth.

• Structural Engineering and Mechanics
• /
• 제87권2호
• /
• pp.107-116
• /
• 2023

The three-dimensional Hashin criterion and user subroutine VUMAT were used to simulate the damage in the composite layer, and the secondary stress criterion was used to simulate the interlayer failure of the cohesive element of the bonding layer and the propagation characteristics under the layer. The results showed that when the shear stress wave (shear wave) propagates on the surface of the laminate, the stress wave attenuation along the fiber strength direction is small, and thus producing a large stress profile. When the compressive stress wave (longitudinal wave) is transmitted between the layers, it is reflected immediately instead of being transmitted immediately. This phenomenon occurs only when the energy has accumulated to a certain degree between the layers. The transmission of longitudinal waves is related to the thickness and the layer orientation. Along the symmetry across the thickness direction, the greater is the stress amplitude along the layer direction. Based on the detailed investigation on the impact on various laminated composites carried out in this paper, the propagation characteristics of stress waves, the damage and the destruction of laminates can be explained from the perspective of stress waves and a reasonable layering sequence of the composite can be designed against damage and failure from low velocity impact.