学习abaqus岩土工程中的应用，一个程序过不去

蛰龙

[Fortran] 纯文本查看 复制代码

查看代码

打印?

01	program test

02	implicit none

03	call FRIC()

04	call FRIC()

05	end program test

06	SUBROUTINE FRIC(LM,TAU,DDTDDG,DDTDDP,DSLIP,SED,SFD,

07	1 DDTDDT,PNEWDT,STATEV,DGAM,TAULM,PRESS,DPRESS,DDPDDH,SLIP,

08	2 KSTEP,KINC,TIME,DTIME,NOEL,CINAME,SLNAME,MSNAME,NPT,NODE,

09	3 NPATCH,COORDS,RCOORD,DROT,TEMP,PREDEF,NFDIR,MCRD,NPRED,

10	4 NSTATV,CHRLNGTH,PROPS,NPROPS)

11

C

12	INCLUDE 'ABA_PARAM.INC'

13

C

14	CHARACTER*80 CINAME,SLNAME,MSNAME

15

C

16	DIMENSION TAU(NFDIR),DDTDDG(NFDIR,NFDIR),DDTDDP(NFDIR),

17	1 DSLIP(NFDIR),DDTDDT(NFDIR,2),STATEV(*),DGAM(NFDIR),

18	2 TAULM(NFDIR),SLIP(NFDIR),TIME(2),COORDS(MCRD),

19	3 RCOORD(MCRD),DROT(2,2),TEMP(2),PREDEF(2,*),PROPS(NPROPS)

20

21

22	if(LM.EQ.2) return

23	FK1=PROPS(1)

24	FK2=PROPS(2)

25	FN=PROPS(3)

26	FRF=PROPS(4)

27	FFAI=TAN(PROPS(5)/180.*3.1415926)

28	FGW=PROPS(6)

29	FPA=PROPS(7)

30	xpress=press-dpress

31	IF(XPRESS.LE.10)THEN

32	XPRESS=10.

33

END IF

34	FEI1=FRF*abs(TAU(1))/(XPRESS*FFAI)

35	FEI2=FRF*abs(TAU(2))/(XPRESS*FFAI)

36	IF(FEI1.GE.0.99)FEI1=0.99

37	IF(FEI2.GE.0.99)FEI2=0.99

38	XK1=(1-FEI1)**2*FK1*FGW*(XPRESS/FPA)**FN

39	XK2=(1-FEI2)**2*FK2*FGW*(XPRESS/FPA)**FN

40

LM=0

41	TAU(1)=TAU(1)+XK1*DGAM(1)

42	TAU(2)=TAU(2)+XK2*DGAM(2)

43	DDTDDG(1,1)=XK1

44	DDTDDG(2,2)=XK2

45	DDTDDG(1,2)=0.0

46	DDTDDG(2,1)=0.0

47

return

48

end

学习abaqus岩土工程中的应用,费康，过不去这个子程序

楚香饭

子程序 FRIC 有十多个参数，而你的主程序里 call FRIC() 一个参数都不给。能通过才怪了。