[Fortran] 纯文本查看 复制代码 subroutine fric_coef (
C Write only -
* fCoef, fCoefDeriv,
C Read only -
* nBlock, nProps, nTemp, nFields,
* jFlags, rData,
* surfInt, surfSlv, surfMst,
* props, slipRate, pressure,
* tempAvg, fieldAvg )
C
include 'aba_param.inc'
real :: a, b, c
C
dimension fCoef(nBlock),
* fCoefDeriv(nBlock,3),
* props(nProps),
* slipRate(nBlock),
* pressure(nBlock),
* tempAvg(nBlock),
* fieldAvg(nBlock,nFields)
C
parameter( iKStep = 1,
* iKInc = 2,
* nFlags = 2,
* a = 0.1,
* b = -1.2,
* c = 0.3)
C
parameter( iTimStep = 1,
* iTimGlb = 2,
* iDTimCur = 3,
* nData = 3)
C
dimension jFlags(nFlags), rData(nData)
C
character*80 surfInt, surfSlv, surfMst
C
xMuk = props(1)
xMus = props(2)
beta = props(3)
DO k = 1, nBlock
fCoef(k)=(a*exp(b*pressure(k))+c)
* *(xMuk+(xMus-xMuk)*exp(-beta*slipRate(k)))
END DO
return
end
|