diff --git a/test/runtests.jl.old b/test/runtests.jl.old
new file mode 100644
index 0000000..eb7a173
--- /dev/null
+++ b/test/runtests.jl.old
@@ -0,0 +1,163 @@
+# This file is a part of JuliaFEM.
+# License is MIT: see https://github.com/JuliaFEM/JuliaFEM.jl/blob/master/LICENSE.md
+
+using JuliaFEM, Test
+
+@testset "JuliaFEM.jl" begin
+    @testset "test_dirichlet.jl" begin
+        include("test_dirichlet.jl")
+    end
+    @testset "test_elasticity_1d.jl" begin
+        include("test_elasticity_1d.jl")
+    end
+    @testset "test_elasticity_2d_linear_with_surface_load.jl" begin
+        include("test_elasticity_2d_linear_with_surface_load.jl")
+    end
+    @testset "test_elasticity_2d_nonhomogeneous_boundary_conditions.jl" begin
+        include("test_elasticity_2d_nonhomogeneous_boundary_conditions.jl")
+    end
+    @testset "test_elasticity_2d_nonlinear_with_surface_load.jl" begin
+        include("test_elasticity_2d_nonlinear_with_surface_load.jl")
+    end
+    @testset "test_elasticity_2d_plane_stress_stiffness_matrix.jl" begin
+        include("test_elasticity_2d_plane_stress_stiffness_matrix.jl")
+    end
+    @testset "test_elasticity_2d_residual.jl" begin
+        include("test_elasticity_2d_residual.jl")
+    end
+    @testset "test_elasticity_3d_linear_with_surface_load.jl" begin
+        include("test_elasticity_3d_linear_with_surface_load.jl")
+    end
+    @testset "test_elasticity_3d_nonlinear_with_surface_load.jl" begin
+        include("test_elasticity_3d_nonlinear_with_surface_load.jl")
+    end
+    @testset "test_elasticity_3d_unit_block.jl" begin
+        include("test_elasticity_3d_unit_block.jl")
+    end
+    @testset "test_elasticity_forwarddiff.jl" begin
+        include("test_elasticity_forwarddiff.jl")
+    end
+    @testset "test_elasticity_hollow_sphere_with_surface_pressure.jl" begin
+        include("test_elasticity_hollow_sphere_with_surface_pressure.jl")
+    end
+    @testset "test_elasticity_med_pyr5_point_load.jl" begin
+        include("test_elasticity_med_pyr5_point_load.jl")
+    end
+    @testset "test_elasticity_plane_strain.jl" begin
+        include("test_elasticity_plane_strain.jl")
+    end
+    @testset "test_elasticity_pyr5_point_load.jl" begin
+        include("test_elasticity_pyr5_point_load.jl")
+    end
+    @testset "test_elasticity_tet4_stiffness_matrix.jl" begin
+        include("test_elasticity_tet4_stiffness_matrix.jl")
+    end
+    @testset "test_elasticity_tet10_mass_matrix.jl" begin
+        include("test_elasticity_tet10_mass_matrix.jl")
+    end
+    @testset "test_elasticity_tet10_stiffness_matrix.jl" begin
+        include("test_elasticity_tet10_stiffness_matrix.jl")
+    end
+    @testset "test_elasticity_tetra.jl" begin
+        include("test_elasticity_tetra.jl")
+    end
+    @testset "test_elasticplastic_2d_nonhomogenious_boundary_conditions.jl" begin
+        include("test_elasticplastic_2d_nonhomogenious_boundary_conditions.jl")
+    end
+    @testset "test_elasticplastic_3d_linear_with_surface_load.jl" begin
+        include("test_elasticplastic_3d_linear_with_surface_load.jl")
+    end
+    @testset "test_heat_2d_one_element.jl" begin
+        include("test_heat_2d_one_element.jl")
+    end
+    @testset "test_heat_3d.jl" begin
+        include("test_heat_3d.jl")
+    end
+    @testset "test_heat_tet10_convection.jl" begin
+        include("test_heat_tet10_convection.jl")
+    end
+    @testset "test_heat_3d_2.jl" begin
+        include("test_heat_3d_2.jl")
+    end
+    @testset "test_heat.jl" begin
+        include("test_heat.jl")
+    end
+    @testset "test_heat_2.jl" begin
+        include("test_heat_2.jl")
+    end
+    @testset "test_heat_3.jl" begin
+        include("test_heat_3.jl")
+    end
+    @testset "test_heat_3d_two_rings.jl" begin
+        include("test_heat_3d_two_rings.jl")
+    end
+    @testset "test_heat_4.jl" begin
+        include("test_heat_4.jl")
+    end
+    @testset "test_modal_analysis.jl" begin
+        include("test_modal_analysis.jl")
+    end
+    @testset "test_modal_analysis_elasticity.jl" begin
+        include("test_modal_analysis_elasticity.jl")
+    end
+    @testset "test_modal_analysis_elasticity_2.jl" begin
+        include("test_modal_analysis_elasticity_2.jl")
+    end
+    @testset "test_modal_analysis_zero_eigenmodes.jl" begin
+        include("test_modal_analysis_zero_eigenmodes.jl")
+    end
+    @testset "test_mortar.jl" begin
+        include("test_mortar.jl")
+    end
+    @testset "test_mortar_2d.jl" begin
+        include("test_mortar_2d.jl")
+    end
+    @testset "test_mortar_2d_assembly.jl" begin
+        include("test_mortar_2d.jl")
+    end
+    @testset "test_mortar_2d_contact.jl" begin
+        include("test_mortar_2d_contact.jl")
+    end
+    @testset "test_mortar_2d_mesh_tie.jl" begin
+        include("test_mortar_2d_mesh_tie.jl")
+    end
+    @testset "test_mortar_2d_weighted_gap.jl" begin
+        include("test_mortar_2d_weighted_gap.jl")
+    end
+    @testset "test_mortar_3d_mesh_tie_modal.jl" begin
+        include("test_mortar_3d_mesh_tie_modal.jl")
+    end
+    @testset "test_mortar_3d_mesh_tie_two_rings.jl" begin
+        include("test_mortar_3d_mesh_tie_two_rings.jl")
+    end
+    @testset "test_mortar_3d_polygon_clip.jl" begin
+        include("test_mortar_3d_polygon_clip.jl")
+    end
+    @testset "test_postprocess.jl" begin
+        include("test_postprocess.jl")
+    end
+    @testset "test_potential_energy.jl" begin
+        include("test_potential_energy.jl")
+    end
+    @testset "test_problems_contact_3d.jl" begin
+        include("test_problems_contact_3d.jl")
+    end
+    @testset "test_problems_elasticity.jl" begin
+        include("test_problems_elasticity.jl")
+    end
+    @testset "test_problems_mortar_3d.jl" begin
+        include("test_problems_mortar_3d.jl")
+    end
+    @testset "test_problems_mortar_3d_lowlevel.jl" begin
+        include("test_problems_mortar_3d_lowlevel.jl")
+    end
+    @testset "test_solvers_postprocess.jl" begin
+        include("test_solvers_postprocess.jl")
+    end
+    @testset "test_virtual_work.jl" begin
+        include("test_virtual_work.jl")
+    end
+    @testset "test_von_mises_material.jl" begin
+        include("test_von_mises_material.jl")
+    end
+end
diff --git a/test/runtests_new.jl b/test/runtests_new.jl
new file mode 100644
index 0000000..aea77c6
--- /dev/null
+++ b/test/runtests_new.jl
@@ -0,0 +1,44 @@
+# JuliaFEM Test Suite - New Structure
+# Educational testing with Literate.jl
+
+using Test, JuliaFEM
+
+# Configuration
+const RUN_TUTORIALS = get(ENV, "JULIAFEM_TEST_TUTORIALS", "true") == "true"
+const RUN_UNIT = get(ENV, "JULIAFEM_TEST_UNIT", "false") == "true"
+const RUN_OLD = get(ENV, "JULIAFEM_TEST_OLD", "false") == "true"
+
+println("="^70)
+println("JuliaFEM Test Suite (New Structure)")
+println("="^70)
+println("Tutorials: ", RUN_TUTORIALS ? "✓" : "✗")
+println("Unit tests: ", RUN_UNIT ? "✓" : "✗")
+println("Old tests: ", RUN_OLD ? "✓" : "✗")
+println("="^70)
+
+# Tutorial Tests
+if RUN_TUTORIALS
+    @testset "Tutorials" begin
+        @testset "01_Fundamentals" begin
+            include("tutorials/01_fundamentals/creating_elements.jl")
+        end
+    end
+end
+
+# Unit Tests  
+if RUN_UNIT
+    @testset "Unit Tests" begin
+        @info "No unit tests yet"
+    end
+end
+
+# Old Tests
+if RUN_OLD
+    @warn "Old tests have 49+ failures - see docs/TEST_FIXES_NEEDED.md"
+    include("runtests.jl")  # Original test suite
+end
+
+println()
+println("="^70)
+println("Complete - See docs/TESTING_PHILOSOPHY.md")
+println("="^70)
diff --git a/test/tutorials/01_fundamentals/creating_elements.jl b/test/tutorials/01_fundamentals/creating_elements.jl
new file mode 100644
index 0000000..2954e6c
--- /dev/null
+++ b/test/tutorials/01_fundamentals/creating_elements.jl
@@ -0,0 +1,159 @@
+# # Tutorial 1: Creating Elements and Fields
+#
+# This tutorial introduces the fundamental concepts of JuliaFEM:
+# - Creating nodes and elements
+# - Updating element fields
+# - Understanding the field concept
+#
+# ## Prerequisites
+# 
+# Basic Julia knowledge. This is the first tutorial - no prior JuliaFEM experience needed!
+
+using JuliaFEM
+using Test
+
+# ## Nodes and Coordinates
+#
+# In finite element analysis, everything starts with **nodes** - points in space
+# where we'll compute displacements, temperatures, or other field values.
+#
+# In JuliaFEM, nodes are simply vectors of coordinates. Let's create a few nodes
+# for a 2D square:
+
+node1 = [0.0, 0.0]  # Bottom-left corner
+node2 = [1.0, 0.0]  # Bottom-right corner  
+node3 = [1.0, 1.0]  # Top-right corner
+node4 = [0.0, 1.0]  # Top-left corner
+
+# We typically store nodes in a dictionary for easy access by node ID:
+
+X = Dict(
+    1 => node1,
+    2 => node2,
+    3 => node3,
+    4 => node4
+)
+
+# ## Elements
+#
+# An **element** connects nodes and defines:
+# - The shape (linear, quadratic, etc.)
+# - The interpolation (how fields vary within the element)
+# - The physics (what equations to solve)
+#
+# Let's create a 4-node quadrilateral element (`Quad4`):
+
+element = Element(Quad4, (1, 2, 3, 4))
+
+# The tuple `(1, 2, 3, 4)` specifies which nodes this element connects.
+# Node numbering follows a counter-clockwise convention:
+#
+# ```
+# 4 ---- 3
+# |      |
+# |      |  
+# 1 ---- 2
+# ```
+
+# ## The Field Concept
+#
+# In JuliaFEM, everything is a **field**. A field is data associated with an element
+# that can vary in space and time:
+#
+# - **Geometry field:** Node coordinates
+# - **Material fields:** Young's modulus, Poisson's ratio
+# - **Load fields:** Body forces, surface tractions
+# - **Solution fields:** Displacements, temperatures
+#
+# Fields are updated using the `update!` function.
+
+# ### Geometry Field
+#
+# First, we tell the element where its nodes are located:
+
+update!(element, "geometry", X)
+
+# Now the element knows its shape and can compute things like its area,
+# jacobian, etc.
+
+# ### Material Fields
+#
+# Let's add material properties (for an elasticity problem):
+
+update!(element, "youngs modulus", 200.0e3)  # 200 GPa steel
+update!(element, "poissons ratio", 0.3)
+
+# ### Load Fields
+#
+# We can apply loads as fields too. For example, a body force in the y-direction:
+
+update!(element, "displacement load 2", 10.0)  # 10 N/m³ in y-direction
+
+# The naming convention:
+# - `"displacement load 1"` = body force in x-direction
+# - `"displacement load 2"` = body force in y-direction  
+# - `"displacement load 3"` = body force in z-direction (3D)
+
+# ## Field Access
+#
+# We can retrieve fields using function call syntax:
+
+geom = element("geometry", 0.0)  # Get geometry at time t=0.0
+E = element("youngs modulus", 0.0)
+ν = element("poissons ratio", 0.0)
+
+# Fields can be time-dependent! The second argument is the time value.
+# For static fields (like geometry), the time doesn't matter.
+
+# ## Testing Our Understanding
+#
+# Let's verify everything works as expected:
+
+@testset "Element Creation" begin
+    # Element should have 4 nodes
+    @test length(element.connectivity) == 4
+    
+    # Geometry field returns a tuple of node coordinates
+    @test length(geom) == 4  # 4 nodes
+    @test geom[1] == [0.0, 0.0]  # First node
+    
+    # Material properties should be retrievable
+    @test E == 200.0e3
+    @test ν == 0.3
+end
+
+# ## Multiple Elements
+#
+# Real problems have many elements. Let's create a vector of elements:
+
+elements = [
+    Element(Quad4, (1, 2, 3, 4)),
+    Element(Quad4, (2, 5, 6, 3))  # Adjacent element (assuming nodes 5, 6 exist)
+]
+
+# We can update fields for all elements at once:
+
+update!(elements, "youngs modulus", 200.0e3)
+update!(elements, "poissons ratio", 0.3)
+
+# This is more efficient than updating each element individually.
+
+# ## What We Learned
+#
+# ✅ Nodes are just coordinate vectors  
+# ✅ Elements connect nodes using connectivity tuples  
+# ✅ Everything is a field (geometry, materials, loads, solutions)  
+# ✅ Fields are updated with `update!(element, "field_name", value)`  
+# ✅ Fields are accessed with `element("field_name", time)`  
+# ✅ Multiple elements can be updated together
+#
+# ## Next Steps
+#
+# - **Tutorial 2:** Reading meshes from files (ABAQUS .inp format)
+# - **Tutorial 3:** Basis functions and shape function evaluation
+# - **Tutorial 4:** Solving your first problem (1D elasticity)
+#
+# ## Further Reading
+#
+# - JuliaFEM paper: [Frondelius & Aho (2017)](https://doi.org/10.23998/rm.64224)
+# - Field concept: Based on Abaqus field definitions