找形时采用很小的弹性模量.docVIP

/prep7 ET,1,LINK10 !找形时采用很小的弹性模量 R,1,0.001468,0.9, MP,EX,1,0.60551e9 MP,NUXY,1,0.3 !定义节点 *do,i,1,5 *do,j,1,2*i-1 k=(i-1)*(i-1)+j n,k,(j-i)*9.15,-36.6+(i-1)*9.15,0 *enddo *enddo *do,i,1,4 *do,j,1,2*i-1 k=(i-1)*(i-1)+j+25 n,k,(j-i)*9.15,36.6-(i-1)*9.15,0 *enddo *enddo !定义单元 *do,i,1,4 *do,j,1,2*i k=i*i+j e,k,k+1 *enddo *enddo *do,i,1,3 *do,j,1,2*i k=i*i+j+25 e,k,k+1 *enddo *enddo *do,i,1,4 *do,j,1,2*i-1 k=(i-1)*(i-1)+j e,k,k+2*i *enddo *enddo *do,i,1,3 *do,j,1,2*i-1 k=(i-1)*(i-1)+j+25 e,k,k+2*i *enddo *enddo *do,i,18,24 e,i,i+17 *enddo !施加位移约束 d,1,ux,0 d,1,uy,0 d,1,uz,3.66 d,2,ux,0 d,2,uy,0 d,2,uz,1.83 d,4,ux,0 d,4,uy,0 d,4,uz,1.83 d,5,ux,0 d,5,uy,0 d,5,uz,0 d,9,ux,0 d,9,uy,0 d,9,uz,0 d,10,ux,0 d,10,uy,0 d,10,uz,-1.83 d,16,ux,0 d,16,uy,0 d,16,uz,-1.83 d,17,ux,0 d,17,uy,0 d,17,uz,-3.66 d,25,ux,0 d,25,uy,0 d,25,uz,-3.66 d,35,ux,0 d,35,uy,0 d,35,uz,-1.83 d,41,ux,0 d,41,uy,0 d,41,uz,-1.83 d,30,ux,0 d,30,uy,0 d,30,uz,0 d,34,ux,0 d,34,uy,0 d,34,uz,0 d,27,ux,0 d,27,uy,0 d,27,uz,1.83 d,29,ux,0 d,29,uy,0 d,29,uz,1.83 d,26,ux,0 d,26,uy,0 d,26,uz,3.66 !设定时间步 TIME,1 AUTOTS,0 NSUBST,20, , ,1 KBC,0 !输出每个子步的结果 OUTPR,BASIC,ALL, OUTRES,ALL,ALL, !设定大变形，应力刚化 NLGEOM,1 SSTIF,ON FINISH !求截 /SOLU /STATUS,SOLU SOLVE FINISH !后处理 /POST1 PRNSOL,DOF, PRESOL,SMISC,1 PLESOL,SMISC,1 PLNSOL,U,X,0,1 finish /PREP7 !更新几何形状 !a.rst为计算结果文件名，最后一个为目录 !这两个参数应根据你的计算情况定 UPGEOM,1,LAST,LAST,a,rst,E:\JZD\1\ !弹性模量恢复为真值 R,1,0.001468,0.0027248, MP,EX,1,2.0e11 !重新施加位移约束 d,1,uz,0 d,2,uz,0 d,4,uz,0 d,5,uz,0 d,9,uz,0 d,10,uz,0 d,16,uz,0 d,17,uz,0 d,25,uz,0 d,35,uz,0 d,41,uz,0 d,30,uz,0 d,34,uz,0 d,27,uz,0 d,29,uz,0 d,26,uz,0 !求截 /SOLU /STATUS,SOLU SOLVE FINISH /POST1 !观察施加荷载之前结构的位移 !可见位移基本为零,预张力基本不变 PLNSOL,U,Z,0,1 PLESOL,SMISC,1 FINISH !写荷载工况文件01 /PREP7 LSWRITE,01, !施加节点荷载 *do,i,1,41 f,i,fz,-167445 *enddo !设置第2荷载步 TIME,2 AUTOTS,0 NSUBST,20, , ,1 KBC,0 LSWR