# Excavation of slope with mesh generated by T3D and hypoplastic material model for unsaturated soils, load and stiffness are in kPa begsec_files slope.t3d mesh_format t3d edge_numbering 0 read_mat_strings no read_mat_kwd yes read_mat_strings no read_crsec_kwd yes endsec_files begsec_mater num_mat_types 2 mattype elisomat num_inst 1 1 e 10e6 nu 0.35 mattype hypoplastusatthermat num_inst 1 # phi, lam_star, kap_star, n_star, nu 1 phi 25 lam_star 0.08 kap_star 0.01 N_star 1.56 nu 0.25 # nu , ns , l , nt , lt , m , n 0.0 l 0.0 nT 0.0 lT 0.0 m 10.0 # alpha_s , kappa_m, sm_star , em_star , csh , alpha_s 0.0 kap_m 0.0 sm_star -140.0e3 em_star 0.0 csh 0.002 # se_ref , em_ref, tref , at , bt , aer , p_t s_airentry0 -200.0 em0 0.50 tr 294.0 at 0.118 bt 0.0 aer 0.5 p_t 1.0 # prescribed suction function presc_suction yes funct_type stat # constant function const_val 0.0 # prescribed temperature presc_tempr yes funct_type stat # constant function const_val 294.0 # stress integration algorithm rkfsra # Runge-Kutta-Fehlsberg # RKF_type=rkf23bst ni, err, h_min, rktype rkf23bst ni 10000 err 1.0e-7 h_min 1.0e-17 endsec_mater begsec_crsec num_crsec_types 1 crstype csplanestr num_inst 1 1 1.0 endsec_crsec #problem description begsec_probdesc Slope stability mespr 1 problemtype growing_mech_structure # straincomp 1 strainpos 2 strainaver 1 # stresscomp 1 stresspos 2 stressaver 1 # othercomp 0 # reactcomp 1 # adaptivity 0 stochasticcalc 0 homogenization 0 noderenumber no_renumbering # time_contr_type adaptivemin start_time 0.0 end_time 11.0 # 11 stages num_imp_times 32 # 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 6.5 7.0 7.5 8.0 8.5 9.0 9.5 10.0 10.5 11.0 11.1 11.2 11.3 11.4 11.5 11.6 11.7 11.8 11.9 12.0 funct_type tab approx_type piecewiseconst ntab_items 4 0.0 0.5 1.0 0.5 11.0 0.1 100.0 0.1 dtmin 0.01 timetypeprin seconds comp_inidispl no smooth_removal yes cpincrnr 0.1 cpminincrnr 0.001 hdbackup nohdb nr_num_iter 20 nr_error 1.0e-3 nr_solvtype modnewton stiffmatstor double_skyline stiffmat_type ijth_tangent_stiff 1 10 #stiffmat_type tangent_stiff typelinsol lu endsec_probdesc begsec_loadcase # # Load case definitions # num_loadcases 1 lc_id 1 num_sublc 1 # time function for the first subloadcase tfunc_lc_id 1 tfunc_slc_id 1 funct_type tab approx_type linear ntab_items 4 0.0 1.0 11.0 1.0 12.0 4.0 100.0 4.0 #temperature load type for the first subloadcase tempr_type_lc_id 1 tempr_type_slc_id 1 temp_load_type 0 # prescribed displacements used in nodal dof time functions for the first subloadcase num_pres_displ_lc_id 1 num_pres_displ_slc_id 1 num_presc_displ 0 # # Definice casovych funkci pro prvky, uzly a podpory # # time functions controlling dofs and elements time_functions # 10 stages of excavation + # 1 time functions for remaining mass + # 1 function for supports num_gfunct 12 # 1st layer (top) gf_id 1 funct_type itab nitab_items 3 0.0 1 1.0 0 100.0 0 # 2nd layer gf_id 2 funct_type itab nitab_items 3 0.0 1 2.0 0 100.0 0 # 3rd layer gf_id 3 funct_type itab nitab_items 3 0.0 1 3.0 0 100.0 0 # 4th layer gf_id 4 funct_type itab nitab_items 3 0.0 1 4.0 0 100.0 0 # 5th layer gf_id 5 funct_type itab nitab_items 3 0.0 1 5.0 0 100.0 0 # 6th layer gf_id 6 funct_type itab nitab_items 3 0.0 1 6.0 0 100.0 0 # 7th layer gf_id 7 funct_type itab nitab_items 3 0.0 1 7.0 0 100.0 0 # 8th layer gf_id 8 funct_type itab nitab_items 3 0.0 1 8.0 0 100.0 0 # 9h layer gf_id 9 funct_type itab nitab_items 3 0.0 1 9.0 0 100.0 0 # 10th layer gf_id 10 funct_type itab nitab_items 3 0.0 1 10.0 0 100.0 0 # remaining mass gf_id 11 funct_type itab nitab_items 3 0.0 1 11.0 1 100.0 1 # supports gf_id 12 funct_type itab nitab_items 2 0.0 0 100.0 0 endsec_loadcase begsec_nodsurfpr ndofn 2 propid 1 ndofn 2 propid 10 ndofn 2 propid 11 ndofn 2 propid 12 ndofn 2 propid 13 ndofn 2 propid 14 ndofn 2 propid 15 ndofn 2 propid 16 ndofn 2 propid 17 ndofn 2 propid 18 ndofn 2 propid 19 nod_inicond propid 1 lc_id 1 cond ini_cd_type inicond nval 2 # e_ini, ascan_ini 0.53 0.5 nod_inicond propid 10 lc_id 1 cond ini_cd_type inicond nval 2 # e_ini, ascan_ini 0.53 0.5 nod_inicond propid 11 lc_id 1 cond ini_cd_type inicond nval 2 # e_ini, ascan_ini 0.53 0.5 nod_inicond propid 12 lc_id 1 cond ini_cd_type inicond nval 2 # e_ini, ascan_ini 0.53 0.5 nod_inicond propid 13 lc_id 1 cond ini_cd_type inicond nval 2 # e_ini, ascan_ini 0.53 0.5 nod_inicond propid 14 lc_id 1 cond ini_cd_type inicond nval 2 # e_ini, ascan_ini 0.53 0.5 nod_inicond propid 15 lc_id 1 cond ini_cd_type inicond nval 2 # e_ini, ascan_ini 0.53 0.5 nod_inicond propid 16 lc_id 1 cond ini_cd_type inicond nval 2 # e_ini, ascan_ini 0.53 0.5 nod_inicond propid 17 lc_id 1 cond ini_cd_type inicond nval 2 # e_ini, ascan_ini 0.53 0.5 nod_inicond propid 18 lc_id 1 cond ini_cd_type inicond nval 2 # e_ini, ascan_ini 0.53 0.5 nod_inicond propid 19 lc_id 1 cond ini_cd_type inicond nval 2 # e_ini, ascan_ini 0.53 0.5 endsec_nodsurfpr begsec_nodedgpr nod_tfunc propid 11 ndir 2 dir 1 tfunc_id 12 dir 2 tfunc_id 12 nod_tfunc propid 10 ndir 1 dir 1 tfunc_id 12 endsec_nodedgpr #element properties begsec_elsurfpr # slope mass el_type propid 1 planeelementlt strastrestate planestrain el_mat propid 1 num_mat 1 type hypoplastusatthermat type_id 1 el_crsec propid 1 type csplanestr type_id 1 el_tfunc propid 1 tfunc_id 11 volume_load propid 1 lc_id 1 slc_id 1 ncomp 2 func_type stat coord_sys 1 load_comp 0.0 -20.0 el_eigstr propid 1 str_type stress ncomp 4 eigstr_comp funct_type pars func_formula 20.0/2.0*(y-40.0) funct_type pars func_formula 20.0*(y-40.0) funct_type stat const_val 0.0 funct_type pars func_formula 20.0/2.0*(y-40.0) # 1st excaved layer el_type propid 10 planeelementlt strastrestate planestrain el_mat propid 10 num_mat 1 type hypoplastusatthermat type_id 1 el_crsec propid 10 type csplanestr type_id 1 el_tfunc propid 10 tfunc_id 1 volume_load propid 10 lc_id 1 slc_id 1 ncomp 2 func_type stat coord_sys 1 load_comp 0.0 -20.0 el_eigstr propid 10 str_type stress ncomp 4 eigstr_comp funct_type pars func_formula 20.0/2.0*(y-40.0) funct_type pars func_formula 20.0*(y-40.0) funct_type stat const_val 0.0 funct_type pars func_formula 20.0/2.0*(y-40.0) # 2nd excaved layer el_type propid 11 planeelementlt strastrestate planestrain el_mat propid 11 num_mat 1 type hypoplastusatthermat type_id 1 el_crsec propid 11 type csplanestr type_id 1 el_tfunc propid 11 tfunc_id 2 volume_load propid 11 lc_id 1 slc_id 1 ncomp 2 func_type stat coord_sys 1 load_comp 0.0 -20.0 el_eigstr propid 11 str_type stress ncomp 4 eigstr_comp funct_type pars func_formula 20.0/2.0*(y-40.0) funct_type pars func_formula 20.0*(y-40.0) funct_type stat const_val 0.0 funct_type pars func_formula 20.0/2.0*(y-40.0) # 3rd excaved layer el_type propid 12 planeelementlt strastrestate planestrain el_mat propid 12 num_mat 1 type hypoplastusatthermat type_id 1 el_crsec propid 12 type csplanestr type_id 1 el_tfunc propid 12 tfunc_id 3 volume_load propid 12 lc_id 1 slc_id 1 ncomp 2 func_type stat coord_sys 1 load_comp 0.0 -20.0 el_eigstr propid 12 str_type stress ncomp 4 eigstr_comp funct_type pars func_formula 20.0/2.0*(y-40.0) funct_type pars func_formula 20.0*(y-40.0) funct_type stat const_val 0.0 funct_type pars func_formula 20.0/2.0*(y-40.0) # 4th excaved layer el_type propid 13 planeelementlt strastrestate planestrain el_mat propid 13 num_mat 1 type hypoplastusatthermat type_id 1 el_crsec propid 13 type csplanestr type_id 1 el_tfunc propid 13 tfunc_id 4 volume_load propid 13 lc_id 1 slc_id 1 ncomp 2 func_type stat coord_sys 1 load_comp 0.0 -20.0 el_eigstr propid 13 str_type stress ncomp 4 eigstr_comp funct_type pars func_formula 20.0/2.0*(y-40.0) funct_type pars func_formula 20.0*(y-40.0) funct_type stat const_val 0.0 funct_type pars func_formula 20.0/2.0*(y-40.0) # 5h excaved layer el_type propid 14 planeelementlt strastrestate planestrain el_mat propid 14 num_mat 1 type hypoplastusatthermat type_id 1 el_crsec propid 14 type csplanestr type_id 1 el_tfunc propid 14 tfunc_id 5 volume_load propid 14 lc_id 1 slc_id 1 ncomp 2 func_type stat coord_sys 1 load_comp 0.0 -20.0 el_eigstr propid 14 str_type stress ncomp 4 eigstr_comp funct_type pars func_formula 20.0/2.0*(y-40.0) funct_type pars func_formula 20.0*(y-40.0) funct_type stat const_val 0.0 funct_type pars func_formula 20.0/2.0*(y-40.0) # 6th excaved layer el_type propid 15 planeelementlt strastrestate planestrain el_mat propid 15 num_mat 1 type hypoplastusatthermat type_id 1 el_crsec propid 15 type csplanestr type_id 1 el_tfunc propid 15 tfunc_id 6 volume_load propid 15 lc_id 1 slc_id 1 ncomp 2 func_type stat coord_sys 1 load_comp 0.0 -20.0 el_eigstr propid 15 str_type stress ncomp 4 eigstr_comp funct_type pars func_formula 20.0/2.0*(y-40.0) funct_type pars func_formula 20.0*(y-40.0) funct_type stat const_val 0.0 funct_type pars func_formula 20.0/2.0*(y-40.0) # 7th excaved layer el_type propid 16 planeelementlt strastrestate planestrain el_mat propid 16 num_mat 1 type hypoplastusatthermat type_id 1 el_crsec propid 16 type csplanestr type_id 1 el_tfunc propid 16 tfunc_id 7 volume_load propid 16 lc_id 1 slc_id 1 ncomp 2 func_type stat coord_sys 1 load_comp 0.0 -20.0 el_eigstr propid 16 str_type stress ncomp 4 eigstr_comp funct_type pars func_formula 20.0/2.0*(y-40.0) funct_type pars func_formula 20.0*(y-40.0) funct_type stat const_val 0.0 funct_type pars func_formula 20.0/2.0*(y-40.0) # 8th excaved layer el_type propid 17 planeelementlt strastrestate planestrain el_mat propid 17 num_mat 1 type hypoplastusatthermat type_id 1 el_crsec propid 17 type csplanestr type_id 1 el_tfunc propid 17 tfunc_id 8 volume_load propid 17 lc_id 1 slc_id 1 ncomp 2 func_type stat coord_sys 1 load_comp 0.0 -20.0 el_eigstr propid 17 str_type stress ncomp 4 eigstr_comp funct_type pars func_formula 20.0/2.0*(y-40.0) funct_type pars func_formula 20.0*(y-40.0) funct_type stat const_val 0.0 funct_type pars func_formula 20.0/2.0*(y-40.0) # 9th excaved layer el_type propid 18 planeelementlt strastrestate planestrain el_mat propid 18 num_mat 1 type hypoplastusatthermat type_id 1 el_crsec propid 18 type csplanestr type_id 1 el_tfunc propid 18 tfunc_id 9 volume_load propid 18 lc_id 1 slc_id 1 ncomp 2 func_type stat coord_sys 1 load_comp 0.0 -20.0 el_eigstr propid 18 str_type stress ncomp 4 eigstr_comp funct_type pars func_formula 20.0/2.0*(y-40.0) funct_type pars func_formula 20.0*(y-40.0) funct_type stat const_val 0.0 funct_type pars func_formula 20.0/2.0*(y-40.0) # 10th excaved layer el_type propid 19 planeelementlt strastrestate planestrain el_mat propid 19 num_mat 1 type hypoplastusatthermat type_id 1 el_crsec propid 19 type csplanestr type_id 1 el_tfunc propid 19 tfunc_id 10 volume_load propid 19 lc_id 1 slc_id 1 ncomp 2 func_type stat coord_sys 1 load_comp 0.0 -20.0 el_eigstr propid 19 str_type stress ncomp 4 eigstr_comp funct_type pars func_formula 20.0/2.0*(y-40.0) funct_type pars func_formula 20.0*(y-40.0) funct_type stat const_val 0.0 funct_type pars func_formula 20.0/2.0*(y-40.0) endsec_elsurfpr begsec_outdrv #------------------------------- # Definition of MEFEL output | #------------------------------- # description of output to the text file textout 0 # # description of output to the text file # 3 = GiD format - one huge file outgr_format 3 # # graphics output file name without extension slope_hypopl_excav # sel_nodstep sel_all sel_nodlc sel_all # displ_nodes sel_all displ_comp sel_all # strain_nodes 0 stress_nodes 0 other_nodes 0 force_nodes sel_all force_comp sel_all # sel_elemstep sel_all sel_elemlc sel_all strain_elems sel_all elemstrain_comp sel_mtx elemstra_transfid 0 stress_elems sel_all elemstress_comp sel_all elemstre_transfid 0 other_elems sel_no # # text output of graphs # number of created files with diagrams numdiag 1 slope_hypopl_excav.dat numunknowns 2 sel_diagstep sel_impvalues # 1. column point atnode # point type = node node 1 # point id = 1. node quant_type pr_time # unknown type = load coefficient/time # 2. column point atnode # point type = node node 5496 # point id = 5496. node = top edge of the slope beginning quant_type pr_displ # unknown type = displacement compid 2 # component id 2 = vertical displacement endsec_outdrv