STLCutters

STL to cell-wise triangulation to solve FE problems in Gridap.jl through GridapEmbedded.jl

Installation

# Type ] to enter package mode
pkg> add STLCutters

Examples

Sub-triangulation examples

Use a test geometry, e.g., 47076.stl (Chichen Itza)

julia> include("examples/SubTriangulation.jl")
julia> filename = "test/data/47076.stl"
julia> SubTriangulation.main(filename,n=50,output="example1")

Download a geometry directly from Thingi10k, e.g, 37384. Please check whether the geometry is solid and manifold in Thingi10k metadata.

julia> include("examples/SubTriangulation.jl")
julia> filename = SubTriangulation.download(37384)
julia> SubTriangulation.main(filename,n=50,output="example2")

Finite Elements examples

Solve a Poisson equation on a test geometry, e.g., 293137.stl (Low-Poly Bunny)

julia> include("examples/Poisson.jl")
julia> filename = "test/data/293137.stl"
julia> Poisson.main(filename,n=20,output="example3")

Solve a Linear Elasticity problem on a test geometry, e.g., 550964.stl (Eiffel Tower in a 5 degree slope)

julia> include("examples/LinearElasticity.jl")
julia> filename = "test/data/550964.stl"
julia> LinearElasticity.main(filename,n=50,force=(tand(5),0,-1),output="example4")

Solve an Incompressible Flow problem on a test geometry, e.g., 47076.stl (Chichen Itza)

julia> # ENV["ENABLE_MKL"] = "" ## Uncomment if GridapPardiso.jl requirements are fulfilled
julia> include("examples/Stokes.jl")
julia> filename = "test/data/47076.stl"
julia> Stokes.main(filename,n=10,output="example5")