Don't double count stiffnesses when using the two part nonlinear stra…

…in stiffening update option
sandialabs · Oct 11, 2024 · 7abd69e · 7abd69e
1 parent e1a73e2
commit 7abd69e
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 8 deletions.
diff --git a/src/modal.jl b/src/modal.jl
@@ -1,5 +1,5 @@
 """
-frequencies = autoCampbellDiagram(FEAinputs,mymesh,myel,system,assembly;
+frequencies = autoCampbellDiagram(FEAinputs,mymesh,myel,system,assembly,sections;
     minRPM = 0.0,
     maxRPM = 40.0,
     NRPM = 9, # int

diff --git a/src/timoshenko.jl b/src/timoshenko.jl
@@ -835,13 +835,13 @@ function calculateTimoshenkoElementNL(input,elStorage;predef=nothing)
     end
 
     if predef == "update"
-        elStorage.K21NLpredef = collect(K21)
-        elStorage.K12NLpredef = collect(K12)
-        elStorage.K31NLpredef = collect(K31)
-        elStorage.K13NLpredef = collect(K13)
-        elStorage.K22NLpredef = collect(K22)
-        elStorage.K23NLpredef = collect(K23)
-        elStorage.K33NLpredef = collect(K33)
+        elStorage.K21NLpredef = collect(2*K12NL')
+        elStorage.K12NLpredef = collect(K12NL)
+        elStorage.K31NLpredef = collect(2*K13NL')
+        elStorage.K13NLpredef = collect(K13NL)
+        elStorage.K22NLpredef = collect(K22NL)
+        elStorage.K23NLpredef = collect(K23NL)
+        elStorage.K33NLpredef = collect(K33NL)
     end
 
     # Only used if (useDisp)