chart corrected

Author: mkraska

NonLinear/Sections/Plastic/README.md

@@ -62,4 +62,4 @@ To create a force-displacement-plot, type
 ```
 Note that this script does not yet account for changes in the input file (maximum displacement is hard-wired).
 
-<img src="df.png" title="Force-displacement-plot">
+<img src="df.png" title="Force-displacement-plot (full model)">

NonLinear/Sections/Plastic/df.gpl

@@ -6,8 +6,10 @@ set ylabel "Force in kN"
 #set nokey
 set key bottom right
 set out "df.png"
-plot "total force fx,fy,fz_NXLSQUARE.txt" using ($1*20):(-1*$4/1000) title "square" w l,\
-	"total force fx,fy,fz_NXLCIRCLE.txt" using ($1*20):(-1*$4/1000) title "circle" w l, \
-	"total force fx,fy,fz_NXLHCIRCLE.txt" using ($1*20):(-1*$4/1000) title "hcircle" w l, \
-	"total force fx,fy,fz_NXLIPROFILE.txt" using ($1*20):(-1*$4/1000) title "Iprofile" w l
+# force factor: half model, remove negative sign
+ff=-2
+plot "total force fx,fy,fz_NXLSQUARE.txt" using ($1*20):(ff*$4/1000) title "square" w l,\
+	"total force fx,fy,fz_NXLCIRCLE.txt" using ($1*20):(ff*$4/1000) title "circle" w l, \
+	"total force fx,fy,fz_NXLHCIRCLE.txt" using ($1*20):(ff*$4/1000) title "hcircle" w l, \
+	"total force fx,fy,fz_NXLIPROFILE.txt" using ($1*20):(ff*$4/1000) title "Iprofile" w l
 set out #