Skip to content

Hypoplasticity + Intergranular Strain Anisotropy (HYPO+ISA):

*Mechanical = Hypo-Isa   
phi [rad], h_s [F/A], n [-], e_d0 [-], e_c0 [-], e_i0 [-], alpha [-], beta [-] 
m_R [-], R [-], beta_h0 [-], chi^0 [-], chi^max [-], epsilon^acc [-], c^z [-], beta_hmax [-]
E_ph [F/A], nu_ph [-], p_{min,ph} [F/A], K^w [F/A], IM

The material parameters of the hypoplastic model with intergranular strain are given in the following.

  • \(\varphi_c\) [rad] is the critical friction angle and has to be provided in radian.
  • \(h_s\) [F/A] and \(n\) [-] describe the decrease of the maximum void ratio \(e_i\) with increasing pressure during isotropic compression. This parameter has no unit.
  • The void ratios \(e_{d0}\) [-], \(e_{c0}\) [-] and \(e_{i0}\) [-] correspond to the densest, the critical and the loosest state at pressure \(p\) = 0 F/A. These parameters have no unit.
  • The material constant \(\alpha\) describes the influence of density on the peak friction angle. This parameter has no unit.
  • \(\beta\) mainly influences the stiffness of dense soil skeletons. This parameter has no unit.
  • \(m_R\) governs the increase of the stiffness due to changes in the direction of the strain path. The parameter has no unit.
  • \(R\) is the radius of intergranular strain. This parameter has no unit.
  • \(\beta_{h0}\) and \(\beta_{hmax}\) are the minimum and maximum hardening parameters of the intergranular strain and describe the evolution of the intergranular strain tensor along a strain path. These parameters have no unit.
  • \(\chi^0\) and \(\chi^{max}\) are the minimum and maximum intergranular strain exponent. These parameters have no unit.
  • \(\epsilon^{acc}\) is the accumulation rate factor. This parameter has no unit.
  • \(c_z\) is the cyclic mobility factor. This parameter has no unit.
  • \(E_{ph}\) and \(\nu_{ph}\) are the parameters of the phantom elasticity. The implementation of W. Fuentes uses a phantom linear elasticity overlying the constitutive model response of the Hypo-ISA model with \(E_{ph}=40\) kPa and \(\nu_{ph}=0.3\). We have also included this in our implementation, but the parameters must be specified in the input for transparency. We recommend turning them off by default (\(E_{ph}=0\)). Deviating from the original implementation from W. Fuentes, the phantom elasticity in the numgeo implementation is only active when the mean effective stress is equal to or drops below \(p_{min,ph}\).
  • \(\bar{K}^w\) is the bulk modulus of pore water for locally undrained simulations (and only for those). If no locally undrained simulation is to be performed, \(K^w=0\) should be set.
  • \(IM\) is a flag controlling which integration method should be used to advance the stress from time \(t_0\) to time \(t=t_0 + \Delta t\). \(IM=1\) corresponds to a Forward Euler (FE) integration scheme, with \(IM=2\) a modified FE is used with (sub)substepping for the integration of the ISA part of the model. A Modified Euler integration is currently being tested.

State variables

In addition to the (effective) stress, the hypoplastic constitutive model takes additional state variables. Performing simulations with this model requires the prescription of those. The following state variables are contained in the model:

  • Void ratio: \(e\), void_ratio. The prescription of the void ratio is mandatory, if the initialization is omitted or wrong values are prescribed the simulation will abort.
  • Intergranular strain anisotropy: \(R_{ij}\), int_strain11, int_strain22, int_strain33, int_strain12, int_strain23, int_strain13. Only in combination with the intergranular strain extension. Prescribing \(R_0\) is not mandatory, if omitted, \(R_{ij}\) is assumed. \(R_{ii} = R_0 / \sqrt{3}\) corresponds to an isotropically consolidated state.

The state variables can be initialized in two different ways:

  • Using the *Initial conditions, type = state variables command in the input file, e.g.

`Initial conditions, type = state variables
element_set_name, void_ratio, <value>
`  
element_set_name, int_strain11, <value>
element_set_name, int_strain22, <value>
element_set_name, int_strain33, <value>
* Using the interface for user-defined initial state variables. In this approach, the vector holding the state variables is filled directly in a fortran subroutine as described in Section User defined state variables. The position of the state variables in the associated vector are as follows:

  • statev(1) = \(e\)
  • statev(3) = \(R_{11}\)
  • statev(4) = \(R_{22}\)
  • statev(5) = \(R_{33}\)
  • statev(6) = \(R_{12}\)
  • statev(7) = \(R_{23}\)
  • statev(8) = \(R_{13}\)

Additional output variables

The following additional output variables are available in the *Output command:

  • Void ratio: void_ratio
  • Intergranular strain: int_strain11, int_strain22, int_strain33, int_strain12, int_strain23, int_strain13.
  • Stress projection: proj. Flag that indicates if the stress correction is active (see Theory Manual). proj takes the value of 0 if the stress projection is inactive.

Implementation

The Hypo+ISA model was developed by W. Fuentes (wiliam.fuentes@uni-norte.de), see Poblete et al. (2016)1 and Fuentes et al. (2020)2 for details on the model. Since version 2023 of numgeo, we use our own implementation of the model.

  1. M. Poblete, W. Fuentes, and given-i=Th family=Triantafyllidis, given=Th. On the simulation of multidimensional cyclic loading with intergranular strain. Acta Geotechnica, 11(6):1263–1285, 2016. URL: http://link.springer.com/10.1007/s11440-016-0492-2 (visited on 2021-08-12), doi:10.1007/s11440-016-0492-2. ↩

  2. William Fuentes, Torsten Wichtmann, Melany Gil, and Carlos Lascarro. ISA-Hypoplasticity accounting for cyclic mobility effects for liquefaction analysis. Acta Geotechnica, 15(6):1513–1531, 2016. URL: http://link.springer.com/10.1007/s11440-019-00846-2 (visited on 2021-08-12), doi:10.1007/s11440-019-00846-2. ↩