extrude_mesh

| main | Tutorials | Functions | website |

Extrudes a 2D mesh to 3D. If the 2D mesh has triangular elements the 3D elements will be prisms. If the 2D mesh has quadrilaterals then the 3D elements will be hexahedral.

Version : 1.0

Author : George Kourakos

email: giorgk@gmail.com

web : http://groundwater.ucdavis.edu/msim

Date 18-Mar-2014

Department of Land Air and Water

University of California Davis

Contents

Usage

[p MSH]=extrude_mesh(p_2d, MSH_2d, top, bot, t, el_order)

Input

p_2d: Coordinates of 2D mesh nodes

MSH_2d: Mesh structure generated by Gmsh.

top: top elevation of the 2D mesh nodes

bot: bottom elevation of the 2D mesh nodes

t: parametric variable that defines how many layers will be generated and their vertical distribution between the top and bottom elevation. The t values must be t(1) = 0 (corresponds to bottom elevation) and t(end) = 1 (coresponds to top elevation) A quick way to plot the layer distribution is plot(ones(length(t),1),t,'.-')

el_order:element order 'linear' or 'quadratic'. At the moment incomplete elements are not supported

Output

p: coordinates of 3D mesh nodes

MSH 3D mesh structure.

see also: Centerfor2points

| main | Tutorials | Functions | website |