﻿<liml version="24">
  <configuration name="ELU" />
  <analysis type="L30" solver="ccx" />
  <timesteps timeperiod="3 s" stepsize="0.001 s" automatictimestepping="true" maxstepsize="0.005 s" />
  <node nid="1" x="-0.5" y="-0.5" z="-0.5" />
  <node nid="2" x="0.5" y="-0.5" z="-0.5" />
  <node nid="3" x="0.5" y="0.5" z="-0.5" />
  <node nid="4" x="-0.5" y="0.5" z="-0.5" />
  <node nid="5" x="-0.5" y="-0.5" z="0.5" />
  <node nid="6" x="0.5" y="-0.5" z="0.5" />
  <node nid="7" x="0.5" y="0.5" z="0.5" />
  <node nid="8" x="-0.5" y="0.5" z="0.5" />
  <node nid="9" x="0" y="-0.5" z="-0.5" />
  <node nid="10" x="0.5" y="0" z="-0.5" />
  <node nid="11" x="0" y="0.5" z="-0.5" />
  <node nid="12" x="-0.5" y="0" z="-0.5" />
  <node nid="13" x="0" y="-0.5" z="0.5" />
  <node nid="14" x="0.5" y="0" z="0.5" />
  <node nid="15" x="0" y="0.5" z="0.5" />
  <node nid="16" x="-0.5" y="0" z="0.5" />
  <node nid="17" x="-0.5" y="-0.5" z="0" />
  <node nid="18" x="0.5" y="-0.5" z="0" />
  <node nid="19" x="0.5" y="0.5" z="0" />
  <node nid="20" x="-0.5" y="0.5" z="0" />
  <node nid="21" x="0" y="-0.5" z="0" />
  <node nid="22" x="0.5" y="0" z="0" />
  <node nid="23" x="0" y="0.5" z="0" />
  <node nid="24" x="-0.5" y="0" z="0" />
  <node nid="25" x="0" y="0" z="-0.5" />
  <node nid="26" x="0" y="0" z="0.5" />
  <node nid="27" x="0" y="0" z="0" />
  <elset name="Mazars" colorrgb="FF8000" expanded="true">
    <elem eid="1" shape="hex8" nodes="24 27 23 20 16 26 15 8" />
    <elem eid="2" shape="hex8" nodes="1 9 25 12 17 21 27 24" />
    <elem eid="3" shape="hex8" nodes="9 2 10 25 21 18 22 27" />
    <elem eid="4" shape="hex8" nodes="25 10 3 11 27 22 19 23" />
    <elem eid="5" shape="hex8" nodes="12 25 11 4 24 27 23 20" />
    <elem eid="6" shape="hex8" nodes="17 21 27 24 5 13 26 16" />
    <elem eid="7" shape="hex8" nodes="21 18 22 27 13 6 14 26" />
    <elem eid="8" shape="hex8" nodes="27 22 19 23 26 14 7 15" />
  </elset>
  <ccxmodeldefinition><![CDATA[**Revisited Mazars
**<YoungModulus>, <PoissonRatio>, <Ac>, <At>, <Bc>, <Bt>, <k>, <ed0>
*Material, name=@CALCULIXBEHAVIOUR_MAZARS
*Depvar
3
*User Material, constants=8
3.2e10, 0.2, 1.15, 0.8, 1391.3, 10000.0 , 0.7 , 0.00009375,
0.0
*DENSITY
1000
***DAMPING,ALPHA=0.03465648,BETA=0.0
*SOLID SECTION, ELSET=Mazars, MATERIAL=@CALCULIXBEHAVIOUR_MAZARS

**INITIAL CONDITIONS,TYPE=VELOCITY
**WALL,1,0.
**WALL,2,0.0
**WALL,3,0.0
**@MaterialProperty<constant> 'YoungModulus'     3.2e10; 
**@MaterialProperty<constant> 'PoissonRatio'     0.2;    
**@MaterialProperty<constant> 'Ac' 1.15;
**@MaterialProperty<constant> 'At' 0.8;
**@MaterialProperty<constant> 'Bc' 1391.3;
**@MaterialProperty<constant> 'Bt' 10000.;
**@MaterialProperty<constant> 'k'  0.7;
**@MaterialProperty<constant> 'ed0' 9.375e-5;
***DV1=@StateVariable real d;
**SDV2=@StateVariable real Y;  Maximum equivalent strain
**SDV3=@StateVariable real eeqcor; Equivalent Strain . It's a Lagrange (Green) Strain
]]></ccxmodeldefinition>
  <ccxstep><![CDATA[*EL FILE
SDV

]]></ccxstep>
  <ccxmodifykeyword action="setparameter" keyword="STATIC" setparameter="SOLVER" value="PARDISO" omitparameter="" />
  <ccxmodifykeyword action="setparameter" keyword="DYNAMIC" setparameter="EXPLICIT" value="" omitparameter="" suppressed="true" />
  <elementtype selection="" shape="hex8" type="C3D8" />
  <displacement selection="displacement_nodes">
    <value>0 mm</value>
    <x>0</x>
    <y>0</y>
    <z>1</z>
  </displacement>
  <group name="Tension-Compression" />
  <surfacepressure selection="pressure_faces" name="S1-TC" group="Tension-Compression">
    <normal unit="MPa" timeunit="s" positionunit="mm">10*t</normal>
  </surfacepressure>
  <surfacepressure selection="pressure_faces(2)" name="S2-TC" group="Tension-Compression">
    <normal unit="MPa" timeunit="s" positionunit="mm">-10*t</normal>
  </surfacepressure>
  <displacement selection="displacement_nodes(2)">
    <value>0 mm</value>
    <x>1</x>
    <y>0</y>
    <z>0</z>
  </displacement>
  <ccxmodifykeyword action="setparameter" keyword="DYNAMIC" setparameter="ALPHA" value="-0.1" omitparameter="" suppressed="true" />
  <ccxmodifykeyword action="setparameter" keyword="STEP" setparameter="Inc" value="50" omitparameter="" suppressed="true" />
  <displacement selection="displacement_nodes(3)">
    <value>0 mm</value>
    <x>0</x>
    <y>1</y>
    <z>0</z>
  </displacement>
  <nodeselection name="WALL">
    <node nid="3" />
    <node nid="15" />
    <node nid="6" />
    <node nid="7" />
    <node nid="4" />
    <node nid="1" />
    <node nid="10" />
    <node nid="8" />
    <node nid="14" />
    <node nid="5" />
    <node nid="2" />
    <node nid="13" />
    <node nid="11" />
    <node nid="16" />
    <node nid="9" />
    <node nid="12" />
    <node nid="17" />
    <node nid="18" />
    <node nid="19" />
    <node nid="20" />
    <node nid="21" />
    <node nid="22" />
    <node nid="23" />
    <node nid="24" />
    <node nid="25" />
    <node nid="26" />
    <node nid="27" />
  </nodeselection>
  <faceselection name="Reaction_Compression" />
  <elementselection name="Mazars">
    <element eid="2" />
    <element eid="3" />
    <element eid="5" />
    <element eid="1" />
    <element eid="4" />
    <element eid="7" />
    <element eid="8" />
    <element eid="6" />
  </elementselection>
  <elementselection name="EALL">
    <element eid="2" />
    <element eid="3" />
    <element eid="5" />
    <element eid="1" />
    <element eid="4" />
    <element eid="7" />
    <element eid="8" />
    <element eid="6" />
  </elementselection>
  <faceselection name="pressure_faces" direct="true">
    <face eid="8" faceid="4" />
    <face eid="4" faceid="4" />
    <face eid="7" faceid="4" />
    <face eid="3" faceid="4" />
    <face eid="6" faceid="6" />
    <face eid="2" faceid="6" />
    <face eid="1" faceid="6" />
    <face eid="5" faceid="6" />
  </faceselection>
  <faceselection name="pressure_faces(2)" direct="true">
    <face eid="6" faceid="3" />
    <face eid="2" faceid="3" />
    <face eid="7" faceid="3" />
    <face eid="3" faceid="3" />
    <face eid="1" faceid="5" />
    <face eid="5" faceid="5" />
    <face eid="8" faceid="5" />
    <face eid="4" faceid="5" />
  </faceselection>
  <nodeselection name="displacement_nodes" direct="true">
    <node nid="27" />
  </nodeselection>
  <nodeselection name="displacement_nodes(2)" direct="true">
    <node nid="26" />
    <node nid="25" />
  </nodeselection>
  <nodeselection name="displacement_nodes(3)" direct="true">
    <node nid="27" />
  </nodeselection>
  <formula name="Eq_Langrange_Strain_SDV3" displayunit="" fieldscalar="formula(8)">
    <unit exponents=""></unit>
    <formula>max(0.00009375,sqrt(heaviside(e1)*e1^2+heaviside(e2)*e2^2+heaviside(e3)*e3^2))</formula>
  </formula>
  <fieldvariable name="displacement" displayunit="mm" />
  <fieldvariable name="principal strain" displayunit="" />
  <fieldvariable name="von Mises stress" displayunit="MPa" />
  <fieldvariable name="strain energy density" displayunit="MPa" />
  <fieldvariable name="principal stress" displayunit="MPa" />
  <fieldvariable name="stress" displayunit="MPa" />
  <fieldvariable name="strain" displayunit="" />
  <fieldvariable name="external force" displayunit="kN" />
  <formula name="Time" displayunit="" fieldscalar="formula">
    <unit exponents="+0+0+1">s</unit>
    <formula>t</formula>
  </formula>
  <solution configuration="ELU">
    <analysis type="L30" solver="ccx" />
    <elset name="Default" colorrgb="DADADA" expanded="true" />
  </solution>
</liml>