[求助]关于摩擦模型的二次开发

leojcf

我做的卷板成形模拟，使用经典库伦摩擦力的话，输出的辊子与板材间的摩擦力曲线不合理。

现在想使用修正的库伦摩擦模型，参考的一篇文献为《基于ABAQUS的板料成形数值模拟摩擦模型的二次开发》，摘一段关于这个修正库伦摩擦模型的内容，大家帮忙看看，我该看些什么东西。我对二次开发不太了解。不知道从哪开始学起。

“修正的库仑摩擦模型在ABAQUS中的实现

　　在ABAQUS/Explicit中如果要选用其他摩擦模型,就必须利用它给用户提供的VFRIC用户子程序接口定义所需的摩擦模型。在VFRIC用户子程序中, ABAQUS/Explicit给用户提供了许多参数,利用这些参数,就可以通过Visual Fortran编写必要的程序,实现与ABAQUS/Explicit之间的数据交流,从而达到摩擦模型的定义。”

有使用过修正的库伦摩擦模型的吗？交流一下啊。我QQ是65204612

在ABAQUS/Explicit中,其用户子程序接口
程序为:
subroutine vfric (
C Write only -
1 fTangential,
C Read/Write -
2 statev,
C Read only -
3 kStep, kInc, nContact, nFacNod, nSlvNod, nMstNod,
4 nFricDir, nDir, nStateVar, nProps, nTemp, nPred,
numDefTfv,
5 jSlvUid, jMstUid, jConSlvid, jConMstid, timStep,
timGlb,
6 dTimCur, surfInt, surfSlv, surfMst, lContType,
7 dSlipFric, fStickForce, fTangPrev, fNormal, frictionWork,
8 shape, coordSlv, coordMst, dirCosSl, dircosN, props,
9 areaSlv, tempSlv, preDefSlv, tempMst, preDefMst)
C
include `vaba _ param·inc'
C
character*8 surfInt, surfSlv, surfMst
C
dimension props (nProps), statev (nStateVar, nSlvNod),
1 dSlipFric (nDir, nContact),
2 fTangential (nFricDir, nContact),
3 fTangPrev (nDir, nContact),
4 fStickForce (nContact), areaSlv (nSlvNod),
5 fNormal (nContact), shape (nFacNod, nContact),
6 coordSlv (nDir, nSlvNod), coordMst (nDir, nMstNod),
7 dirCosSl (nDir, nContact), dircosN (nDir, nContact),
8 jSlvUid (nSlvNod), jMstUid (nMstNod),
9 jConSlvid (nContact), jConMstid (nFacNod, nContact)
1 tempSlv (nContact), preDefSlv (nContact, nPred),
2 tempMst (numDefTfv), preDefMst (numDefTfv, nPred)
user coding to define fTangential

and, optionally, statev
return
end
程序中各变量的定义限于篇幅,不一一说明。
用户子程序在command窗口运行,其格式为:

abq651 job=<inp文件名> user=<子程序文件名>。