graininear 发表于 2013-6-4 23:59
除了FKN,还有单元的各种关键字定义,都可能会产生影响的
你好版主,我把命令流贴上,你能帮我看看哪里有问题吗
/PREP7
ET,1,solid45 !3-D
MP,EX,1,2.7e10
MP,nuxy,1,0.2
MP,dens,1,2500
BLC4,,,0.25,0.25,5
BLOCK,0,0.25,0,0.25,5,10
esize,0.05 !将边长0.25分为5等分
vmesh,all
ET,2,TARGE170 !3-D目标单元(指压杆)
ET,3,CONTA173 !3-D试件上节点表面接触单元
r,2 !建立一个空的实常数(硬接触),以用于建立接触
r,2,,,1.0,0.1,0
rmore,,,1.0e20,0.0,1.0
keyopt,3,4,0
keyopt,3,5,0
keyopt,3,7,0
keyopt,3,8,0
keyopt,3,9,1
keyopt,3,10,1
keyopt,3,11,0
keyopt,3,12,4
keyopt,3,2,0
ASEL,S,,,7 !选择目标表面 A7
NSLA,S,1 !选择面上的所有节点,必须是这样
ESLN,S,0
TYPE,2
real,2
ESURF,all
ESEL,S,TYPE,,2 !查看单元类型为2的单元,目的是确认是否为一个面
eplot
ASEL,S,,,2 !选择目标表面 A2
NSLA,S,1 !选择面上的所有节点,必须是这样
ESLN,S,0
TYPE,3
real,2
ESURF,all
ESEL,S,TYPE,,3 !查看单元类型为3的单元,目的是确认是否为一个面
eplot
allsel
*dim,force,table,5,1,1 !5行1列
*tread,force,'book1','txt','',1 !跳过第一行,因为文件force.txt的第一行为标识符号
*vplot,force(1,0),force(1,1)
/SOLU
ANTYPE,TRANS
TRNOPT,FULL
!NLGEOM,1 !打开“大变形”开关,此时是进行非线性分析
NLGEOM,0 !关闭“大变形”开关,此时是进行线性分析
OUTRES,ALL,ALL
AUTOTS,0 !关闭自动时间步长 由于下面*do命令有时间增量
KBC,0 !斜坡荷载 若采用stepping load结果也是一致的。
tm_start=1e-12 !开始时间必须大于0
tm_end=40e-3 !瞬态结束时间
tm_incr=2e-3 !时间增量
/PSF,PRES,NORM,2,1,1
NSEL,ALL
NSEL,S,LOC,Z,40-0.001,40+0.001
NPLOT
CM,nodes,NODE !选择桩顶节点,建立节点组元nodes
*do,tm,tm_start,tm_end,tm_incr ! tm循环从tm_start到tm_end,增量是tm_incr
time,tm !时间值
SF,nodes,PRES,force(tm) !施加 时间~荷载 并查看施加的压力方向是否正确(压力是指向单元面为正)
!最好通过节点组元来施加荷载,不能使用在area上施加的方法
allsel
solve !开始计算
*enddo |