Create a source term from a given "Strain/Stress Load"

Now the use of these keywords is independent of the keywords Pre Strain/Stress
ElmerCSC · Apr 15, 2024 · 48ade7a · 48ade7a
1 parent 0c46e1d
commit 48ade7a
Show file tree

Hide file tree

Showing 2 changed files with 7 additions and 4 deletions.
diff --git a/fem/src/SOLVER.KEYWORDS b/fem/src/SOLVER.KEYWORDS
@@ -469,11 +469,13 @@ BodyForce:Real:     'Sound Source'
 BodyForce:Real:     'Source'
 BodyForce:Real:     'Source Field'
 BodyForce:Real:     'Strain'
+BodyForce:Real:     'Strain Load'
+BodyForce:Real:     'Stress'
 BodyForce:Real:     'Stress Bodyforce 1'
 BodyForce:Real:     'Stress Bodyforce 2'
 BodyForce:Real:     'Stress Bodyforce 3'
+BodyForce:Real:     'Stress Load'
 BodyForce:Real:     'Stress Pressure'
-BodyForce:Real:     'Stress'
 bodyforce:real:     'Temperature'
 BodyForce:Real:     'Temperature Condition'
 BodyForce:Real:     'Temperature Load'

diff --git a/fem/src/Stress.F90 b/fem/src/Stress.F90
@@ -175,9 +175,10 @@ SUBROUTINE StressCompose( MASS, DAMP, STIFF, FORCE, FORCE_im, LOAD, LOAD_im, Ela
 
      NeedHeat = ANY( NodalTemperature(1:ntot) /= 0.0d0 )
      IF (.NOT.EvaluateLoadAtIP) NeedHarmonic = ANY( LOAD_im(:,1:n) /= 0.0d0 ) 
-     NeedPreStress = ANY( NodalPreStrain(1:6,1:n) /= 0.0d0 ) 
-     NeedPreStress = NeedPreStress .OR. ANY( NodalPreStress(1:6,1:n) /= 0.0d0 ) 
-
+     NeedPreStress = ANY( NodalPreStrain(1:6,1:n) /= 0.0d0 ) .OR. &
+         ANY( NodalStrainLoad(1:6,1:n) /= 0.0d0 )
+     NeedPreStress = NeedPreStress .OR. ANY( NodalPreStress(1:6,1:n) /= 0.0d0 ) .OR. &
+         ANY( NodalStressLoad(1:6,1:n) /= 0.0d0 )
 
      BF => GetBodyForce()
      GPA = .FALSE.