force-displacement chart added

Author: mkraska

NonLinear/Sections/Plastic/README.md

@@ -7,14 +7,15 @@ Tested with CGX/CCX 2.10
 + Elasto-plastic material without hardening (ideal plastic)
 + Objective:
   + Evolution of the plastic zone
-  + Force-displacement curves (TBD)
+  + Force-displacement curves
   + Ultimate bending load - plastic section modulus (TBD)
 
 
 File                           | Contents    
  :-------------                | :-------------
  [par.pre.fbl](par.pre.fbl)    | CGX script, pre-processing, parametrized with param.py
  [movie.fbl](movie.fbl)        | CGX script, movie generation
+ [chart.fbl](chart.fbl)        | CGX script, chart generation
  [Biegung.inp](Biegung.inp)    | CCX input
 
 

NonLinear/Sections/Plastic/chart.fbl

@@ -0,0 +1,2 @@
+sys dat2txt.py
+sys gnuplot df.gpl

NonLinear/Sections/Plastic/df.gpl

@@ -0,0 +1,13 @@
+set term png
+set grid
+set xlabel "Displacement in mm"
+set ylabel "Force in kN"
+#set yrange [-100:100]
+#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
+set out #