Question

我这段ANSYS命令流，建模的问题在哪啊，求解时总出现节点位移过大的问题

如题，以下是命令流
/PREP7
!衬砌内轮廓
K,1,0,4. 7 !圆1的圆心,半径为2.9m
K,14,0,7.6 !圆1的圆弧定点
K,15,-2.05,6.75 !圆1的圆弧左起点
K,16,2.05,6.75 !圆1的圆弧右终点
Larc,15,16,1,2.9
K,2,0.05,4.69 !左圆2的圆心，半径3.61
K,21,-3.04,4.87 !左圆2的圆弧左起点
Larc,15,21,2,3.61

K,20,-0.05, 4.69 !右圆形2，半径3.61
K,22,3.04,4.87 !右圆2的圆弧右终点
Larc,16,22,20,3.61

K,3,5.71,3.2 !左圆形3的圆心，半径8.91
K,31,-2.5,-0.25 ! 左圆3的圆弧左起点
Larc,21,31,3,8.91

K,30,-5.71,3.2 !右圆形3的圆心，半径8.91
K,32,2.5,-0.25 !右圆3的圆弧右终点
Larc,22,32,30,8.91

K,4,0,2.86 !圆心4，半径4.3m
K,41,-2.05,-0.92 ! 圆4的圆弧左起点
K,42,2.05,-0.92 ! 圆4的圆弧左起点
Larc,41,42,4,4.3
K,51,-2.05,-0.25 !
K,52,2.05,-0.25 !
L,31,51
L,51,41
L,32,52
L,52,42

!衬砌外轮廓 衬砌厚度为0.3m，
K,114,0,7.9 !圆1的圆弧定点，半径增加0.3为3.2
K,115,-2.26,6.96 !圆1的圆弧左起点
K,116,2.26,6.96 !圆1的圆弧右终点
Larc,115,116,1,3.2
!左圆2的圆心，半径3.61增加至3.91
K,121,-3.34,4.93 !左圆2的圆弧左起点，
Larc,115,121,2,3.91
!右圆2的圆心，半径3.61增加至3.91
K,122,3.34,4.93 !右圆2的圆弧右起点，
Larc,116,122,20,3.91
!左圆3的圆心，半径8.91增加至9.21
K,131,-3.48,3.14 ! 左圆3的圆弧左起点
Larc,121,131,3,9.21
!右圆3的圆心，半径8.91增加至9.21
K,132,3.48,3.14 !左圆3的圆弧左起点
Larc,122,132,30,9.21
!圆心4，半径4.3增加至4.6m
K,141,-2.19,-1.18 ! 圆4的圆弧左起点
K,142,2.19,-1.18 ! 圆4的圆弧左起点
Larc,141,142,4,4.6
K,151,-3.23,-1.18 !
K,152, 3.23,-1.18 !
L,131,151
L,151,141
L,132,152
L,152,142

K,100,32.4,49.5
k,101,-32.4,49.5
k,102,-32.4,-49.5
k,103,32.4,-49.5
l,100,101
l,101,102
l,102,103
l,103,100
K,104,32.4,0.76
K,105,3.34,0.76
K,106,-3.34,0.76
K,107,-32.4,0.76

L,104,105
L,106,107
A,101,100,104,105,132,122,116,115,121,131,106,107
A,107,106,151,141,142,152,105,104,103,102
LDELE,17
LDELE,19
AL,1,2,4,7,8,6,10,9,5,3
A,121,115,116,122,132,105,152,142,141,151,106,131
AOVLAP,3,4

ET,1,PLANE42 !定义围岩1单元类型(500KPa)
ET,2,PLANE42 !定义围岩2单元类型（800KPa）
ET,3,PLANE42

MPTEMP,1,0
MPDATA,EX,1,,2E9 !围岩1弹性模量
MPDATA,PRXY,1,,0.32 !围岩1泊松比
MP,DENS,1,2200
TB,DP,1 !D-P模型
TBDATA,1,1.4E5,32 !围岩1粘结力0.14MPA 内摩擦角32度

MPDATA,EX,2,,2E9
MPDATA,PRXY,2,,0.32
MP,DENS,2,2200
TB,DP,2
TBDATA,2,2.1E6,35 !围岩2粘结力0.21MPA 内摩擦角35度

MPDATA,EX,3,,2E9
MPDATA,PRXY,3,,0.32
MP,DENS,2,2500
TB,DP,3
TBDATA,1,3.5E6,55 !衬砌

ASEL,S,AREA,,1 !选择面1
AATT,1 ! 赋予材料属性1
MSHKEY,0 !自由划分网格
SMRTSIZE,3 !单元大小为3
AMESH,1 !划分面1

ASEL,S,AREA,,2
AATT,2
MSHKEY,0
SMRTSIZE,3
AMESH,2

ASEL,S,AREA,,5
AATT,3
MSHKEY,0
SMRTSIZE,1
AMESH,5

下面是求解的命令流
/SOLU
ALPHAD,3
BETAD,0.001
ANTYPE,TRANS
TIMINT,OFF
ACEL,,-9.8
TIME,1E-6
SOLVE
这时就出问题，告诉我某某点位移超极限了，到底建模问题在哪啊
约束部分在这里
DL,16,2,UX,0
DL,16,2,UY,0
DL,18,2,UY,0
DL,18,2,UX,0
DL,20,2,UX,0
DL,20,2,UY,0
DL,31,2,UX,0
DL,32,2,UX,0
DL,29,1,UX,0
DL,14,1,UX,0
DL,12,1,UX,0
DL,13,1,UX,0
DL,15,1,UX,0
DL,28,1,UX,0 !两侧水平约束，底边固支
DL,6,3,UX,0
DL,6,3,UY,0
DL,8,3,UY,0
DL,8,3,UX,0
DL,7,3,UX,0
DL,7,3,UY,0
DL,10,3,UX,0
DL,10,3,UY,0
DL,9,3,UX,0
DL,9,3,UY,0
DL,2,3,UX,0
DL,4,3,UX,0
DL,3,3,UX,0
DL,5,3,UX,0

Answer 1

选好节点后 *get,l5ncou,node,,count *get,l5nnum,node,,num,min *dim,l5sy,array,l5ncou,2 *do,i,1,l5ncou *get,l5sy(i,1),node,l5nnum,loc,x *get,l5sy(i,2),node,l5nnum,locy l5nnum=ndnext(l5nnum) !知道坐标了，自己算位移就行 *enddo

追问

....根本不是选节点的问题吧，，