byb 发表于 2007-1-26 21:35

一个施加载荷的例子,请大家帮忙修改

编写的如下命令流,却无法收敛,请版主,大侠们帮忙修改修改~~~~谢谢~~~


finsih
/clear
BREADTH=341
LENGTH=-1000
B1=83
DEGREE=8
/prep7
BLC4,0,0,BREADTH,LENGTH
MP,DENS,1,7800
MP,EX,1,210e9   
MP,NUXY,1,0.3   
TB,BKIN,1,,,,   
TBDAT,1,235e6   
TBDAT,2,105e7
ET,1,SHELL181   
KEYOPT,1,1,0
KEYOPT,1,3,0
KEYOPT,1,8,0
KEYOPT,1,9,0
KEYOPT,1,10,0   
R,1,8, , , , , ,
!定义材料
LSEL,S,LOC,X,0
LSEL,A,LOC,X,BREADTH
LESIZE,ALL,20
ALLSEL,ALL
LSEL,S,LOC,z,0
LSEL,A,LOC,z,LENGTH
LESIZE,ALL,20
AMESH,ALL
!划分网格
LSEL,S,LOC,X,0
DL,ALL,,SYMM
!施加对称载荷
ALLSEL,ALL   
wpoff,B1
NSEL,S,LOC,Y,0
NSEL,R,LOC,X,0,B1
D,ALL,,,,,,ALL
!在全局坐标系下0~B1段施加约束
ALLSEL,ALL
NSEL,S,LOC,Y,0
NSEL,R,LOC,X,B1,BREADTH
CSWPLA,11,0,1,1,
CSYS,11,
*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNCCSYS   
*SET,_FNCNAME,'XX'
*SET,_FNCCSYS,11
! /INPUT,XX.func
*DIM,%_FNCNAME%,TABLE,6,5,2,,,,%_FNCCSYS%   
!   
! Begin of equation: {X}
*SET,%_FNCNAME%(0,0,1), 0.0, -999   
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), 0.0
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, 99, 0, 1, 2, 0, 0
*SET,%_FNCNAME%(0,2,1),   0
*SET,%_FNCNAME%(0,3,1),   1
*SET,%_FNCNAME%(0,4,1),   0
*SET,%_FNCNAME%(0,5,1),   0
! End of equation: {X}
!   
! Begin of equation: (1-cos(DEGREE))*{X}
*SET,%_FNCNAME%(0,0,2), 1000, -999
*SET,%_FNCNAME%(2,0,2), 0.0
*SET,%_FNCNAME%(3,0,2), 0.0
*SET,%_FNCNAME%(4,0,2), 0.0
*SET,%_FNCNAME%(5,0,2), 0.0
*SET,%_FNCNAME%(6,0,2), 0.0
*SET,%_FNCNAME%(0,1,2), 1.0, -1, 10, 8, 0, 0, 0
*SET,%_FNCNAME%(0,2,2), 0.0, -2, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,3,2),   1, -3, 0, 1, -2, 2, -1
*SET,%_FNCNAME%(0,4,2), 0.0, -1, 0, 1, -3, 3, 2
*SET,%_FNCNAME%(0,5,2), 0.0, 99, 0, 1, -1, 0, 0
! End of equation: (1-cos(DEGREE))*{X}   
!-->
!建立(1-cos(DEGREE))*{X} 数组
*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNCCSYS   
*SET,_FNCNAME,'ZZ'
*SET,_FNCCSYS,11
! /INPUT,ZZ.func
*DIM,%_FNCNAME%,TABLE,6,3,2,,,,%_FNCCSYS%   
!   
! Begin of equation: {X}
*SET,%_FNCNAME%(0,0,1), 0.0, -999   
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), 0.0
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, 99, 0, 1, 2, 0, 0
*SET,%_FNCNAME%(0,2,1),   0
*SET,%_FNCNAME%(0,3,1),   1
! End of equation: {X}
!   
! Begin of equation: sin(DEGREE)*{X}
*SET,%_FNCNAME%(0,0,2), 1000, -999
*SET,%_FNCNAME%(2,0,2), 0.0
*SET,%_FNCNAME%(3,0,2), 0.0
*SET,%_FNCNAME%(4,0,2), 0.0
*SET,%_FNCNAME%(5,0,2), 0.0
*SET,%_FNCNAME%(6,0,2), 0.0
*SET,%_FNCNAME%(0,1,2), 1.0, -1, 9, 8, 0, 0, 0
*SET,%_FNCNAME%(0,2,2), 0.0, -2, 0, 1, -1, 3, 2
*SET,%_FNCNAME%(0,3,2),   1, 99, 0, 1, -2, 0, 0
! End of equation: sin(DEGREE)*{X}   
!-->
!建立sin(DEGREE)*{X}数组
D,ALL, , %XX% , , , ,UX, , , , ,   
D,ALL, , %ZZ% , , , ,UZ, , , , ,
!在B1~BREADTH段施加X方向的(1-cos(DEGREE))*{X} 位移载荷,施加Z方向的sin(DEGREE)*{X} 位移载荷
ALLSEL,ALL
/REPLOT
FINISH
/SOLU
ANTYPE,0
NLGEOM,1
NSUBST,30,50,0
AUTOTS,1
LNSRCH,1
PRED,ON,,ON
TIME,1
!设置求解非线性选项
SOLVE
页: [1]
查看完整版本: 一个施加载荷的例子,请大家帮忙修改