2005-5-13
根据你的建议,修改原程序如下:
记为1程序
!******************************
!单元生死,分层计算受力场
! J积分法计算KI:
! 读入crack-J-path, 定义路径
! 读入crack-J-cacul, 完成KI计算
!*******************************
N=3 !层数
M=15 !层单元数
E=30e6 !弹性常数
v=0.3
/SOL
ANTYPE,4
TRNOPT,FULL
LUMPM,0
/SOLU
ESEL,ALL
*DO,J,1,N,1
*IF,J,EQ,1,THEN
ANTYPE,STATIC,NEW
*ELSE
ANTYPE,,REST
PARRES, NEW, PARAFILE, DAT !读取数据文件
*ENDIF
TIME,J
NSUBST,20,100,1
ESEL,S,ELEM,,CENG1(1,J)
*DO,I,1,M,1
ESEL,A,ELEM,,CENG1(I,J)
*ENDDO
CM,C1, ELEM
ESEL,ALL
ESEL,S,ELEM,,CENG2(1,J)
*DO,I,1,M,1
ESEL,A,ELEM,,CENG2(I,J)
*ENDDO
CM,C2, ELEM !层选择
/SOLU
CMSEL,,C1
SFE, C1, 3, PRES, ,-0.5641895, , ,
CMSEL, , C2
SFE, C2, 3, PRES, ,-0.5641895, , , !!层加载
ESEL, , LIVE
/PSF, PRES, NORM, 2, 0, 1 !
/PBF, DEFA, ,1 !
/PIC, DEFA, ,1 !显示载荷标量和边界条件
/PBC, ALL, ,1
/REP
ESEL, ALL
/SOLU
SOLVE
J=J+1
PARSAV, ALL, PARAFILE, DAT !存储数据文件
/INPUT,';CRACK-J-PATH';,';LOG';,';G:\ANSYS WORK FILE\';,, 0 !
/INPUT,';CRACK-J-CACUL';,';LOG';,';G:\ANSYS WORK FILE\';,, 0 !读入文件,进入/post1计算
VALUEKI (1,J)=KI
VALUEKI (2,J)=(KI-1.0249)/1.0249
/SOLU
SFEDELE, C1, ALL, PRES !层卸载
SFEDELE, C2, ALL, PRES
EKILL, C1 !杀死层单元
EKILL, C2
*ENDDO
!************************************************
|