求面面接触分析实例，桩土分析最好，谢啦！

cuiyibin-1 · 发表于 2006-3-28 16:32

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

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

x

求面面接触分析实例

VibInfo · 发表于 2006-3-28 23:54

帖子http://forum.vibunion.com/forum/viewthread.php?tid=8426中给出了一个

AaronSpark · 发表于 2006-4-24 07:11

桩土相互作用接触模型 /filname,pile-soil contact /title，pile /config,nres,2000000 /units,si !国际单位制 fe=0.3 !摩擦系数 fk1=10 !法向接触刚度 fk2=1 !法向接触刚度 fd=3 !切向接触刚度 ftoln=1 !初始渗透因子 A=3.1415926*0.6*0.6 !桩的横截面积 /prep7 et,1,solid45 mp,mu,1,0.2 !桩体材料属性 mp,ex,1,2.5e10 mp,nuxy,1,0.2 mp,dens,1,2500 mp,ex,2,2.0e8 !土层2材料属性 mp,nuxy,2,0.4 mp,dens,2,2000 mp,ex,3,2.5e8 !土层3材料属性 mp,nuxy,3,0.4 mp,dens,3,2000 mp,ex,4,3.0e8 !土层4材料属性 mp,nuxy,4,0.4 mp,dens,4,2000 mp,ex,5,1.8e9 !岩层5材料属性 mp,nuxy,5,0.29 mp,dens,5,2600 tb,concr,1 !桩体的参数 tbdata,,0.7,0.9,2.6,-1 tb,dp,2 !土层2的D-P参数 tbdata,1,19,32,30 tb,dp,3 !土层3的D-P参数 tbdata,1,19,32,30 tb,dp,4 !土层4的D-P参数 tbdata,1,19,32,30 tb,dp,5 !土层5的D-P参数 tbdata,1,0,35,0 cylind,0,0.6,0,21.0,0,90 cylind,0.6,3.0,0,3.0,0,90 cylind,0.6,3.0,3.0,8.0,0,90 cylind,0.6,3.0,8.0,14.0,0,90 cylind,0.6,3.0,14.0,21.0,0,90 cylind,0,0.6,21.0,24.0,0,90 !建立桩尖岩层几何体 cylind,0.6,3.0,21.0,24.0,0,90 !布尔运算粘结体所有的除桩之外的部分 vsel,s,,,2,7,1 vglue,2,3,4,5,6,7 allsel,all numcmp,all /pnum,volu,1 vplot vlist save /pnum,line,1 allsel,all vsel,s,,,1 !桩体附属性，划分网格 aslv,s lsla,s lplot lesize,1,0.1 lesize,2,0.1 lesize,3,,,6 lesize,4,,,6 lesize,5,0.1 lesize,6,0.1 lesize,7,0.5 lesize,8,0.5 lesize,9,0.5 type,1 mat,1 mshape,0,3d mshkey,1 vmesh,all vsel,s,,,2 !土层2附属性，划分网格 vplot aslv,s lsla,s lplot lesize,10,,,6 lesize,11,,,6 lesize,12,,,6 lesize,13,,,6 lesize,14,,,6 lesize,15,,,6 lesize,16,,,6 lesize,17,,,6 lesize,18,0.5 lesize,19,0.5 lesize,20,0.5 lesize,21,0.5 type,1 mat,2 mshape,0,3d mshkey,1 vmesh,all vsel,s,,,4 !土层3附属性，并划分网格 vplot aslv,s lsla,s lplot lesize,22,,,6 lesize,23,,,6 lesize,24,,,6 lesize,25,,,6 lesize,43,0.5 lesize,44,0.5 lesize,45,0.5 lesize,46,0.5 type,1 mat,3 mshape,0,3d mshkey,1 vmesh,all vsel,s,,,5 !土层4附属性，并划分网格 vplot aslv,s lsla,s lplot lesize,26,,,6 lesize,27,,,6 lesize,28,,,6 lesize,29,,,6 lesize,47,0.5 lesize,48,0.5 lesize,49,0.5 lesize,50,0.5 type,1 mat,4 mshape,0,3d mshkey,1 vmesh,all vsel,s,,,6 !岩层5附属性，并划分网格 vplot aslv,s lsla,s lplot lesize,30,,,6 lesize,31,,,6 lesize,32,,,6 lesize,33,,,6 lesize,51,0.5 lesize,52,0.5 lesize,53,0.5 lesize,54,0.5 type,1 mat,5 mshape,0,3d mshkey,1 vmesh,all vsel,s,,,3 !桩底岩层6附属性，并划分网格 vplot aslv,s lsla,s lplot lesize,34,,,6 lesize,35,,,6 lesize,36,,,6 lesize,37,,,6 lesize,41,,,6 lesize,42,,,6 type,1 mat,5 mshape,0,3d mshkey,1 vmesh,all vsel,s,,,7 !桩底岩层6附属性，并划分网格 vplot aslv,s lsla,s lplot lesize,38,,,6 lesize,57,,,6 lesize,58,,,6 lesize,41,0.5 lesize,42,0.5 lesize,55,0.5 lesize,56,0.5 type,1 mat,5 mshape,0,3d mshkey,1 vmesh,all !用接触向导建立接触对 et,2,170 et,3,173 keyopt,3,4,1 keyopt,3,9,0 !消除初始渗透 keyopt,3,12,2 !不分开接触，允许滑动，法向不分开 !桩侧，建立目标面 r,1,,,fk1,ftoln, , , rmore,,,,,,5 r,2,,,fk2,ftoln, , , rmore,,,,,,5 mp,mu,6,0.2 allsel,all vsel,s,,,1 aslv,s asel,s,,,3 /pnum,area,1 aplot cm,target1,area type,2 real,1 mat,6 nsla,s,1 esurf,all ! 桩侧，建立接触面 allsel vsel,s,,,2 vsel,a,,,4,6,1 aslv,s asel,s,,,9 asel,a,,,22 asel,a,,,26 asel,a,,,30 aplot cm,contact1,area type,3 real,1 mat,6 nsla,s,1 esurf,all !桩尖，建立目标面 allsel,all vsel,s,,,3 aslv,s asel,r,,,15 aplot nsla,s,1 type,2 real,2 mat,6 cm,target2,elem esurf,all esel,s,type,,2 eplot !桩尖，建立接触面 allsel,all vsel,s,,,1 aslv,s asel,s,,,2 aplot nsla,s,1 type,3 real,2 mat,6 cm,contact2,area esurf,all esel,s,type,,3 eplot /psymb,esys,1 allsel gplot save finish !求解过程 /solu !定义面约束条件 asel,s,loc,x,0 !侧面施加对称约束 da,all,symm asel,s,loc,y,0 da,all,symm asel,s,loc,z,24 !底面试加约束 da,all,all allsel,all !土外层施加约束 asel,s,,,8 asel,a,,,20 asel,a,,,24 asel,a,,,28 asel,a,,,32 aplot da,all,all allsel,all !设置非线性选项 acel,0,0,9.8 !定义重力加速度 neqit,50, !迭代次数 nropt,Modi ! 用改变的牛顿-拉普拉斯方程求解 nlgeom,on ! 打开大变形选项 autot,on ! 打开自动时间步 lnsrch,on !自动线性搜索 pred,on !打开预测求解器 outpr,all,all OUTRES,all,all allsel,all cncheck,detail !检查接触状态 cncheck,summary /solu CNCHECK,POST FINISH !载荷步,加载阶段 /solu FirstF=2E+06/A LastF=3E+06/A Increment=2E+05/A *do,Force,FirstF,LastF,Increment asel,s,loc,z,0 sfa,all,,pres,Force allsel,all nsubst,10,100,5 solve *enddo finish !进入POST1后处理器 /post1 set,last pldisp,3 plnsol,u,x pinsol,s,x pinsol,s,y pinsol,s,z pinsol,s,xy pinsol,s,eqv pinsol,epto,eqv pinsol,cont,pres pinsol,cont,sfric finish !退出 /exit,all

AaronSpark · 发表于 2006-4-24 07:15

桩土接触二维问题 这是二维平面应变问题，比较简单。 桩为混凝土，采用线弹性，土采用DP模型。 二维接触分别采用TARGE169和CONTA171单元。 /PREP7 ET,1,PLANE42 KEYOPT,1,3,2 ET,2,TARGE169 ET,3,CONTA171 MP,EX,1,2E10 MP,PRXY,1,0.2 MP,DENS,1,2500 MP,EX,2,20E6 MP,PRXY,2,0.45 TB,DP,2 TBDATA,,10,30,20 MP,DENS,2,1700 R,1 BLC4,,,5,10 BLC4,5,15,10,-15 BLC4,5,15,2.6,-0.5 BLC4,5,14.5,0.6,-12 BLC4,7,14.5,0.6,-12 A***A,2,3 A***A,6,4 A***A,2,5 AGLUE,ALL NUMCMP,ALL BLC4,5,15,2.6,-0.5 BLC4,5,14.5,0.6,-12 BLC4,7,14.5,0.6,-12 AGLUE,3,4,5 NUMCMP,ALL WPOFFS,,2.5 WPROTA,,90 ASEL,S,,,1,2 A***W,ALL NUMCMP,ALL ALLSEL,ALL ESIZE,0.2 TYPE,1 MAT,1 AMESH,1,3 MAT,2 AMESH,4,8 LSEL,S,,,14 LSEL,A,,,16,20 LSEL,A,,,21,25,2 NSLL,S,1 TYPE,2 REAL,1 ESURF,ALL LSEL,S,,,4,9 LSEL,A,,,11,12 NSLL,S,1 TYPE,3 REAL,1 ESURF,ALL ALLSEL,ALL FINI /SOLU ANTYPE,STATIC LSEL,S,LOC,X DL,ALL,,SYMM LSEL,S,LOC,X,15 DL,ALL,,UX, LSEL,S,LOC,Y, DL,ALL,,UY, ALLSEL,ALL SFL,13,PRES,10000 ACEL,,9.8 TIME,1 NSUBST,50,200,20 AUTOTS,ON NLGEOM,ON PRED,ON LNSRCH,ON OUTRES,ALL SOLVE FINI /POST1 PLNSOL,S,EQV ETABLE,CONX,NMISC,21 ESEL,U,TYPE,,1 PLETAB,CONX,NOAV ETABLE,CONP,SMISC,5 PLETAB,CONP,NOAV

zhaofuc · 发表于 2006-4-24 14:30

太感谢了 我正打算用二维接触模拟压力分散型抗浮锚杆 可以参考

zhaofuc · 发表于 2006-4-25 14:42

我运行了一下你的程序，怎么不行？ 下面是我用二维接触做的压力型抗浮锚干的一段命令流，请指点一下哪里有问题？非常感谢！ /prep7 et,1,plane42 et,2,solid45 mp,ex,1,2.5e10 !桩的弹性模量 mp,nuxy,1,0.2 !桩的泊松比 mp,dens,1,2500 !桩的密度 mp,ex,2,2.5e8 !土的弹性模量 mp,nuxy,2,0.4 !土的泊松比 mp,dens,2,2000 !土的密度 tb,dp,2 tbdata,1,19,32,30 !粘聚力c为19,摩擦角为32度,膨胀角为30
RECTNG,0,1,0,8, !面1 RECTNG,1,10,0,7, !面2 RECTNG,0,10,0,-8, !面3 /pnum,area,1 /pnum,line,1 asel,s,,,2,3,1 aglue,all numcmp,all allsel aplot
lsel,s,,,1,3,2 !划分面1 lesize,all,,,2 lsel,s,,,2,4,2 lesize,all,,,16 type,2 mat,1 amesh,1 
lsel,s,,,6,8,2 !划分面2 lesize,all,0.5 lsel,s,,,5 lesize,all,0.5 lsel,s,,,7 lesize,all,0.5 type,2 mat,2 amesh,2 
lsel,s,,,11 !划分面3 lesize,all,,,2 lsel,s,,,5 lesize,all,0.5 lsel,s,,,12 lesize,all,0.5 lsel,s,,,10 lesize,all,0.5 lsel,s,,,9 lesize,all,0.5 type,2 mat,2 amesh,3

!接触单元设置 allsel et,3,169 et,4,171 keyopt,4,9,0 keyopt,4,12,2 !r,1 mp,mu,2,0.2
asel,s,,,1 nsla,s,1 nsel,r,loc,x,1 cm,target,line r,1 type,3 esurf
allsel asel,s,,,2 nsla,s,1 nsel,r,loc,x,1 cm,target,line !r,1 type,4 esurf allsel
 !底端接触 allsel asel,s,,,1 nsla,s,1 nsel,r,loc,y,0 r,1 type,3 esurf allsel asel,s,,,3
nsla,s,1 nsel,r,loc,y,0 type,4 esurf allsel,all
 !设置边界条件 /solu asel,s,,,2 nsla,s,1 dl,6,2,ux asel,s,,,3 nsla,s,1 dl,12,3,ux dl,10,3,ux !dL,7,2,UX dl,9,3,all,all allsel
 !施加重力加速度及荷载 !acel,0,9.8,0, asel,s,,,1 !lplot sfl,3,pres,600

!设置分析选项并求解 antype,static nlgeom,on time,1 auto,on nropt,full nsub,100,10000,100 outres,all,all allsel solve finish

yanwei_69 · 发表于 2008-11-10 11:24

我也是刚学习这部分，运行了一下你的程序，感觉是不是没考虑对称约束，我见你在左侧也是加的固定约束，你看看是不是啊，是不是应该加对称约束，你不知道单元设置里设置了对称没有。不知道你解决了没有？要是解决了和小弟说说，

博博 · 发表于 2010-9-3 10:13

不错的东西，正在学习

ymcheng123 · 发表于 2012-4-16 09:04

我也在做桩土分析

ymcheng123 · 发表于 2012-4-17 16:15

遇到问题了，好难解决啊

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

[材料与接触] 求面面接触分析实例，桩土分析最好，谢啦！

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

回复：（cuiyibin-1）求面面接触分析实例，桩土分析最...

回复：（cuiyibin-1）求面面接触分析实例，桩土分析最...

回复：（cuiyibin-1）求面面接触分析实例，桩土分析最...

回复

菜鸟回答