ANSYS应用——m法计算桩基时模拟桩与土作用地命令流_第1页
ANSYS应用——m法计算桩基时模拟桩与土作用地命令流_第2页
ANSYS应用——m法计算桩基时模拟桩与土作用地命令流_第3页
ANSYS应用——m法计算桩基时模拟桩与土作用地命令流_第4页
ANSYS应用——m法计算桩基时模拟桩与土作用地命令流_第5页
已阅读5页，还剩3页未读， 继续免费阅读

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1、实用文档ANSYS用一一m法计算桩基时模拟桩与土作用的命令流!z1spring1.mac/PREP7WPCSYS,-1 ! 工作平面恢复到默认状态CSYS,0dpile1=1.5! 桩直径 *allsel,allcmsel,s,znsll,s,1NSEL,r,LOC,Z,-5-5-5.36,-6-5-5.36! 建立泥面以下桩的节点集合CM,SPNODE1,NODEallsel,all,allcmsel,s,znsll,s,1NSEL,r,LOC,Z,-6-5-5.36,-20-5-5.36CM,SPNODE2,NODEallsel,all,allcmsel,s,znsll,s,1NSEL,r

2、,LOC,Z,-20-5-5.36,-25-5-5.36CM,SPNODE3,NODEallsel,all,allcmsel,s,znsll,s,1NSEL,r,LOC,Z,-25-5-5.36,-27-5-5.36CM,SPNODE4,NODEallsel,all,allcmsel,s,znsll,s,1NSEL,r,LOC,Z,-27-5-5.36,-31-5-5.36CM,SPNODE5,NODEallsel,all,all!spring1/PREP7CMSEL,S,SPNODE1 !选择泥面以下, 桩的节点集合* GET,enum,ELEM,COUNT, , , ,* GET, NNU

3、M, NODE, 0, COUNT,* GET, nmin, NODE, 0, num, min,*DiM,AA1,ARRAY,NNUM*DiM,AA2,ARRAY,NNUM*DiM,AA3,ARRAY,NNUMCMSEL,S,SPNODE1*do,i,1,nnumAA1(i)=NMINnmin=ndnext(nmin)*enddo ! 给集合 AA1 赋值*do,i,1,nnumNSEL,S,AA1(i)ESLN,S*GET,emin,ELEM,0,num,min ! 获取最小的单元编号n1=NELEM(emin,1)n2=NELEM(emin,2)Z12=NZ(N1)-NZ(N2)X12=

4、NX(N1)-NX(N2)*IF,(Z12*X12),GE,0,THENSIG=1*ELSESIG=-1*ENDIFZ12=ABS(Z12)L12=DISTND(N1,N2)CLI=TAN(ACOS(Z12/L12)NGEN,2,100000,AA1(I), , ,-1,1*CLI*SIGNGEN,2,200000,AA1(I), , ,1AA2(i)=AA1(i)+100000AA3(i)=AA1(i)+200000TYPE,3 ! 将单元赋为弹簧单元, 为 etsp 单元类型, 见 begin!MAT,mpsp !mpsp 为材料号, 见 beginESYS,0R,enum+i,8000*

5、(-10-5.36-NZ(AA1(I)*(L12)*dpile1,!根据 M法给弹簧赋予刚度, 其中 zlsize 为桩单元长!REAL,enum+ie,AA1(I),AA2(I)e,AA1(I),AA3(I)*enddo/SOLALLSEL,ALL* do,i,1,NNUMD,AA2(i),allD,AA3(i),all*ENDDO!spring1!spring2/PREP7CMSEL,S,SPNODE2 !选择泥面以下, 桩的节点集合* GET,enum,ELEM,COUNT, , , ,* GET, NNUM, NODE, 0, COUNT,* GET, nmin, NODE, 0, n

6、um, min,* DiM,AA1,ARRAY,NNUM* DiM,AA2,ARRAY,NNUM* DiM,AA3,ARRAY,NNUMCMSEL,S,SPNODE2*do,i,1,nnumAA1(i)=NMINnmin=ndnext(nmin)* enddo ! 给集合 AA1 赋值*do,i,1,nnumNSEL,S,AA1(i)ESLN,S* GET,emin,ELEM,0,num,min ! 获取最小的单元编号n1=NELEM(emin,1)n2=NELEM(emin,2)Z12=NZ(N1)-NZ(N2)X12=NX(N1)-NX(N2)*IF,(Z12*X12),GE,0,THEN

7、SIG=1*ELSESIG=-1*ENDIFZ12=ABS(Z12)L12=DISTND(N1,N2)CLI=TAN(ACOS(Z12/L12)NGEN,2,100000,AA1(I), , ,-1,1*CLI*SIGNGEN,2,200000,AA1(I), , ,1AA2(i)=AA1(i)+100000AA3(i)=AA1(i)+200000TYPE,3 ! 将单元赋为弹簧单元, 为 etsp 单元类型, 见 begin!MAT,mpsp !mpsp 为材料号, 见 beginESYS,0R,enum+i,4500*(-10-5.36-NZ(AA1(I)*(L12)*dpile1,!根据

8、 M法给弹簧赋予刚度, 其中 zlsize 为桩单元长!REAL,enum+ie,AA1(I),AA2(I)e,AA1(I),AA3(I)*enddo/SOLALLSEL,ALL* do,i,1,NNUMD,AA2(i),allD,AA3(i),all*ENDDO!spring2!spring3/PREP7CMSEL,S,SPNODE3 !选择泥面以下, 桩的节点集合* GET,enum,ELEM,COUNT, , , ,* GET, NNUM, NODE, 0, COUNT,* GET, nmin, NODE, 0, num, min,*DiM,AA1,ARRAY,NNUM*DiM,AA2,

9、ARRAY,NNUM*DiM,AA3,ARRAY,NNUMCMSEL,S,SPNODE3*do,i,1,nnumAA1(i)=NMINnmin=ndnext(nmin)*enddo ! 给集合 AA1 赋值*do,i,1,nnumNSEL,S,AA1(i)ESLN,S*GET,emin,ELEM,0,num,min ! 获取最小的单元编号n1=NELEM(emin,1)n2=NELEM(emin,2)Z12=NZ(N1)-NZ(N2)X12=NX(N1)-NX(N2)*IF,(Z12*X12),GE,0,THENSIG=1*ELSESIG=-1*ENDIFZ12=ABS(Z12)L12=DIS

10、TND(N1,N2)CLI=TAN(ACOS(Z12/L12)NGEN,2,100000,AA1(I), , ,-1,1*CLI*SIGNGEN,2,200000,AA1(I), , ,1AA2(i)=AA1(i)+100000AA3(i)=AA1(i)+200000TYPE,3 ! 将单元赋为弹簧单元, 为 etsp 单元类型, 见 begin!MAT,mpsp !mpsp 为材料号, 见 beginESYS,0R,enum+i,10000*(-10-5.36-NZ(AA1(I)*(L12)*dpile1,!根据 M法给弹簧赋予刚度, 其中 zlsize 为桩单元长!REAL,enum+ie

11、,AA1(I),AA2(I)e,AA1(I),AA3(I)*enddo/SOLALLSEL,ALL*do,i,1,NNUMD,AA2(i),allD,AA3(i),all*ENDDO!spring3!spring4/PREP7CMSEL,S,SPNODE4 !选择泥面以下, 桩的节点集合* GET,enum,ELEM,COUNT, , , ,* GET, NNUM, NODE, 0, COUNT,* GET, nmin, NODE, 0, num, min,* DiM,AA1,ARRAY,NNUM* DiM,AA2,ARRAY,NNUM* DiM,AA3,ARRAY,NNUMCMSEL,S,S

12、PNODE3*do,i,1,nnumAA1(i)=NMINnmin=ndnext(nmin)* enddo ! 给集合AA1 赋值*do,i,1,nnumNSEL,S,AA1(i)ESLN,S* GET,emin,ELEM,0,num,min ! 获取最小的单元编号n1=NELEM(emin,1)n2=NELEM(emin,2)Z12=NZ(N1)-NZ(N2)X12=NX(N1)-NX(N2)*IF,(Z12*X12),GE,0,THENSIG=1*ELSESIG=-1*ENDIFZ12=ABS(Z12)L12=DISTND(N1,N2)CLI=TAN(ACOS(Z12/L12)NGEN,2

13、,100000,AA1(I), , ,-1,1*CLI*SIGNGEN,2,200000,AA1(I), , ,1AA2(i)=AA1(i)+100000AA3(i)=AA1(i)+200000TYPE,3 ! 将单元赋为弹簧单元, 为 etsp 单元类型, 见 begin!MAT,mpsp !mpsp 为材料号, 见 beginESYS,0R,enum+i,20000*(-10-5.36-NZ(AA1(I)*(L12)*dpile1,!根据 M法给弹簧赋予刚度, 其中 zlsize 为桩单元长!REAL,enum+ie,AA1(I),AA2(I)e,AA1(I),AA3(I)*enddo/S

14、OLALLSEL,ALL* do,i,1,NNUMD,AA2(i),allD,AA3(i),all*ENDDO!spring4!spring5/PREP7CMSEL,S,SPNODE5 !选择泥面以下, 桩的节点集合* GET,enum,ELEM,COUNT, , , ,* GET, NNUM, NODE, 0, COUNT,* GET, nmin, NODE, 0, num, min,*DiM,AA1,ARRAY,NNUM*DiM,AA2,ARRAY,NNUM*DiM,AA3,ARRAY,NNUMCMSEL,S,SPNODE3*do,i,1,nnumAA1(i)=NMINnmin=ndnex

15、t(nmin)*enddo ! 给集合 AA1 赋值*do,i,1,nnumNSEL,S,AA1(i)ESLN,S*GET,emin,ELEM,0,num,min ! 获取最小的单元编号n1=NELEM(emin,1)n2=NELEM(emin,2)Z12=NZ(N1)-NZ(N2)X12=NX(N1)-NX(N2)*IF,(Z12*X12),GE,0,THENSIG=1*ELSESIG=-1*ENDIFZ12=ABS(Z12)L12=DISTND(N1,N2)CLI=TAN(ACOS(Z12/L12)NGEN,2,100000,AA1(I), , ,-1,1*CLI*SIGNGEN,2,200000,AA1(I), , ,1AA2(i)=AA1(i)+100000AA3(i)=AA1(i)+200000TYPE,3 ! 将单元赋为弹簧单元, 为 etsp

人人文库> 全部分类> 教育资料 > 备课教案

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论