From a22963d2e579438a7b8a02bdc8af7716710ba734 Mon Sep 17 00:00:00 2001 From: Olli Vainola Date: Sat, 16 May 2015 09:19:22 +0300 Subject: [PATCH 1/2] shape functions added --- src/shape_functions.jl | 1 + 1 file changed, 1 insertion(+) create mode 100644 src/shape_functions.jl diff --git a/src/shape_functions.jl b/src/shape_functions.jl new file mode 100644 index 0000000..c40ebaf --- /dev/null +++ b/src/shape_functions.jl @@ -0,0 +1 @@ +function rectangle_linear(\eta, \zeta) From 25f5fe541e5264133c1d6a1a542a03223eef49bb Mon Sep 17 00:00:00 2001 From: Olli Vainola Date: Sat, 16 May 2015 10:04:27 +0300 Subject: [PATCH 2/2] hexa shape functions added --- src/shape_functions.jl | 20 +++++++++++++++++++- 1 file changed, 19 insertions(+), 1 deletion(-) diff --git a/src/shape_functions.jl b/src/shape_functions.jl index c40ebaf..f035904 100644 --- a/src/shape_functions.jl +++ b/src/shape_functions.jl @@ -1 +1,19 @@ -function rectangle_linear(\eta, \zeta) +function C2D4(ξ::Float64, η::Float64) + values::Array{Float64, 1} = [(1 - ξ) * (1 - η), + ξ * (1 - η), + ξ * η, + (1 - ξ) * η] + return values +end + +function C3D8(ξ::Float64, η::Float64, μ::Float64) + values::Array{Float64, 1} = 1 / 8 * [(1 - ξ) * (1 - η) * (1 - μ), + (1 + ξ) * (1 - η) * (1 - μ), + (1 + ξ) * (1 + η) * (1 - μ), + (1 - ξ) * (1 + η) * (1 - μ), + (1 - ξ) * (1 - η) * (1 + μ), + (1 + ξ) * (1 - η) * (1 + μ), + (1 + ξ) * (1 + η) * (1 + μ), + (1 - ξ) * (1 + η) * (1 + μ)] + return values +end