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