SIFEL Home page

Exacavtion problem with hypoplastic model for unsaturated expansive solis

All problems has simple geometry involving only one element 1x1 m with axis of symmetry in vertical direction at the origin of coordinate sytem x-y. Therefor input files were prepared directly without the preprocesor usage.

1. Preparation of mesh file:
   1. Generation of the finite element mesh is carried out with the help of input file slope.t3d.in by the command
      T3d -p 264 -r 1 -i slope.t3d.in -o slope.t3d -d 2.0 -X -$
   2. Resulting mesh file can be downloaded here and thus this step can be skipped
2. Generation of MEFEL input file:
   1. Compile the MEFEL preprocessor in the directory SIFEL/MEFEL/PREP by statement make.
   2. Copy the preprocessor mechprep from the directory SIFEL/MEFEL/PREP/_DBG to the working directory.
   3. Generate the input file slope_hypopl_excav.in by statement mechprep slope_hypopl_excav.pr slope_hypopl_excav.in
   4. Resulting MEFEL input file can be downloaded here and thus this step can be skipped
3. Solution of the problem:
   1. Compile the MEFEL in the directory SIFEL/MEFEL/SRC by statement make.
   2. Copy the executable file mefel from the directory SIFEL/MEFEL/SRC/_DBG to the working one.
   3. Solve problem by the execution of the following commands:
      mefel slope_hypopl_excav.in

Results

Resulting distribution of vertical stress sig_y (left) vertical displacement component (right) at initial stage.

Resulting distribution of vertical stress sig_y (left) vertical displacement component (right) at 5. stage.

Resulting distribution of vertical stress sig_y (left) vertical displacement component (right) at final stage.