From 3bd152db2d700145cb9f34b2859c3b6df27bd654 Mon Sep 17 00:00:00 2001 From: Pratheek Shanthraj Date: Mon, 13 Oct 2014 18:08:54 +0000 Subject: [PATCH] temporary fix if marc/abaqus is solving a thermal problem and providing a temperature you want DAMASK to use. --- code/CPFEM.f90 | 9 ++++++++- 1 file changed, 8 insertions(+), 1 deletion(-) diff --git a/code/CPFEM.f90 b/code/CPFEM.f90 index cc156638a..820a29c66 100644 --- a/code/CPFEM.f90 +++ b/code/CPFEM.f90 @@ -327,7 +327,8 @@ subroutine CPFEM_general(mode, ffn, ffn1, temperature, dt, elFE, ip) materialpoint_results, & materialpoint_sizeResults, & materialpoint_stressAndItsTangent, & - materialpoint_postResults + materialpoint_postResults, & + field_putFieldTemperature use IO, only: & IO_write_jobRealFile, & IO_warning @@ -485,6 +486,9 @@ subroutine CPFEM_general(mode, ffn, ffn1, temperature, dt, elFE, ip) !* If no parallel execution is required, there is no need to collect FEM input if (.not. parallelExecution) then +#if defined(Marc4DAMASK) || defined(Abaqus) + call field_putFieldTemperature(ip,elCP,temperature) +#endif materialpoint_F0(1:3,1:3,ip,elCP) = ffn materialpoint_F(1:3,1:3,ip,elCP) = ffn1 @@ -494,6 +498,9 @@ subroutine CPFEM_general(mode, ffn, ffn1, temperature, dt, elFE, ip) if (rnd < 0.5_pReal) rnd = rnd - 1.0_pReal CPFEM_cs(1:6,ip,elCP) = rnd * CPFEM_odd_stress CPFEM_dcsde(1:6,1:6,ip,elCP) = CPFEM_odd_jacobian * math_identity2nd(6) +#endif +#if defined(Marc4DAMASK) || defined(Abaqus) + call field_putFieldTemperature(ip,elCP,temperature) #endif materialpoint_F0(1:3,1:3,ip,elCP) = ffn materialpoint_F(1:3,1:3,ip,elCP) = ffn1