added tests and improved latticeslip. Small changes

Author: githubgrasp

src/sm/Materials/LatticeMaterials/latticeframeconcreteplastic.C

@@ -891,6 +891,7 @@ namespace oofem {
         //fprintf(file, "kappad %.8e ", this->kappaD);
         // fprintf(file, "damage %.8e ", this->damage);
         fprintf(file, ", kappaD %.8e, damage %.8e \n", this->kappaD, this->damage);
+	fprintf(file, "\n");
     }
 
     void

src/sm/Materials/LatticeMaterials/latticeframesteelplastic.C

@@ -669,6 +669,7 @@ namespace oofem {
             fprintf(file, "% .8e ", s);
         }
         fprintf(file, "kappa  % .8e ", this->kappa);
+	fprintf(file, "\n");
     }
 
 

src/sm/Materials/LatticeMaterials/latticeslip.C

@@ -67,13 +67,12 @@ namespace oofem {
     {
         LatticeLinearElastic::initializeFrom(ir);
 
-        //Parameter which relates the shear stiffness to the normal stiffness. Default is 1000. Reread this again, because the
-        alphaOne = 1000.;
-        IR_GIVE_OPTIONAL_FIELD(ir, alphaOne, _IFT_LatticeSlip_a1); // Macro
 
-        //Parameter which is used for the definition of bending stiffness. Default is 1000.
-        alphaTwo = 1000.;
-        IR_GIVE_OPTIONAL_FIELD(ir, alphaTwo, _IFT_LatticeSlip_a2); // Macro
+        //Parameter which relates the shear stiffness to the normal stiffness. Reread this again because we use the stiffness matrix from latticelinearelastic material where the default is 1.
+        IR_GIVE_FIELD(ir, this->alphaOne, _IFT_LatticeSlip_a1); // Macro
+
+        //Parameter which is used for the definition of bending stiffness. Reread this again because we use the stiffness matrix from latticelinearelastic material where the default is 0. It cannot be zero for latticeslip.
+        IR_GIVE_FIELD(ir, this->alphaTwo, _IFT_LatticeSlip_a2); // Macro
 
         //Parameter which limits the stress in slip direction.
         IR_GIVE_FIELD(ir, tauMax, _IFT_LatticeSlip_t0); // Macro

tests/sm/latticelink1.in

@@ -9,23 +9,23 @@ node 1 coords 3 0. 0. 0.
 node 2 coords 3 0.1 0. 0.
 node 3 coords 3 0.2 0. 0. 
 #Lattice nodes
-node 4 coords 3 0.005 0.01 0.01
-node 5 coords 3 -0.001 -0.008 0.012
-node 6 coords 3 0.01 0.02 -0.008
+node 4 coords 3 0. 0. 0.
+node 5 coords 3 0.1 0. 0.
+node 6 coords 3 0.2 0.0 0.
 #BeamElements
 latticeframe3d 1 nodes 2 1 2 crossSect 1 mat 1
 latticeframe3d 2 nodes 2 2 3 crossSect 1 mat 1
 #links for first node
 latticelink3d 3 nodes 2 1 4 crossSect 2 mat 2 length 0.05 diameter 1.e-2 dirvector 3 1.000000e+00 0.000000e+00 0.000000e+00
 latticelink3d 4 nodes 2 2 5 crossSect 2 mat 2 length 0.1 diameter 1.e-02 dirvector 3 1.000000e+00 0.000000e+00 0.000000e+00
 latticelink3d 5 nodes 2 3 6 crossSect 2 mat 2 length 0.05 diameter 1.e-02 dirvector 3 1.000000e+00 0.000000e+00 0.000000e+00
-latticecs 1 material 1 iy 7.857e-9 iz 7.857e-9 ik 1.6e-8 shearCoeff 0.83 area 3.1416e-4
+latticecs 1 material 1 iy 7.857e-9 iz 7.857e-9 ik 1.6e-8 shearCoeff 0.83 area 3.1416e-2
 latticecs 2 material 2
-latticeframeelastic 1 d 0 e 2.1e21 n 0.3 talpha 0
-latticeslip 2 talpha 0. d 0. e 30.e9 a1 1000 type 2 alpha 0.4 s1 1.e-3 s2 2e-3 s3 6.5e-3 t0 1.541e7 tf 6.164e6
+latticeframeelastic 1 d 0 e 2.1e11 n 0.3 talpha 0
+latticeslip 2 talpha 0. d 0. e 30.e9 type 2 alpha 0.4 s1 1.e-3 s2 2e-3 s3 6.5e-3 t0 1.541e7 tf 6.164e6 a1 100. a2 100.
 BoundaryCondition 1 loadTimeFunction 1 dofs 6 1 2 3 4 5 6 values 6 0. 0. 0. 0. 0. 0. set 4
 BoundaryCondition 2 loadTimeFunction 1 dofs 1 1 values 1 0.01 set 2
-PiecewiseLinFunction 1 nPoints 2 t 3 0. 7. 10. f(t) 3 0. 1. 0.9 
+PiecewiseLinFunction 1 nPoints 3 t 3 0. 9. 10. f(t) 3 0. 1. 0.93 
 set 1 elements 2 1 2
 set 2 nodes 1 3 
 set 3 elements 3 3 4 5

tests/sm/linearhardeningframesteel.in

@@ -0,0 +1,20 @@
+linearhardeningframesteel.out
+Cantilever test with latticeframesteelplastic material model
+NonLinearStatic nmsteps 1 nsteps 20  contextOutputStep 100
+nsteps 4 rtolv 1.e-0 maxiter 2000 stiffMode 1 controllmode 1 refloadmode 1 ddm 2 3 1 ddv 1 1 ddltf 1
+domain 3dLattice
+OutputManager tstep_all element_all dofman_all
+ndofman 3 nelem 2 ncrosssect 1 nmat 1 nbc 1 nic 0 nltf 1
+node 1 coords 3 0.000000e+00 0.000000e+00 0.000000e+00 bc 6 1 1 1 1 1 1 
+node 2 coords 3 0.500000e+00 0.000000e+00 0.000000e+00  
+node 3 coords 3 1.000000e+00 0.000000e+00 0.000000e+00 bc 6 0 0 0 0 0 0
+latticeframe3d 1 nodes 2 1 2 zaxis 3 0 0 1 crossSect 1 
+latticeframe3d 2 nodes 2 2 3 zaxis 3 0 0 1 crossSect 1
+latticecs 1 material 1 iy 0.083 iz 0.083 ik 0.167 shearCoeff 1. area 1.
+latticeframesteelplastic 1 d 1. E 1. n 0.15 talpha 0 nx0 1 mx0 1 my0 1 mz0 1 tol 1.e-6 iter 100 sub 10 h 0.3
+BoundaryCondition 1 loadTimeFunction 1 prescribedvalue 0.0
+ConstantFunction 1 f(t) 1.
+#%BEGIN_CHECK%
+#ELEMENT number 2 gp 1 keyword latticestrain component 1
+#ELEMENT number 2 gp 1 keyword latticestress component 1
+#%END_CHECK%

tests/sm/multilinearhardeningframesteel.in

@@ -0,0 +1,21 @@
+multilinearhardeningframesteel.out
+Cantilever test with latticeframesteelplastic material model
+NonLinearStatic nmsteps 1 nsteps 20  contextOutputStep 100 nmodules 0
+nsteps 4 rtolv 1.e-0 maxiter 2000 stiffMode 2 controllmode 1 refloadmode 1 ddm 2 3 1 ddv 1 2 ddltf 1 
+#vtkxml primvars 1 1 cellvars 6 138 139 140 156 157 158 tstep_step 1 domain_all
+domain 3dLattice
+OutputManager tstep_all element_all dofman_all
+ndofman 3 nelem 2 ncrosssect 1 nmat 1 nbc 1 nic 0 nltf 1
+node 1 coords 3 0.000000e+00 0.000000e+00 0.000000e+00 bc 6 1 1 1 1 1 1 
+node 2 coords 3 0.500000e+00 0.000000e+00 0.000000e+00  
+node 3 coords 3 1.000000e+00 0.000000e+00 0.000000e+00 bc 6 0 0 0 0 0 0
+latticeframe3d 1 nodes 2 1 2 zaxis 3 0 0 1 crossSect 1 
+latticeframe3d 2 nodes 2 2 3 zaxis 3 0 0 1 crossSect 1
+latticecs 1 material 1 iy 0.083 iz 0.083 ik 0.167 shearCoeff 1. area 1.
+latticeframesteelplastic 1 d 1. E 1. n 0.15 talpha 0 nx0 1 mx0 1 my0 1 mz0 1 tol 1.e-6 iter 100 sub 10 htype 1 h_eps 10 0. 0.3 0.6 1. 1.2 2.5 4. 4.5 10 100 h(eps) 10 1. 1. 1.143 1.214 1.281 1.319 1.333 0.666 0.1 0.1
+BoundaryCondition 1 loadTimeFunction 1 prescribedvalue 0.0
+ConstantFunction 1 f(t) 1.
+#%BEGIN_CHECK% 
+#ELEMENT number 2 gp 1 keyword latticestrain component 1
+#ELEMENT number 2 gp 1 keyword latticestress component 1
+#%END_CHECK%