#
# run with  'mefel drained.in -kwd=1'
#
#------------------------
#  Problem description  |
#------------------------
Simple axisymmetric triaxial test with constant radial stress, temperature and suction on hypoplastic model from D. Masin, time dependent problem
mespr 1
problemtype mech_timedependent_prob

straincomp 1 
strainpos  2 
strainaver 0 

stresscomp 1 
stresspos  2 
stressaver 0 

othercomp 1 
otherpos  2 
otheraver 0 

reactcomp  1

adaptivity      0 
stochasticcalc  0 
homogenization  0 
noderenumber    0

time_contr_type  adaptivemin
start_time       0.0
end_time         2000.0
num_imp_times    0
funct_type       stat
const_val        1.0
dtmin            1.0e-20
timetypeprin     seconds

hdbackup         nohdb


#solver data
nr_num_iter      200 
nr_error         1.0e-6 
check_div        off

stiffmatstor   double_skyline
stiffmat_type  tangent_stiff
typelinsol     lu



#------------------------
#  Definition of nodes  |
#------------------------
4 #number of nodes

# node_id, x,  y,  z, num_dof, cross_sec_type=nocrsec, loc_coord_sys
1   0.0  1.0  0.0  2  0 0
2   0.0  0.0  0.0  2  0 0
3   1.0  0.0  0.0  2  0 0
4   1.0  1.0  0.0  2  0 0



#----------------------------------------------
#  Definition of boundary conditons at nodes  |
#----------------------------------------------

4 #number of boundary conditions
# nod_id cond_dof1 cond_dof2
# cond_dof < 0 represents prescribed displacements defined later in the load case section
1  0 -1
2  0  0
3  1  0
4  1 -1 



#----------------------------
#  Definition of elements   |
#----------------------------

1 #number of elements

# elem_id, elem_type=axisymlq=63, node_1, node_2, node_3, elem_code_num_flag, cross_sec_type=nocrsec, num_mat_types, mat_type1=hypoplastusatthermmat=421, mat_id1
1   63      1 2 3 4      0 0    1 421 1



#----------------------------
#  Definition of materials  |
#----------------------------

1 #number of different material types

# hypoplastusattherma_mat num_param_sets 
421 1
# param_set_id, 
1
# phi, lam_star, kap_star, n_star,  
  27.0   0.08     0.008     1.05
# nu    ,   ns  ,  ls  ,   nt  ,   lt  ,   m  , 
  0.25    0.025   0.0    -0.01    0.0    2.0
# alpha_s , kappa_m,  sm_star ,  em_star , csh , se_ref ,  em_ref,  
  3.5e-5     0.04     -2400.0     0.38     0.0   -200.0    0.18     
# tref  ,   at   ,    bt   ,  aer  , p_t
  294.0   0.118   -1.54e-4    0.25   0.0
# prescribed suction function
1 # yes
0 # constant function
-1900.0

# prescribed temperature
1 # yes
0 # constant function
278.0

# stress integration algorithm
11 # Runge-Kutta-Fehlberg
# RKF_type ni, err, h_min, 
4 10000 1.0e-7 1.0e-17



#--------------------------------
#  Definition of cross-section  |
#--------------------------------

0 # number of cross_sec_types



#--------------
#    Load     |
#--------------

1 # type of load = time independent load in subloadcases
1 # number of time dependent load cases
2 # number of subloadcases

#
# subloadcase 1 - proportional load on the top edge and confining pressure
#
0 # number of loaded nodes
0 # number of loaded elements
1 # number of prescribed displacements
-0.0001
0 # temperature changes will not be assumed
1 t # time function defined by parser -> time equals magnitude of the total strain

#
# subloadcase 2 - constant load due to initial stresses
#
0 # number of loaded nodes
1 # number of loaded elements
1  2   0 0 1 -200.0 0.0 -200.0 0.0   0
1 # number of prescribed displacements
0.0
0 # temperature changes will not be assumed
0 1.0 # constant time function


#---------------------------------------------
#  Definition of initial conditons at nodes  |
#---------------------------------------------

4 # number of initial conditions at nodes
# node_id, initial_cond_type=inicond, num_of_init_values, e_0, ascan_0
1 16 2  0.6  0.5
2 16 2  0.6  0.5
3 16 2  0.6  0.5
4 16 2  0.6  0.5



4 # eigenstresses will be assumed
# folowing lines defines hydrostatic pressure stress state -200.0 kPa
20 # axisymmetric strain/stress state
1 1 1 2 # indeces of general functions for each eigenstress component
2 # total number of eigenstress general functions used
0 -200.0 # the first function is constant
0    0.0 # the second function is constant
# the above setup results in assignment of stress vector sig={sig_x, sig_z, sig_r, tau_xz}^T = {-200.0, -200.0, -200.0, 0.0}^T to all integration points




#outdriver
#text output
0

#GiD output
3
drained
4 10 1
1 1
0
0
0
0

4 10 1
1 1 0
1 1 0
0 




# table text output
1
drained.dat
16 1
2 1 1 2 1  # eps_x
2 1 1 2 2  # eps_y
2 1 1 2 3  # eps_r
2 1 1 2 4  # gamma_xy
2 1 1 3 1  # sig_x
2 1 1 3 2  # sig_y
2 1 1 3 3  # sig_r
2 1 1 3 4  # tau_xy
2 1 1 8 12 # porosity e
2 1 1 8 13 # suction
2 1 1 8 14 # degree of saturation
2 1 1 8 23 # epsv
2 1 1 8 24 # e_ax
2 1 1 8 25 # q_ax
2 1 1 8 21 # dtsub
1 1 7
