Glassvision/JuliaFEM.jl
1
0
Fork 0
mirror of https://github.com/JuliaFEM/JuliaFEM.jl.git synced 2026-09-22 02:40:51 +00:00
Code Issues Packages Projects Releases Wiki Activity

more integration rules

Browse Source
This commit is contained in:
Jukka Aho
2016-02-05 22:22:40 +02:00
parent baf6074de9
commit 69ebba5e11
1 changed files with 8 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files
src/integrate.jl
+8 -1
View File
@@ -119,6 +119,13 @@ function get_integration_points(::Type{Quad4}, ::Type{Val{2}})
]
end
function get_integration_points(::Type{Quad4}, ::Type{Val{3}})
p = [-sqrt(3/5), 0.0, sqrt(3/5)]
w = [5/9, 8/9, 5/9]
pts = vec([IntegrationPoint([p[i], p[j]], w[i]*w[j]) for i=1:3, j=1:3])
return pts
end
function get_integration_points(::Type{Quad4})
return get_integration_points(Quad4, Val{2})
end
@@ -128,7 +135,7 @@ end
function get_integration_points(::Type{Hex8}, ::Type{Val{2}})
p = 1.0/sqrt(3.0)*[-1.0, 1.0]
w = 1.0
w = [1.0, 1.0]
return vec([IntegrationPoint([p[i], p[j], p[k]], w) for i=1:2, j=1:2, k=1:2])
end