【讨论】断裂（VM１４３）特请chenhu 和curtain二位大虾帮我

fshj

我是一致菜菜鸟，我做VM１４３时有个问题要请教大家：
　　那段算J积分的宏该怎样调用，我是这样做的，可行不通，
/CLEAR, NOSTART ! CLEAR DATABASE FOR 2ND SOLUTION
/PREP7
SMRT,OFF
/TITLE, VM143, FRACTURE MECHANIC STRESS INTENSITY - CRACK IN A FINITE WIDTH PLATE
/COM,   ****** CRACK IN 2-DIMENSIONS USING 2-D PLANE82 ******
ET,1,PLANE82,,,2              ! PLANE82 (PLANE STRAIN)
MP,EX,1,30E6
MP,NUXY,1,0.3
K,1                           ! DEFINE KEYPOINTS AND LINE SEGMENTS   
K,2,4
K,3,4,5
K,4,-1,5
K,5,-1
L,1,2
L,2,3
LESIZE,2,,,4
L,3,4
LESIZE,3,,,4
L,4,5,
LESIZE,4,,,6,.2
L,5,1
ESIZE,,5
KSCON,1,.15,1,8               ! DEFINE CRACK TIP ELEMENT SIZE
AL,1,2,3,4,5
DL,1,1,SYMM                   ! APPLY SOLID MODEL BOUNDARY CONDITIONS
DL,4,1,SYMM
SFL,3,PRES,-.5641895
AMESH,1
OUTPR,ALL     
FINISH
/COM         
/OUTPUT,SCRATCH
/SOLU     
SOLVE
FINISH
/OUTPUT
/POST1
ETABLE,SENE,SENE              ! RETRIEVE STRAIN ENERGY PER ELEMENT
ETABLE,VOLU,VOLU              ! RETRIEVE VOLUME PER ELEMENT
  
（下面这段宏该怎么用呢，前面两行和最后一行还要吗，）
!*CREATE,JIN1
!STINFC                        ! DATA BLOCK NAME
SEXP,W,SENE,VOLU,1,-1         ! CALCULATE STRAIN ENERGY DENSITY
PATH,JINT,4,50,48             ! DEFINE PATH WITH NAME = "JINT"
PPATH,1,ARG1                  ! DEFINE PATH POINTS BY NODE
PPATH,2,ARG2
PPATH,3,ARG3
PPATH,4,ARG4
PDEF,W,ETAB,W                 ! PUT STRAIN ENERGY DENSITY ON THE PATH
PCALC,INTG,J,W,YG             ! INTEGRATE ENERGY W.R.T. GLOBAL Y
*GET,JA,PATH,,LAST,J          ! GET FINAL VALUE OF INTEGRAL FOR 1ST TERM OF J
PDEF,CLEAR                    ! CLEAR OLD PATH VARIABLES
PVECT,NORM,NX,NY,NZ           ! DEFINE THE PATH UNIT NORMAL VECTOR
PDEF,INTR,SX,SX               ! PUT STRESS SX ON THE PATH
PDEF,INTR,SY,SY               ! PUT STRESS SY ON THE PATH
PDEF,INTR,SXY,SXY             ! PUT STRESS SXY ON THE PATH
PCALC,MULT,TX,SX,NX           ! CALCULATE TRACTION TX
PCALC,MULT,C1,SXY,NY          !   TX = SX*NX + SXY*NY
PCALC,ADD,TX,TX,C1
PCALC,MULT,TY,SXY,NX          ! CALCULATE TRACTION TY
PCALC,MULT,C1,SY,NY           !   TY = SXY*NX + SY*NY
PCALC,ADD,TY,TY,C1
*GET,DX,PATH,,LAST,S          ! DEFINE PATH SHIFT AS 1% OF PATH LENGTH
DX=DX/100
PCALC,ADD,XG,XG,,,,-DX/2      ! SHIFT PATH FROM X TO X-DX/2 (GLOBAL X DIR.)
PDEF,INTR,UX1,UX              ! DEFINE UX AT X-DX
PDEF,INTR,UY1,UY              ! DEFINE UY AT X-DX
PCALC,ADD,XG,XG,,,,DX         ! SHIFT PATH FROM X-DX/2 TO X+DX/2
PDEF,INTR,UX2,UX              ! DEFINE UX AT X+DX
PDEF,INTR,UY2,UY              ! DEFINE UY AT X+DX
PCALC,ADD,XG,XG,,,,-DX/2      ! SHIFT PATH BACK TO ORIGINAL POSITION
C=(1/DX)
PCALC,ADD,C1,UX2,UX1,C,-C     ! CALCULATE DERIVATIVE DUX/DX
PCALC,ADD,C2,UY2,UY1,C,-C     ! CALCULATE DERIVATIVE DUY/DX
PCALC,MULT,C1,TX,C1           ! DEFINE INTEGRAND
PCALC,MULT,C2,TY,C2           !   = TX*DUX/DX + TY*DUY/DX
PCALC,ADD,C1,C1,C2
PCALC,INTG,J,C1,S             ! FORM SECOND INTEGRAL (W.R.T. PATH LENGTH S)
*GET,JB,PATH,,LAST,J          ! GET FINAL VALUE OF INTEGRAL FOR 2ND TERM OF J
JINT=2*(JA-JB)                ! ADD BOTH TERMS AND DOUBLE FOR HALF MODELS
PDEF,CLEAR                    ! CLEAR PATH VARIABLES
!*END
万分感激！！！

yujb025

将那段文字拷到记事本中并保存为*.mac文件，把文件放到你的工作目录中，然后在ansys中运行文件名就可以了。
好象是抢了二位大虾的生意了，真是不好意思！

fshj

我明白了，anyway,谢谢上面兄弟！！