Glassvision/JuliaFEM.jl
1
0
Fork 0
mirror of https://github.com/JuliaFEM/JuliaFEM.jl.git synced 2026-08-06 04:21:33 +00:00
Code Issues Packages Projects Releases Wiki Activity

test(solvers): add modal elasticity test 2

Browse Source 
New 59-line test file for modal elasticity analysis:
- Tests modal analysis of two cylinder meshes (CYLINDER1 and CYLINDER20)
- Compares eigenfrequencies with Code Aster results (Tet4 elements)
- Tests 3 lowest eigenfrequencies for each cylinder
- Uses legacy Element/Problem API

Validates modal analysis accuracy for different cylinder geometries
by comparing with reference Code Aster results.
This commit is contained in:
Jukka Aho
2025-12-15 13:45:46 +02:00
parent f3f42fd115
commit 782c252719
1 changed files with 59 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files
test/solvers/test_modal_elasticity_2.jl
+59
View File
@@ -0,0 +1,59 @@
# This file is a part of JuliaFEM.
# License is MIT: see https://github.com/JuliaFEM/JuliaFEM.jl/blob/master/LICENSE.md
using JuliaFEM, Test
# eigenvalues of CYLINDER1
meshfile = @__DIR__() * "/testdata/primitives.med"
mesh = aster_read_mesh(meshfile, "CYLINDER_1_TET4")
cylinder = Problem(Elasticity, "CYLINDER", 3)
cylinder_elements = create_elements(mesh, "CYLINDER")
update!(cylinder_elements, "youngs modulus", 10000.0)
update!(cylinder_elements, "poissons ratio", 0.3)
update!(cylinder_elements, "density", 10.0)
add_elements!(cylinder, cylinder_elements)
bc = Problem(Dirichlet, "bc", 3, "displacement")
bc_elements = create_elements(mesh, "FACE_YZ1")
update!(bc_elements, "displacement 1", 0.0)
update!(bc_elements, "displacement 2", 0.0)
update!(bc_elements, "displacement 3", 0.0)
add_elements!(bc, bc_elements)
analysis = Analysis(Modal)
analysis.properties.nev = 3
add_problems!(analysis, cylinder, bc)
run!(analysis)
freqs_jf = sqrt.(analysis.properties.eigvals)/(2.0*pi)
freqs_ca = [4.84532E+00, 4.90698E+00, 8.33813E+00]
@test isapprox(freqs_jf, freqs_ca; rtol=1.0e-5)
# eigenvalues of CYLINDER20
meshfile = @__DIR__() * "/testdata/primitives.med"
mesh = aster_read_mesh(meshfile, "CYLINDER_20_TET4")
cylinder = Problem(Elasticity, "CYLINDER", 3)
cylinder_elements = create_elements(mesh, "CYLINDER")
#update!(cylinder_elements, "youngs modulus", 10.0e6)
update!(cylinder_elements, "youngs modulus", 50475.5)
update!(cylinder_elements, "poissons ratio", 0.3)
#update!(cylinder_elements, "density", 10.0)
update!(cylinder_elements, "density", 1.0)
add_elements!(cylinder, cylinder_elements)
bc = Problem(Dirichlet, "bc", 3, "displacement")
bc_elements = create_elements(mesh, "FACE1", "FACE2")
update!(bc_elements, "displacement 1", 0.0)
update!(bc_elements, "displacement 2", 0.0)
update!(bc_elements, "displacement 3", 0.0)
add_elements!(bc, bc_elements)
analysis = Analysis(Modal)
analysis.properties.nev = 3
add_problems!(analysis, cylinder, bc)
run!(analysis)
freqs_jf = sqrt.(analysis.properties.eigvals)/(2.0*pi)
freqs_ca = [1.19789E+00, 1.20179E+00, 3.07391E+00]
@test isapprox(freqs_jf, freqs_ca; rtol=1.0e-5)