mesh transition question with python API
hello see the attachment for my question. I want to have the 4 bottom volumic parts meshed by extrusion with hexaeders. But the extrude method:
gmsh.model.geo.extrude( [(2,11),(2,12),(2,13),(2,14),], length,0,0, numElements=[1,1,1], heights=[ 0.2, 0.5, 1.], recombine=True)
does not work here (meshing fails with "Exception: Could not find extruded node (1.6, -1.732050807568877, -1) in surface 87" ) and would not create coincident nodes along the lateral surfaces.
it seems that it fails because gmsh "wants" to avoid non-coincident nodes on the lateral surfaces. I would have to calculate by hand the values in heights=[] ( the above [ 0.2, 0.5, 1.] is wrong because only estimates). Calculating these values seems impossible because the mesh biasing of the hexahedral part has been on other lines set with something like
setTransfiniteCurve( l1, nb_of_nodes, meshType="Progression", coef=2)
and this methods has a "built-in" rounding which is not transparent to me.
My supposition is that for the extrusion instead of
numElements=[1,1,1], heights=[ 0.2, 0.5, 1.]
there should be a method like the above "nb_of_nodes, meshType="Progression", coef=2"
Does this exist? i didn't find it in the doc (gmsh.py). Or how can this task be resolved?
Many thanks
Johannes_ACKVmesh_transition_problem.zip