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

下载本文档

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

文档简介

1、ANSYS应用m法计算桩基时模拟桩与土作用的命令流!z1spring1.mac/PREP7 WPCSYS,-1 !工作平面恢复到默认状态CSYS,0 dpile1=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,LOC,Z,

2、-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, NNUM, NODE, 0

3、, COUNT, *GET, nmin, NODE, 0, num, min, *DiM,AA1,ARRAY,NNUM *DiM,AA2,ARRAY,NNUM *DiM,AA3,ARRAY,NNUM CMSEL,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=NX(N1)

4、-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*(-10-5.36-NZ(AA

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

6、AY,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,THENSIG=1*ELSESIG=-1*ENDIFZ12=ABS

7、(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, !根据M法给弹簧赋予刚度,其中zlsize为桩单元长!REAL,enum+ie,

8、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/PREP7 CMSEL,S,SPNODE3 !选择泥面以下,桩的节点集合*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*

9、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,100000,AA1(

10、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,AA1(I),AA2(I)e,AA1(I),AA3(I)*enddo/SOLALLSEL,ALL*do,i,1,NNUM

11、D,AA2(i),allD,AA3(i),all*ENDDO!spring3!spring4/PREP7 CMSEL,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,SPNODE3*do,i,1,nnumAA1(i)=NMINnmin=ndnext(nmin)*enddo !给集合AA1赋值*do,i

12、,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)+1

13、00000AA3(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/SOLALLSEL,ALL*do,i,1,NNUMD,AA2(i),allD,AA3(i),all*ENDDO!spring4!spring5/PREP7 CMSEL,S

14、,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=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,100000,AA1(I), , ,-1,1*CLI*SIGNGEN,2,200000,AA1(I), , ,1AA2(i)=AA1(i)+100000AA3(i)=AA1(i)+200000TYPE,3 !将单元赋为弹簧单元,

人人文库> 全部分类> 教育资料 > 课件下载

温馨提示

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

最新文档

评论

0/150

提交评论