223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

removed Koen's loading subroutine from abaqus interface

This commit is contained in:
Franz Roters 2009-10-13 06:55:15 +00:00
parent 0a34da292e
commit 2c1a610b14
1 changed files with 1 additions and 59 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

60
code/mpie_cpfem_abaqus.f
View File

@ -214,62 +214,4 @@ subroutine UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,&
integer(pInt) mpie_error
call xit
end subroutine
!############################################################################
!
! include "KJ_Disp.f"
subroutine disp(u,kstep,kinc,time,node,noel,jdof,coords)
! hardwired aba_param.inc
implicit real*8(a-h,o-z)
parameter (nprecd=2)
!
dimension u(3),time(2), coords(3)
real ktime,ktcl,ktmax,ktmin,ktdeltaup,ktdeltadown
real klmin,klmax,kldelta,klbegincycle
real kdeltau, kru, kret
! When using stdb_abqerr for debugging
! dimension intv(2),realv(4)
! character*8 charv(1)
! ratchet displacement per cycle
kru = 0.1
! ratcheting ends at time kret
kret = 400.
! displacement amplitude
kdeltau = 0.5
! time cycle length:
ktcl = 4.
ktmax = ktcl/4
ktmin = 3.*ktmax
ktdeltadown = ktmin - ktmax
ktdeltaup = ktcl - ktdeltadown
! load minimum & maximum:
klmin = -kdeltau
klmax = kdeltau
kldelta = klmax - klmin
klbegincycle = klmin + kldelta * (ktcl-ktmin) / ktdeltaup
! load as a function of (total time); trianglar loading cycle
ktime = time(2)
kru = kru * MIN(kret, ktime) / ktcl
if ( ktime .lt. ktmax ) then
! special case for path to first maximum
u(1) = kru + (ktime/ktmax) * klmax
else
do while ( ktime .ge. ktcl )
ktime = ktime - ktcl
end do
if ( ktime .le. ktmax ) then
u(1) = kru + klbegincycle + ktime * (klmax-klbegincycle) / ktmax
else if ( ktime .lt. ktmin ) then
u(1) = kru + klmax - (ktime-ktmax) * kldelta / (ktmin-ktmax)
else
u(1) = kru + klmin + (ktime-ktmin) * (klbegincycle-klmin)/(ktcl-ktmin)
end if
endif
return
end
!
end subroutine