加载中…单元SOLID285与单元SOLID187的对比
(2015-11-10 10:44:34)
标签：
ansys |
分类： 软件（ANSYS、ADINA、SAP2000 |
T=0.01时，模型SOLID285土体的竖向位移和应力如下图所示。
T=0.01时，模型SOLID187土体的竖向位移和应力如下图所示。
附件:
1 、samenode02.txt
/filname,samenode02,1
/PREP7
SELTOL,1.0E-8
!-----------geometry parameters--------------- !mm
pilelength
= 18000
!斜向长度
freelength
= 6000
!斜向长度
bondlength
= 12000
!斜向长度
hori_pilelength =
16914.47
!cosa
水平长度
hori_freelength =
5638.16
hori_bondlength =
hori_pilelength-hori_freelength
verti_pilelength =
6156.36
!sina 竖向长度
verti_freelength =
2052.12
verti_bondlength =
verti_pilelength-verti_freelength
d_plate =200
D_pile =600
d_shaft =45.6
BOUDARY_SIZE =4000
DIS_PLATE =3000
!---------inital meshing parameters------------
soil_esize = 800
shaft_esize = 200
!------------loading parameters----------------
shaft_disp = 600
!------------NLAD parameters------------------
skw_ = 0.9 !skewness threshold
freq = 1 !NLAD frequency (1 means check every substep)
start_time = 0 !beginning time for NLAD interval
end_time = 1 !end time for NLAD interval
!------------remeshing parameters---------------
sculpt_layers = 2 !number of sculpting layers
conc_ = 10 !concave patch splitting angle
conv_ = 10 !convex patch splitting angle
edge_ = 10 !edge splitting angle
!---------element types & settings-----------
et,1,285
! for
soil
et,2,285
! for
pile
et,7,285
! for
pile
!******************material properties*****************
MPTEMP,,,,,,,,
MPTEMP,1,0
MP,EX,1,11.2 !soil material properties
MP,NUXY,1,0.4
!poisson's
ratio
MPDATA,DENS,1,,1.9e-9
TB,DP,1,,,
TBMODIF,1,1,0.012
TBMODIF,1,2,13.1
TBMODIF,1,3,0
MPTEMP,,,,,,,,
MPTEMP,1,0
MP,EX,2,279 !pile material properties
MP,NUXY,2,0.25
!poisson's
ratio
MPDATA,DENS,2,,1.9e-9
MPTEMP,,,,,,,,
MPTEMP,1,0
MP,EX,7,2e5 !shaft material properties
MP,NUXY,7,0.2
!poisson's
ratio
MPDATA,DENS,7,,7.85E-009
!------------geometry creation-------------------
/PREP7
CYL4, , ,d_shaft/2,0, ,180,-(pilelength+BOUDARY_SIZE)
!建立杆轴几何体，便于布尔运算
CYL4, , ,d_plate/2,0, ,180,-(pilelength+BOUDARY_SIZE)
!分散板范围内的整个桩体
CYL4, , ,D_pile/2,0, ,180,-(pilelength+BOUDARY_SIZE)
!整个桩体
BLOCK,-1000,1000,0,1000,-(pilelength+BOUDARY_SIZE),0,
BLOCK,-BOUDARY_SIZE,BOUDARY_SIZE,0,BOUDARY_SIZE,-(pilelength+BOUDARY_SIZE),0,
VOVLAP,ALL
wpoff,0,0,-freelength
VSBW,ALL
!用工作面切割
wpoff,0,0,-(bondlength-2*DIS_PLATE)
VSBW,ALL
!生成第一块分散板
wpoff,0,0,-DIS_PLATE
VSBW,ALL
wpoff,0,0,-DIS_PLATE
VSBW,ALL
!生成第三块分散板
CSYS,1
VSEL,R,LOC,X,0,d_shaft/2
VSEL,R,LOC,Z,-freelength,0
VDELE,2, , ,1
!选择不上，只好手工删除钢绞线自由段
ALLS
CSYS,0
WPCSYS,-1,0
!将工作面与整体坐标对齐
wpro,,90.000000,
wpro,,,30
VSBW,ALL
wpro,,,30
VSBW,ALL
wpro,,,30
VSBW,ALL
wpro,,,30
VSBW,ALL
wpro,,,30
VSBW,ALL
ALLS
CSYS,1
VSEL,S,LOC,X,0,d_shaft/2
VSEL,R,LOC,Z,-freelength,-pilelength
CM,SHAFT,VOLU
!定义钢绞线锚固段杆轴组合体
CSYS,1
VSEL,S,LOC,X,D_pile/2,d_shaft/2
VSEL,U,LOC,Z,-pilelength,-(pilelength+BOUDARY_SIZE)
CM,PILE,VOLU
!定义桩体组合体
ALLS
CMSEL,U,SHAFT
CMSEL,U,PILE
CM,SOIL,VOLU
CMSEL,S,SOIL
CSYS,0
VSEL,R,LOC,X,-1000,1000
VSEL,R,LOC,Y,0,1000
CM,INSOIL,VOLU
CMSEL,S,SOIL
CMSEL,U,INSOIL
CM,OUTSOIL,VOLU
ALLS
!-------------soil meshing----------------
type,1
mat,1
!real,1
CMSEL,S,INSOIL
ASLV,S
LSLA,S
CSYS,1
LSEL,R,LOC,X,0,400
!选择范围放大了些
!LSEL,U,LOC,Z,-freelength+100,-100
!LSEL,U,LOC,Z,100,bondlength-100
!LSEL,U,LOC,Z,bondlength+10,bondlength+1000
CM,SMALLSOIL,LINE
LESIZE,ALL,100, , , , , , ,1
CMSEL,S,INSOIL
ASLV,S
LSLA,S
CMSEL,U,SMALLSOIL
CM,MODERNSOIL,LINE
LESIZE,ALL,400, , , , , , ,1
CMSEL,S,SOIL
ASLV,S
LSLA,S
CMSEL,U,SMALLSOIL
CMSEL,U,MODERNSOIL
CM,LONGSOIL,LINE
LESIZE,ALL,1200, , , , , , ,1
ALLSEL
VMESH,SOIL
!-------------pile meshing----------------
type,2
mat,2
!real,2
CMSEL,S,PILE
ASLV,S
LSLA,S
CSYS,1
LSEL,R,LOC,X,0,d_plate/2
CM,INPILE,LINE
LSEL,R,LOC,X,d_shaft/2
CM,SMALLPILE,LINE
LESIZE,ALL,25, , , , , , ,1
CMSEL,S,INPILE
CMSEL,U,SMALLPILE
CM,MIDDLEPILE,LINE
LESIZE,ALL,50, , , , , , ,1
CMSEL,S,PILE
ASLV,S
LSLA,S
CMSEL,U,SMALLPILE
CMSEL,U,MIDDLEPILE
CM,LONGPILE,LINE
LESIZE,ALL,100, , , , , , ,1
ALLS
VMESH,PILE
!-------------shaft meshing----------------
type,7
mat,7
CMSEL,S,SHAFT
ASLV,S
LSLA,S
LESIZE,ALL,100, , , , , , ,1
VMESH,SHAFT
ALLS
!--------------create components-----------
!-----boundary conditions & loading-------
/SOL
CSYS,0
CMSEL,S,SOIL
CMSEL,A,PILE
CMSEL,A,SHAFT
ASLV,S
ASEL,R,LOC,Y,0
NSLA,S,1
DSYM,SYMM,Y, , !施加对称约束
alls
ASEL,S,LOC,Y,BOUDARY_SIZE
NSLA,S,1
D,ALL, , , , , ,UY, , , , ,
alls
ASEL,S,LOC,X,BOUDARY_SIZE
NSLA,S,1
D,ALL, , , , , ,UX, , , , ,
alls
ASEL,S,LOC,X,-BOUDARY_SIZE
NSLA,S,1
D,ALL, , , , , ,UX, , , , ,
alls
ASEL,S,LOC,Z,-(pilelength+BOUDARY_SIZE)
NSLA,S,1
D,ALL, , , , , ,UZ, , , , ,
alls
/SOL
CSYS,0
CMSEL,S,SHAFT
ASLV,S
nsla,s,1
NSEL,R,LOC,Z,-freelength
D,ALL, ,100, , , ,UZ, , , , ,
!位移约束的参考坐标是CSYS0
ALLS
/SOLU
time,1
nsubst,100,10000,10
OUTRES,ERASE
OUTRES,ALL,ALL
CSYS,0
ACEL,0,0,9800,
alls
CNVTOL,F, ,0.08,2, , !2范数
solve
finish
2、samenode03.txt
/filname,samenode03,1
/PREP7
SELTOL,1.0E-8
!-----------geometry parameters--------------- !mm
pilelength
= 18000
!斜向长度
freelength
= 6000
!斜向长度
bondlength
= 12000
!斜向长度
hori_pilelength =
16914.47
!cosa
水平长度
hori_freelength =
5638.16
hori_bondlength =
hori_pilelength-hori_freelength
verti_pilelength =
6156.36
!sina 竖向长度
verti_freelength =
2052.12
verti_bondlength =
verti_pilelength-verti_freelength
d_plate =200
D_pile =600
d_shaft =45.6
BOUDARY_SIZE =4000
DIS_PLATE =3000
!---------inital meshing parameters------------
soil_esize = 800
shaft_esize = 200
!------------loading parameters----------------
shaft_disp = 600
!------------NLAD parameters------------------
skw_ = 0.9 !skewness threshold
freq = 1 !NLAD frequency (1 means check every substep)
start_time = 0 !beginning time for NLAD interval
end_time = 1 !end time for NLAD interval
!------------remeshing parameters---------------
sculpt_layers = 2 !number of sculpting layers
conc_ = 10 !concave patch splitting angle
conv_ = 10 !convex patch splitting angle
edge_ = 10 !edge splitting angle
!---------element types & settings-----------
et,1,187
! for
soil
et,2,187
! for
pile
et,7,187
! for
pile
!******************material properties*****************
MPTEMP,,,,,,,,
MPTEMP,1,0
MP,EX,1,11.2 !soil material properties
MP,NUXY,1,0.4
!poisson's
ratio
MPDATA,DENS,1,,1.9e-9
TB,DP,1,,,
TBMODIF,1,1,0.012
TBMODIF,1,2,13.1
TBMODIF,1,3,0
MPTEMP,,,,,,,,
MPTEMP,1,0
MP,EX,2,279 !pile material properties
MP,NUXY,2,0.25
!poisson's
ratio
MPDATA,DENS,2,,1.9e-9
MPTEMP,,,,,,,,
MPTEMP,1,0
MP,EX,7,2e5 !shaft material properties
MP,NUXY,7,0.2
!poisson's
ratio
MPDATA,DENS,7,,7.85E-009
!------------geometry creation-------------------
/PREP7
CYL4, , ,d_shaft/2,0, ,180,-(pilelength+BOUDARY_SIZE)
!建立杆轴几何体，便于布尔运算
CYL4, , ,d_plate/2,0, ,180,-(pilelength+BOUDARY_SIZE)
!分散板范围内的整个桩体
CYL4, , ,D_pile/2,0, ,180,-(pilelength+BOUDARY_SIZE)
!整个桩体
BLOCK,-1000,1000,0,1000,-(pilelength+BOUDARY_SIZE),0,
BLOCK,-BOUDARY_SIZE,BOUDARY_SIZE,0,BOUDARY_SIZE,-(pilelength+BOUDARY_SIZE),0,
VOVLAP,ALL
wpoff,0,0,-freelength
VSBW,ALL
!用工作面切割
wpoff,0,0,-(bondlength-2*DIS_PLATE)
VSBW,ALL
!生成第一块分散板
wpoff,0,0,-DIS_PLATE
VSBW,ALL
wpoff,0,0,-DIS_PLATE
VSBW,ALL
!生成第三块分散板
CSYS,1
VSEL,R,LOC,X,0,d_shaft/2
VSEL,R,LOC,Z,-freelength,0
VDELE,2, , ,1
!选择不上，只好手工删除钢绞线自由段
ALLS
CSYS,0
WPCSYS,-1,0
!将工作面与整体坐标对齐
wpro,,90.000000,
wpro,,,30
VSBW,ALL
wpro,,,30
VSBW,ALL
wpro,,,30
VSBW,ALL
wpro,,,30
VSBW,ALL
wpro,,,30
VSBW,ALL
ALLS
CSYS,1
VSEL,S,LOC,X,0,d_shaft/2
VSEL,R,LOC,Z,-freelength,-pilelength
CM,SHAFT,VOLU
!定义钢绞线锚固段杆轴组合体
CSYS,1
VSEL,S,LOC,X,D_pile/2,d_shaft/2
VSEL,U,LOC,Z,-pilelength,-(pilelength+BOUDARY_SIZE)
CM,PILE,VOLU
!定义桩体组合体
ALLS
CMSEL,U,SHAFT
CMSEL,U,PILE
CM,SOIL,VOLU
CMSEL,S,SOIL
CSYS,0
VSEL,R,LOC,X,-1000,1000
VSEL,R,LOC,Y,0,1000
CM,INSOIL,VOLU
CMSEL,S,SOIL
CMSEL,U,INSOIL
CM,OUTSOIL,VOLU
ALLS
!模型中杆轴是建3D实体，还是3d梁单元？若是建3D梁单元，则杆轴与桩体之间的接触是否在1/2模型中是否可以成功算出？？
!-------------soil meshing----------------
type,1
mat,1
!real,1
CMSEL,S,INSOIL
ASLV,S
LSLA,S
CSYS,1
LSEL,R,LOC,X,0,400
!选择范围放大了些
!LSEL,U,LOC,Z,-freelength+100,-100
!LSEL,U,LOC,Z,100,bondlength-100
!LSEL,U,LOC,Z,bondlength+10,bondlength+1000
CM,SMALLSOIL,LINE
LESIZE,ALL,100, , , , , , ,1
CMSEL,S,INSOIL
ASLV,S
LSLA,S
CMSEL,U,SMALLSOIL
CM,MODERNSOIL,LINE
LESIZE,ALL,400, , , , , , ,1
CMSEL,S,SOIL
ASLV,S
LSLA,S
CMSEL,U,SMALLSOIL
CMSEL,U,MODERNSOIL
CM,LONGSOIL,LINE
LESIZE,ALL,1200, , , , , , ,1
ALLSEL
VMESH,SOIL
!-------------pile meshing----------------
type,2
mat,2
!real,2
CMSEL,S,PILE
ASLV,S
LSLA,S
CSYS,1
LSEL,R,LOC,X,0,d_plate/2
CM,INPILE,LINE
LSEL,R,LOC,X,d_shaft/2
CM,SMALLPILE,LINE
LESIZE,ALL,25, , , , , , ,1
CMSEL,S,INPILE
CMSEL,U,SMALLPILE
CM,MIDDLEPILE,LINE
LESIZE,ALL,50, , , , , , ,1
CMSEL,S,PILE
ASLV,S
LSLA,S
CMSEL,U,SMALLPILE
CMSEL,U,MIDDLEPILE
CM,LONGPILE,LINE
LESIZE,ALL,100, , , , , , ,1
ALLS
VMESH,PILE
!-------------shaft meshing----------------
type,7
mat,7
CMSEL,S,SHAFT
ASLV,S
LSLA,S
LESIZE,ALL,100, , , , , , ,1
VMESH,SHAFT
ALLS
!--------------create components-----------
!-----boundary conditions & loading-------
/SOL
CSYS,0
CMSEL,S,SOIL
CMSEL,A,PILE
CMSEL,A,SHAFT
ASLV,S
ASEL,R,LOC,Y,0
NSLA,S,1
DSYM,SYMM,Y, , !施加对称约束
alls
ASEL,S,LOC,Y,BOUDARY_SIZE
NSLA,S,1
D,ALL, , , , , ,UY, , , , ,
alls
ASEL,S,LOC,X,BOUDARY_SIZE
NSLA,S,1
D,ALL, , , , , ,UX, , , , ,
alls
ASEL,S,LOC,X,-BOUDARY_SIZE
NSLA,S,1
D,ALL, , , , , ,UX, , , , ,
alls
ASEL,S,LOC,Z,-(pilelength+BOUDARY_SIZE)
NSLA,S,1
D,ALL, , , , , ,UZ, , , , ,
alls
/SOL
CSYS,0
CMSEL,S,SHAFT
ASLV,S
nsla,s,1
NSEL,R,LOC,Z,-freelength
D,ALL, ,100, , , ,UZ, , , , ,
!位移约束的参考坐标是CSYS0
ALLS
/SOLU
time,1
nsubst,100,10000,10
OUTRES,ERASE
OUTRES,ALL,ALL
CSYS,0
ACEL,0,0,9800,
alls
CNVTOL,F, ,0.08,2, , !2范数
solve
finish
喜欢
0
赠金笔
前一篇：三维抗拔桩模型为何无剪应力？