加筋挡土墙 ADINA 大牛看看哪儿有错误，里面有详细说明

peter888

====================================================
*定义点坐标，建立地基土，加筋体，面板三部分边角9个点
====================================================
COORDINATES POINT SYSTEM=0
@CLEAR
1 0.00000000000000 30.00000000000000 18.00000000000000 0
2 0.00000000000000 48.0000000000000 18.00000000000000 0
3 0.00000000000000 48.0000000000000 0.0000000000000 0
4 0.00000000000000 0.0000000000000 0.0000000000000 0
5 0.00000000000000 0.0000000000000 18.00000000000000 0
6 0.00000000000000 0.0000000000000 24.00000000000000 0
7 0.00000000000000 30.00000000000000 24.00000000000000 0
8 0.00000000000000 30.100000000000000 24.00000000000000 0
9 0.00000000000000 30.100000000000000 18.00000000000000 0
====================================================
*定义面，1.地基土体；2.加筋体；3.面板。由点生成面时，拾取点按照逆时针方向
====================================================
SURFACE VERTEX NAME=1 P1=2 P2=5 P3=4 P4=3
SURFACE VERTEX NAME=2 P1=7 P2=6 P3=5 P4=1
SURFACE VERTEX NAME=3 P1=8 P2=7 P3=1 P4=9
====================================================
*定义并施加约束，TRANSLATION（平移）
====================================================
FIXITY NAME=Y-FIX
'Y-TRANSLATION'
'OVALIZATION'
FIXITY NAME=Z-FIX
'Z-TRANSLATION'
'OVALIZATION'
FIXBOUNDARY LINES FIXITY=ALL
2  'Y-FIX'
4  'Y-FIX'
6  'Y-FIX'
3  'Z-FIX'
====================================================
*定义质量力，并施加荷载
====================================================
LOAD MASS-PROPORTIONAL NAME=1 MAGNITUD=9.80000000000000,
     AX=0.00000000000000 AY=0.00000000000000 AZ=-1.00000000000000,
     INTERPRE=BODY-FORCE
APPLY-LOAD BODY=0
1  'MASS-PROPORTIONAL' 1  'MODEL' 0 0 1 0.00000000000000 0 -1 0 0 0  ,
'NO' 0.00000000000000 0.00000000000000 1 0
====================================================
*定义材料参数，
MATERIAL（材料）；NU（泊松比）；PHI（内摩擦角）；
PSI（剪胀角：剪胀角是用来表示材料在剪切过程中体积变化率的一个物理量。剪切过程中产生的位移分为法向位移和切向位移，剪胀角的正切值为法向位移同切向位移的比值。）；
COH（粘聚力）；TCUT（抗剪强度）；ALPHA（热膨胀系数－只有弹性材料有）
====================================================
MATERIAL MOHR-COULOMB NAME=1 E=2.40000000000000E+07,
     NU=0.270000000000000 PHI=31.0000000000000 PSI=20.0000000000000,
     COH=28000.0000000000 TCUT=1.00000000000000E+10,
     DENSITY=2000.00000000000 DILATION=YES MDESCRIP='dijitu'
*
MATERIAL MOHR-COULOMB NAME=2 E=2.20000000000000E+07,
     NU=0.2300000000000000 PHI=25.0000000000000 PSI=20.0000000000000,
     COH=12000.0000000000000 TCUT=1.00000000000000E+10,
     DENSITY=1800.00000000000 DILATION=YES MDESCRIP='tiantu'
*
MATERIAL ELASTIC NAME=3 E=3.45000000000000E+11 NU=0.210000000000000,
     DENSITY=2400.00000000000 ALPHA=0.00000000000000 MDESCRIP='mianban'
*
MATERIAL ELASTIC NAME=4 E=2.05000000000000E+11 NU=0.20000000000000,
     DENSITY=7800.00000000000 ALPHA=0.00000000000000 MDESCRIP='jiajin'
====================================================
*划分网格－设定划分边界线，划分密度，网格间距0.5米
====================================================
SUBDIVIDE LINE NAME=1 MODE=DIVISIONS NDIV=96 RATIO=1.00000000000000,
     PROGRESS=GEOMETRIC CBIAS=NO
3
*
SUBDIVIDE LINE NAME=2 MODE=DIVISIONS NDIV=32 RATIO=1.00000000000000,
     PROGRESS=GEOMETRIC CBIAS=NO
4
*
SUBDIVIDE LINE NAME=6 MODE=DIVISIONS NDIV=12 RATIO=1.00000000000000,
     PROGRESS=GEOMETRIC CBIAS=NO
11
8
*
SUBDIVIDE LINE NAME=5 MODE=DIVISIONS NDIV=60 RATIO=1.00000000000000,
     PROGRESS=GEOMETRIC CBIAS=NO
7
*
SUBDIVIDE LINE NAME=10 MODE=DIVISIONS NDIV=1 RATIO=1.00000000000000,
     PROGRESS=GEOMETRIC CBIAS=NO
9
*
====================================================
*定义单元群1：地基土；2：填土；3：面板；4：筋材；
====================================================
EGROUP TWODSOLID NAME=1 SUBTYPE=STRAIN DISPLACE=DEFAULT,
     STRAINS=DEFAULT MATERIAL=1 INT=DEFAULT RESULTS=STRESSES,
     DEGEN=YES FORMULAT=0 STRESSRE=GLOBAL INITIALS=NONE FRACTUR=NO,
     CMASS=DEFAULT STRAIN-F=0 UL-FORMU=DEFAULT PNTGPS=0 NODGPS=0,
     LVUS1=0 LVUS2=0 SED=NO RUPTURE=ADINA INCOMPAT=DEFAULT,
     TIME-OFF=0.00000000000000 POROUS=NO WTMC=1.00000000000000,
     OPTION=NONE DESCRIPT='dijitu' THICKNES=1.00000000000000,
     PRINT=DEFAULT SAVE=DEFAULT TBIRTH=0.00000000000000,
     TDEATH=0.00000000000000 TMC-MATE=1
*
EGROUP TWODSOLID NAME=2 SUBTYPE=STRAIN DISPLACE=DEFAULT,
     STRAINS=DEFAULT MATERIAL=2 INT=DEFAULT RESULTS=STRESSES,
     DEGEN=YES FORMULAT=0 STRESSRE=GLOBAL INITIALS=NONE FRACTUR=NO,
     CMASS=DEFAULT STRAIN-F=0 UL-FORMU=DEFAULT PNTGPS=0 NODGPS=0,
     LVUS1=0 LVUS2=0 SED=NO RUPTURE=ADINA INCOMPAT=DEFAULT,
     TIME-OFF=0.00000000000000 POROUS=NO WTMC=1.00000000000000,
     OPTION=NONE DESCRIPT='tiantu' THICKNES=1.00000000000000,
     PRINT=DEFAULT SAVE=DEFAULT TBIRTH=0.00000000000000,
     TDEATH=0.00000000000000 TMC-MATE=1
*
EGROUP TWODSOLID NAME=3 SUBTYPE=STRAIN DISPLACE=DEFAULT,
     STRAINS=DEFAULT MATERIAL=3 INT=DEFAULT RESULTS=STRESSES,
     DEGEN=YES FORMULAT=0 STRESSRE=GLOBAL INITIALS=NONE FRACTUR=NO,
     CMASS=DEFAULT STRAIN-F=0 UL-FORMU=DEFAULT PNTGPS=0 NODGPS=0,
     LVUS1=0 LVUS2=0 SED=NO RUPTURE=ADINA INCOMPAT=DEFAULT,
     TIME-OFF=0.00000000000000 POROUS=NO WTMC=1.00000000000000,
     OPTION=NONE DESCRIPT='mianban' THICKNES=1.00000000000000,
     PRINT=DEFAULT SAVE=DEFAULT TBIRTH=0.00000000000000,
     TDEATH=0.00000000000000 TMC-MATE=1
*
EGROUP TRUSS NAME=4 SUBTYPE=GENERAL DISPLACE=DEFAULT MATERIAL=4,
     INT=DEFAULT GAPS=NO INITIALS=NONE CMASS=DEFAULT,
     TIME-OFF=0.00000000000000 OPTION=REBAR RB-LINE=1 DESCRIPT='NONE',
     AREA=0.000270000000000000 PRINT=DEFAULT SAVE=DEFAULT,
     TBIRTH=0.00000000000000 TDEATH=0.00000000000000 TMC-MATE=1

*
EGROUP TRUSS NAME=5 SUBTYPE=GENERAL DISPLACE=DEFAULT MATERIAL=4,
     INT=DEFAULT GAPS=NO INITIALS=NONE CMASS=DEFAULT,
     TIME-OFF=0.00000000000000 OPTION=REBAR RB-LINE=2 DESCRIPT='NONE',
     AREA=0.000270000000000000 PRINT=DEFAULT SAVE=DEFAULT,
     TBIRTH=0.00000000000000 TDEATH=0.00000000000000 TMC-MATE=1

*
EGROUP TRUSS NAME=6 SUBTYPE=GENERAL DISPLACE=DEFAULT MATERIAL=4,
     INT=DEFAULT GAPS=NO INITIALS=NONE CMASS=DEFAULT,
     TIME-OFF=0.00000000000000 OPTION=REBAR RB-LINE=3 DESCRIPT='NONE',
     AREA=0.000270000000000000 PRINT=DEFAULT SAVE=DEFAULT,
     TBIRTH=0.00000000000000 TDEATH=0.00000000000000 TMC-MATE=1

*
EGROUP TRUSS NAME=7 SUBTYPE=GENERAL DISPLACE=DEFAULT MATERIAL=4,
     INT=DEFAULT GAPS=NO INITIALS=NONE CMASS=DEFAULT,
     TIME-OFF=0.00000000000000 OPTION=REBAR RB-LINE=4 DESCRIPT='NONE',
     AREA=0.000270000000000000 PRINT=DEFAULT SAVE=DEFAULT,
     TBIRTH=0.00000000000000 TDEATH=0.00000000000000 TMC-MATE=1
*
EGROUP TRUSS NAME=8 SUBTYPE=GENERAL DISPLACE=DEFAULT MATERIAL=4,
     INT=DEFAULT GAPS=NO INITIALS=NONE CMASS=DEFAULT,
     TIME-OFF=0.00000000000000 OPTION=REBAR RB-LINE=5 DESCRIPT='NONE',
     AREA=0.000270000000000000 PRINT=DEFAULT SAVE=DEFAULT,
     TBIRTH=0.00000000000000 TDEATH=0.00000000000000 TMC-MATE=1

*
EGROUP TRUSS NAME=9 SUBTYPE=GENERAL DISPLACE=DEFAULT MATERIAL=4,
     INT=DEFAULT GAPS=NO INITIALS=NONE CMASS=DEFAULT,
     TIME-OFF=0.00000000000000 OPTION=REBAR RB-LINE=6 DESCRIPT='NONE',
     AREA=0.000270000000000000 PRINT=DEFAULT SAVE=DEFAULT,
     TBIRTH=0.00000000000000 TDEATH=0.00000000000000 TMC-MATE=1

*
EGROUP TRUSS NAME=10 SUBTYPE=GENERAL DISPLACE=DEFAULT MATERIAL=4,
     INT=DEFAULT GAPS=NO INITIALS=NONE CMASS=DEFAULT,
     TIME-OFF=0.00000000000000 OPTION=REBAR RB-LINE=7 DESCRIPT='NONE',
     AREA=0.000270000000000000 PRINT=DEFAULT SAVE=DEFAULT,
     TBIRTH=0.00000000000000 TDEATH=0.00000000000000 TMC-MATE=1

*
EGROUP TRUSS NAME=11 SUBTYPE=GENERAL DISPLACE=DEFAULT MATERIAL=4,
     INT=DEFAULT GAPS=NO INITIALS=NONE CMASS=DEFAULT,
     TIME-OFF=0.00000000000000 OPTION=REBAR RB-LINE=8 DESCRIPT='NONE',
     AREA=0.000270000000000000 PRINT=DEFAULT SAVE=DEFAULT,
     TBIRTH=0.00000000000000 TDEATH=0.00000000000000 TMC-MATE=1

*
EGROUP TRUSS NAME=12 SUBTYPE=GENERAL DISPLACE=DEFAULT MATERIAL=4,
     INT=DEFAULT GAPS=NO INITIALS=NONE CMASS=DEFAULT,
     TIME-OFF=0.00000000000000 OPTION=REBAR RB-LINE=9 DESCRIPT='NONE',
     AREA=0.000270000000000000 PRINT=DEFAULT SAVE=DEFAULT,
     TBIRTH=0.00000000000000 TDEATH=0.00000000000000 TMC-MATE=1

*
EGROUP TRUSS NAME=13 SUBTYPE=GENERAL DISPLACE=DEFAULT MATERIAL=4,
     INT=DEFAULT GAPS=NO INITIALS=NONE CMASS=DEFAULT,
     TIME-OFF=0.00000000000000 OPTION=REBAR RB-LINE=10 DESCRIPT='NONE',
     AREA=0.000270000000000000 PRINT=DEFAULT SAVE=DEFAULT,
     TBIRTH=0.00000000000000 TDEATH=0.00000000000000 TMC-MATE=1

*
EGROUP TRUSS NAME=14 SUBTYPE=GENERAL DISPLACE=DEFAULT MATERIAL=4,
     INT=DEFAULT GAPS=NO INITIALS=NONE CMASS=DEFAULT,
     TIME-OFF=0.00000000000000 OPTION=REBAR RB-LINE=11 DESCRIPT='NONE',
     AREA=0.000270000000000000 PRINT=DEFAULT SAVE=DEFAULT,
     TBIRTH=0.00000000000000 TDEATH=0.00000000000000 TMC-MATE=1

*
EGROUP TRUSS NAME=15 SUBTYPE=GENERAL DISPLACE=DEFAULT MATERIAL=4,
     INT=DEFAULT GAPS=NO INITIALS=NONE CMASS=DEFAULT,
     TIME-OFF=0.00000000000000 OPTION=REBAR RB-LINE=12 DESCRIPT='NONE',
     AREA=0.000270000000000000 PRINT=DEFAULT SAVE=DEFAULT,
     TBIRTH=0.00000000000000 TDEATH=0.00000000000000 TMC-MATE=1

====================================================
*定义点坐标，建立了筋带点
====================================================
COORDINATES POINT SYSTEM=0
10 0.00000000000000 30.100000000000000 18.2500000000000000 0
11 0.00000000000000 30.100000000000000 18.7500000000000000 0
12 0.00000000000000 30.100000000000000 19.250000000000000 0
13 0.00000000000000 30.100000000000000 19.7500000000000000 0
14 0.00000000000000 30.100000000000000 20.250000000000000 0
15 0.00000000000000 30.100000000000000 20.750000000000000 0
16 0.00000000000000 30.100000000000000 21.250000000000000 0
17 0.00000000000000 30.100000000000000 21.750000000000000 0
18 0.00000000000000 30.100000000000000 22.250000000000000 0
19 0.00000000000000 30.100000000000000 22.750000000000000 0
20 0.00000000000000 30.100000000000000 23.250000000000000 0
21 0.00000000000000 30.100000000000000 23.750000000000000 0
22 0.00000000000000  21.00000000000000   18.2500000000000000 0
23 0.00000000000000  21.00000000000000   18.750000000000000 0
24 0.00000000000000  21.00000000000000   19.250000000000000 0
25 0.00000000000000  21.00000000000000   19.750000000000000 0
26 0.00000000000000  20.00000000000000   20.250000000000000 0
27 0.00000000000000  20.00000000000000   20.750000000000000 0
28 0.00000000000000  20.0000000000000   21.250000000000000 0
29 0.00000000000000  20.0000000000000   21.750000000000000 0
30 0.00000000000000  20.0000000000000   22.250000000000000 0
31 0.00000000000000  19.0000000000000   22.750000000000000 0
32 0.00000000000000  19.0000000000000   23.250000000000000 0
33 0.00000000000000  19.0000000000000   23.750000000000000 0
====================================================
*定义线，定义筋带
====================================================
LINE STRAIGHT NAME=12 P1=10 P2=22
LINE STRAIGHT NAME=13 P1=11 P2=23
LINE STRAIGHT NAME=14 P1=12 P2=24
LINE STRAIGHT NAME=15 P1=13 P2=25
LINE STRAIGHT NAME=16 P1=14 P2=26
LINE STRAIGHT NAME=17 P1=15 P2=27
LINE STRAIGHT NAME=18 P1=16 P2=28
LINE STRAIGHT NAME=19 P1=17 P2=29
LINE STRAIGHT NAME=20 P1=18 P2=30
LINE STRAIGHT NAME=21 P1=19 P2=31
LINE STRAIGHT NAME=22 P1=20 P2=32
LINE STRAIGHT NAME=23 P1=21 P2=33
====================================================
*筋带划分网格－设定划分边界线，划分密度
====================================================
SUBDIVIDE LINE NAME=12 MODE=DIVISIONS NDIV=18 RATIO=1.00000000000000,
     PROGRESS=GEOMETRIC CBIAS=NO
@CLEAR
12
13
14
15
*
SUBDIVIDE LINE NAME=16 MODE=DIVISIONS NDIV=20 RATIO=1.00000000000000,
     PROGRESS=GEOMETRIC CBIAS=NO
@CLEAR
16
17
18
19
20
*
SUBDIVIDE LINE NAME=21 MODE=DIVISIONS NDIV=22 RATIO=1.00000000000000,
     PROGRESS=GEOMETRIC CBIAS=NO
@CLEAR
21
22
23
====================================================
*定义筋带为rebar单元
====================================================
REBAR-LINE NAME=1
12
REBAR-LINE NAME=2
13
REBAR-LINE NAME=3
14
REBAR-LINE NAME=4
15
REBAR-LINE NAME=5
16
REBAR-LINE NAME=6
17
REBAR-LINE NAME=7
18
REBAR-LINE NAME=8
19
REBAR-LINE NAME=9
20
REBAR-LINE NAME=10
21
REBAR-LINE NAME=11
22
REBAR-LINE NAME=12
23
====================================================
*根据定义的网格密度，划分网格，采用自由网格划分的方法
====================================================
*
GSURFACE NODES=4 PATTERN=AUTOMATIC NCOINCID=BOUNDARIES NCEDGE=1234,
     NCVERTEX=1234 NCTOLERA=1.00000000000000E-05 SUBSTRUC=0 GROUP=1,
     PREFSHAP=AUTOMATIC MESHING=MAPPED SMOOTHIN=NO DEGENERA=NO,
     COLLAPSE=NO MIDNODES=CURVED METHOD=ADVFRONT FLIP=NO
@CLEAR
1
@
*
GSURFACE NODES=4 PATTERN=AUTOMATIC NCOINCID=BOUNDARIES NCEDGE=1234,
     NCVERTEX=1234 NCTOLERA=1.00000000000000E-05 SUBSTRUC=0 GROUP=2,
     PREFSHAP=AUTOMATIC MESHING=MAPPED SMOOTHIN=NO DEGENERA=NO,
     COLLAPSE=NO MIDNODES=CURVED METHOD=ADVFRONT FLIP=NO
@CLEAR
2
*
GSURFACE NODES=4 PATTERN=AUTOMATIC NCOINCID=BOUNDARIES NCEDGE=1234,
     NCVERTEX=1234 NCTOLERA=1.00000000000000E-05 SUBSTRUC=0 GROUP=3,
     PREFSHAP=AUTOMATIC MESHING=MAPPED SMOOTHIN=NO DEGENERA=NO,
     COLLAPSE=NO MIDNODES=CURVED METHOD=ADVFRONT FLIP=NO
@CLEAR
3
*
GSURFACE NODES=4 PATTERN=AUTOMATIC NCOINCID=BOUNDARIES NCEDGE=1234,
     NCVERTEX=1234 NCTOLERA=1.00000000000000E-05 SUBSTRUC=0 GROUP=3,
     PREFSHAP=AUTOMATIC MESHING=MAPPED SMOOTHIN=NO DEGENERA=NO,
     COLLAPSE=NO MIDNODES=CURVED METHOD=ADVFRONT FLIP=NO
====================================================
*施加初始地应力
====================================================
INITIAL STRAINS
....初始地应力太多省去
====================================================
*打开大位移开关
====================================================
KINEMATICS DISPLACE=LARGE STRAINS=SMALL UL-FORMU=DEFAULT PRESSURE=NO,
     INCOMPAT=AUTOMATIC RIGIDLIN=NO
====================================================
*由于初始应力场是加载在单元坐标系下，但给定的应力是基于整体坐标系下，所以要专门定义一个
初始应变轴方向的正交坐标系轴， 步骤：model-orthotropic axes systems-define
====================================================
AXES CONSTANT NAME=1 AX=1.00000000000000 AY=0.00000000000000,
     AZ=0.00000000000000 BX=0.00000000000000 BY=1.00000000000000,
     BZ=0.00000000000000
====================================================
*（接上面定义的正交坐标系轴），制定局部坐标系的关系，施加在零时刻的模型上面
====================================================
SET-AXES-STR SURFACES
1 1 2 3
2 1 2 3
====================================================
*每个单元节点施加地应力－－－此部分可参考地应力施加部分解释
====================================================
EGROUP TWODSOLID NAME=1 SUBTYPE=STRAIN DISPLACE=DEFAULT,
     STRAINS=DEFAULT MATERIAL=1 INT=DEFAULT RESULTS=STRESSES,
     DEGEN=YES FORMULAT=0 STRESSRE=GLOBAL INITIALS=NODAL FRACTUR=NO,
     CMASS=DEFAULT STRAIN-F=0 UL-FORMU=DEFAULT PNTGPS=0 NODGPS=0,
     LVUS1=0 LVUS2=0 SED=NO RUPTURE=ADINA INCOMPAT=DEFAULT,
     TIME-OFF=0.00000000000000 POROUS=NO WTMC=1.00000000000000,
     OPTION=NONE DESCRIPT='dijitu' THICKNES=1.00000000000000,
     PRINT=DEFAULT SAVE=DEFAULT TBIRTH=0.00000000000000,
     TDEATH=0.00000000000000 TMC-MATE=1
*
EGROUP TWODSOLID NAME=2 SUBTYPE=STRAIN DISPLACE=DEFAULT,
     STRAINS=DEFAULT MATERIAL=2 INT=DEFAULT RESULTS=STRESSES,
     DEGEN=YES FORMULAT=0 STRESSRE=GLOBAL INITIALS=NODAL FRACTUR=NO,
     CMASS=DEFAULT STRAIN-F=0 UL-FORMU=DEFAULT PNTGPS=0 NODGPS=0,
     LVUS1=0 LVUS2=0 SED=NO RUPTURE=ADINA INCOMPAT=DEFAULT,
     TIME-OFF=0.00000000000000 POROUS=NO WTMC=1.00000000000000,
     OPTION=NONE DESCRIPT='tiantu' THICKNES=1.00000000000000,
     PRINT=DEFAULT SAVE=DEFAULT TBIRTH=0.00000000000000,
     TDEATH=0.00000000000000 TMC-MATE=1
*
EGROUP TWODSOLID NAME=3 SUBTYPE=STRAIN DISPLACE=DEFAULT,
     STRAINS=DEFAULT MATERIAL=3 INT=DEFAULT RESULTS=STRESSES,
     DEGEN=YES FORMULAT=0 STRESSRE=GLOBAL INITIALS=NONE FRACTUR=NO,
     CMASS=DEFAULT STRAIN-F=0 UL-FORMU=DEFAULT PNTGPS=0 NODGPS=0,
     LVUS1=0 LVUS2=0 SED=NO RUPTURE=ADINA INCOMPAT=DEFAULT,
     TIME-OFF=0.00000000000000 POROUS=NO WTMC=1.00000000000000,
     OPTION=NONE DESCRIPT='mianban' THICKNES=1.00000000000000,
     PRINT=DEFAULT SAVE=DEFAULT TBIRTH=0.00000000000000,
     TDEATH=0.00000000000000 TMC-MATE=1
*
EGROUP TRUSS NAME=4 SUBTYPE=GENERAL DISPLACE=DEFAULT MATERIAL=4,
     INT=DEFAULT GAPS=NO INITIALS=NONE CMASS=DEFAULT,
     TIME-OFF=0.00000000000000 OPTION=REBAR RB-LINE=1 DESCRIPT='NONE',
     AREA=0.000270000000000000 PRINT=DEFAULT SAVE=DEFAULT,
     TBIRTH=0.00000000000000 TDEATH=0.00000000000000 TMC-MATE=1

peter888

我计算出来的强顶最大位移0.75mm，显然不对，大家看看是哪儿的问题

peter888

整个计算过程没问题

houandy

楼主，你可以把你要解决问题的情况介绍下吗？

houandy

筋带不能用rebar单元，楼主也没有定义筋带和面板，以及面板和土体的接触问题