From c0964dfd88772b77814f0ac5069ffcf631abeab2 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Thu, 6 Jun 2019 13:13:22 +0200 Subject: [PATCH] disable coordinate update for the moment --- src/DAMASK_marc.f90 | 13 ++++++------- 1 file changed, 6 insertions(+), 7 deletions(-) diff --git a/src/DAMASK_marc.f90 b/src/DAMASK_marc.f90 index 8c83b96de..822b7abc4 100644 --- a/src/DAMASK_marc.f90 +++ b/src/DAMASK_marc.f90 @@ -29,9 +29,11 @@ #include "prec.f90" module DAMASK_interface + use prec implicit none private + character(len=4), parameter, public :: InputFileExtension = '.dat' character(len=4), parameter, public :: LogFileExtension = '.log' @@ -95,8 +97,7 @@ end subroutine DAMASK_interface_init !-------------------------------------------------------------------------------------------------- !> @brief solver job name (no extension) as combination of geometry and load case name !-------------------------------------------------------------------------------------------------- -function getSolverJobName() - use prec +function getSolverJobName implicit none character(1024) :: getSolverJobName, inputName @@ -115,8 +116,6 @@ end function getSolverJobName end module DAMASK_interface - - #include "commercialFEM_fileList.f90" !-------------------------------------------------------------------------------------------------- @@ -320,7 +319,7 @@ subroutine hypela2(d,g,e,de,s,t,dt,ngens,m,nn,kcus,matus,ndi,nshear,disp, & endif do node = 1,theMesh%elem%nNodes CPnodeID = mesh_element(4+node,cp_en) - mesh_node(1:ndeg,CPnodeID) = mesh_node0(1:ndeg,CPnodeID) + numerics_unitlength * dispt(1:ndeg,node) + !mesh_node(1:ndeg,CPnodeID) = mesh_node0(1:ndeg,CPnodeID) + numerics_unitlength * dispt(1:ndeg,node) enddo endif @@ -332,8 +331,8 @@ subroutine hypela2(d,g,e,de,s,t,dt,ngens,m,nn,kcus,matus,ndi,nshear,disp, & call debug_reset() ! and resets debugging outdatedFFN1 = .false. cycleCounter = cycleCounter + 1 - mesh_cellnode = mesh_build_cellnodes() ! update cell node coordinates - call mesh_build_ipCoordinates() ! update ip coordinates + !mesh_cellnode = mesh_build_cellnodes() ! update cell node coordinates + !call mesh_build_ipCoordinates() ! update ip coordinates endif if (outdatedByNewInc) then computationMode = ior(computationMode,CPFEM_AGERESULTS)