本帖最后由 senja 于 2012-4-11 20:38 编辑
我是大笨蛋 发表于 2012-4-11 16:32
先谢谢斑竹了。不过我试了一下,每次运算都提示我数组CC没有定义,不知道该怎么解决。恳请斑竹帮忙呀。
...
之前没写好,应该再重启动之后再读入cc,就没问题了。另外发现命令D处有点小小问题,都改过来了。运行没有问题。
finish
/clear
/PREP7
ET,1,VISCO106
ET,2,PLANE42
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,3.2e11
MPDATA,PRXY,1,,0.38
TB,ANAN,1,,,0
TBMODIF,1,1,5.633e7
TBMODIF,2,1,10830
TBMODIF,3,1,1.49e7
TBMODIF,4,1,11
TBMODIF,5,1,0.303
TBMODIF,6,1,2.6408e9
TBMODIF,7,1,8.0415e7
TBMODIF,8,1,0.0231
TBMODIF,9,1,1.34
mpdata,ex,2,,2e9
mpdata,prxy,2,,0.3
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,8400
mpdata,dens,2,,7800
rect,0,112,0,0.2
rect,0,160,0.2,12.8
rect,0,112,0,-0.2
rect,0,160,-0.2,-12.8
aglue,1,2
aglue,3,4
CYL4,152.5,6.5,3
CYL4,152.5,-6.5,3
ASBA,5,4
ASBA,2,6
TYPE,1
MAT,1
REAL,
ESYS,0
SECNUM,
LESIZE,2,0.1, , , , , , ,1
LESIZE,10,0.1, , , , , , ,1
LESIZE,4,0.1, , , , , , ,1
LESIZE,12,0.1, , , , , , ,1
LESIZE,1,0.2, , , , , , ,1
LESIZE,11,0.2, , , , , , ,1
LESIZE,9,0.2, , , , , , ,1
LESIZE,3,0.2, , , , , , ,1
TYPE,2
MAT,2
REAL,
ESYS,0
SECNUM,
LESIZE,17,1, , , , , , ,1
LESIZE,6,1, , , , , , ,1
LESIZE,7,1, , , , , , ,1
LESIZE,18,1, , , , , , ,1
LESIZE,5,1, , , , , , ,11
LESIZE,13,1, , , , , , ,1
LESIZE,14,1, , , , , , ,1
LESIZE,15,1, , , , , , ,1
LESIZE,16,1, , , , , , ,1
LESIZE,19,1, , , , , , ,1
LESIZE,20,1, , , , , , ,1
LESIZE,21,1, , , , , , ,1
LESIZE,22,1, , , , , , ,1
LESIZE,23,1, , , , , , ,1
LESIZE,24,1, , , , , , ,1
LESIZE,8,1, , , , , , ,1
amesh,all
finish
/SOLU
ANTYPE,4
TRNOPT,FULL
LUMPM,0
ANTYPE,4
NLGEOM,1
DELTIM,0.25,0,0
OUTRES,ERASE
OUTRES,ALL,ALL
KBC,0
RESCONTRL,DEFINE,ALL,ALL,9
TIME,1
DL,18,,ALL,
DL,4,,ALL,
DL,12,,ALL,
DL,5,,ALL,
D,3,ALL,,,560,1
D,2249,ALL,,,2806,1
BFA,1,TEMP,300
BFA,3,TEMP,300
TIME,0.25
AUTOTS,-1
DELTIM,0.2, , ,1
KBC,0
TSRES,ERASE
allsel
solve
TIME,1
AUTOTS,-1
DELTIM,0.2, , ,1
KBC,0
TSRES,ERASE
nsel,all
DL,15, ,UY,3
DL,16, ,UY,3
DL,23, ,UY,-3
DL,24, ,UY,-3
allsel
solve
/post1
*get,u11,node,561,u,x
*get,u12,node,2248,u,x
*get,w11,node,561,u,y
*get,w12,node,2248,u,y
*get,x1,node,560,rf,fx
*get,z1,node,560,rf,fy
E=1/2*(x1*(u11-u12)+z1*(w11-w12))
G=E/0.2
*status
*dim,kk,array,1,2
kk(1,1)=1
kk(1,2)=1
*cfopen,cycle,txt
*vwrite,kk(1,1),kk(1,2)
(2f5.2)
*cfclos
*do,i,1,5
*if,G,lt,1e30,then
/solu
antype,trans,rest,2,4
*dim,cc,array,1,2
*vread,cc,cycle,txt,,jik,2,1
(2f5.2)
ddele,560-cc(1,1),uy,560,1
ddele,560-cc(1,1),ux,560,1
ddele,2249+cc(1,1),ux,2249,1
ddele,2249+cc(1,1),uy,2249,1
allsel
solve
/post1
*get,u1,node,560-cc(1,1),u,x
*get,u2,node,2249+cc(1,1),u,x
*get,w1,node,560-cc(1,1),u,y
*get,w2,node,2249+cc(1,1),u,y
*get,x1,node,559-cc(1,1),rf,fx
*get,z1,node,559-cc(1,1),rf,fy
E=1/2*(x1*(u1-u2)+z1*(w1-w2))
G=-E/0.2
*status
cc(1,1)=cc(1,1)+1
*cfopen,cycle,txt
*vwrite,cc(1,1),cc(1,2)
(2f5.2)
*cfclos
*enddo
|