声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1132|回复: 1

塑性土体+弹性桩+粘结 在荷载下的情况

[复制链接]
发表于 2005-10-26 20:07 | 显示全部楼层 |阅读模式

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

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

x
<DIV >!塑性土体+弹性桩+粘结 <BR>!集中荷载:216KN,瞬态分析 <BR>/com,structural <BR>fini <BR>/clear <BR>/prep7 <BR>et,1,solid45 <BR>mp,ex,1,2e7 <BR>mp,prxy,1,0.45 <BR>mp,dens,1,1180 <BR>mp,damp,1,0.05 <BR>TB,dp,1 <BR>TBDATA,1,34000,16.7,16.7 <BR>et,2,solid45 <BR>mp,ex,2,2e10 <BR>mp,prxy,2,0.3 <BR>mp,dens,2,2300 <BR>!!!!kelvin单元 <BR>ET,3,COMBIN14 <BR>R,3,3.4e6,1.7e6 <BR><BR>BLOCK,-6,6,0,3,0,-10, <BR>CYL4,0,0,0,0,0.25,180,-10 <BR>VSBV, 1, 2 <BR>CYL4,0,0,0,0,0.25,180,-10 <BR>LSTR, 18, 7 <BR>LSTR, 18, 8 <BR>numcmp,all <BR>FLST,2,2,4,ORDE,2 <BR>FITEM,2,29 <BR>FITEM,2,28 <BR>ADRAG,P51X, , , , , ,27 <BR>FITEM,2,1 <BR>FITEM,2,2 <BR>FLST,3,2,5,ORDE,2 <BR>FITEM,3,14 <BR>FITEM,3,15 <BR>VSBA,P51X,P51X <BR>numcmp,all <BR><BR>vglue,all <BR>numcmp,all <BR><BR>lesize,13,,,1 <BR>lesize,14,,,1 <BR>lesize,26,,,1 <BR>lesize,21,,,1 <BR>lesize,4,,,1 <BR>lesize,5,,,1 <BR>lesize,18,,,1 <BR>lesize,23,,,1 <BR><BR>LESIZE,1, , ,6 <BR>LESIZE,3, , ,6 <BR>LESIZE,6, , ,6 <BR>LESIZE,8, , ,6 <BR>LESIZE,19, , ,6 <BR>LESIZE,20, , ,6 <BR>LESIZE,24, , ,6 <BR>LESIZE,25, , ,6 <BR><BR>LESIZE,2, , ,24 <BR>LESIZE,7, , ,24 <BR>LESIZE,29, , ,24 <BR>LESIZE,28, , ,24 <BR><BR>LESIZE,9, , ,20 <BR>LESIZE,10, , ,20 <BR>LESIZE,11, , ,20 <BR>LESIZE,12, , ,20 <BR>LESIZE,27, , ,20 <BR>LESIZE,15, , ,20 <BR>LESIZE,17, , ,20 <BR>LESIZE,16, , ,20 <BR>LESIZE,22, , ,20 <BR><BR>LESIZE,30, , ,10,0.5 <BR>LESIZE,31, , ,10,0.5 <BR>LESIZE,32, , ,10,0.5 <BR>LESIZE,33, , ,10,0.5 <BR>LESIZE,36, , ,10,0.5 <BR>LESIZE,37, , ,10,0.5 <BR>LESIZE,34, , ,10,0.5 <BR>LESIZE,35, , ,10,0.5 <BR><BR>type,1 <BR>mat,1 <BR>vsweep,5,23,20 <BR>vsweep,6,26,19 <BR>vsweep,4,20,19 <BR>type,2 <BR>mat,2 <BR>vsweep,3,13,9 <BR>vsweep,2,5,13 <BR>vsweep,1,4,9 <BR>numcmp,all <BR>!!!!生成KELVIN单元 <BR>ASEL,S, , , 2 <BR>nsla,r,1 <BR>cm,xl,node <BR>ngen,2,9000,xl,,,-0.5,,,1 <BR>allsel,all <BR>nsla,r,1 <BR>cm,xr,node <BR>ngen,2,9000,xr,,,0.5,,,1 <BR>ALLSEL,ALL <BR>NUMCMP,NODE <BR>Real,1 <BR>*do,i,361,381 <BR>e,i,i+8208 <BR>*enddo <BR>*do,i,529,654 <BR>e,i,i+8061 <BR>*enddo <BR>*do,i,1978,1998 <BR>e,i,i+6738 <BR>*enddo <BR>*do,i,2146,2271 <BR>e,i,i+6591 <BR>*enddo <BR>numcmp,all <BR>!加底面固端约束条件 <BR>FLST,2,6,5,ORDE,6 <BR>FITEM,2,6 <BR>FITEM,2,10 <BR>FITEM,2,14 <BR>FITEM,2,17 <BR>FITEM,2,21 <BR>FITEM,2,24 <BR>DA,P51X,all, <BR>!对称面Y方向约束 <BR>FLST,2,5,5,ORDE,5 <BR>FITEM,2,1 <BR>FITEM,2,4 <BR>FITEM,2,5 <BR>FITEM,2,23 <BR>FITEM,2,26 <BR>DA,P51X,uy, <BR><BR>!Kelvin单元约束 <BR>nsel,s,loc,x,6.5 <BR>d,all,all <BR>nsel,s,loc,x,-6.5 <BR>d,all,all <BR>allsel,all <BR>gplot <BR>!施加惯性力和侧向压力 <BR>ACEL,0,0,9.8 <BR>SFGRAD,PRES, ,Z,0,-11800*0.65 <BR>SFA,3,1,PRES,0 <BR>SFA,2,1,PRES,0 <BR>sftran <BR>/pbc,all,,1 <BR>/psf,pres,norm,2,0,1 <BR>/SOLU <BR><BR>*DO,I,1,10 !做10次循环 <BR>FLST,2,1,3,ORDE,1 <BR>FITEM,2,14 <BR>FK,P51X,FX,108e3/10*i <BR>time,i <BR>LSWRITE,i, <BR>*ENDDO <BR><BR>*DO,I,1,21 <BR>FLST,2,1,3,ORDE,1 <BR>FITEM,2,14 <BR>FK,P51X,FX,108e3/10*(11-i) <BR>time,(10+i) <BR>LSWRITE,(10+i), <BR>*ENDDO <BR><BR>*DO,I,1,11 <BR>FLST,2,1,3,ORDE,1 <BR>FITEM,2,14 <BR>FK,P51X,FX,-108e3/10*(11-i) <BR>time,(31+i) <BR>LSWRITE,(31+i), <BR>*ENDDO <BR><BR>lssolve,1,42,1 <BR>fini <BR><BR>/post1 <BR>!桩中轴结点X向位移 <BR>*do,i,1,42 <BR><BR>*if,i,lt,10,then <BR>*cfopen,dis0%i%,txt <BR>*endif <BR><BR>*if,i,ge,10,then <BR>*cfopen,dis%i%,txt <BR>*endif <BR><BR>set,i <BR>*GET,zc,NODE,8568,LOC,Z <BR>*GET,ux,NODE,8568,U,X <BR>*vwrite,ux,zc <BR>(2f20.10) <BR>*do,j,1,20 <BR>*GET,zc,NODE,8547+j,LOC,Z <BR>*GET,ux,NODE,8547+j,U,X <BR>*vwrite,ux,zc <BR>(2f20.10) <BR><BR>*enddo <BR>*cfclose <BR><BR>!桩右(桩、土)结点X向位移 <BR><BR>*if,i,lt,10,then <BR>*cfopen,RTdis0%i%,txt <BR>*endif <BR><BR>*if,i,ge,10,then <BR>*cfopen,RTdis%i%,txt <BR>*endif <BR><BR>set,i <BR>*GET,zc,NODE,2145,LOC,Z <BR>*GET,ux,NODE,2145,U,X <BR>*vwrite,ux,zc <BR>(2f20.10) <BR>*do,j,1,19 <BR>*GET,zc,NODE,2114+j,LOC,Z <BR>*GET,ux,NODE,2114+j,U,X <BR>*vwrite,ux,zc <BR>(2f20.10) <BR>*enddo <BR>*GET,zc,NODE,2139,LOC,Z <BR>*GET,ux,NODE,2139,U,X <BR>*vwrite,ux,zc <BR>(2f20.10) <BR>*cfclose <BR><BR>!桩左(桩、土)结点X向位移 <BR><BR>*if,i,lt,10,then <BR>*cfopen,LTdis0%i%,txt <BR>*endif <BR><BR>*if,i,ge,10,then <BR>*cfopen,LTdis%i%,txt <BR>*endif <BR><BR>set,i <BR>*GET,zc,NODE,528,LOC,Z <BR>*GET,ux,NODE,528,U,X <BR>*vwrite,ux,zc <BR>(2f20.10) <BR>*do,j,1,19 <BR>*GET,zc,NODE,497+j,LOC,Z <BR>*GET,ux,NODE,497+j,U,X <BR>*vwrite,ux,zc <BR>(2f20.10) <BR>*enddo <BR>*GET,zc,NODE,522,LOC,Z <BR>*GET,ux,NODE,522,U,X <BR>*vwrite,ux,zc <BR>(2f20.10) <BR>*cfclose <BR>*enddo <BR><BR>!/post1 <BR>!FLST,2,6,1 <BR>!FITEM,2,8567 <BR>!FITEM,2,8563 <BR>!FITEM,2,8559 <BR>!FITEM,2,8555 <BR>!FITEM,2,8551 <BR>!FITEM,2,8568 <BR>!PATH,bh,6,30,20, <BR>!PPATH,P51X,1 <BR>!PDEF,ux,U,X,AVG</DIV>
回复
分享到:

使用道具 举报

发表于 2006-4-25 14:34 | 显示全部楼层
<P>请问楼主你这个命令流哪里弄来的??<BR>怎么跟我们这边一个人的一模一样啊?<BR>是否熟人???</P>
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-5-15 02:35 , Processed in 0.118587 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表