Continue to Site

Eng-Tips is the largest engineering community on the Internet

Intelligent Work Forums for Engineering Professionals

  • Congratulations MintJulep on being selected by the Eng-Tips community for having the most helpful posts in the forums last week. Way to Go!

Simulation problem which using user subroutine

Status
Not open for further replies.

Epik

Mechanical
Jun 8, 2016
1
Hello.

I'm a ABAQUS user. I have some problem about my analysis. So I ask about my simulation.

First, my PC system is below.

- ABAQUS 6.12-1
- Intel fortran compiler
- visual studio 2008

My simulation and user subroutine is below.

- Soil step, considering permeability with heat transfer
- UMAT
- UMATHT
- USDFLD

I combine the all subroutine in one file. So analysis was successfully finish every time.

But here is the problem.

I use the USDFLD for change the permeability following the STATEV(2) value.

But permeability didn't change during analysis.

I check the Field variable value (FV1) and FV1 change following the STATEV(2) value.

URL]


In the analysis, only first permeability(9.741e-08, 0.8, 1) was applied during all analysis.

It was compared using none USDFLD used analysis. (only first permeability value was used)

When I change the first permeability value(9.741e-08, 0.8, 1), then result was changed.

And if I change the second permeability value(4.8705e-08, 0.8, 2), then result was not changed.

Below is the my analysis data part.

I'll wait your's professional many comment.

Thank you.

Materials.

*Permeability, Type=iso, dependencies = 1, specific=10000
9.741e-08, 0.8, 1
4.8705e-08, 0.8, 2
9.741e-07, 0.8, 3
3.6042e-06, 0.8, 4
*User Defined Field

Subroutine of USDFLD

subroutine usdfld(field,statev,pnewdt,direct,t,celent,time,dtime,
1 cmname,orname,nfield,nstatv,noel,npt,layer,kspt,kstep,kinc,
2 ndi,nshr,coord,jmac,jmatyp,matlayo,laccfla)
c
include 'aba_param.inc'
c
character*80 cmname,orname
character*8 flgray(15)
dimension field(nfield),statev(nstatv),direct(3,3),t(3,3),time(2)
dimension array(15),jarray(15),jmac(*),jmatyp(*),coord(*)
if (KSTEP == 1) THEN
Field(1)=1.D0
else if (KSTEP == 2) THEN
If (STATEV(2) .LE. 20) Then
Field(1)=1.D0
else if (STATEV(2) .GT. 20 .and. STATEV(2) .LE. 1000) THEN
Field(1)=2.D0
else if (STATEV(2) .GT. 1000 .and. STATEV(2) .LE. 1200) THEN
Field(1)=3.D0
else if (STATEV(2) .GT. 1200) THEN
Field(1)=4.D0
end if
end if
return
end
 
Status
Not open for further replies.

Part and Inventory Search

Sponsor