abaqus的fric子程序实例，网上找.docxVIP

subroutine fric(lm,tau,ddtddg,ddtddp,dslip,sed,spd, 1 ddtddt,pnewdt,statev,dgam,taulm,press,dpress,ddpddh, 2 slip,kstep,kinc,time,dtime,noel,ciname,slname, 3 msname,npt,node,npatch,coords,rcoord,drot,temp, 4 predef,nfdir,mcrd,npred,nstatv,chrlngth,props,nprops) c include aba_param.inc c character*80 ciname,slname,msname dimension tau(nfdir),ddtddg(nfdir,nfdir),ddtddp(nfdir), 1 dslip(nfdir),ddtddt(nfdir,2),statev(*), 2 dgam(nfdir),taulm(nfdir),slip(nfdir),time(2), 3 coords(mcrd),rcoord(mcrd),drot(2,2),temp(2), 4 predef(2,*),props(nprops) c 摩擦系数 F = 0.7 C 对对对对对对对对对对对对对对对由相滑移量判断滑移与粘着的界 CRI = 1E-7 c 没有用到以下量直接清零 DDTDDG( 1 , 1 ) = 0.0 DDTDDG( 1 , 2 ) = 0.0 DDTDDG( 2 , 1 ) = 0.0 DDTDDG( 2 , 2 ) = 0.0 c 对对对对对对对对如果上一迭代分离直接返回 IF( LM == 2 ) THEN return END IF c 对对对对对对对对对对对对对对对对对对对对对 如果接触力下次迭代的接触状更新分离 IF( PRESS = 0 ) THEN LM = 2 RETURN END IF IF( LM == 1 ) THEN c 对对对对对对对对对对对对对对对判断算果是否和假接触状一致 IF( DABS( TAULM( 1 ) ) = DABS( F * PRESS ) ) THEN c 对不一致重新确定接触状 LM = 0 IF( TAULM( 1 ) 0 ) THEN TAU( 1 ) = -F * PRESS DDTDDP( 1 ) = -F ELSE TAU( 1 ) = F * PRESS DDTDDP( 1 ) = F END IF RETURN END IF c 对对对对对对对对 一致且粘着直接返回 RETURN END IF IF( LM == 0 ) THEN TT = DGAM( 1 ) c 对对对对对对对对对对对对对对对 判断算果是否和假接触状一致 IF( DABS( TT ) CRI ) THEN c 对对对对对对对对对对对对对对对对对对对对不一致改接触状定接触状粘着直接返回 LM = 1 RETURN END IF c 对对对对对对对对对对对对对对对对对对对 一致界面状滑移根据摩擦改界面切向力 IF( DGAM( 1 ) 0 ) THEN TAU( 1 ) = -F * PRESS DDTDDP( 1 ) = -F ELSE TAU( 1 ) = F * PRESS DDTDDP( 1 ) = F END IF RETURN END IF r