共同分析ANSYS网格划分 致命错误

hahahuoshen

大家来分析下

hahahuoshen

大家来说说这是怎么回事？

hahahuoshen

!  ——————————————————————————————
!  |         =====某铜矿采空区稳定性数值模拟研究 =====      |
!  |                       矿体与围岩模型                      |
!  ——————————————————————————————
!**********************************************************
!此模型的规格是：300*200*100
!坐标原点
!(0,0,0)为模型原点,x在33线,y为矿石的走向线上,z的0点位海平面+0m
!这个坐标的0点在19线和3129KM的交点处
!布尔运算容差为0.1
!**********************************************************
! 矿体模型
/prep7
/VIEW,1,1,1,1   
/ANG,1  
/REP,FAST  

! 1号开挖模型，最原始的矿块一共三个几何模型add得来的——-280中段
B=21000
M=2
N=12.5
P=0
K, B+1011,152.9+M ,56.1+N , 58.1+P
K, B+1012,168.9+M, 56.1+N, 58.1+P
K, B+1013, 168.9 +M, 59.1+N, 58.1+P
K, B+1014, 152.9 +M, 59.1+N, 58.1+P
K, B+1015, 152.9 +M, 56.1 +N, 60.8+P
K, B+1016, 168.9+M, 56.1+N, 60.8+P
K, B+1017, 168.9+M, 59.1 +N, 60.8+P
K, B+1018, 152.9+M, 59.1+N, 60.8+P
V , B+1011, B+ 1012, B+1013, B+1014, B+ 1015, B+ 1016, B+1017, B+1018

K, B+1035, 152.9+M , 56.1 +N, 60.8+P
K, B+1036, 168.9+M, 56.1+N, 60.8+P
K, B+1037, 168.9 +M, 59.1+N , 60.8+P
K, B+1038, 152.9+M, 59.1+N, 60.8+P
K, B+1031,152.9+M,50.5+N, 67.575+P
K, B+1032,168.9+M, 50.5+N, 67.575+P
K, B+1033, 168.9 +M, 63+N, 64.575+P
K, B+1034, 152.9 +M, 63+N, 64.575+P
V , B+1035, B+ 1036, B+1037, B+1038, B+1031, B+1032, B+ 1033, B+ 1034

K, B+1021,152.9+M ,50.5 +N, 67.575+P
K, B+1022,168.9+M, 50.5+N, 67.575+P
K, B+1023, 168.9+M , 63+N, 64.575+P
K, B+1024, 152.9 +M, 63+N, 64.575+P
K, B+1025, 152.9+M , 50.5 +N, 98.1+P
K, B+1026, 168.9+M, 50.5+N, 98.1+P
K, B+1027, 168.9 +M, 63+N , 98.1+P
K, B+1028, 152.9+M, 63+N, 98.1+P
V , B+1021, B+1022, B+1023, B+1024, B+1025, B+1026, B+1027, B+1028

Vadd,1,2,3

! 2号开挖模型一共三个几何模型add得来的，它于1号相连——-280中段
B=22000
M=2
N=0
P=0
K, B+1011,152.9+M ,53.4+N , 58.1+P
K, B+1012,168.9+M, 53.4+N, 58.1+P
K, B+1013, 168.9 +M, 56.4+N, 58.1+P
K, B+1014, 152.9 +M, 56.4+N, 58.1+P
K, B+1015, 152.9 +M, 53.4 +N, 60.8+P
K, B+1016, 168.9+M, 53.4+N, 60.8+P
K, B+1017, 168.9+M, 56.4 +N, 60.8+P
K, B+1018, 152.9+M, 56.4+N, 60.8+P
V , B+1011, B+ 1012, B+1013, B+1014, B+ 1015, B+ 1016, B+1017, B+1018,

K, B+1021, 152.9+M , 53.4 +N, 60.8+P
K, B+1022, 168.9+M, 53.4+N, 60.8+P
K, B+1023, 168.9 +M, 56.4+N , 60.8+P
K, B+1024, 152.9+M, 56.4+N, 60.8+P
K, B+1025,152.9+M,50.5+N, 64.575+P
K, B+1026,168.9+M, 50.5+N, 64.575+P
K, B+1027, 168.9 +M, 63+N, 67.575+P
K, B+1028, 152.9 +M, 63+N, 67.575+P
V , B+1021, B+ 1022, B+1023, B+1024, B+1025, B+1026, B+ 1027, B+ 1028

K, B+1031,152.9+M ,50.5 +N, 64.575+P
K, B+1032,168.9+M, 50.5+N, 64.575+P
K, B+1033, 168.9+M , 63+N, 67.575+P
K, B+1034, 152.9 +M, 63+N, 67.575+P
K, B+1035, 152.9+M , 50.5 +N, 98.1+P
K, B+1036, 168.9+M, 50.5+N, 98.1+P
K, B+1037, 168.9 +M, 63+N , 98.1+P
K, B+1038, 152.9+M, 63+N, 98.1+P
V , B+1031, B+1032, B+1033, B+1034, B+1035, B+1036, B+1037, B+1038

Vadd,1,2,3


! 3号开挖模型，最原始的矿块一共三个几何模型add得来的——-280中段
B=23000
M=21
N=12.5
P=0
K, B+1011,152.9+M ,56.1+N , 58.1+P
K, B+1012,168.9+M, 56.1+N, 58.1+P
K, B+1013, 168.9 +M, 59.1+N, 58.1+P
K, B+1014, 152.9 +M, 59.1+N, 58.1+P
K, B+1015, 152.9 +M, 56.1 +N, 60.8+P
K, B+1016, 168.9+M, 56.1+N, 60.8+P
K, B+1017, 168.9+M, 59.1 +N, 60.8+P
K, B+1018, 152.9+M, 59.1+N, 60.8+P
V , B+1011, B+ 1012, B+1013, B+1014, B+ 1015, B+ 1016, B+1017, B+1018

K, B+1035, 152.9+M , 56.1 +N, 60.8+P
K, B+1036, 168.9+M, 56.1+N, 60.8+P
K, B+1037, 168.9 +M, 59.1+N , 60.8+P
K, B+1038, 152.9+M, 59.1+N, 60.8+P
K, B+1031,152.9+M,50.5+N, 67.575+P
K, B+1032,168.9+M, 50.5+N, 67.575+P
K, B+1033, 168.9 +M, 63+N, 64.575+P
K, B+1034, 152.9 +M, 63+N, 64.575+P
V , B+1035, B+ 1036, B+1037, B+1038, B+1031, B+1032, B+ 1033, B+ 1034

K, B+1021,152.9+M ,50.5 +N, 67.575+P
K, B+1022,168.9+M, 50.5+N, 67.575+P
K, B+1023, 168.9+M , 63+N, 64.575+P
K, B+1024, 152.9 +M, 63+N, 64.575+P
K, B+1025, 152.9+M , 50.5 +N, 98.1+P
K, B+1026, 168.9+M, 50.5+N, 98.1+P
K, B+1027, 168.9 +M, 63+N , 98.1+P
K, B+1028, 152.9+M, 63+N, 98.1+P
V , B+1021, B+1022, B+1023, B+1024, B+1025, B+1026, B+1027, B+1028



Vadd,1,2,3
! 4号开挖模型一共三个几何模型add得来的，它于1号相连——-280中段
B=24000
M=21
N=0
P=0
K, B+1011,152.9+M ,53.4+N , 58.1+P
K, B+1012,168.9+M, 53.4+N, 58.1+P
K, B+1013, 168.9 +M, 56.4+N, 58.1+P
K, B+1014, 152.9 +M, 56.4+N, 58.1+P
K, B+1015, 152.9 +M, 53.4 +N, 60.8+P
K, B+1016, 168.9+M, 53.4+N, 60.8+P
K, B+1017, 168.9+M, 56.4 +N, 60.8+P
K, B+1018, 152.9+M, 56.4+N, 60.8+P
V , B+1011, B+ 1012, B+1013, B+1014, B+ 1015, B+ 1016, B+1017, B+1018,

K, B+1021, 152.9+M , 53.4 +N, 60.8+P
K, B+1022, 168.9+M, 53.4+N, 60.8+P
K, B+1023, 168.9 +M, 56.4+N , 60.8+P
K, B+1024, 152.9+M, 56.4+N, 60.8+P
K, B+1025,152.9+M,50.5+N, 64.575+P
K, B+1026,168.9+M, 50.5+N, 64.575+P
K, B+1027, 168.9 +M, 63+N, 67.575+P
K, B+1028, 152.9 +M, 63+N, 67.575+P
V , B+1021, B+ 1022, B+1023, B+1024, B+1025, B+1026, B+ 1027, B+ 1028

K, B+1031,152.9+M ,50.5 +N, 64.575+P
K, B+1032,168.9+M, 50.5+N, 64.575+P
K, B+1033, 168.9+M , 63+N, 67.575+P
K, B+1034, 152.9 +M, 63+N, 67.575+P
K, B+1035, 152.9+M , 50.5 +N, 98.1+P
K, B+1036, 168.9+M, 50.5+N, 98.1+P
K, B+1037, 168.9 +M, 63+N , 98.1+P
K, B+1038, 152.9+M, 63+N, 98.1+P
V , B+1031, B+1032, B+1033, B+1034, B+1035, B+1036, B+1037, B+1038


Vadd,1,2,3


! 1号开挖模型，——-235中段
B=121000
M=-48
N=12.5
P=50

K, B+1041,152.9+M ,50.5 +N, 58.1+P    !新加入的体
K, B+1042,168.9+M, 50.5+N, 58.1+P
K, B+1043, 168.9+M , 63+N, 58.1+P
K, B+1044, 152.9 +M, 63+N, 58.1+P
K, B+1045, 152.9+M , 50.5 +N, 98.1+P
K, B+1046, 168.9+M, 50.5+N, 98.1+P
K, B+1047, 168.9 +M, 63+N , 98.1+P
K, B+1048, 152.9+M, 63+N, 98.1+P
V , B+1041, B+1042, B+1043, B+1044, B+1045, B+1046, B+1047, B+1048

Vadd,1,2,3
! 2号开挖模型一共三个几何模型add得来的，它于1号相连——-235中段
B=122000
M=-48
N=0
P=50
K, B+1041,152.9+M ,50.5 +N, 58.1+P    !新加入的体
K, B+1042,168.9+M, 50.5+N, 58.1+P
K, B+1043, 168.9+M , 63+N, 58.1+P
K, B+1044, 152.9 +M, 63+N, 58.1+P
K, B+1045, 152.9+M , 50.5 +N, 98.1+P
K, B+1046, 168.9+M, 50.5+N, 98.1+P
K, B+1047, 168.9 +M, 63+N , 98.1+P
K, B+1048, 152.9+M, 63+N, 98.1+P
V , B+1041, B+1042, B+1043, B+1044, B+1045, B+1046, B+1047, B+1048

Vadd,1,2,3


! 3号开挖模型，最原始的矿块一共三个几何模型add得来的——-235中段
B=123000
M=-28
N=12.5
P=50
K, B+1041,152.9+M ,50.5 +N, 58.1+P    !新加入的体
K, B+1042,168.9+M, 50.5+N, 58.1+P
K, B+1043, 168.9+M , 63+N, 58.1+P
K, B+1044, 152.9 +M, 63+N, 58.1+P
K, B+1045, 152.9+M , 50.5 +N, 98.1+P
K, B+1046, 168.9+M, 50.5+N, 98.1+P
K, B+1047, 168.9 +M, 63+N , 98.1+P
K, B+1048, 152.9+M, 63+N, 98.1+P
V , B+1041, B+1042, B+1043, B+1044, B+1045, B+1046, B+1047, B+1048

Vadd,1,2,3

! 4号开挖模型一共三个几何模型add得来的，它于1号相连——-235中段
B=124000
M=-28
N=0
P=50
K, B+1041,152.9+M ,50.5 +N, 58.1+P    !新加入的体
K, B+1042,168.9+M, 50.5+N, 58.1+P
K, B+1043, 168.9+M , 63+N, 58.1+P
K, B+1044, 152.9 +M, 63+N, 58.1+P
K, B+1045, 152.9+M , 50.5 +N, 98.1+P
K, B+1046, 168.9+M, 50.5+N, 98.1+P
K, B+1047, 168.9 +M, 63+N , 98.1+P
K, B+1048, 152.9+M, 63+N, 98.1+P
V , B+1041, B+1042, B+1043, B+1044, B+1045, B+1046, B+1047, B+1048

Vadd,1,2,3





!模型 ————上盘围岩  V1
K,2031,275 ,-50 ,0
K,2032,300, -50,0
K,2033,300, 150,0
K,2034,275 ,150 ,0
K,2035,75, -50 ,200
K,2036,300, -50,200
K,2037,300, 150,200
K,2038,75,150,200
V , 2031, 2032, 2033, 2034, 2035, 2036, 2037, 2038

!模型————矿体 V12
K,2021,225 ,-50 ,0
K,2022,275, -50,0
K,2023,275, 150, 0
K,2024,225 ,150 , 0
K,2025,25 , -50 ,200
K,2026,75 , -50,200
K,2027,75 , 150,200
K,2028,25 ,150,200
V , 2021, 2022, 2023, 2024, 2025, 2026, 2027, 2028

!模型————下盘围岩V15
K,2011,0, -50, 0
K,2012,225, -50, 0
K,2013,225, 150, 0
K,2014,0,150 ,0
K,2015,0, -50 ,200
K,2016,25 , -50,200
K,2017,25, 150,200
K,2018,0,150,200
V , 2011, 2012, 2013, 2014, 2015, 2016, 2017, 2018



!模型 ————外围大模型 1
K,2041,-100 ,-150 ,-100
K,2042,400, -150, -100
K,2043,400, 250, -100
K,2044,-100 ,250 , -100
K,2045,-100, -150 , 300
K,2046,400, -150, 300
K,2047,400, 250, 300
K,2048,-100,250, 300
V , 2041, 2042, 2043, 2044, 2045, 2046, 2047, 2048

Voverlap,all






et,1,45
!上盘
MP,EX,1,32E9
MP,PRXY,1,0.26
MP,DENS,1,271
!矿体
MP,EX,2,40E9
MP,PRXY,2,0.3
MP,DENS,2,3500
TB,DP,2
TBDATA,1,22e6,36
!下盘
MP,EX,3,53E9
MP,PRXY,3,0.3
MP,DENS,3,2700
TB,DP,3
TBDATA,1,20e6,40
!1比4充填体
MP,EX,4,0.69E9
MP,PRXY,4,0.14
MP,DENS,4,2000
TB,DP,4
TBDATA,1,0.6e6,18
!1比8充填体
MP,EX,5,0.5E9
MP,PRXY,5,0.07
MP,DENS,5,2000
TB,DP,5
TBDATA,1,0.5e6,18

!1比16充填体
MP,EX,6,0.55E9
MP,PRXY,6,0.07
MP,DENS,6,2000
TB,DP,6
TBDATA,1,0.5e6,18

!1比24充填体
MP,EX,7,0.65E9
MP,PRXY,7,0.07
MP,DENS,7,2000
TB,DP,7
TBDATA,1,0.5e6,18


!1比24充填体
MP,EX,8,0.655E9
MP,PRXY,8,0.075
MP,DENS,8,2000
TB,DP,8
TBDATA,1,0.5e6,18

!1比24充填体
MP,EX,9,0.6555E9
MP,PRXY,9,0.0765
MP,DENS,9,2000
TB,DP,9
TBDATA,1,0.5e6,18

!1比24充填体
MP,EX,10,0.6525E9
MP,PRXY,10,0.0725
MP,DENS,10,2000
TB,DP,10
TBDATA,1,0.5e6,18

eborn

划分网格时雅克比矩阵值为负值，你可以尝试两种解决方法，
（1）把你的几何体剖分的比较规则点，划分；
（2）更换单元类型再试试。

sijieshenhua

没遇到过啊

tanhao365

单元划分的不好,不规则的形状经常出现这问题,最好的结果是在切分一下让模型过度 平滑一点,不要有切角之类的地方,还有一个方法就是单元质量检查,然后选取质量不好的单元,删除.虽然这么搞不科学,但是据我做的几个例子这些小单元没了在岩土中对结果影响也不怎么大,放心删除就行

heanjiang

出现负网格的了吧