سابروتين نويسی usdfld در آباکوس
با سلام...
من برای آناليز flow که دارم انجام ميدم ..نياز دارم که ضريب نفوذپذيری رو تابع فشار جانبی
در نظر بگيرم..که بايد سابروتين USDFLD رو بنويسم..
سابروتينی که من استفاده ميکنم وقت آناليز..هيچ error اي وقت آناليز نشون نميده و
مدل run ميشه....ولی جوابها وقت استفاده از سابروتين و بدون استفاده از اون يکی
در ميآن....
که نميدونم مشکل از سابروتين هست يا در محيط CAE اشتباه فراخونی ميکنم..
ممنون ميشم راهنمايی و کمکم کنيد
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(*)
C Constant values for uncrack and loading concrete
C a0 = 0
C a1 = 382.5651
C b = 0.7808
C T = 25
call getvrm('SP',array,jarray,flgray,jrcd,
&jmac, jmatyp, matlayo, laccfla)
SPS = abs(array(1))
C Minimum value of principal stress ( Confining Pressure)
call getvrm('SDV',array,jarray,flgray,jrcd,
&jmac, jmatyp, matlayo, laccfla)
SPSMIN = array(1)
C Use the min stress as a field variable
field(1)=Min(SPS, SPSMIN)
C Store the min stress as a solution dependent state
C variable
statev(1)=field(1)
K = (40.151*( 10*EXP6))*EXP( -statev(1) / 0.7807)
C If error, write comment to .DAT file:
IF(JRCD.NE.0)THEN
WRITE(6,*) 'REQUEST ERROR IN USDFLD FOR ELEMENT NUMBER ',
& NOEL,'INTEGRATION POINT NUMBER ',NPT
ENDIF
C
return
end