![]() |
MechSys
1.0
Computing library for simulations in continuum and discrete mechanics
|
#include <unsatflow.h>
Public Types | |
enum | WRC_t { BC_t, ZI_t, HZ_t, PW_t } |
Public Member Functions | |
UnsatFlow (int NDim, SDPair const &Prms, Model const *EquilibMdl=NULL) | |
void | InitIvs (SDPair const &Ini, State *Sta) const |
Initialize internal values. | |
void | Update (double Dpw, double DEv, UnsatFlowState *Sta) |
Update state. | |
void | TgVars (UnsatFlowState const *Sta) const |
Calculate c, C, chi, and kwb. | |
void | TgVars (UnsatFlowState const *Sta, Vec_t const &Ww, double &Cpw, double &Cvs, Vec_t &fwd) const |
void | TgIncs (UnsatFlowState const *Sta, double Dpw, double DEv, double &DSw, double &Dchi) const |
double | FindSw (double pc) |
Find Sw corresponding to pc by integrating from (pc,Sw)=(0,1) to pc (disregarding Dev) | |
double | Findpc (double Sw) |
Find pc corresponding to Sw by integrating from (pc,Sw)=(0,1) to Sw (disregarding Dev) | |
double | rw (double Sw) const |
kwsat multiplier | |
void | GenCurve (Array< double > pcs, char const *Filekey, size_t Np=100) |
Generate WRC. | |
void | Stiffness (State const *Sta, UnsatFlowState const *FSta, Mat_t &D) const |
Coupled "effective" stiffness. | |
void | RateAndUpdate (size_t Idx, Vec_t const &d, double dpwdt, double dt, UnsatFlowState *FSta, EquilibState *Sta) const |
Public Attributes | |
Model const * | EMdl |
Mat_t | kwsatI |
Inverse of saturated kw. | |
Mat_t | kwsatb |
Saturated kw bar. | |
double | akw |
Exponent of kw model. | |
WRC_t | WRC |
Water retention curve model. | |
String | Name |
Model name. | |
double | Cw |
Compressibility of water. | |
double | bc_lam |
double | bc_sb |
double | bc_wr |
Brooks & Corey model parameters. | |
double | zi_del |
double | zi_bet |
double | zi_gam |
double | zi_a |
double | zi_b |
double | zi_alp |
Zienkiewicz et al. 1990 (zi_bet and zi_a => [m]) | |
double | hz_a |
double | hz_b |
double | hz_A |
double | hz_B |
Huang and Zienkiewicz (Gawin, Schrefler ...) | |
double | ld |
double | xRd |
double | yR |
double | xRw |
double | bd |
double | bw |
double | b1 |
Pedroso and Williams (PW) | |
double | y0 |
double | c1d |
double | c2d |
double | c3d |
double | c1w |
double | c2w |
double | c3w |
double | c |
double | C |
double | chi |
double | Cpc |
double | Ceps |
HMStressUpdate | HMSUp |
enum UnsatFlow::WRC_t |
UnsatFlow::UnsatFlow | ( | int | NDim, |
SDPair const & | Prms, | ||
Model const * | EquilibMdl = NULL |
||
) | [inline] |
double UnsatFlow::Findpc | ( | double | Sw | ) | [inline] |
Find pc corresponding to Sw by integrating from (pc,Sw)=(0,1) to Sw (disregarding Dev)
double UnsatFlow::FindSw | ( | double | pc | ) | [inline] |
Find Sw corresponding to pc by integrating from (pc,Sw)=(0,1) to pc (disregarding Dev)
void UnsatFlow::GenCurve | ( | Array< double > | pcs, |
char const * | Filekey, | ||
size_t | Np = 100 |
||
) | [inline] |
Generate WRC.
void UnsatFlow::InitIvs | ( | SDPair const & | Ini, |
State * | Sta | ||
) | const [inline, virtual] |
Initialize internal values.
Implements Model.
void UnsatFlow::RateAndUpdate | ( | size_t | Idx, |
Vec_t const & | d, | ||
double | dpwdt, | ||
double | dt, | ||
UnsatFlowState * | FSta, | ||
EquilibState * | Sta | ||
) | const [inline] |
double UnsatFlow::rw | ( | double | Sw | ) | const [inline] |
kwsat multiplier
void UnsatFlow::Stiffness | ( | State const * | Sta, |
UnsatFlowState const * | FSta, | ||
Mat_t & | D | ||
) | const [inline] |
Coupled "effective" stiffness.
void UnsatFlow::TgIncs | ( | UnsatFlowState const * | Sta, |
double | Dpw, | ||
double | DEv, | ||
double & | DSw, | ||
double & | Dchi | ||
) | const [inline] |
void UnsatFlow::TgVars | ( | UnsatFlowState const * | Sta | ) | const [inline] |
Calculate c, C, chi, and kwb.
void UnsatFlow::TgVars | ( | UnsatFlowState const * | Sta, |
Vec_t const & | Ww, | ||
double & | Cpw, | ||
double & | Cvs, | ||
Vec_t & | fwd | ||
) | const [inline] |
void UnsatFlow::Update | ( | double | Dpw, |
double | DEv, | ||
UnsatFlowState * | Sta | ||
) | [inline] |
Update state.
double UnsatFlow::akw |
Exponent of kw model.
double UnsatFlow::b1 |
Pedroso and Williams (PW)
double UnsatFlow::bc_lam |
double UnsatFlow::bc_sb |
double UnsatFlow::bc_wr |
Brooks & Corey model parameters.
double UnsatFlow::bd |
double UnsatFlow::bw |
double UnsatFlow::c [mutable] |
double UnsatFlow::C [mutable] |
double UnsatFlow::c1d |
double UnsatFlow::c1w |
double UnsatFlow::c2d |
double UnsatFlow::c2w |
double UnsatFlow::c3d |
double UnsatFlow::c3w |
double UnsatFlow::Ceps [mutable] |
double UnsatFlow::chi [mutable] |
double UnsatFlow::Cpc [mutable] |
double UnsatFlow::Cw |
Compressibility of water.
Model const* UnsatFlow::EMdl |
HMStressUpdate UnsatFlow::HMSUp [mutable] |
double UnsatFlow::hz_a |
double UnsatFlow::hz_A |
double UnsatFlow::hz_B |
Huang and Zienkiewicz (Gawin, Schrefler ...)
double UnsatFlow::hz_b |
Saturated kw bar.
Inverse of saturated kw.
double UnsatFlow::ld |
Water retention curve model.
double UnsatFlow::xRd |
double UnsatFlow::xRw |
double UnsatFlow::y0 |
double UnsatFlow::yR |
double UnsatFlow::zi_a |
double UnsatFlow::zi_alp |
Zienkiewicz et al. 1990 (zi_bet and zi_a => [m])
double UnsatFlow::zi_b |
double UnsatFlow::zi_bet |
double UnsatFlow::zi_del |
double UnsatFlow::zi_gam |