Failed to divide the grid for intersecting faults

Coreform Cubit
#1

When I try to mesh a volume for an intersecting fault, an error occurs:

Command: mesh volume all
ERROR: D:/Trelis/Trelis 16.5/ziduanc.jou (383)
Bounding mesh has connectivity problems.
(Bounding mesh has odd Euler characteristic 3.)
ERROR: D:/Trelis/Trelis 16.5/ziduanc.jou (383)
Inconsistent edge to triangle count relationship in bounding mesh.
ERROR: D:/Trelis/Trelis 16.5/ziduanc.jou (383)
Mesh of Volume 5 is invalid.
ERROR: D:/Trelis/Trelis 16.5/ziduanc.jou (383)
Volume 5 (Volume 5) meshing unsuccessful.
ERROR: D:/Trelis/Trelis 16.5/ziduanc.jou (383)
1 volume(s) did not mesh.

I don’t know the reason, but from my analysis it seems that the error occurs at the intersecting part. How should I modify my commands to mesh correctly?
Any help is greatly appreciated.
ziduanc.jou (16.0 KB)

#2

Hello @xiongfen,
i am sorry for the late reply.

Unfortunately i can’t get a mesh with the spline surface that you created for the volume.

Where did you get your data for the fault creation?

#3

I am currently using a simpler fault model
I obtained the fault data from this paper https://pubs.geoscienceworld.org/ssa/bssa/article/110/4/1688/587386/Kinematics-of-Fault-Slip-Associated-with-the-4-6
And I want to use a mesh size of 200 on the main fault surfaces (24, 25, 26, 27)
An error still occurs

mesh volume all
Error: D:/Trelis/Trelis 16.5/ziduanc.jou (64)
Inconsistency between edges and number of triangles in the boundary mesh.
Error: D:/Trelis/Trelis 16.5/ziduanc.jou (64)
Mesh invalid for body 6.
Error: D:/Trelis/Trelis 16.5/ziduanc.jou (64)
Failed to generate mesh for body 6.
Error: D:/Trelis/Trelis 16.5/ziduanc.jou (64)
1 body failed to be meshed.

The code is as follows:

reset
#top

Contour

create vertex x 0 y 0 z 0.000000e+03
create vertex x -21547.563 y 20951.061 z 0.000000e+03
create vertex x-29150.425 y 33754.393 z 0.000000e+03
create vertex x -44488.000 y 46605.212 z 0.000000e+03
create vertex x 0 y 0 z -17e+03
create vertex x -21630 y 20880 z -17e+03
create vertex x-29280 y 33670 z -17e+03
create vertex x -44660 y 46470 z -17e+03
create vertex x -10787 y 20118.15 z 0.000000e+03
create vertex x -24169.593 y 5255.259 z 0.000000e+03
create vertex x -10786.981 y 20118.15 z -17e+03
create vertex x -24169.593 y 5255.259 z -17e+03
create vertex x -11623.726 y 15352.847 z 0.000000e+03
create vertex x -25516.893 y 29739.643 z 0.000000e+03
create vertex x -11623.726 y 15352.847 z -17e+03
create vertex x -25516.893 y 29739.643 z -17e+03
create curve vertex 1 5
create curve vertex 5 6
create curve vertex 6 7
create curve vertex 7 8
create curve vertex 8 4
create curve vertex 4 3
create curve vertex 3 2
create curve vertex 2 1
create curve vertex 2 6
create curve vertex 3 7
create surface curve 1 2 9 8
create surface curve 9 3 10 7
create surface curve 10 4 5 6
create curve vertex 9 10
create curve vertex 10 12
create curve vertex 12 11
create curve vertex 11 9
create curve vertex 13 14
create curve vertex 14 16
create curve vertex 16 15
create curve vertex 15 13
create surface curve 17 18 19 20
create surface curve 13 14 15 16
delete curve all
unit volume all
brick x 150000 y 150000 z 40000
move volume 6 x -20000 y 20000 z -20000 include_merged
imprint volume all
merge volume all
subtract body 1 from volume 6
surface 24 25 26 27 size 200
volume all scheme tetmesh
mesh volume all

#4

Hello @xiongfen,
this fault is much simpler, we can now decompose the geometry.
You can now do a tet or hex mesh.

#!cubit
reset
create vertex x 0 y 0 z 0.000000e+03
create vertex x -21547.563 y 20951.061 z 0.000000e+03
create vertex x-29150.425 y 33754.393 z 0.000000e+03
create vertex x -44488.000 y 46605.212 z 0.000000e+03
create vertex x 0 y 0 z -17e+03
create vertex x -21630 y 20880 z -17e+03
create vertex x-29280 y 33670 z -17e+03
create vertex x -44660 y 46470 z -17e+03
create vertex x -10787 y 20118.15 z 0.000000e+03
create vertex x -24169.593 y 5255.259 z 0.000000e+03
create vertex x -10786.981 y 20118.15 z -17e+03
create vertex x -24169.593 y 5255.259 z -17e+03
create vertex x -11623.726 y 15352.847 z 0.000000e+03
create vertex x -25516.893 y 29739.643 z 0.000000e+03
create vertex x -11623.726 y 15352.847 z -17e+03
create vertex x -25516.893 y 29739.643 z -17e+03
create curve vertex 1 5
create curve vertex 5 6
create curve vertex 6 7
create curve vertex 7 8
create curve vertex 8 4
create curve vertex 4 3
create curve vertex 3 2
create curve vertex 2 1
create curve vertex 2 6
create curve vertex 3 7
create surface curve 1 2 9 8
create surface curve 9 3 10 7
create surface curve 10 4 5 6
create curve vertex 9 10
create curve vertex 10 12
create curve vertex 12 11
create curve vertex 11 9
create curve vertex 13 14
create curve vertex 14 16
create curve vertex 16 15
create curve vertex 15 13
create surface curve 17 18 19 20
create surface curve 13 14 15 16
delete curve all
brick x 150000 y 150000 z 40000
move volume 6 x -20000 y 20000 z -20000 include_merged
brick x 80000 y 30000 z 40000
rotate volume 7 angle -45 about z
move volume 7 x -20000 y 20000 z -20000 include_merged
subtract volume 7 from volume 6 keep_tool 
block 1 add vol 7
webcut volume 7 with plane zplane offset -17e3 
webcut volume 7 with sheet extended from surface 5
webcut volume 7 with sheet extended from surface 1
webcut volume 7 with sheet extended from surface 4
webcut volume 9  with plane normal to curve 17  fraction 0 from end 
webcut volume 12 with sheet extended from surface 4 
webcut volume 13  with sheet extended from surface 1 2 
webcut volume 9  with sheet extended from surface 2 3 
webcut volume 6 with plane zplane offset -17e3 
imprint body all
delete body all with is_sheet
merge vol all
surface all with z_coord=0 scheme pave
surface all with z_coord=0 size 1000
mesh surface all with z_coord=0

grafik
grafik853×506 48.4 KB 

reset
create vertex x 0 y 0 z 0.000000e+03
create vertex x -21547.563 y 20951.061 z 0.000000e+03
create vertex x-29150.425 y 33754.393 z 0.000000e+03
create vertex x -44488.000 y 46605.212 z 0.000000e+03
create vertex x 0 y 0 z -17e+03
create vertex x -21630 y 20880 z -17e+03
create vertex x-29280 y 33670 z -17e+03
create vertex x -44660 y 46470 z -17e+03
create vertex x -10787 y 20118.15 z 0.000000e+03
create vertex x -24169.593 y 5255.259 z 0.000000e+03
create vertex x -10786.981 y 20118.15 z -17e+03
create vertex x -24169.593 y 5255.259 z -17e+03
create vertex x -11623.726 y 15352.847 z 0.000000e+03
create vertex x -25516.893 y 29739.643 z 0.000000e+03
create vertex x -11623.726 y 15352.847 z -17e+03
create vertex x -25516.893 y 29739.643 z -17e+03
create curve vertex 1 5
create curve vertex 5 6
create curve vertex 6 7
create curve vertex 7 8
create curve vertex 8 4
create curve vertex 4 3
create curve vertex 3 2
create curve vertex 2 1
create curve vertex 2 6
create curve vertex 3 7
create surface curve 1 2 9 8
create surface curve 9 3 10 7
create surface curve 10 4 5 6
create curve vertex 9 10
create curve vertex 10 12
create curve vertex 12 11
create curve vertex 11 9
create curve vertex 13 14
create curve vertex 14 16
create curve vertex 16 15
create curve vertex 15 13
create surface curve 17 18 19 20
create surface curve 13 14 15 16
delete curve all
brick x 150000 y 150000 z 40000
move volume 6 x -20000 y 20000 z -20000 include_merged
brick x 80000 y 30000 z 40000
rotate volume 7 angle -45 about z
move volume 7 x -20000 y 20000 z -20000 include_merged
subtract volume 7 from volume 6 keep_tool 
block 1 add vol 7
webcut volume 7 with plane zplane offset -17e3 
webcut volume 7 with sheet extended from surface 5
webcut volume 7 with sheet extended from surface 1
webcut volume 7 with sheet extended from surface 4
webcut volume 9  with plane normal to curve 17  fraction 0 from end 
webcut volume 12 with sheet extended from surface 4 
webcut volume 13  with sheet extended from surface 1 2 
webcut volume 9  with sheet extended from surface 2 3 
webcut volume 6 with plane zplane offset -17e3 
imprint body all
delete body all with is_sheet
merge vol all
vol all scheme tetmesh
vol all size 3000
vol all in block 1 size 1500
Surface 162 164 148 60 85 84 74 97 160 size 200
mesh Surface 162 164 148 60 85 84 74 97 160 
mesh vol all

grafik
grafik826×526 44.1 KB