T21: Strain induced precipitates

This tutorial was tested on
MatCalc version 6.00 rel 0.282
license: free
database: mc_fe.tdb; mc_fe.ddb

Complimentary files

Click here to view the script for this tutorial


  • Effect of volumetric misfit on precipitation kinetics
  • Introduction of deformation into the thermo-mechanical treatment
  • Effects of the applied deformation on the precipitation kinetics

Deformation process influences the microstructure of the alloy. One of the effects is the occurence of the strain induced precipitates. This tutorial demonstrates how the deformation processes can be included in the simulation of the processed material. Moreover, the effects of the deformation on the precipitation kinetics are demonstrated and discussed.

Setting up the system

Create a new workspace and open the 'mc_fe.tdb' database. Select the elements 'Fe', 'Nb' and 'C', together with the phases 'BCC_A2', 'FCC_A1' and 'Cementite'. Click on 'Read' instead of 'Read & Close', as the subsequent step is to read the diffusion database in this window. Select 'Diffusion data' on the left side and read the 'mc_fe.ddb' database. Enter the composition of 0.2 wt.% C and 0.04 wt.% Nb. Click on 'Set automatic startvalues' and calculate an equilibrium at the initial 1000°C.

Phase solubilities

In order to have some insight into the system equilibria, perform a stepped equilibrium calculation. In 'Stepped euilibrium ' window select the temperature range of 400-1400°C with the step of 10°C. Click on 'OK'

 MatCalc plot

The phase fractions can be plotted by calling the appropriate user-defined window. In 'View''Create new window', select the 'user defined' tab. There, select '01_all_phase_fractions over T_celsius_logY' and click on 'OK'.

 MatCalc plot

This will result in the plot looking like the one shown below. From this diagram it can be concluded that NbC precipitates below 1300°C and the ferrite will form below 850°C.  MatCalc plot

Precipiation domains and phases

Next, create a precipitation domain called 'austenite' in the 'Precipitation domains …' window. Select FCC_A1 as the thermodynamic matrix phase.

Now, define the precipitate phase settings in 'Phase status' window. Select 'FCC_A1#01' and click on 'Create' and 'precipitate (_Pnn'). In 'Nucleation' → 'Sites' set 'dislocations' as nucleation sites (and remove the checkmark from 'bulk').

 MatCalc plot

Thermo-mechanical treatment

The last step before the initial kinetics simulation is the definition of the thermo-mehanical treatment. In 'Global''Thermo-mech. treatments …' create a new treatment with the name 'cooling'. Next, create a segment in which the austenite domain will be cooled from 1300°C to 750°C with 1°C/s cooling rate

 MatCalc plot

 MatCalc plot

Kinetics simulation of simple cooling process

With all the setup procedures done, perform the kinetics simulation. Click on 'Calc' → 'Precipitation kinetics'. Select the 'cooling' treatment in the 'Temperature control …' area and click on 'Go'.

 MatCalc plot

Once the calculation is completed, use the 'user defined plots' for visualizing the phase fraction, number density and mean radius evolution of the precipitates, as well as the temperature profile. Click on 'View' → 'Create new window', select the 'user-defined' tab and choose the last entry '03_kinetics_4_frames_T_f_n_r_linX'. Click on 'OK'. Set the start of the scaling range to '1'. The analysis of the created plots reveals a somewhat unspectacular precipitaition of NbC phase after 200 s (below 1100°C) that reaches the number density of 1e18 m-3 and stays in the nanometer range.

 MatCalc plot

 MatCalc plot

 MatCalc plot

 MatCalc plot

By clicking on 'Global' → 'Buffers' → 'Rename', rename the current buffer to 'Cooling'.

Effect of volumetric misfit

In general, the precipitates show usually different crystallographic features from the matrix phase. This will be the source of the stress acting on the precipitates within the grain which needs to be considered especially during the nucleation stage. As for now, all the simulations were performed in the previous tutorial were neglecting this factor. Now it will be demonstrated how this effect can be intorduced into the simulation and what is the effect of it.

Click on 'Global' → 'Phase status' and select 'Structure' tab. In the field 'vol. misfit (dV/V)' type in the values of '0.1' for the precipitate phase. This represents the stress coming from the precipitate volume differing by 10% from the matrix volume.

 MatCalc plot

Next, in 'Nucleation' → 'Controls' tab put a checkmark at the 'account for coherent misfit stress' for the precipitate phase.

 MatCalc plot

Click on 'OK' to close the window. Create a new buffer with the name 'Vol_misfit'. With this buffer selected, perform again the kinetic simulation. Once the simulation is completed, duplicate and lock all series in the plots, then switch the displayed buffer to 'Vol_misfit'. As it can be noticed, the definition of the volumetric misfit decreased amount of the precipitated phase.

 MatCalc plot

 MatCalc plot

 MatCalc plot

Simulation of deformation process

When deformation process is to be simulated, there are basically three things to be modified in the treatment:

  • Activation of the dislocation density evolution model.
  • Definition of the deformation rates and temperatures.
  • Deactivation of the volumetric misfit effect for the nucleation during the deformation.

Evolution of dislocation density

The deformation will change the microstructure of the alloy. Here, the change of the dislocation density by the deformation is presented.

Click on 'Global' → 'Precipitation domains'. The simulated deformation will be performed in the austenite region, so select austenite as a precipitation domain and switch to the 'MS Evolution' tab. Inside, select 'Substructure' tab and select '1-param - Sherstnev-Lang-Kozeschnik - 'ABC'' as the model for the substructure evolution.

You might notice some parameters in the 'Dislocation generation and anihilation …' section which describe the generation of the dislocations ('A'-parameter) and their anihilation during the static ('B'-parameter) and dynamic recovery ('C'-parameter). Leave all the parameters there on the default value.

Thermo-mechanical treatmens settings

In this simulation, three deformation runs will be represented. Click on 'Global' → 'Thermo-mech. treatments' and create a new treatment with the name 'deformation'. The first segment will have the starting temperature of 1300°C but the end temperature will be 1000°C with the cooling rate of 1 K/s. In the next segment, change the ramp control setting to 'Accumulated strain' and set the 'Accumulated strain' value to '0.7'.

Afterwards, switch to the 'MS Evolution' tab and type in the strain_rate value of '0.15'.

The definition of the next segments should follow the settings below:

  1. End temperature of 950°C with 1 K/s rate.
  2. Deformation at 950°C with strain rate of 0.15 increasing the strain by 0.7.
  3. End temperature of 900°C with 1 K/s rate.
  4. Deformation at 900°C with strain rate of 0.15 increasing the strain by 0.7.
  5. End temperature of 750°C with 1 K/s rate.

The complete treatment data are shown below.

Deactivation of the volumetric misfit

As previously mentioned, some volumetric misfit might be present related to the precipitates appearing in the grains, decreasing thus their nucleation rate. However, the stress introduced during the applied deformation can cancel this effect. To account for this in MatCalc, clicking on 'Global' → 'Phase status' and select the 'Nucleation' → 'Controls' tab. Put a checkmark at 'ignore misfit stress during deformation' for both precipitate phases. Click on 'OK' to close the window.

Kinetic simulation

Click on 'Global' → 'Buffers' → 'Create' to create a new buffer with the name 'Deformation'. Afterwards, click on 'Calc' → 'Precipitation kinetics' and select the deformation treatment. Start the calculation by clicking on 'OK'. The results of the precipitation kinetics calculation are shown now with the blue curve.

 MatCalc plot  MatCalc plot  MatCalc plot

A considerable difference of three order of magnitudes is observed in the phase fraction of the precipitates. The same difference could be noticed for the number density. In both diagrams, three steps for the increase can be distinguished. The curve for the mean radius has a wavy form showing minima at the time frames of the sudden increases of the phase fraction and number density.

tutorials/t21.txt · Last modified: 2017/08/04 15:23 by pwarczok