将索网划分为link10单元，约束不够啊！

guxing2001 · 发表于 2006-5-15 15:18

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

在索网固有频率分析时将其划分为link10单元，但是前80多阶都是为零。静力分析时，提示说节点位移过大，不能进行。求大家帮帮忙，怎么解决。如何给link10单元加上足够的约束！

genghp · 发表于 2006-5-15 18:11

呵呵！建模的问题

guxing2001 · 发表于 2006-5-15 20:11

怎么解决阿，我都困了两个星期了，高手们给指点一二吧！

AaronSpark · 发表于 2006-5-17 07:51

!!!!!!说明 mj1=0.00014 !!!!!单索的截面积 mj2=0.00056 !!!!!凸索上的截面积 mj3=0.00056 !!!!!凹索上的截面积 mj4=0.00084 !!!!!凸索下的截面积 mj5=0.00028 !!!!!凹索下的截面积 tugao=1.8 !!!!!凸形失高 aogao=1.5 !!!!!凹形失高 longth=30.55 !单元长度 width=23.4 !单元宽度 /prep7 /title,element of cable-edged net structure ET,1,LINK10 !定义索单元 ET,2,LINK10 ET,3,LINK10 ET,4,LINK10 ET,5,LINK10 ET,6,PIPE16 !定义撑杆单元 R,1,mj1,0.9 !!定义单索的实常数 MP,EX,1,1.9e8 MP,NUXY,1,0.3 r,2,mj2,0.9 !凸的上索 MP,EX,2,1.9e8 MP,NUXY,2,0.3 r,3,mj3,0.9 !凹的上索 MP,EX,3,1.9e8 MP,NUXY,3,0.3 r,4,mj4,0.9 !凸的下索 MP,EX,4,1.9e8 MP,NUXY,4,0.3 r,5,mianj5,0.9 !凹的下索 MP,EX,5,1.9e8 MP,NUXY,5,0.3 r,6,0.081,0.005 !撑杆实常数 MP,EX,6,1.9e8 MP,NUXY,6,0.3 a1=tugao/16*7 a2=tugao/16*12 a3=tugao/16*15 a4=tugao b1=-aogao/9*5 b2=-aogao/9*8 b3=-aogao !定义节点 n,1, n,9,longth fill,1,9 ngen,7,9,1,9,1,,width/6 NPLOT !显示节点 !!!内索单索 type,1 real,1 mat,1 *do,i,1,8 *do,j,2,6 e,(j-1)*9+i,(j-1)*9+i+1 *enddo *enddo *do,i,2,8 *do,j,1,6 e,(j-1)*9+i,j*9+i *enddo *enddo esel,s,type,,1 cm,neisuo,elem !!凸形索---上 type,2 real,2 mat,2 *do,i,1,8,1 *do,j,1,7,6 e,(j-1)*9+i,(j-1)*9+i+1 *enddo *enddo esel,s,type,,2 cm,tuxingsuoshang,elem !!aoxingsuo---shang type,3 real,3 mat,3 *do,i,1,9,8 *do,j,1,6 e,(j-1)*9+i,j*9+i *enddo *enddo esel,s,type,,3 cm,aoxingsuoshang,elem ALLSEL,ALL GPLOT /solu !施加位移约束凸 d,2,ux,0 d,2,uy,0 d,2,uz,a1 d,8,ux,0 d,8,uy,0 d,8,uz,a1 d,56,ux,0 d,56,uy,0 d,56,uz,a1 d,62,ux,0 d,62,uy,0 d,62,uz,a1 d,3,ux,0 d,3,uy,0 d,3,uz,a2 d,7,ux,0 d,7,uy,0 d,7,uz,a2 d,57,ux,0 d,57,uy,0 d,57,uz,a2 d,61,ux,0 d,61,uy,0 d,61,uz,a2 d,4,ux,0 d,4,uy,0 d,4,uz,a3 d,6,ux,0 d,6,uy,0 d,6,uz,a3 d,58,ux,0 d,58,uy,0 d,58,uz,a3 d,60,ux,0 d,60,uy,0 d,60,uz,a3 d,5,ux,0 d,5,uy,0 d,5,uz,a4 d,59,ux,0 d,59,uy,0 d,59,uz,a4 !施加位移约束凹 d,10,ux,0 d,10,uy,0 d,10,uz,b1 d,18,ux,0 d,18,uy,0 d,18,uz,b1 d,46,ux,0 d,46,uy,0 d,46,uz,b1 d,54,ux,0 d,54,uy,0 d,54,uz,b1 d,19,ux,0 d,19,uy,0 d,19,uz,b2 d,37,ux,0 d,37,uy,0 d,37,uz,b2 d,27,ux,0 d,27,uy,0 d,27,uz,b2 d,45,ux,0 d,45,uy,0 d,45,uz,b2 d,28,ux,0 d,28,uy,0 d,28,uz,b3 d,36,ux,0 d,36,uy,0 d,36,uz,b3 !施加位移约束 4个角点 d,1,ux,0 d,1,uy,0 d,1,uz,0 d,9,ux,0 d,9,uy,0 d,9,uz,0 d,55,ux,0 d,55,uy,0 d,55,uz,0 d,63,ux,0 d,63,uy,0 d,63,uz,0 !设定时间步 TIME,1 AUTOTS,0 NSUBST,30, , ,1 KBC,0 !输出每个子步的结果 OUTPR,BASIC,ALL, OUTRES,ALL,ALL, !设定大变形，应力刚化 NLGEOM,1 SSTIF,ON finish /solu /STATUS,SOLU SOLVE EPLOT finish /prep7 UPGEOM,1,LAST,LAST,file,rst,
!施加位移约束 allsel,all nsel,all ddele,all,all d,2,ux,0 d,2,uy,0 d,2,uz,0 d,8,ux,0 d,8,uy,0 d,8,uz,0 d,56,ux,0 d,56,uy,0 d,56,uz,0 d,62,ux,0 d,62,uy,0 d,62,uz,0 
d,3,ux,0 d,3,uy,0 d,3,uz,0 d,7,ux,0 d,7,uy,0 d,7,uz,0 d,57,ux,0 d,57,uy,0 d,57,uz,0 d,61,ux,0 d,61,uy,0 d,61,uz,0 d,4,ux,0 d,4,uy,0 d,4,uz,0 d,6,ux,0 d,6,uy,0 d,6,uz,0 d,58,ux,0 d,58,uy,0 d,58,uz,0 d,60,ux,0 d,60,uy,0 d,60,uz,0 d,5,ux,0 d,5,uy,0 d,5,uz,0 d,59,ux,0 d,59,uy,0 d,59,uz,0 !施加位移约束凹 d,10,ux,0 d,10,uy,0 d,10,uz,0 d,18,ux,0 d,18,uy,0 d,18,uz,0 d,46,ux,0 d,46,uy,0 d,46,uz,0 d,54,ux,0 d,54,uy,0 d,54,uz,0 d,19,ux,0 d,19,uy,0 d,19,uz,0 d,37,ux,0 d,37,uy,0 d,37,uz,0 d,27,ux,0 d,27,uy,0 d,27,uz,0 d,45,ux,0 d,45,uy,0 d,45,uz,0 d,28,ux,0 d,28,uy,0 d,28,uz,0 d,36,ux,0 d,36,uy,0 d,36,uz,0 !施加位移约束 4个角点 d,1,ux,0 d,1,uy,0 d,1,uz,0 d,9,ux,0 d,9,uy,0 d,9,uz,0 d,55,ux,0 d,55,uy,0 d,55,uz,0 d,63,ux,0 d,63,uy,0 d,63,uz,0 TIME,1 AUTOTS,0 NSUBST,30, , ,1 KBC,0 !输出每个子步的结果 OUTPR,BASIC,ALL, OUTRES,ALL,ALL, !设定大变形，应力刚化 NLGEOM,1 SSTIF,ON FINISH /SOLU /STATUS,SOLU SOLVE EPLOT finish /prep7 UPGEOM,1,LAST,LAST, !!!第一次循环更新坐标，释放不平衡应力 n=2 *do,i,1,n /solu antype,static NLGEOM,1 SSTIF,ON /solu solve finish /prep7 UPGEOM,1,LAST,LAST,file,rst, *enddo finish !!!施加真实应力 /prep7 R,1,mj1,0.02 !!定义单索的实常数 r,2,mj2,0.012 !凸的上索 r,3,mj3,0.01 !凹的上索 /solu /STATUS,SOLU NLGEOM,1 SSTIF,ON SOLVE FINISH !!!第二次循环更新坐标，释放不平衡应力 n=2 *do,i,1,n /solu antype,static NLGEOM,1 SSTIF,ON /solu solve finish /prep7 UPGEOM,1,LAST,LAST,file,rst, *enddo !!!还原索的真实弹性模量 /PREP7 R,1,mj1,0.02 !!定义单索的实常数 MP,EX,1,1.9e11 MP,NUXY,1,0.3 r,2,mj2,0.12 !凸的上索 MP,EX,2,1.9e11 MP,NUXY,2,0.3 r,3,mj3,0.01 !凹的上索 MP,EX,3,1.9e11 MP,NUXY,3,0.3 /SOL /STATUS,SOLU NLGEOM,1 SSTIF,ON SOLVE FINISH !!!第三次循环更新坐标，释放不平衡应力,找形结束 n=2 *do,i,1,n /solu antype,static NLGEOM,1 SSTIF,ON /solu solve finish /prep7 UPGEOM,1,LAST,LAST,file,rst, *enddo finish !!!!加索桁架 /prep7 !!!!凸桁架下索 wpoffs,0,0,-0.50 csys,4 !激活已经定义的坐标 nsym,z,1000,1,9,1 !映射节点 nsym,z,1000,55,63,1 type,4 real,4 mat,4 *do,i,1001,1008,1 e,i,i+1 *enddo *do,i,1055,1062 e,i,i+1 *enddo esel,s,type,,4 cm,tusuoxia,elem !!!!!!!凹桁架下索 wpoffs,0,0,b3+0.5 csys,4 !激活已定义的坐标 nsym,z,2000,1,55,9 !映射节点 nsym,z,2000,9,63,9 type,5 real,5 mat,5 *do,i,2001,2046,9 e,i,i+9 *enddo *do,i,2009,2054,9 e,i,i+9 *enddo nummrg,all,all esel,s,type,,5 cm,aosuoxia,elem !!!!连接撑杆 type,6 real,6 mat,6 !!!凸－－撑杆 *do,i,2,8,1 e,i,i+1000 *enddo *do,i,56,62,1 e,i,i+1000 *enddo !!!!凹－－撑杆 *do,i,10,19,9 e,i,i+2000 *enddo *do,i,37,46,9 e,i,i+2000 *enddo *do,i,18,27,9 e,i,i+2000 *enddo *do,i,45,54,9 e,i,i+2000 *enddo esel,s,type,,6 cm,chenggan,elem nummrg,all,all !合并凹的重合点 wpoffs,0,0,1.5 csys,4 !激活已定义的坐标 csys,0 finish /prep7 ddele,all,all !!!!还原真实弹性模量 R,1,mj1,0.002 !!定义单索的实常数 MP,EX,1,1.9e11 MP,NUXY,1,0.3 mp,dens,1,7850 r,2,mj2,0.0012 !凸的上索 MP,EX,2,1.9e11 MP,NUXY,2,0.3 mp,dens,2,7850 r,3,mj3,0.001 !凹的上索 MP,EX,3,1.9e11 MP,NUXY,3,0.3 mp,dens,3,7850 r,4,mj4,0.0018 !凸的下索 MP,EX,4,1.9e11 MP,NUXY,4,0.3 mp,dens,4,7850 r,5,mj5,0.002 !凹的下索 MP,EX,5,1.9e11 MP,NUXY,5,0.3 mp,dens,5,7850 r,6,0.081,0.005 !撑杆实常数 MP,EX,6,1.9e11 MP,NUXY,6,0.3 mp,dens,6,7850 !!!!约束 allsel,all nsel,all ddele,all,all !!!重新约束 *do,i,2,8,1 d,i,uy,0 *enddo *do,i,56,62,1 d,i,uy,0 *enddo *do,i,10,46,9 d,i,ux,0 *enddo *do,i,18,54,9 d,i,ux,0 *enddo
*do,i,1010,1046,9 d,i,ux,0 *enddo *do,i,1018,1054,9 d,i,ux,0 *enddo
*do,i,1,2001,1000 d,i,ux,0 d,i,uy,0 d,i,uz,0 *enddo *do,i,9,2009,1000 d,i,ux,0 d,i,uy,0 d,i,uz,0 *enddo *do,i,55,2055,1000 d,i,ux,0 d,i,uy,0 d,i,uz,0 *enddo *do,i,63,2063,1000 d,i,ux,0 d,i,uy,0 d,i,uz,0 *enddo
NSUBST,20, , ,1 AUTOTS,1 !输出每个子步的结果 OUTPR,BASIC,ALL, OUTRES,ALL,ALL, !设定大变形，应力刚化 NLGEOM,1 SSTIF,ON finish solve finish /solu /STATUS,SOLU d,all,all ANTYPE,0 NLGEOM,off PSTRES,ON FINISH /SOLU SOLVE finish /solu allsel,all ddele,all *do,i,2,8,1 d,i,uy,0 *enddo *do,i,56,62,1 d,i,uy,0 *enddo *do,i,10,46,9 d,i,ux,0 *enddo *do,i,18,54,9 d,i,ux,0 *enddo
*do,i,1010,1046,9 d,i,ux,0 *enddo *do,i,1018,1054,9 d,i,ux,0 *enddo
*do,i,1,2001,1000 d,i,ux,0 d,i,uy,0 d,i,uz,0 *enddo *do,i,9,2009,1000 d,i,ux,0 d,i,uy,0 d,i,uz,0 *enddo *do,i,55,2055,1000 d,i,ux,0 d,i,uy,0 d,i,uz,0 *enddo *do,i,63,2063,1000 d,i,ux,0 d,i,uy,0 d,i,uz,0 *enddo ANTYPE,2 !模态分析 PSTRES,ON !考虑Prestress效应 MODOPT,lanb,50,0,0, !取前50阶振型,子空间迭代法 
FINISH /SOLU SOLVE finish

账号		自动登录	找回密码
密码			我要加入

将索网划分为link10单元，约束不够啊！

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（guxing2001）将索网划分为link10单元，约束不...