From 3064c14ff688be4351cbd46c30f123af02b0ade0 Mon Sep 17 00:00:00 2001 From: Haiming Zhang Date: Tue, 9 Dec 2014 18:56:19 +0000 Subject: [PATCH] fig a bug of energy barrier of vacancy generation. --- code/vacancy_generation.f90 | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/code/vacancy_generation.f90 b/code/vacancy_generation.f90 index efd87eac3..a77bc074a 100644 --- a/code/vacancy_generation.f90 +++ b/code/vacancy_generation.f90 @@ -320,9 +320,9 @@ subroutine vacancy_generation_dotState(nSlip, accumulatedSlip, Tstar_v, Temperat constituent = mappingConstitutive(1,ipc,ip,el) instance = phase_vacancyInstance(phase) pressure = math_trace33(math_Mandel6to33(Tstar_v)) - energyBarrier = vacancy_generation_formationEnergy(instance) - & - pressure*vacancy_generation_atomicVol(instance) - & - sum(accumulatedSlip)*vacancy_generation_plasticityCoeff(instance) + energyBarrier = (vacancy_generation_formationEnergy(instance) - & + pressure)*vacancy_generation_atomicVol(instance) - & + sum(accumulatedSlip)*vacancy_generation_plasticityCoeff(instance) vacancyState(phase)%dotState(1,constituent) = & vacancy_generation_freq(instance)* &