今天让我们学习vuviscosity子程序。这是一个 Abaqus/Explicit 中用于定义**粘性行为(Viscosity)**的用户子程序,适用于 Viscous 材料模型(例如模拟高分子材料、流体或非牛顿流体)。
subroutine vuviscosity (C Read only -* nblock,* jElem, kIntPt, kLayer, kSecPt,* stepTime, totalTime, dt, cmname,* nstatev, nfieldv, nprops,* props, tempOld, tempNew, fieldOld, fieldNew,* stateOld,* shrRate,C Write only -* viscosity,* stateNew )Cinclude 'vaba_param.inc'Cdimension props(nprops),* tempOld(nblock),* fieldOld(nblock,nfieldv),* stateOld(nblock,nstatev),* shrRate(nblock),* tempNew(nblock),* fieldNew(nblock,nfieldv),* viscosity(nblock),* stateNew(nblock,nstatev)Ccharacter*80 cmnameCparameter ( one = 1.d0 )CC Cross viscosityCeta0 = props(1)rlambda = props(2)rn = props(3)Cdo k = 1, nblockviscosity(k) = eta0/(one+(rlambda*shrRate(k))**(one-rn))end doCreturnend






