ABAQUS中fric_coef子程序

KHX

想用Fortran语言把这个公式编入.for文档中，但总是出现错误，会的朋友能否指点一下

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

查看代码

打印?

01	subroutine fric_coef (

02	C Write only -

03	*   fCoef, fCoefDeriv,

04	C Read only -

05	*   nBlock, nProps, nTemp, nFields,

06	*   jFlags, rData,

07	*   surfInt, surfSlv, surfMst,

08	*   props, slipRate, pressure,

09	*   tempAvg, fieldAvg

10	*   real a, b, c )

11

C

12	include 'aba_param.inc'

13

C

14	dimension fCoef(nBlock),

15	*   fCoefDeriv(nBlock,3),

16	*   props(nProps),

17	*   slipRate(nBlock),

18	*   pressure(nBlock),

19	*   tempAvg(nBlock),

20	*   fieldAvg(nBlock,nFields)

21

C

22	parameter( iKStep   = 1,

23	*           iKInc    = 2,

24	*           nFlags   = 2

25	*           a        = 0.1,

26	*           b        = -1.2,

27	*           c        = 0.3)

28

C

29	parameter( iTimStep = 1,

30	*           iTimGlb  = 2,

31	*           iDTimCur = 3,

32	*           nData    = 3,)

33

C

34	dimension jFlags(nFlags), rData(nData)

35

C

36	character*80 surfInt, surfSlv, surfMst

37

C

38	xMuk = props(1)

39	xMus = props(2)

40	beta = props(3)

41	DO  k = 1, nBlock

42	fCoef(k)=(a*exp(b*pressure(k))+c)

43	*[xMuk+(xMus-xMuk)*exp(-beta*slipRate(k))]

44

END DO

45

46

return

47

end

青衣巷

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

查看代码

打印?

01	subroutine fric_coef (

02	C Write only -

03	*   fCoef, fCoefDeriv,

04	C Read only -

05	*   nBlock, nProps, nTemp, nFields,

06	*   jFlags, rData,

07	*   surfInt, surfSlv, surfMst,

08	*   props, slipRate, pressure,

09	*   tempAvg, fieldAvg )

10

C

11	include 'aba_param.inc'

12	real :: a, b, c

13

C

14	dimension fCoef(nBlock),

15	*   fCoefDeriv(nBlock,3),

16	*   props(nProps),

17	*   slipRate(nBlock),

18	*   pressure(nBlock),

19	*   tempAvg(nBlock),

20	*   fieldAvg(nBlock,nFields)

21

C

22	parameter( iKStep   = 1,

23	*           iKInc    = 2,

24	*           nFlags   = 2,

25	*           a        = 0.1,

26	*           b        = -1.2,

27	*           c        = 0.3)

28

C

29	parameter( iTimStep = 1,

30	*           iTimGlb  = 2,

31	*           iDTimCur = 3,

32	*           nData    = 3)

33

C

34	dimension jFlags(nFlags), rData(nData)

35

C

36	character*80 surfInt, surfSlv, surfMst

37

C

38	xMuk = props(1)

39	xMus = props(2)

40	beta = props(3)

41	DO  k = 1, nBlock

42	fCoef(k)=(a*exp(b*pressure(k))+c)

43	*             *(xMuk+(xMus-xMuk)*exp(-beta*slipRate(k)))

44

END DO

45

46

return

47

end

KHX

青衣巷 发表于 2021-11-19 15:58
[mw_shl_code=fortran,true]      subroutine fric_coef (
C Write only -
     *   fCoef, fCoefDeriv,

很感谢您的回复，我试着运行了一下，又出现了这个错误the implicit statement is not position corrently within the scoping unit，我查了一下好像是implicit语句位置不对，刚接触Fortran语言不是很懂，想再请教一下

青衣巷

       include 'aba_param.inc'
         real :: a, b, c

交换一下位置

KHX

青衣巷 发表于 2021-11-20 10:33
include 'aba_param.inc'
         real :: a, b, c

能运行了，谢谢