10 Useful Abaqus Subroutine (Part 2)

engineering applications of finite element

Before reading the article “engineering applications of finite element method”, you can read the first part of article: “10 Useful Abaqus Subroutine part-1

Introduction of ABAQUS subroutine | Engineering applications of finite element method

When you want to simulate an analysis in the ABAQUS, you may encounter some complications that require user-defined coding. So, ABAQUS has left you a coding platform for this matter. It’s called the ABAQUS subroutine. Here, we have presented some Abaqus user subroutine examples. We explained some of them with engineering applications of finite element method.


When the extended versions of the classical Coulomb friction model provided in the ABAQUS can no longer help you with your related problem, and when you need a complex definition of shear transmission between contacting surfaces, these two could come in handy. In other words, if you need contact with user-defined friction, the FRIC (standard solver) or VFRIC (explicit solver) subroutine can be used.

1.1. Some FRIC & VFRIC instances

1.1.1. Pile foundation

In 2020, the Geotechnical and Structural Engineering Research Center at Shandong University of China investigated on a finite element analysis with the ABAQUS subroutine that predicts the pile groups’ response to individual piles’ different layouts in pile groups. The FRIC subroutine was developed as the secondary platform, a hyperbolic model of end resistance, and a softening model of skin friction for the contact pair calculation of the ABAQUS.

contour plot of deformation

Figure 1: Response of a pile group under the total load of 6400KN [Ref]

engineering applications of finite element ⭐⭐⭐ Abaqus Course |10 hours Video  👩‍🎓+1000 Students   ♾️ Lifetime Access

✅ Module by Module Training                                  ✅ Standard/Explicit Analyses Tutorial

✅ Subroutines (UMAT) Training                    …           ✅ Python Scripting Lesson & Examples

…………………………                 …………………….. ……………   …………………………………..

1.1.2.   Subsea pipelines

The submarine pipelines are mostly used to transport Oil and Gas.

In 2019, the School of Civil Engineering University of Tehran developed research about the local buckling of subsea pipelines. A model was simulated by the ABAQUS software and FRIC subroutine that can eliminate the walking occurrence. The pipe-soil interaction model was implemented by the FRIC subroutine, which controls the lateral and axial movement of the pipeline.

subsea pipelines

Figure 2: Subsea pipelines

Pipe analysis during lateral displacement | engineering applications of finite element

Figure 3: Cyclic pipe-soil interaction response during lateral displacements [Ref]

1.1.3. Tire- asphalt pavement contact behavior

In 2020, the National Science Foundation of China funded an investigation about tire-pavement contact behaviors under dry conditions. The researchers used the ABAQUS software to simulate a 3D finite element model to analyze these behaviors. Also, they used the FRIC subroutine to define the contact between the tire and the asphalt pavement.

Tire- asphalt contact simulation

Figure 4: Simulation of the Tire-asphalt pavement contact [Ref]

Tire tracks on asphalt

Figure 5: Tire tracks on asphalt

1.1.4. Metal cutting modeling of titanium alloy

Titanium alloys are hard metals. They are used in many industries, such as aerospace, automotive, etc. Machining these alloys is a challenging process.

In 2011, China Scholarship Council (CSC) supported research to evaluate stress and damage during the cutting process. The researchers simulated a titanium alloy cutting with the help of the ABAQUS subroutine. The VFRIC subroutine developed an improved friction model and an orthogonal cutting model.

contour plot of cutting process

Figure 6: Cutting process simulation [Ref]

Note: If you want to learn more about this useful subroutine and see some examples, click on this picture:

Introduction to VFRICTION and VFRIC Subroutines in ABAQUS

2. UMESHMOTION subroutine

If you want to define motion for mesh nodes, you must use the UMESHMOTION subroutine. This subroutine helps you to specify Mesh Motion constraints during adaptive meshing. This subroutine is called at the end of any increment where adaptive meshing is performed. Engineering applications of finite element method in the UMESHMOTION subroutine are described below.

2.1. Engineering applications of finite element method for the UMESHMOTION subroutine

2.1.1. Artificial knee joints

Department of Mechanical Engineering University of Tokyo in 2020 worked on wear simulation of artificial knee joints. The investigators used the ABAQUS software and UMESHMOTION subroutine to evaluate the wear of artificial joints. The UMESHMOTION subroutine was used to update the geometry of the contact area. This is the important Abaqus user subroutine example.

Artificial knee joint

Figure 7: Total knee replacement

Simulation of knee joint

Figure 8: Simulation of total knee replacement [Ref]

2.1.2. Fretting wear

In 2020, China Scholarship Council funded an investigation about Fretting wear finite element analysis under variable coefficient of friction. The fretting wear may cause losing fasteners or other problems; moreover, it is challenging to measure micromotion in experiments; So, FEM simulation and the ABAQUS subroutine can analyze its behavior. The UMESHMOTION subroutine executed the continuous change of the contact surface geometry.

Fretting wear modeling

Figure 9: Dimensions and geometry of the basic model [Ref]

Note: You’re interested in this complicated yet very useful subroutine?! Want to see some examples? Click on this figure:

UMESHMOTION Subroutine in ABAQUS-package

3. HETVAL subroutine

This subroutine provides internal heat generation during heat transfer analysis. For instance, it could be associated with phase changes occurring during the solution. This subroutine could depend on state variables such as the fraction of material transformed.

3.1. HETVAL applications

3.1.1. Phase-field fracture method

The Phase-field method is a common technique to simulate crack initiation, propagation, and coalescence without the need to trace the fracture surface.

Department of Construction and Manufacturing Engineering University of Oviedo in 2021 presented a robust and simple application of the phase-field fracture method in the ABAQUS software with the HETVAL subroutine. With this method, you can use the ABAQUS’ in-built features and bypass the need for defining user elements. The HETVAL subroutine defined the internal heat flux and its derivative with respect to the temperature.

contour plot of screw tension | engineering applications of finite element

Figure 10: Final Phase-field contours for screw tension tests [Ref]

3.1.2. Rolling tire

In 2018, Korea Automotive Technology Institute developed a simulation method on a dynamic rolling tire to predict temperature rise. The investigators used the ABAQUS subroutine for the simulation. The HETVAL subroutine calculated the heat generation rate using total dissipated strain energy. With this Abaqus user subroutine example, you can make better designs.

contour plot of a tire temperature

Figure 11: Rolling tire temperature distribution [Ref]

3.1.3. Carbon-Carbon composites

ASME congress in 2015 disseminated research that predicts the influence of thermal shock conditions on C/C composites. The researchers worked on this investigation with the help of the ABAQUS subroutines, the HETVAL and UMAT. The HETVAL subroutine implemented the generated internal heat caused by the carbon material degradation.

Heat flux distribution in c-c composite

Figure 12: Heat flux distribution within the C/C composite [Ref]

3.1.4. P91 steel welding

National Key Research and Development Program of China in 2021 supported an investigation that found a method to calculate the welding residual Stress of P91 steel. They used ABAQUS subroutines such as the HETVAL to achieve this method. The HETVAL subroutine developed the phase transformation latent heat generated by the martensite transformation.

contour plot of welding process

Figure 13: Stress distribution of one of the cases in the welding process [Ref]

I hope you enjoyed engineering applications of finite element method in the mentioned subroutines.

4. DFLUX subroutine

This subroutine can be used to define nonuniformly distributed flux in a mass diffusion or heat transfer analysis. It can be defined as a function of time, position, temperature, integration point number, element number, etc.

engineering applications of finite element

Figure 14: Arc welding simulation using DFLUX

4.1. Several applications of the DFLUX subroutine

4.1.1. Asphalt pavement

In 2018, the Dalian University of Technology in China developed a simulation to observe the Asphalt pavement behavior under cyclic temperature when a reflective crack is applied. The researchers used the ABAQUS subroutines such as the DFLUX and FILM. The DFLUX subroutine calculated the temperature distribution in the pavement structure. This is just one instance of this Abaqus user subroutine example. You can search and find many of them in the web.

Temperature distribution in asphalt

Figure 15: Temperature distribution in asphalt pavement structure [Ref]

4.1.2. Laser cladding

At the 46th SME conference held in 2018, a numerical analysis was presented to determine the hardness and residual stresses of multi-layered cladding for applications of die repair. This analysis was implemented by the ABAQUS subroutine. The DFLUX subroutine applied the heat load.

Laser cladding

Figure 16: Laser cladding

Boundary condition of the analysis

Figure 17: Boundary conditions and computational domains of the
coupled metallo-thermomechanical analysis [Ref]

4.1.3. Additive manufacturing (3D printing)

The 3D printing process uses a high-power laser to melt compositions of metallic feedstock from a powder bed.

In 2018, the Department of Mechanical Engineering University of Sheffield published an article that developed a model to assess residual stresses in selective laser-melted Ti6Al4V. The investigators simulated the model with the ABAQUS software and the DFLUX subroutine. They used the subroutine to write a moving volumetric heat source program to simulate the laser. An Abaqus user subroutine example like this is one of the most important areas of research nowadays.

contour plot of the analysis

Figure 18: Location for temperature estimation [Ref]

4.1.4. Friction stir welding

In 2019, Pandit Deendayal Petroleum University of India developed a 3D finite element model with the ABAQUS subroutine to predict the thermal gradient of cooling-assisted friction stir welding. The generated heat during friction stir welding was applied as the surface flux over time. This was done by the DFLUX subroutine.

Friction stir welding

Figure 19: Friction stir welding

Thermal distribution in welding

Figure 20: Thermal profile during weld stage [Ref]

Note: Tutorial and more examples of this fascinating subroutine:

DFLUX subroutine in ABAQUS-package

5. DISP & VDISP subroutine

The last ABAQUS subroutines that we are going to explain are DISP and VDISP. These are for specifying prescribed boundary conditions. The DISP is used in the ABAQUS/Standard solver, while the VDISP is used for the ABAQUS/Explicit solver. When you need to define time-dependent and location-dependent boundary conditions simultaneously, while the ABAQUS features cannot support these problems, the DISP and VDISP subroutines will come to your aid.

5.1. DISP and VDISP applications

5.1.1. Groundwater

In 2018, China Postdoctoral Science Foundation funded an article that simulated a finite element model by the ABAQUS and DISP subroutine to investigate the groundwater level fluctuation over time.

meshed model

Figure 21: Calculation model [Ref]

5.1.2. Soft tissue

In 2016, the Department of Mechanical Engineering of Sharif University developed a constitutive and micromechanics modeling of connective soft Tissues. The researchers simulated the model with the ABAQUS subroutines such as UMAT and DISP. They used the DISP subroutine to model the periodic boundary conditions.

Rat tail stress distribution

Figure 22: Shear stress distribution of rat tail tendon [Ref]

5.1.3. Sheet metal forming

In 2021, Shenzhen Key Laboratory of Cross-scale Manufacturing Mechanics investigated on sheet metal forming assisted by Ultrasonic vibration. The investigators simulated an analysis with the ABAQUS subroutines such as VUMAT and VDISP. The VDISP subroutine established the forming tool’s predefined path.

Boundary condition of the model

Figure 23: Boundary conditions and predefined forming tool path [Ref]

Note: Need more info about this subroutine? Check this out:


Get this post as a PDF file: 10-useful-abaqus-subroutinepart2

One thought on “10 Useful Abaqus Subroutine (Part 2)

  1. User Avatar arun says:

    This article was great. I read part 1 of this article and I was eager to read part 2. Now, in your opinion, which educational package should I buy to start learning?

Leave a Reply