1、/config,nres,2000000 /units,si !国际单位制fk1=1!法向接触刚度 ftoln=0.1!初始渗透因子/prep7et,1,solid45 A=3.14*0.3*0.3cylind,0,0.192,0,11,0,90 !旋喷桩上部内侧 v1cylind,0.192,0. 2,0,11,0,90!钢管 v2cylind,0.2,0.3,0,11,0,90 !旋喷桩上部外侧 v3 cylind,0,0.3,11,22,0,90 !旋喷桩下部 v4cylind,0.3,1.5,0,2.6,0,90 !2.6m 淤泥质粉质粘土 v5cylind,0.3,1.5,2.6,7
2、.92,0,90 !5.32m 粉土 v6cylind,0.3,1.5,7.92,16.3,0,90 !8.38m 淤泥 v7cylind,0.3,1.5,16.3,18.67,0,90!2.37m 粉质粘土 v8cylind,0.3,1.5, 18.67,20.47 ,0,90!1.8m 粉质粘土 v9cylind,0.3,1.5, 20.47,25.3 ,0,90!4.83m 粉土 v10cylind,0,0.3,22,42,0,90 !桩底内侧 20m 粉质粘土 v11cylind,0.3,1.5,25.3,42,0,90 !桩底外侧 20m 粉质粘土 v12allsel,allvsel
3、,s,2mp,ex,2,2e11 !钢管mp,nuxy,2,0.3mp,dens,2,7850allsel,allvsel,s,1,4vsel,u,2mp,ex,3,1.5e8mp,nuxy,3,0.2mp,dens,3,2131!tb,concr,3 !水泥土破坏模型，材料属性!tbdata,0.7,0.9,0.3e6,2e6!张开裂缝的剪力传递系数 ,闭合裂缝的剪力传递系数 ,单轴抗拉强度, 单轴抗压强度allsel,allvsel,s,5vplotmp,ex,5,3e6 !淤泥质粉质黏土mp,nuxy,5,0.35mp,dens,5,1950 !mp,mu,5,0.633 !设置摩擦系数
4、为 0.633tb,dp,5 !土层 5 的 D-P 参数 tbdata,1,13.2e3,10.1,0 !土的粘聚力，摩擦角，膨胀角allsel,allvsel,s,6vplotmp,ex,6,9e6 !粉土mp,nuxy,6,0.3mp,dens,6,1970 !mp,mu,6,0.633 !设置摩擦系数为 0.633tb,dp,6 !土层 6 的 D-P 参数 tbdata,1,15.2e3,21.2,30 !土的粘聚力，摩擦角，膨胀角allsel,allvsel,s,7vplotmp,ex,7,1.5e6 !淤泥mp,nuxy,7,0.4mp,dens,7,1670 !mp,mu,6,
5、0.633 !设置摩擦系数为 0.633tb,dp,7 !土层 7 的 D-P 参数 tbdata,1,12.9e3,9,0 !土的粘聚力，摩擦角，膨胀角allsel,allvsel,s,8vplotmp,ex,8,5e6 !粉质粘土mp,nuxy,8,0.35mp,dens,8,1880 !mp,mu,8,0.633 !设置摩擦系数为 0.633tb,dp,8 !土层 8 的 D-P 参数 tbdata,1,15.8e3,15.2,0 !土的粘聚力，摩擦角，膨胀角allsel,allvsel,s,9vplotmp,ex,9,8e6 !粉质粘土mp,nuxy,9,0.35mp,dens,9,2
6、000 !mp,mu,9,0.633 !设置摩擦系数为 0.633tb,dp,9 !土层 9 的 D-P 参数 tbdata,1,37.1e3,19.8,0 !土的粘聚力，摩擦角，膨胀角allsel,allvsel,s,10vplotmp,ex,10,10e6 !粉土mp,nuxy,10,0.3mp,dens,10,1980 !mp,mu,10,0.633 !设置摩擦系数为 0.633tb,dp,10 !土层 10 的 D-P 参数 tbdata,1,14.3e3,20.1,10 !土的粘聚力，摩擦角，膨胀角allsel,allvsel,s,11,12vplotmp,ex,11,6e6 !粉质
7、粘土mp,nuxy,11,0.35mp,dens,11,1910 !mp,mu,10,0.633 !设置摩擦系数为 0.633tb,dp,11 !土层 11 的 D-P 参数 tbdata,1,23.4e3,15.9,0 !土的粘聚力，摩擦角，膨胀角!布尔运算粘结体所有的除桩之外的部分ALLSEL,ALLvsel,s,5,12,1vplotvglue,allvsel,s,1,4vplotvglue,allallsel,all/pnum,line,1saveallsel,allallsel,alllsel,s,3,4lsel,a,13,14lesize,all,6lsel,s,7,8lsel,a
8、,18,19lesize,all,0.1lsel,s,55,58lesize,all,1type,1mat,2 !钢管vmesh,6allsel,allvsel,s,1,7,6vsel,a,8vplottype,1mat,3 !水泥土lsel,s,1,2lesize,all,6lsel,s,55,56lesize,all,1lsel,s,22,24,2lesize,all,10,0.8,1lsel,s,3lesize,all,6lsel,s,13lesize,13,6lsel,s,25lesize,25,6lsel,s,7,9lsel,a,18,19lsel,a,30,31lsel,a,67,
9、69lesize,all,0.1lsel,s,4,6 !体 8lsel,a,14,26,12lesize,all,6lsel,s,57,58lesize,all,1lsel,s,65,66lesize,all,10,0.8,1aadd,35,31,2aslv,s lsla,s FLST,2,3,4,ORDE,3 FITEM,2,6 FITEM,2,57 FITEM,2,65 LCCAT,P51X FLST,2,3,4,ORDE,3 FITEM,2,5 FITEM,2,58 FITEM,2,66 LCCAT,P51X vsweep,1vsweep,7vsweep,8allsel,allvsel,
10、u,1,6,5vsel,u,7,8vsel,u,18vplotlsel,s,44,46,2lesize,all,6lsel,s,43,45,2lesize,all,5,0.8,1lsel,s,51,54,lesize,all,0.1lsel,s,47,49,2lesize,all,6lsel,s,48,50,2lesize,all,5,0.8,1type,1mat,5vsweep,5vsel,s,13lsel,s,136,139lesize,all,0.1lsel,s,59,61,2lesize,all,6lsel,s,60,62,2lesize,all,5,0.8,1type,1mat,6v
11、sweep,13vsel,s,14lsel,s,140,143lesize,all,0.1lsel,s,71,73,2lesize,all,6lsel,s,72,74,2lesize,all,5,0.8,1type,1mat,7vsweep,14vsel,s,15lsel,s,144,147lesize,all,0.1lsel,s,83,85,2lesize,all,6lsel,s,84,86,2lesize,all,5,0.8,1type,1mat,8vsweep,15vsel,s,16lsel,s,148,151lesize,all,0.1lsel,s,95,97,2lesize,all,
12、6lsel,s,96,98,2lesize,all,5,0.8,1type,1mat,9vsweep,16vsel,s,19lsel,s,158,163lesize,all,0.1lsel,s,107,109,2lesize,all,6lsel,s,108,110,2lesize,all,5,0.8,1type,1mat,10vsweep,19vsel,s,17lsel,s,152,156lesize,all,0.1lsel,s,118,128,10lesize,all,6lsel,s,155,157,2lesize,all,5,0.8,1type,1mat,11vsweep,17vsel,s
13、,18lsel,s,154,158,2lsel,a,159lsel,a,123lesize,all,0.1lsel,s,118,120,lesize,all,6lsel,s,115,117lesize,all,6 type,1mat,11vsweep,18!用接触向导建立接触对et,2,170et,3,173keyopt,3,9,0! 包括由几何形状和接触表面偏移所引起的穿透keyopt,3,12,2 !不分开接触，允许滑动，法向不分开!桩建立目标面r,20 ,fk1,ftoln, , ,mp,mu,100,0.85allsel,allvsel,s,1,6,5vsel,a,7,8aslv,sa
14、plotasel,s,14asel,a,19,39,20aplotcm,target,areatype,2real,20mat,100nsla,s,1esurf,allesel,s,type,2 eplot r,21,fk1,ftoln, , ,mp,mu,101,0.3allsel,allvsel,s,1,6,5vsel,a,7,8aslv,saplotasel,s,14asel,a,19,39,20aplotcm,target,areatype,2real,21mat,101nsla,s,1esurf,allesel,s,type,2 eplot r,22,fk1,ftoln, , ,mp
15、,mu,102,0.15allsel,allvsel,s,1,6,5vsel,a,7,8aslv,saplotasel,s,14asel,a,19,39,20aplotcm,target,areatype,2real,22mat,102nsla,s,1esurf,allesel,s,type,2 eplot r,23,fk1,ftoln, , ,mp,mu,103,0.3allsel,allvsel,s,1,6,5vsel,a,7,8aslv,saplotasel,s,14asel,a,19,39,20aplotcm,target,areatype,2real,23mat,103nsla,s,
16、1esurf,allesel,s,type,2 eplot r,24,fk1,ftoln, , ,mp,mu,104,0.35allsel,allvsel,s,1,6,5vsel,a,7,8aslv,saplotasel,s,14asel,a,19,39,20aplotcm,target,areatype,2real,24mat,104nsla,s,1esurf,allesel,s,type,2 eplot r,25,fk1,ftoln, , ,mp,mu,105,0.3allsel,allvsel,s,1,6,5vsel,a,7,8aslv,saplotasel,s,14asel,a,19,
17、39,20aplotcm,target,areatype,2real,25mat,105nsla,s,1esurf,allesel,s,type,2 eplot !桩建立接触面 allsel,allvsel,u,1,6,5vsel,u,7,8vplotasel,s,26asel,a,72asel,a,76asel,a,80asel,a,84asel,a,95asel,a,59aplottype,3real,20mat,100nsla,s,1esurf,allallsel,allvsel,u,1,6,5vsel,u,7,8vplotasel,s,72aplottype,3real,21mp,mu
18、,101,0.3mat,101nsla,s,1esurf,allallsel,allvsel,u,1,6,5vsel,u,7,8vplotasel,s,76aplottype,3real,22mp,mu,102,0.15mat,102nsla,s,1esurf,allallsel,allvsel,u,1,6,5vsel,u,7,8vplotasel,s,80aplottype,3real,23mp,mu,103,0.3mat,103nsla,s,1esurf,allallsel,allvsel,u,1,6,5vsel,u,7,8vplotasel,s,84aplottype,3real,24m
19、p,mu,104,0.35mat,104nsla,s,1esurf,allallsel,allvsel,u,1,6,5vsel,u,7,8vplotasel,s,59,95,36aplottype,3real,25mp,mu,105,0.3mat,105nsla,s,1esurf,all/psymb,esys,1 allsel gplotsavefinish!求解过程/solu!定义面约束条件asel,s,loc,x,0 !侧面施加对称约束 da,all,symmasel,s,loc,y,0 da,all,symm asel,s,loc,z,42 !底面试加约束 da,all,allallse
20、l,all !土外层施加约束asel,s,86asel,a,94asel,a,82asel,a,78asel,a,74asel,a,70asel,a,25aplot da,all,allALLSEL,ALLantype,static nlgeom,on solcontrol,onautot,onpred,onneqit,50!每个子步允许的最大迭代次数outpr,all,allOUTRES,all,allFirstF=840e3/ALastF=840e3/A!考虑注浆会增加 30的部分Increment=60e3/A*do,Force,FirstF,LastF,Incrementallsel,allvsel,s,1 vsel,a,6 vsel,a,7 aslv,s, asel,s,1, asel,a,29 asel,a,34 aplot sfa,all,pres,Forceallsel,allnsubst,10,50,5solve*enddofinish
