################################################################################
#section la plus generale
################################################################################
Section MultiScale RealUnits

UNITCODE RealUnits
DIMENSION 3

LET stress = 624.15096e-6 # = 100 MPa in eV per Å^3
LET shearstrain = $ELSTRAINO*2


# model declaration
ELAST_CODE AKANTU fem_model

COMPUTE disp EXTRACT INPUT fem_model FIELD position
COMPUTE shear_disp PYTHON INPUT disp FILENAME stimulate_shear
STIMULATION stim_shear FIELD INPUT dofs=fem_model INPUT field=shear_disp FIELD displacement  ADDITIVE ONESHOT 0
STIMULATION damp ZERO INPUT fem_model FIELD velocity COEF .997

COMPUTE strain EXTRACT INPUT fem_model FIELD strain
COMPUTE mean_strain REDUCE INPUT strain OPERATOR AVERAGE

DUMPER strain_dump TEXT INPUT mean_strain PRECISION 12 FREQ 10

endSection



################################################################################
## FEM section
Section AKANTU:fem_model RealUnits
DOMAIN_GEOMETRY 1
MESH_FILENAME cube_3d.msh
MATERIAL_FILENAME material.dat
TIMESTEP 10
PBC 1 0 1
PRESTRESSED_PERIODICITY
TRACTION_BNDY_TAG shear
TRACTION_BNDY_TENSOR 0 stress 0 stress 0 0 0 0 0
endSection