EvgenAbaqus
Mechanical
- Jan 27, 2015
- 17
Hello,
I have a txt-file which looks like this:
1 2.
2 2.
3 2.
........
5758 4.
5759 4.
5760 4.
where 1st column - elem number, 2nd -some parameter
I'd like read data from text file into subroutine SDVINI
and I need use this data for further purposes in UMAT subroutine.
SUBROUTINE UMAT whitout SDVINI works ok, but with SDVINI gives an error in Abaqus....
I think it's because I'm not correctly written code SDVINI subroutine.
please check my code
SUBROUTINE SDVINI(STATEV,COORDS,NSTATV,NCRDS,NOEL,NPT,LAYER,KSPT)
C
INCLUDE 'aba_param.inc'
C
DIMENSION STATEV(NSTATV),COORDS(NCRDS)
real, dimension(5760,1) :: RHON
STATEV(7) = 0
IF (STATEV(7).eq.0) then
fname='D:\d113.txt'
Open(66,file=fname,form='formatted',status='old')
READ(66,*) RHON
STATEV(1)= RHON
CLOSE(66)
ENDIF
RETURN
END
and in UMAT i need use next eq.:
EMOD=STATEV(1)**2.5
so i need use in umat parametr from txt-file
but i dont know how to read correctly data from file
I would really happy if you can help me in this.
I have a txt-file which looks like this:
1 2.
2 2.
3 2.
........
5758 4.
5759 4.
5760 4.
where 1st column - elem number, 2nd -some parameter
I'd like read data from text file into subroutine SDVINI
and I need use this data for further purposes in UMAT subroutine.
SUBROUTINE UMAT whitout SDVINI works ok, but with SDVINI gives an error in Abaqus....
I think it's because I'm not correctly written code SDVINI subroutine.
please check my code
SUBROUTINE SDVINI(STATEV,COORDS,NSTATV,NCRDS,NOEL,NPT,LAYER,KSPT)
C
INCLUDE 'aba_param.inc'
C
DIMENSION STATEV(NSTATV),COORDS(NCRDS)
real, dimension(5760,1) :: RHON
STATEV(7) = 0
IF (STATEV(7).eq.0) then
fname='D:\d113.txt'
Open(66,file=fname,form='formatted',status='old')
READ(66,*) RHON
STATEV(1)= RHON
CLOSE(66)
ENDIF
RETURN
END
and in UMAT i need use next eq.:
EMOD=STATEV(1)**2.5
so i need use in umat parametr from txt-file
but i dont know how to read correctly data from file
I would really happy if you can help me in this.