# # run with 'mefel constvolume.in -kwd=1' # #------------------------ # Problem description | #------------------------ Simple axisymmetric specimen with constant suction pressure and temparature subjected to constant increment of volume - hypoplastic model from D. Masin, time dependent problem mespr 1 problemtype mech_timedependent_prob straincomp 1 strainpos 2 strainaver 1 stresscomp 1 stresspos 2 stressaver 1 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 20000.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 incr_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 -2 0 4 -2 -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=hypoplastusattherma_matt=421, mat_id1 1 63 1 2 3 4 0 0 1 421 1 #---------------------------- # Definition of materials | #---------------------------- 2 #number of different material types # elastic material 1 1 1 50.0 0.3 # 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, 1 10000 1.0e-3 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 1 # number of subloadcases # # subloadcase 1 - load will be induced by the change of suction # 0 # number of loaded nodes 0 # number of loaded elements 2 # number of prescribed displacements -0.00001 0.000005 0 # temperature changes will not be assumed 1 t # time function defined by parser -> time equals magnitude of the total strain #--------------------------------------------- # 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 0 #constvolume #1 1 #1 1 #0 #0 #0 #1 1 #1 1 #1 1 0 #1 1 0 #0 # table text output 1 constvolume.dat 13 4 10 # strains 2 1 1 2 1 2 1 1 2 2 2 1 1 2 3 2 1 1 2 4 # stresses 2 1 1 3 1 2 1 1 3 2 2 1 1 3 3 2 1 1 3 4 # state variables 2 1 1 8 12 # e porosity 2 1 1 8 13 # suction 2 1 1 8 14 # degree of saturation Sr 2 1 1 8 21 # dtsub - substep length # time 1 1 7