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

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

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

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


Image