【求助】一个支挡结构的结构与土的接触问题,试过了本论给的几乎全部参考意见,如调子步数,收敛准则,非线性搜索,迭代次数等都不行。本人已经改无可改,黔驴技穷了,真心求教各位版主,高手。
命令如下:
!准备工作,定义单元类型,材料属性 fk1=10 !法向接触刚度 fk1=1 !法向接触刚度 ftoln=1 !初始渗透因子 finish /clear /filname,ZJTLDTQ,1 /prep7 et,1,solid45 r,1 et,2,170 et,3,174 keyopt,3,9,0 keyopt,3,12,2 r,2,,,fk1,ftoln, , , rmore,,,,,,5 r,3,,,fk2,ftoln, , , rmore,,,,,,5 mp,ex,1,1e9 mp,prxy,1,0.3 mp,dens,1,2400 mp,ex,2,1e8 mp,prxy,2,0.3 mp,dens,2,2200 mp,ex,3,3e7 mp,prxy,3,0.36 mp,dens,3,2100 mp,ex,4,3e7 mp,prxy,4,0.36 mp,dens,4,2100 mp,ex,5,2.5e10 mp,prxy,5,0.2 mp,dens,5,2500 mp,ex,6,3e9 mp,prxy,6,0.23 mp,dens,6,2300 mp,mu,7,0.4 tb,dp,1 tbdata,,1e6,35,0 tb,dp,2 tbdata,,1e8,40,0 tb,dp,3 tbdata,,2e4,30,28 tb,dp,4 tbdata,,2e4,30,27 !定义桩长 Hpile=14 !建模,基岩 k,,0,0,0 k,,0,20 k,,20,26.5 k,,28.6,32.1 k,,50.6,52.8 k,,56.2,59.3 k,,62.0,59.5 k,,64.2,59.9 k,,66.0,64.0 k,,69.8,65.9 k,,73.4,66.0 k,,77.5,64.3 k,,97.5,71.9 k,,97.5,0,0 *do,i,1,13 l,i,i+1 *enddo l,1,14 al,all vext,1,,,,,-5 !创建承台 block,41.2,45.2,48.2,49.8,0,-5 numcmp,all !创建第二层土体 k,,0,24.3 k,,20,30.8 k,,31,39.3 k,,46,52.5 k,,51.7,57.5 k,,56.2,61.7 k,,62.8,61.7 k,,64.2,65.9 k,,66.8,66.3 k,,69.7,68.5 k,,73.6,68.4 k,,76.7,66.5 k,,97.5,73 Vplot L,2,37 L,37,38 L,38,39 L,39,33 L,35,40 *do,i,40,48 l,i,i+1 *enddo L,49,13 Lsel,s,,,2,12 Lsel,a,,,55,69 Lsel,a,,,47,48 al,all vext,23,,,,,-5 numcmp,all !将实心土变成空心,创建桩 block,41.7,44.7,48.2-hpile,48.2,-2,-4 vsbv,1,4 block,41.7,44.7,48.2-hpile,48.2,-2,-4 vsbv,3,1 block,41.7,44.7,48.2-hpile,48.2,-2,-4 vplot numcmp,all !创建挡土墙 allsel Vplot k,,41.2,50.3,0 k,,41.5,50.3,0 k,,41.5,50.8,0 k,,41.8,50.8,0 k,,42.3,61.8,0 k,,42.9,61.8,0 k,,45.1,57,0 k,,47.1,57,0 a,36,102,103,104,105,106,107,108,109,40,35 vext,68,,,,,-5 numcmp,all !创建第三层土体 k,,52.7,61.3 k,,50.6,60.2 a,40,41,42,121,122,109 vext,81,,,,,-5 k,,0,28.5 k,,27.6,40 a,129,130,36,33,39,38,37 vext,89,,,,,-5 numcmp,all !创建挡土墙墙背土体 a,107,108,109,122,121,42 vext,98,,,,,-5 numcmp,all !粘结为一体 vglue,1,2 vglue,3,6 vglue,3,7 vglue,2,8 Allsel numcmp,all !网格划分准备工作,切分 !切割基岩 wpoff,35 wprota,,,90 vsbw,all wpoff,,,21.2 vsbw,all wpoff,,48.2-hpile wprota,,90 vsbw,4 wpcsys,-1 wpoff,41.7,48.2-hpile wprota,,,90 vsbw,6 vsbw,14 wpoff,,,3 vsbw,16 vsbw,17 numcmp,all aadd,53,95 numcmp,all !划分网格 !基岩 vsel,s,,,13 vatt,1,1,1 esize,0.5 vsweep,all vsel,s,,,20 vatt,1,1,1 lesize,254,,,50 lesize,292,,,26 lesize,282,,,16 lesize,5,,,16 vsweep,all !桩 vsel,s,,,1 vatt,5,1,1 esize,0.5 vsweep,all !基岩 vsel,s,,,16 vsel,a,,,19 vatt,1,1,1 lesize,138,,,20 lesize,271,,,20 vsweep,all vsel,s,,,6 vatt,1,1,1 lesize,256,,,8 lesize,271,,,21 vsweep,all !承台 vsel,s,,,8 vatt,5,1,1 esize,0.5 vsweep,all !第二层 vsel,s,,,4 vatt,2,1,1 lesize,260,,,1 lesize,259,,,7 lesize,241,,,7 lesize,219,,,19 vsweep,all vsel,s,,,18 vatt,2,1,1 lesize,284,,,1 lesize,280,,,2 lesize,48,,,3 lesize,59,,,6 lesize,60,,,14 lesize,61,,,12 vsweep,all vsel,s,,,14 vsel,a,,,17 vatt,2,1,1 vsweep,all !第三层 vsel,s,,,9 vatt,3,1,1 lesize,220,,,19 vsweep,all !回填土 vsel,s,,,7 vatt,4,1,1 vsweep,all !第三层 vsel,s,,,2 vatt,3,1,1 vsweep,all !挡土墙 vsel,s,,,5 vatt,6,1,1 vsweep,all !外围 !基岩 vsel,s,,,3 vatt,1,1,1 esize,1 vsweep,all vsel,s,,,10 vatt,1,1,1 esize,1 vsweep,all !二层 vsel,s,,,12 vatt,2,1,1 esize,0.5 vsweep,all vsel,s,,,15 vatt,2,1,1 esize,0.8 vsweep,all !三层 vsel,s,,,11 vatt,3,1,1 esize,0.5 vsweep,all !设置接触面 !桩 allsel asel,s,,,6 asel,a,,,54,56,2 asel,a,,,27 asel,a,,,48 nsla,s,1 type,2 real,2 mat,7 esurf,all !承台 Allsel asel,s,,,20,22 asel,a,,,151 nsla,s,1 type,2 real,2 mat,7 esurf,all !挡墙 Allsel asel,s,,,76,79 nsla,s,1 type,2 real,2 mat,7 esurf,all !设置目标面 !基岩 allsel asel,s,,,57,60 asel,a,,,55 nsla,s,1 type,3 real,3 esurf,all !第二层 Allsel asel,s,,,62,65 asel,a,,,46,47 asel,a,,,118 asel,a,,,127 asel,a,,,135 asel,a,,,140 nsla,s,1 type,3 real,3 esurf,all !第三层 Allsel asel,s,,,82 asel,s,,,95 nsla,s,1 type,3 real,3 esurf,all !杂填土 Allsel asel,s,,,88 asel,a,,,94 nsla,s,1 type,3 real,3 esurf,all !挡墙 Allsel asel,s,,,80 nsla,s,1 type,3 real,3 esurf,all !施加约束 asel,a,loc,z,0 da,all,symm asel,s,loc,z,-5 da,all,uz asel,s,loc,x,0 asel,a,loc,y,0 asel,a,loc,x,97.5 da,all,all !施加荷载 nsel,s,loc,x,44,47.5 nsel,r,loc,y,61.8 nsel,r,loc,z,-5,0 sf,all,pres,57600 allsel !历程后处理 /solu antype,static nlgeom,on autots,on acel,,9.8 nropt,full Lnsrch,on time,1 Eqslv,-1 nsubst,10,1000,10 Neqit,50 outres,all,all solve
|