#include <consol_awf1.h>
Public Member Functions | |
| con_awf1mat () | |
| double | get_alpha () |
| double | get_capww (double pw) |
| double | get_dsw_dpw (double pw) |
| double | get_fw1 (double pw, long ipp) |
| double | get_kintr () |
| double | get_krw (double pw) |
| double | get_kw () |
| double | get_kww (double pw, long ipp) |
| double | get_muw () |
| double | get_othervalue (long compother, double pw, long ipp) |
| double | get_phi () |
| double | get_sw (double pw) |
| void | give_dof_names (namevart *dofname, long ntm) |
| returns ordered dof names | |
| void | give_reqntq (long *antq) |
| marks required non-transport quantities | |
| double | give_saturation_degree (long ipp) |
| double | give_suction (long ipp) |
| double | give_water_pressure (long ipp) |
| void | initval (long ipp) |
| void | matcap (double &cc, long ri, long ci, long ipp) |
| void | matcond (matrix &d, long ri, long ci, long ipp) |
| void | matcond1d (matrix &d, long ri, long ci, long ipp) |
| void | matcond2d (matrix &d, long ri, long ci, long ipp) |
| void | matcond3d (matrix &d, long ri, long ci, long ipp) |
| void | print_othervalue_name (FILE *out, long compother) |
| void | read (XFILE *in) |
| void | rhs1d1 (matrix &d, long ri, long ci, long ipp) |
| void | rhs2d1 (matrix &d, long ri, long ci, long ipp) |
| void | rhs3d1 (matrix &d, long ri, long ci, long ipp) |
| void | rhs_volume (matrix &d, long ri, long ci, long ipp) |
| void | updateval (long ipp) |
| void | values_correction (vector &nv, long ipp) |
| void | water_pressure_check (double &pw, long ipp) |
| ~con_awf1mat () | |
Private Attributes | |
| double | alfa |
| double | alpha |
| double | delta |
| double | e0 |
| double | expn |
| double | gamaw |
| double | hp_min |
| double | htz |
| double | k0 |
| double | kappa |
| double | kintr |
| double | ks |
| double | ksat |
| double | kw |
| double | lam |
| double | lam0 |
| double | m |
| waterflowtype | model_type |
| double | muw0 |
| double | phi0 |
| double | r |
| double | rhow |
| double | sigma_m_bar |
| double | sigma_m_eff_bar |
| double | sigma_m_eff_init |
| double | sigma_m_tot_bar |
| double | sirr |
| double | ssat |
This class defines model for water flow in soils.
Definition at line 11 of file consol_awf1.h.
| con_awf1mat | ( | ) |
| ~con_awf1mat | ( | ) |
Definition at line 83 of file consol_awf1.cpp.
| double get_alpha | ( | ) |
function returns Biot's constant
| alpha | - Biot's constant |
12/9/2008, TKr
Definition at line 554 of file consol_awf1.cpp.
References alpha.
| double get_capww | ( | double | pw | ) |
function creates capacity coefficient of the general material
| pw | - water pressure |
| capww | - capacity coefficient |
12/9/2008, TKr
Definition at line 719 of file consol_awf1.cpp.
References alpha, gardner_exponential, get_dsw_dpw(), get_sw(), ks, kuklik_camclay, kuklik_camclay_mech, kw, lewis_and_schrefler, model_type, phi0, potts_log_linear, and van_genuchten.
Referenced by matcap().
| double get_dsw_dpw | ( | double | pw | ) |
function computes partial derivative of degree of saturation with respect to pw, specific wter content
| pw | - water pressure |
| dsw_dpw | - partial derivative of degree of saturation with respect to pw |
12/9/2008, TKr
Definition at line 425 of file consol_awf1.cpp.
References alfa, delta, expn, gamaw, gardner_exponential, hp_min, htz, lewis_and_schrefler, model_type, potts_log_linear, r, sirr, ssat, and van_genuchten.
Referenced by get_capww().
| double get_fw1 | ( | double | pw, | |
| long | ipp | |||
| ) |
function returns coefficient for righ-hand side of the general material
| pw | - water pressure |
| fw | - first part for right-hand side for continutiy equation |
12/9/2008, TKr
Definition at line 775 of file consol_awf1.cpp.
References gardner_exponential, get_krw(), get_kww(), kintr, kuklik_camclay, kuklik_camclay_mech, lewis_and_schrefler, model_type, muw0, potts_log_linear, rhow, and van_genuchten.
| double get_kintr | ( | ) |
function computes intrinsic permeability
| kintr | - intrinsic permeability |
12/9/2008, TKr
Definition at line 541 of file consol_awf1.cpp.
References kintr.
| double get_krw | ( | double | pw | ) |
function computes water relative permeability
| pw | - water pressure |
| krw | - water relative permeability |
12/9/2008, TKr
Definition at line 476 of file consol_awf1.cpp.
References delta, expn, gamaw, gardner_exponential, get_sw(), hp_min, htz, lewis_and_schrefler, model_type, potts_log_linear, r, and van_genuchten.
| double get_kw | ( | ) |
function returns kw bulk modulus of water [Pa]
| kw | - bulk modulus of water [Pa] |
12/9/2008, TKr
Definition at line 567 of file consol_awf1.cpp.
References kw.
| double get_kww | ( | double | pw, | |
| long | ipp | |||
| ) |
function creates conductivity coefficient of the general material
| pw | - water pressure |
| kww | - conductivity coefficient |
12/9/2008, TKr
Definition at line 594 of file consol_awf1.cpp.
References e0, gamaw, gardner_exponential, get_krw(), transmat::givenontransq(), k0, kappa, kintr, ks, ksat, kuklik_camclay, kuklik_camclay_mech, lam, lewis_and_schrefler, m, mean_stress_eff, model_type, muw0, transmat::nontransq, potts_log_linear, precons_press, sigma_m_eff_bar, sigma_m_eff_init, sigma_m_tot_bar, Tm, and van_genuchten.
Referenced by get_fw1(), matcond1d(), matcond2d(), and matcond3d().
| double get_muw | ( | ) |
function returns muw water viscosity [Pa.s]
| muw | - water viscosity [Pa.s] |
12/9/2008, TKr
Definition at line 580 of file consol_awf1.cpp.
References muw0.
| double get_othervalue | ( | long | compother, | |
| double | pw, | |||
| long | ipp | |||
| ) |
function computes all other variables at nodes
| compother | - number of other components | |
| pw | - water capillary pressure on actual node | |
| ipp | - first integration point on element |
| other | - other variable |
03/03/2011, TKr
Definition at line 829 of file consol_awf1.cpp.
References get_sw(), and other.
Referenced by med1::compute_othervalues().
| double get_phi | ( | ) |
function returns porosity
| phi | - porosity |
12/9/2008, TKr
Definition at line 524 of file consol_awf1.cpp.
References phi0.
| double get_sw | ( | double | pw | ) |
function computes degree of saturation(sorption curve)
| pw | - water pressure |
| sw | - degree of saturation |
12/9/2008, TKr
Definition at line 371 of file consol_awf1.cpp.
References delta, expn, gamaw, gardner_exponential, get_krw(), kuklik_camclay, kuklik_camclay_mech, lewis_and_schrefler, model_type, potts_log_linear, sirr, ssat, and van_genuchten.
Referenced by get_capww(), get_krw(), get_othervalue(), and give_saturation_degree().
| void give_dof_names | ( | namevart * | dofname, | |
| long | ntm | |||
| ) |
returns ordered dof names
The function returns ordered dof names of primary unknowns required by the model.
| dofname | - array of uknown name for particular nodal dofs (output) dofname[i] = name of i-th nodal unknown (for names see aliast.h - enum namevart) | |
| ntm | - number of transported media = number of nodal dof = length of array dofname |
Created by Tomas Krejci according to Tomas Koudelka, 28/04/2014
Definition at line 957 of file consol_awf1.cpp.
References print_err(), and trf_water_press.
Referenced by transmat::give_dof_names().
| void give_reqntq | ( | long * | antq | ) |
marks required non-transport quantities
The funtion marks required non-transport quantities in the array antq.
| antq | - array with flags for used material types antq[i] = 1 => quantity type nontransquant(i+1) is required antq[i] = 0 => quantity type nontransquant(i+1) is not required |
Definition at line 1063 of file consol_awf1.cpp.
References gardner_exponential, kuklik_camclay, kuklik_camclay_mech, lewis_and_schrefler, mean_stress_eff, model_type, potts_log_linear, precons_press, print_err(), and van_genuchten.
Referenced by transmat::give_reqntq().
| double give_saturation_degree | ( | long | ipp | ) |
function computes suction stress (s = -pw = pc);
| ipp | - integration point number |
| suction | - suction stress [Pa] |
25/9/2012, TKr
Definition at line 1044 of file consol_awf1.cpp.
References intpointst::av, get_sw(), transmat::ip, and Tm.
Referenced by transmat::give_saturation_degree().
| double give_suction | ( | long | ipp | ) |
function computes suction stress (s = -pw = pc);
| ipp | - integration point number |
| suction | - suction stress [Pa] |
25/9/2012, TKr
Definition at line 996 of file consol_awf1.cpp.
References intpointst::av, gardner_exponential, transmat::ip, kuklik_camclay, kuklik_camclay_mech, lewis_and_schrefler, model_type, potts_log_linear, suction, Tm, and van_genuchten.
Referenced by transmat::give_suction().
| double give_water_pressure | ( | long | ipp | ) |
function returns water pressure
| ipp | - integration point number |
| eigstr | - eigenstrain for mechanical part |
Definition at line 978 of file consol_awf1.cpp.
References intpointst::av, transmat::ip, and Tm.
Referenced by transmat::give_pore_pressure(), and transmat::give_water_pressure().
| void initval | ( | long | ipp | ) |
This function initializes material model data with respect of consistency parameter gamma.
| ipp | - integration point number |
24/08/2012, TKr
Definition at line 941 of file consol_awf1.cpp.
Referenced by transmat::initvalues().
| void matcap | ( | double & | cc, | |
| long | ri, | |||
| long | ci, | |||
| long | ipp | |||
| ) |
function creates capacity matrix of the material
| c | - capacity matrix of the material | |
| ri | - row index | |
| ci | - column index | |
| ipp | - number of integration point |
Definition at line 198 of file consol_awf1.cpp.
References intpointst::av, get_capww(), transmat::ip, and Tm.
Referenced by med1::matcap().
| void matcond | ( | matrix & | d, | |
| long | ri, | |||
| long | ci, | |||
| long | ipp | |||
| ) |
function computes conductivity matrix of the material in the required integration point
| d | - conductivity matrix of material | |
| ri | - row index | |
| ci | - column index | |
| ipp | - number of integration point |
Definition at line 97 of file consol_awf1.cpp.
References matcond1d(), matcond2d(), matcond3d(), matrix::n, and print_err().
Referenced by med1::matcond().
| void matcond1d | ( | matrix & | d, | |
| long | ri, | |||
| long | ci, | |||
| long | ipp | |||
| ) |
function creates conductivity matrix of the material for 1D problems
| d | - conductivity matrix of the material | |
| ri | - row index | |
| ci | - column index | |
| ipp | - number of integration point |
Definition at line 130 of file consol_awf1.cpp.
References intpointst::av, fillm(), get_kww(), transmat::ip, and Tm.
Referenced by matcond().
| void matcond2d | ( | matrix & | d, | |
| long | ri, | |||
| long | ci, | |||
| long | ipp | |||
| ) |
function creates conductivity matrix of the material for 2D problems
| d | - conductivity matrix of the material | |
| ri | - row index | |
| ci | - column index | |
| ipp | - number of integration point |
Definition at line 151 of file consol_awf1.cpp.
References intpointst::av, fillm(), get_kww(), transmat::ip, and Tm.
Referenced by matcond().
| void matcond3d | ( | matrix & | d, | |
| long | ri, | |||
| long | ci, | |||
| long | ipp | |||
| ) |
function creates conductivity matrix of the material for 3D problems
| d | - conductivity matrix of the material | |
| ri | - row index | |
| ci | - column index | |
| ipp | - number of integration point |
Definition at line 174 of file consol_awf1.cpp.
References intpointst::av, fillm(), get_kww(), transmat::ip, and Tm.
Referenced by matcond().
| void print_othervalue_name | ( | FILE * | out, | |
| long | compother | |||
| ) |
function prints names of all other variables at nodes
| out | - output rhle | |
| compother | - number of other components |
03/03/2011, TKr
Definition at line 864 of file consol_awf1.cpp.
Referenced by med1::print_othervaluesnames().
| void read | ( | XFILE * | in | ) |
function reads parameters
| in | - input file |
12/9/2008, TKr
Definition at line 213 of file consol_awf1.cpp.
References alfa, alpha, delta, e0, expn, gamaw, gardner_exponential, hp_min, htz, k0, kappa, kintr, ks, ksat, kuklik_camclay, kuklik_camclay_mech, kw, lam, lewis_and_schrefler, m, model_type, muw0, phi0, potts_log_linear, r, rhow, sigma_m_eff_bar, sigma_m_eff_init, sigma_m_tot_bar, sirr, ssat, van_genuchten, waterflowtype_kwdset(), and xfscanf().
Referenced by transmat::readmattype().
| void rhs1d1 | ( | matrix & | d, | |
| long | ri, | |||
| long | ci, | |||
| long | ipp | |||
| ) |
function creates volume right-hand side matrix of the material for 1D problems
| d | - right-hand matrix of the material | |
| ri | - row index | |
| ci | - column index | |
| ipp | - number of integration point |
Definition at line 296 of file consol_awf1.cpp.
References intpointst::av, f, fillm(), g, get_fw1(), probdesct::gr1, transmat::ip, Tm, and Tp.
Referenced by rhs_volume().
| void rhs2d1 | ( | matrix & | d, | |
| long | ri, | |||
| long | ci, | |||
| long | ipp | |||
| ) |
function creates volume right-hand side matrix of the material for 2D problems
| d | - right-hand matrix of the material | |
| ri | - row index | |
| ci | - column index | |
| ipp | - number of integration point |
Definition at line 315 of file consol_awf1.cpp.
References intpointst::av, f, fillm(), g, get_fw1(), probdesct::gr1, probdesct::gr2, transmat::ip, Tm, and Tp.
Referenced by rhs_volume().
| void rhs3d1 | ( | matrix & | d, | |
| long | ri, | |||
| long | ci, | |||
| long | ipp | |||
| ) |
function creates volume right-hand side matrix of the material for 3D problems
| d | - right-hand matrix of the material | |
| ri | - row index | |
| ci | - column index | |
| ipp | - number of integration point |
Definition at line 340 of file consol_awf1.cpp.
References intpointst::av, f, fillm(), g, get_fw1(), probdesct::gr1, probdesct::gr2, probdesct::gr3, transmat::ip, Tm, and Tp.
Referenced by rhs_volume().
| void rhs_volume | ( | matrix & | d, | |
| long | ri, | |||
| long | ci, | |||
| long | ipp | |||
| ) |
function computes volume part of right-hand side matrix in the required integration point
| d | - right-hand side matrix of material | |
| ri | - row index | |
| ci | - column index | |
| ipp | - number of integration point |
Definition at line 262 of file consol_awf1.cpp.
References matrix::m, m, rhs1d1(), rhs2d1(), and rhs3d1().
Referenced by med1::rhs_volume().
| void updateval | ( | long | ipp | ) |
This function updates values in the other array reached in the previous equlibrium state to values reached in the new actual equilibrium state.
| ipp | - integration point number in the mechmat ip array. |
24/08/2012, TKr
Definition at line 928 of file consol_awf1.cpp.
Referenced by transmat::updateipvalmat().
| void values_correction | ( | vector & | nv, | |
| long | ipp | |||
| ) |
function checks if computed unknowns are physically reasonable
| nv | - vector of unknowns | |
| ipp | - number of integration point |
28/02/2011, TKr
Definition at line 895 of file consol_awf1.cpp.
References water_pressure_check().
Referenced by transmat::values_correction().
| void water_pressure_check | ( | double & | pw, | |
| long | ipp | |||
| ) |
function checks if water pressure is non-positive
| pw | - pore water pressure | |
| ipp | - number of integration point |
28/02/2011, TKr
Definition at line 910 of file consol_awf1.cpp.
Referenced by values_correction().
double alfa [private] |
Definition at line 70 of file consol_awf1.h.
Referenced by con_awf1mat(), get_dsw_dpw(), and read().
double alpha [private] |
Definition at line 67 of file consol_awf1.h.
Referenced by con_awf1mat(), get_alpha(), get_capww(), and read().
double delta [private] |
Definition at line 72 of file consol_awf1.h.
Referenced by con_awf1mat(), get_dsw_dpw(), get_krw(), get_sw(), and read().
double e0 [private] |
Definition at line 73 of file consol_awf1.h.
Referenced by con_awf1mat(), get_kww(), and read().
double expn [private] |
Definition at line 72 of file consol_awf1.h.
Referenced by con_awf1mat(), get_dsw_dpw(), get_krw(), get_sw(), and read().
double gamaw [private] |
Definition at line 70 of file consol_awf1.h.
Referenced by con_awf1mat(), get_dsw_dpw(), get_krw(), get_kww(), get_sw(), and read().
double hp_min [private] |
Definition at line 71 of file consol_awf1.h.
Referenced by con_awf1mat(), get_dsw_dpw(), get_krw(), and read().
double htz [private] |
Definition at line 71 of file consol_awf1.h.
Referenced by con_awf1mat(), get_dsw_dpw(), get_krw(), and read().
double k0 [private] |
Definition at line 73 of file consol_awf1.h.
Referenced by con_awf1mat(), get_kww(), and read().
double kappa [private] |
Definition at line 74 of file consol_awf1.h.
Referenced by con_awf1mat(), get_kww(), and read().
double kintr [private] |
Definition at line 69 of file consol_awf1.h.
Referenced by con_awf1mat(), get_fw1(), get_kintr(), get_kww(), and read().
double ks [private] |
Definition at line 67 of file consol_awf1.h.
Referenced by con_awf1mat(), get_capww(), get_kww(), and read().
double ksat [private] |
Definition at line 70 of file consol_awf1.h.
Referenced by con_awf1mat(), get_kww(), and read().
double kw [private] |
Definition at line 68 of file consol_awf1.h.
Referenced by con_awf1mat(), get_capww(), get_kw(), and read().
double lam [private] |
Definition at line 73 of file consol_awf1.h.
Referenced by con_awf1mat(), get_kww(), and read().
double lam0 [private] |
Definition at line 73 of file consol_awf1.h.
Referenced by con_awf1mat().
double m [private] |
Definition at line 73 of file consol_awf1.h.
Referenced by con_awf1mat(), get_kww(), read(), and rhs_volume().
waterflowtype model_type [private] |
Definition at line 66 of file consol_awf1.h.
Referenced by get_capww(), get_dsw_dpw(), get_fw1(), get_krw(), get_kww(), get_sw(), give_reqntq(), give_suction(), and read().
double muw0 [private] |
Definition at line 68 of file consol_awf1.h.
Referenced by con_awf1mat(), get_fw1(), get_kww(), get_muw(), and read().
double phi0 [private] |
Definition at line 67 of file consol_awf1.h.
Referenced by con_awf1mat(), get_capww(), get_phi(), and read().
double r [private] |
Definition at line 71 of file consol_awf1.h.
Referenced by con_awf1mat(), get_dsw_dpw(), get_krw(), and read().
double rhow [private] |
Definition at line 68 of file consol_awf1.h.
Referenced by con_awf1mat(), get_fw1(), and read().
double sigma_m_bar [private] |
Definition at line 73 of file consol_awf1.h.
Referenced by con_awf1mat().
double sigma_m_eff_bar [private] |
Definition at line 74 of file consol_awf1.h.
Referenced by con_awf1mat(), get_kww(), and read().
double sigma_m_eff_init [private] |
Definition at line 74 of file consol_awf1.h.
Referenced by con_awf1mat(), get_kww(), and read().
double sigma_m_tot_bar [private] |
Definition at line 74 of file consol_awf1.h.
Referenced by con_awf1mat(), get_kww(), and read().
double sirr [private] |
Definition at line 70 of file consol_awf1.h.
Referenced by con_awf1mat(), get_dsw_dpw(), get_sw(), and read().
double ssat [private] |
Definition at line 70 of file consol_awf1.h.
Referenced by con_awf1mat(), get_dsw_dpw(), get_sw(), and read().
1.6.2