ansys中出现单元过度扭曲,怎么解决？
ansys中出现单元过度扭曲,怎么解决？我的命令流
fini
/clear
/PLOPTS,date,0
/prep7
/UNITS,SI
!设置单元类型
ET,1,SOLID45
!设置参数
mp,ex,1,1.5e7 $mp,nuxy,1,0.31 $mp,dens,1,1950 !粉质粘土1参数
mp,ex,2,3e7 $mp,nuxy,2,0.28 $mp,dens,2,1900 !细砂参数
mp,ex,3,2.1e7 $mp,nuxy,3,0.35 $mp,dens,3,1950 !粉质粘土2参数
mp,ex,4,8e7 $mp,nuxy,4,0.3 $mp,dens,4,2180 !卵石参数
mp,ex,5,3e10$mp,nuxy,5,0.2 $mp,dens,5,3140 !PHC桩参数
mp,ex,6,3e10 $mp,nuxy,6,0.2 $mp,dens,6,2500 !钢筋混凝土板参数
mp,ex,7,8e7 $mp,nuxy,7,0.25 $mp,dens,7,2200 !碎石垫层参数
tb,dp,1
tbdata,1,45e3,18.5,0 !粉质粘土1粘聚力，内摩擦角
tb,dp,2
tbdata,1,15e3,32,0 !细砂粘聚力，内摩擦角
tb,dp,3
tbdata,1,40e3,21,0 !粉质粘土2粘聚力，内摩擦角
tb,dp,4
tbdata,1,10e3,45,0 !卵石粘聚力，内摩擦角
K,1,0,0,0
K,2,2,0,0
K,3,2,0,-1.5
K,4,0,0,-1.5
A,1,2,3,4
VOFFST,1,-0.5,, !混凝土板
VOFFST,2,-0.15,, !碎石垫层
VOFFST,7,-12,, !土层1
VOFFST,12,-3,, !土层2
VOFFST,17,-20,, !土层3
VOFFST,22,-9,, !土层3
VOFFST,27,-26,, !土层4
wprota,,-90 !创建与管桩相同直径的实心桩
CYL4,1,0.75,0,0,0.2,360,-70.65
vsel,s,,,1,7,1
vplot
CM,CM_1,volu
CMSEL,s,CM_1
allsel
vsel,s,,,8
vplot
CM,CM_2,volu
CMSEL,s,CM_2
allsel
VSBV,CM_1,CM_2
!布尔运算
CYL4,1,0.75,0.12,0,0.2,360,-0.5 !创建管桩体
wpoff,0,0,-0.5
CYL4,1,0.75,0.12,0,0.2,360,-0.15
wpoff,0,0,-0.15
CYL4,1,0.75,0.12,0,0.2,360,-35
wpoff,0,0,-35
CYL4,1,0.75,0.12,0,0.2,360,-9
wpoff,0,0,-9
CYL4,1,0.75,0.12,0,0.2,360,-26
VGLUE,1,2,14,9,10,11,12 !布尔运算粘结体所有的除桩之外的部分，即土层
VGLUE,12,13,15,4,5
NUMCMP,ALL !压缩编号消除空号
/PNUM,VOLU,1!显示体号码
VPLOT !显示所选体
SAVE
!采用扫掠方法划分网格
VSEL,S,,,7
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,1,0
LSEL,,,,1,4,1
LESIZE,all,,,4
LSEL,,,,81,84,1
LESIZE,all,,,2
VSEL,,,,7
TYPE,1
MAT,6
VSWEEP,7,70,69
allsel
vplot
VSEL,S,,,12
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,1,0
LSEL,,,,13,16,1
LESIZE,all,,,4
LSEL,,,,93,96,1
LESIZE,all,,,2
VSEL,,,,12
TYPE,1
MAT,7
VSWEEP,12,68,70
allsel
vplot
VSEL,S,,,6
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,12,0
LSEL,,,,13,16,1
LESIZE,all,,,4
LSEL,,,,93,96,1
LESIZE,all,,,2
VSEL,,,,6
TYPE,1
MAT,1
VSWEEP,6,68,45
allsel
vplot
VSEL,S,,,10
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,3,0
LSEL,,,,21,24,1
LESIZE,all,,,4
LSEL,,,,69,70,1
LESIZE,all,,,4
VSEL,,,,10
TYPE,1
MAT,2
VSWEEP,10,45,48
allsel
vplot
VSEL,S,,,8
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,20,0
LSEL,,,,29,32,1
LESIZE,all,,,4
LSEL,,,,73,74,1
LESIZE,all,,,4
VSEL,,,,8
TYPE,1
MAT,3
VSWEEP,8,48,42 !划分第5层土层
allsel
vplot
VSEL,S,,,9
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,9,0
LSEL,,,,37,40,1
LESIZE,all,,,4
LSEL,,,,129,132,1
LESIZE,all,,,2
VSEL,,,,9
TYPE,1
MAT,3
VSWEEP,9,42,43 !划分第5层土层
allsel
vplot
VSEL,S,,,11
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,26,0
LSEL,,,,45,48,1
LESIZE,all,,,4
LSEL,,,,121,124,1
LESIZE,all,,,2
VSEL,,,,11
TYPE,1
MAT,4
VSWEEP,11,43,44 !划分第5层土层
allsel
vplot
VSEL,S,,,1
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,1,0
LSEL,,,,81,84,1
LESIZE,all,,,2
LSEL,,,,85,88,1
LESIZE,all,,,2
VSEL,,,,1
TYPE,1
MAT,6
VSWEEP,1,2,1
allsel
vplot
VSEL,S,,,5
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,1,0
LSEL,,,,93,96,1
LESIZE,all,,,2
LSEL,,,,97,100,1
LESIZE,all,,,2
VSEL,,,,5
TYPE,1
MAT,7
VSWEEP,5,27,1
allsel
vplot
VSEL,S,,,3
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,100,0
LSEL,,,,113,116,1
LESIZE,all,,,2
LSEL,,,,109,112,1
LESIZE,all,,,2
VSEL,,,,3
TYPE,1
MAT,5
VSWEEP,3,50,49
allsel
vplot
VSEL,S,,,4
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,9,0
LSEL,,,,121,124,1
LESIZE,all,,,2
LSEL,,,,125,128,1
LESIZE,all,,,2
VSEL,,,,4
TYPE,1
MAT,3
VSWEEP,4,55,56
allsel
vplot
VSEL,S,,,2
EXTOPT,ACLEAR,1
EXTOPT,VSWE,AUTO,0
EXTOPT,ESIZE,26,0
LSEL,,,,121,124,1
LESIZE,all,,,2
LSEL,,,,125,128,1
LESIZE,all,,,2
VSEL,,,,2
TYPE,1
MAT,4
VSWEEP,2,55,61
allsel
vplot
!接触命令流
allsel
ET,2,170! 定义目标单元类型
ET,3,173! 定义接触单元类型
keyopt,3,4,1
keyopt,3,5,1
keyopt,3,9,1 !消除初始渗透
keyopt,3,12,2 !不分开接触，允许滑动，法向不分开
mp,mu,8,0.2 !接触单元的摩擦系数
allsel,all
r,1,,,10,1,0.00001,, !定义单元实常数
rmore,,,,-0.00001,,5 !切向接触刚度
! 设置管桩桩侧目标面
VSEL,S,,,3
ASLV,S
ASEL,S,,,51,52,1
cm,target,area
real,1
TYPE,2
mat,8
NSLA,S,1
ESURF,ALL
esel,s,type,,2
eplot
/psymb,esys,1
allsel
! 创建土层1，建立接触面
ALLSEL,ALL
VSEL,S,,,6
ASLV,S
ASEL,S,,,66,67,1
cm,contact1,area
real,1
TYPE,3
mat,8
NSLA,S,1
ESURF,ALL
esel,s,type,,3
eplot
/psymb,esys,1
allsel
! 创建土层2接触面
ALLSEL,ALL
VSEL,s,,,10
ASLV,S
ASEL,S,,,46,47,1
cm,contact2,area
real,1
TYPE,3
mat,8
NSLA,S,1
ESURF,ALL
esel,s,type,,3
eplot
/psymb,esys,1
allsel
! 创建土层3接触面
allsel,all
VSEL,s,,,8
ASLV,S
ASEL,S,,,40,41,1
cm,contact3,area
real,1
TYPE,3
mat,8
NSLA,S,1
ESURF,ALL
esel,s,type,,3
eplot
/psymb,esys,1
allsel
!（桩端与土的接触不知道用不用设置）
ET,4,170! 定义目标单元类型
ET,5,173! 定义接触单元类型
keyopt,5,4,1
keyopt,5,5,1
keyopt,5,9,1 !消除初始渗透
keyopt,5,12,2 !不分开接触，允许滑动，法向不分开
mp,mu,9,0.2 !接触单元的摩擦系数
r,2,,,1,1,,, !定义单元实常数
rmore,,,,-0.00001,,5 !切向接触刚度
!桩端，建立目标面
allsel,all
VSEL,S,,,3
ASLV,S
ASEL,S,,,49,50,1
cm,target2,area
real,2
TYPE,4
mat,9
NSLA,S,1
ESURF,ALL
esel,s,type,,4
eplot
/psymb,esys,1
allsel
!桩端（上部），建立接触面
ALLSEL,ALL
VSEL,S,,,5
ASLV,S
ASEL,S,,,27
cm,contact4,area
real,2
TYPE,5
mat,9
NSLA,S,1
ESURF,ALL
esel,s,type,,5
eplot
/psymb,esys,1
allsel
!桩端（下部），建立接触面
ALLSEL,ALL
VSEL,S,,,4
ASLV,S
ASEL,S,,,56
cm,contact5,area
real,2
TYPE,5
mat,9
NSLA,S,1
ESURF,ALL
esel,s,type,,5
eplot
/psymb,esys,1
allsel
vgen,2,1,12,1,2
vgen,2,1,24,1,,,-1.5
VGLUE,7,19,43,31,12,24,48,36
VGLUE,12,24,48,36,6,18,42,30
VGLUE,6,18,42,30,10,22,46,34
VGLUE,10,22,46,34,8,20,44,32
VGLUE,8,20,44,32,9,21,45,33
VGLUE,9,21,45,33,11,23,47,35
ALLSEL,ALL !选择所有实体
NUMCMP,ALL !压缩编号消除空号
/solu
asel,s,loc,x,0 ！模型左边界
da,all,all,0
asel,s,loc,x,4 ！模型右边界
da,all,all,0
asel,s,loc,z,0 ！模型前边界
da,all,all,0
asel,s,loc,z,-3 ！模型后边界
da,all,all,0
asel,s,loc,y,-70.65 !模型下边界
da,all,all,0
allsel
!!!施加重力和预加力
/solu
antype,static
nlgeom,on
autots,on ! 打开自动时间步
nropt,full
pred,on !打开时间步长预测器
lnsrch,on
outres,all,all
cncheck,detail !检查接触状态
cncheck,summary
time,1
acel,0,9.8,0,
lswrite,1
time,2
asel,s,loc,y,0
sfa,all,,pres,90000
lswrite,2
lssolve,1,2,1
求高手指点 希望早点整出来啊 无回复，无回复啊
页:
[1]