Removed dependency on NBodySimulator

Author: axla-io

Project.toml

@@ -4,27 +4,13 @@ author = ["Axel Larsson"]
 version = "0.2.0"
 
 [deps]
+StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
+LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
+Graphs = "86223c79-3864-5bf0-83f7-82e725a168b6"
+StaticGraphs = "4c8beaf5-199b-59a0-a7f2-21d17de635b6"
 DiffEqBase = "2b5f629d-d688-5b77-993f-72d75c75574e"
 DiffEqCallbacks = "459566f4-90b8-5000-8ac3-15dfb0a30def"
-FiniteDiff = "6a86dc24-6348-571c-b903-95158fe2bd41"
-Flux = "587475ba-b771-5e3f-ad9e-33799f191a9c"
 ForwardDiff = "f6369f11-7733-5829-9624-2563aa707210"
-GraphRecipes = "bd48cda9-67a9-57be-86fa-5b3c104eda73"
-Graphs = "86223c79-3864-5bf0-83f7-82e725a168b6"
-LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
-NBodySimulator = "0e6f8da7-a7fc-5c8b-a220-74e902c310f9"
-Optimization = "7f7a1694-90dd-40f0-9382-eb1efda571ba"
-OptimizationNLopt = "4e6fcdb7-1186-4e1f-a706-475e75c168bb"
-OptimizationOptimJL = "36348300-93cb-4f02-beb5-3c3902f8871e"
-OptimizationOptimisers = "42dfb2eb-d2b4-4451-abcd-913932933ac1"
-OptimizationPolyalgorithms = "500b13db-7e66-49ce-bda4-eed966be6282"
-OrdinaryDiffEq = "1dea7af3-3e70-54e6-95c3-0bf5283fa5ed"
-Plots = "91a5bcdd-55d7-5caf-9e0b-520d859cae80"
-SciMLSensitivity = "1ed8b502-d754-442c-8d5d-10ac956f44a1"
-StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
-StaticGraphs = "4c8beaf5-199b-59a0-a7f2-21d17de635b6"
-SteadyStateDiffEq = "9672c7b4-1e72-59bd-8a11-6ac3964bc41f"
-Zygote = "e88e6eb3-aa80-5325-afca-941959d7151f"
 
 [extras]
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"

src/DynamicRelaxation.jl

@@ -11,8 +11,6 @@ using StaticGraphs
 # Differential equation solving
 using DiffEqBase
 using DiffEqCallbacks
-using NBodySimulator
-using SteadyStateDiffEq
 using ForwardDiff
 
 include("elem.jl")
@@ -30,18 +28,18 @@ include("optimization/load_finding.jl")
 
 # Elements and nodes
 export ElementProperties, CoordinateSystem, Node3DOF, Node6DOF
-
+export Bar, Beam, AbstractElement
 # Callbacks
 export velocitydecay!, velocityreset!, ke_condition, ke_termination_cond
 
 # Constraints
-export clamped, free, pinned, roller
+export BigonTorqueCondition, clamped, free, pinned, roller
 
 # Loads
 export Px, Py, Pz, Mx, My, Mz, uniform_load, point_loads
 
 # System 
-export StructuralGraphSystem, default_system
+export StructuralGraphSystem, default_system, get_cs
 
 # Simulation
 export LoadScaleRodSimulation, RodSimulation, BigonRodSimulation, get_u0, get_vel_ids, get_state

src/node.jl

@@ -1,5 +1,7 @@
 # Struct for storing nodal properties
-struct Node3DOF <: NBodySimulator.Body
+
+abstract type AbstractNode end
+struct Node3DOF <: AbstractNode
     # Initial conditions
     r::SVector{3,Float64}
     v::SVector{3,Float64}
@@ -13,7 +15,7 @@ function Node3DOF(pos, constrained, constraints)
     Node3DOF(pos, @SVector(zeros(Float64, 3)), constrained, constraints)
 end
 
-struct Node6DOF <: NBodySimulator.Body
+struct Node6DOF <: AbstractNode
     # Initial conditions
     r::SVector{3,Float64}
     q::SVector{4,Float64}