abaqus磨损分析
模型可简化为两个方块压在一起的周期性的相互滑动(尺寸不一样),材料可设为钢,压紧力自己给。目的:1 磨损率;2 预测较小滑块在一定滑动一定周期后的磨损量,参数不全可以自己给主要是作为参考。
目前建模中的难点:将archard损伤模型写入自适应子程序
提示:帮助文档有个轮胎模型的自适应子程序,本模型的子程序算法应该比帮助文档容易的多,因为只有平东,磨损节点一直处于连续接触过程中,没有转动时节点接触状态不连续的情况——可惜本人才,abaqus接触时间短,对fortran更是完全陌生。
关于奖励:除了仿真币外,另外可以追加200——500RMB的报酬。
邮箱yun417182775@yeah.net
模型可简化为两个方块压在一起的周期性的相互滑动(尺寸不一样),材料可设为钢,压紧力自己给。
目的:1 磨损率;2 预测较小滑块在一定滑动一定周期后的磨损量,参数不全可以自己给,主要是作为参考。
目前建模中的难点:将archard损伤模型写入自适应子程序
提示:帮助文档有个轮胎模型的自适应子程序,本模型的子程序算法应该比帮助文档容易的多,因为只有平动,磨损节点一直处于连续接触过程中,没有转动时节点接触状态不连续的情况——可惜本人不才,abaqus接触时间短,对fortran更是完全陌生。
关于奖励:除了仿真币外,另外可以追加200——500RMB的报酬。
邮箱yun417182775@yeah.net
稍后给出我的INP文件。
附件INP文件,跪求ABAQUS先知指导 417182775 发表于 2012-3-25 22:31
附件INP文件,跪求ABAQUS先知指导
哥们儿弄出来了吗?我现在也是做这个 钢丝刷洗属于这个问题么 :lol,路过路过 帮顶。。。。。。。。
,路过路过 高大上啊 在线等大神:) 路过,楼主成功了吗,可以分享一下程序吗 这是Fortran语言的一个磨损的例子,涉及到的archard公式比较简单,我是做模具磨损;你可以参考一下
C USER INPUT FOR ADAPTIVE MESH CONSTRAINT
C
SUBROUTINE UMESHMOTION(UREF,ULOCAL,NODE,NNDOF,
$ LNODETYPE,ALOCAL,NDIM,TIME,DTIME,PNEWDT,
$ KSTEP,KINC,KMESHSWEEP,JMATYP,JGVBLOCK,LSMOOTH)
C
include 'ABA_PARAM.INC'
CHARACTER*80 PARTNAME
DIMENSION ARRAY(3), JPOS(15)
DIMENSION ULOCAL(NDIM)
DIMENSION ALOCAL(NDIM,*), TIME(2)
DIMENSION JMATYP(*), JGVBLOCK(*)
PARAMETER (NELEMMAX=100)
DIMENSION JELEMLIST(NELEMMAX), JELEMTYPE(NELEMMAX)
PARAMETER (k = 0.00000251)
CALL GETNODETOELEMCONN(NODE,NELEMS,JELEMLIST,JELEMTYPE,
$ JRCD,JGVBLOCK)
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
CALL GETVRMAVGATNODE(LOCNUM,JTYP,'CDISP',ARRAY,JRCD,
$ JELEMLIST,NELEMS,JMATYP,JGVBLOCK)
CSLIP = SQRT(ARRAY(2)**2+ARRAY(3)**2)
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
CALL GETVRMAVGATNODE(LOCNUM,JTYP,'CSTRESS',ARRAY,JRCD,
$ JELEMLIST,NELEMS,JMATYP,JGVBLOCK)
CPRESS = ARRAY(1)
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
if (CPRESS.GT.0) then
ULOCAL(NDIM) = -1 * k * CSLIP * CPRESS
end if
RETURN
END
另外,如果往复式(反复)磨损,需要考虑磨损之后的表面节点的变化,这就需要用到python提取表面节点坐标,然后修改模型的inp文件里的节点坐标,重构表面,再进行第二次磨损分析,依次下去。。。:)我也正在做这个 同求啊,有偿的
页:
[1]