!让MSC.MARC根据总应变杀死单元
!by 陆新征,清华大学

subroutine uactive(m,n,mode,irststr,irststn,inc,time,timinc)
include '../common/implicit'
include "../common/ELMCOM"
dimension m(2)
integer ElemNo, IntPNo, KC, ICODE
real*8 var1(3,8, 6), Strainmax, Strainmin
var1=0.
ElemNo=M(1)


do IntPNo=1,8
KC=1
ICODE=301 !total strain
CALL ELMVAR (ICODE,ElemNo,IntPNo,KC,VAR1(1,IntPNo, :)) ! 得到总应变
end do

Strainmax=maxval((var1(3,:,:)))

if(Strainmax>5.E-2) then !Kill the element
MODE=-1; IRSTSTR=1; IRSTSTN=1;
end if

RETURN
END

个人信息
研究工作
实际工程
论文工作
教学工作
资料下载
专题
其他

我们的实验室

抗倒塌专业委员会