1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138
|
# ------------------------------------------------------------------------------
#
# Gmsh Julia tutorial 3
#
# Extruded meshes, ONELAB parameters, options
#
# ------------------------------------------------------------------------------
import gmsh
gmsh.initialize()
function createGeometryAndMesh()
# Clear all models and create a new one
gmsh.clear()
gmsh.model.add("t3")
# Copied from `t1.jl'...
lc = 1e-2
gmsh.model.geo.addPoint(0, 0, 0, lc, 1)
gmsh.model.geo.addPoint(.1, 0, 0, lc, 2)
gmsh.model.geo.addPoint(.1, .3, 0, lc, 3)
gmsh.model.geo.addPoint(0, .3, 0, lc, 4)
gmsh.model.geo.addLine(1, 2, 1)
gmsh.model.geo.addLine(3, 2, 2)
gmsh.model.geo.addLine(3, 4, 3)
gmsh.model.geo.addLine(4, 1, 4)
gmsh.model.geo.addCurveLoop([4, 1, -2, 3], 1)
gmsh.model.geo.addPlaneSurface([1], 1)
gmsh.model.geo.synchronize()
gmsh.model.addPhysicalGroup(1, [1, 2, 4], 5)
gmsh.model.addPhysicalGroup(2, [1], -1, "My surface")
# As in `t2.jl', we plan to perform an extrusion along the z axis. But
# here, instead of only extruding the geometry, we also want to extrude the
# 2D mesh. This is done with the same `extrude()' function, but by
# specifying element 'Layers' (2 layers in this case, the first one with 8
# subdivisions and the second one with 2 subdivisions, both with a height of
# h/2). The number of elements for each layer and the (end) height of each
# layer are specified in two vectors:
h = 0.1
angle = 90.
ov = gmsh.model.geo.extrude([(2, 1)], 0, 0, h, [8, 2], [0.5, 1])
# The extrusion can also be performed with a rotation instead of a
# translation, and the resulting mesh can be recombined into prisms (we use
# only one layer here, with 7 subdivisions). All rotations are specified by
# an an axis point (-0.1, 0, 0.1), an axis direction (0, 1, 0), and a
# rotation angle (-Pi/2):
ov = gmsh.model.geo.revolve([(2, 28)], -0.1, 0, 0.1, 0, 1, 0, -pi / 2, [7])
# Using the built-in geometry kernel, only rotations with angles < Pi are
# supported. To do a full turn, you will thus need to apply at least 3
# rotations. The OpenCASCADE geometry kernel does not have this limitation.
# A translation (-2 * h, 0, 0) and a rotation ((0, 0.15, 0.25), (1, 0, 0),
# angle * Pi / 180) can also be combined to form a "twist". The last
# (optional) argument for the extrude() and twist() functions specifies
# whether the extruded mesh should be recombined or not. The `angle'
# parameter is retrieved from the ONELAB database (it can be set
# interactively in the GUI -- see below):
angle = gmsh.onelab.getNumber("Parameters/Twisting angle")[1]
ov = gmsh.model.geo.twist([(2, 50)], 0, 0.15, 0.25, -2 * h, 0, 0, 1, 0, 0,
angle * pi / 180., [10], [], true)
gmsh.model.geo.synchronize()
# All the extrusion functions return a vector of extruded entities: the
# "top" of the extruded surface (in `ov[1]'), the newly created volume (in
# `ov[2]') and the tags of the lateral surfaces (in `ov[3]', `ov[4]', ...).
# We can then define a new physical volume (with tag 101) to group all the
# elementary volumes:
gmsh.model.addPhysicalGroup(3, [1, 2, ov[2][2]], 101)
gmsh.model.mesh.generate(3)
gmsh.write("t3.msh")
end
# Let us now change some options... Since all interactive options are accessible
# through the API, we can for example make point tags visible or redefine some
# colors:
gmsh.option.setNumber("Geometry.PointNumbers", 1)
gmsh.option.setColor("Geometry.Color.Points", 255, 165, 0)
gmsh.option.setColor("General.Color.Text", 255, 255, 255)
gmsh.option.setColor("Mesh.Color.Points", 255, 0, 0)
# Note that for conciseness "Color." can be ommitted in color options:
r, g, b, a = gmsh.option.getColor("Geometry.Points")
gmsh.option.setColor("Geometry.Surfaces", r, g, b, a)
# We create a ONELAB parameter to define the angle of the twist. ONELAB
# parameters can be modified interactively in the GUI, and can be exchanged with
# other codes connected to the same ONELAB database. The database can be
# accessed through the Gmsh Julia API using JSON-formatted strings (see
# https://gitlab.onelab.info/doc/tutorials/-/wikis/ONELAB-JSON-interface for
# more information):
gmsh.onelab.set("""[
{
"type":"number",
"name":"Parameters/Twisting angle",
"values":[90],
"min":0,
"max":120,
"step":1
}
]""")
# Create the geometry and mesh it:
createGeometryAndMesh()
# Launch the GUI and handle the "check" event (recorded in the "ONELAB/Action"
# parameter) to recreate the geometry with a new twisting angle if necessary:
function checkForEvent()
action = gmsh.onelab.getString("ONELAB/Action")
if length(action) > 0 && action[1] == "check"
gmsh.onelab.setString("ONELAB/Action", [""])
createGeometryAndMesh()
gmsh.graphics.draw()
end
return true
end
if !("-nopopup" in ARGS)
gmsh.fltk.initialize()
while gmsh.fltk.isAvailable() == 1 && checkForEvent()
gmsh.fltk.wait()
end
end
# When the GUI is launched, you can use the `Help->Current Options and
# Workspace' menu to see the current values of all options. To save the options
# in a file, use `File->Export->Gmsh Options', or through the api:
# gmsh.write("t3.opt");
gmsh.finalize()
|