SIFEL/METR/SRC/pcsolver.h File Reference

#include <stdio.h>

Go to the source code of this file.

Functions

void newton_raphson_parcoupl (long lcid)
void newton_raphson_parcoupl_common_dt (long lcid)
void newton_raphson_parcoupl_comp (long lcid)
void newton_raphson_parcoupl_lin (long lcid)
void newton_raphson_parcoupl_nonlin (long lcid)
void newton_raphson_parcoupl_nonlin_new (long lcid)
void newton_raphson_parcoupl_nonlin_old (long lcid)
void solve_pcouplprob ()
void vector_assemb (double *c, double *m, double *t)
void vector_decomp (double *c, double *m, double *t)

Function Documentation

void newton_raphson_parcoupl ( long  lcid  ) 

function solves partially coupled thermo-hydro-mechanical time-dependent problem by the Newton-Raphson method

Parameters:
lcid - load case id

11.6.2008

Definition at line 849 of file pcsolver.cpp.

References newton_raphson_parcoupl_lin(), newton_raphson_parcoupl_nonlin(), nonlinear_nonstationary_problem, nonstationary_problem, print_err(), Tp, and probdesct::tprob.

void newton_raphson_parcoupl_common_dt ( long  lcid  ) 

function solves partially coupled thermo-hydro-mechanical time-dependent problem by the Newton-Raphson method

Parameters:
lcid - load case id

TKo, 27/08/2014

Definition at line 52 of file pcsolver.cpp.

References approximation(), clean_yes, probdescc::cleanmatrix, Cmat, compute_req_val(), compute_req_valt(), copy_data(), copyv(), Cp, discont_nonlin_nonstat_problem, discont_nonstat_problem, timecontr::dtmin, timecontr::endtime(), mt_glob_vec::f, np_glob_vec::f, mt_glob_vec::fb, mt_glob_vec::fl, mt_glob_vec::fp, probdesc::hdbcont, probdesct::hdbcont, init_trfel_mefel(), timecontr::initialtimeincr(), timecontr::isitimptime(), probdesc::istep, probdesct::istep, mt_glob_vec::istep, np_glob_vec::istep, Kmat, np_glob_vec::lhs, np_glob_vec::lhsb, mefel_trfel(), Mespr, Mesprt, Mm, Mp, Ndofm, Ndoft, timecontr::newtime(), nonlinear_nonstationary_problem, nonstat_solver_dform_init(), nonstat_solver_init(), nonstationary_problem, nullv(), timecontr::oldtime(), one_step(), one_step_linear(), one_step_nonlinear(), one_step_nonlinear_dform(), print_close(), print_closet(), print_err(), print_flush(), print_flusht(), print_step(), print_step_forced(), print_stept(), print_stept_forced(), mt_glob_vec::r, np_glob_vec::rhs, hdbcontr::save_stat(), Smat, solver_save(), solvert_save(), timecontr::starttime(), timecontr::tct, np_glob_vec::tdlhs, np_glob_vec::tdlhsb, timecontr::time, probdesc::time, probdesct::time, probdescc::time, probdesc::timecon, probdescc::timecon, probdesct::timecont, Tp, probdesct::tprob, trfel_mefel(), mechmat::updateipval(), and visco_solver_init().

Referenced by solve_pcouplprob().

void newton_raphson_parcoupl_comp ( long  lcid  ) 
void newton_raphson_parcoupl_lin ( long  lcid  ) 

Function solves system of linear TRFEL algebraic and non-linear MEFEL equations by Newton-Raphson method for time-dependent problems

Parameters:
lcid - load case id

15.7.2005 JK, corrected by TKr 5.4.2007 12/06/2012 TKr updated

Definition at line 879 of file pcsolver.cpp.

References actual_previous_nodval(), timecontr::actualforwtimeincr(), gmatrix::addgm(), probdesct::alpha, approximation(), capacity_matrix(), Cmat, mechbclc::comp_sum(), compute_req_val(), compute_req_valt(), conductivity_matrix(), copy_data(), Cp, nonlinman::divc_err, nonlinman::divc_step, timecontr::dtmin, probdesc::eigstrains, timecontr::endtime(), probdesc::errnr, f, probdescc::fcsolv, fillm(), fillv(), fullnewtonc, lhsrhs::give_lhs(), lhsrhst::give_lhs(), lhsrhs::give_rhs(), lhsrhst::give_rhs(), lhsrhst::give_tdlhs(), gmatrix::gmxv(), Gtm, Gtt, probdesct::hdbcont, probdesc::hdbcont, hdbcontr::hdbtype, incr_internal_forces(), init_trfel_mefel(), timecontr::initialtimeincr(), transmat::initmaterialmodels(), mechmat::initmaterialmodels(), internal_forces(), timecontr::isitimptime(), Kmat, lsm_quad(), Lsrs, Lsrst, Mb, mefel_right_hand_side(), Mespr, Mesprc, Mesprt, Mm, Mp, Ndofm, Ndoft, timecontr::newtime(), probdesc::niilnr, probdesc::nlman, nullv(), timecontr::oldtime(), probdesc::otherstate, Out, Outt, p, print_close(), print_closet(), print_flush(), print_flusht(), print_init(), print_initt(), print_step(), print_step_forced(), print_stept(), hdbcontr::restore_stat(), gmatrix::scalgm(), Smat, solution_correction(), slesolv::solve_system(), solver_restore(), solver_save(), solvert_restore(), solvert_save(), ss(), probdesc::ssle, probdesct::ssle, timecontr::starttime(), stiffness_matrix(), nonlinman::stmat, probdesc::strainstate, probdesc::stressstate, mechbclc::sumcomp, timecontr::take_values(), tangent_stiff, probdesct::time, probdesc::time, probdescc::time, probdesc::timecon, probdescc::timecon, probdesct::timecont, Tm, Tp, trfel_right_hand_side(), mechmat::updateipval(), transmat::updateipval(), probdesc::zero, and probdesct::zero.

Referenced by newton_raphson_parcoupl().

void newton_raphson_parcoupl_nonlin ( long  lcid  ) 

Function solves system of non-linear TRFEL algebraic and non-linear MEFEL equations by Newton-Raphson method for time-dependent problems.

Parameters:
lcid - load case id

TKr 5.4.2007, pozn.: 31/08/2012 TKr stare aktualizovat??!!!

Definition at line 1401 of file pcsolver.cpp.

References actual_previous_nodval(), timecontr::actualforwtimeincr(), gmatrix::addgm(), probdesct::alpha, approximation(), approximation_humid(), approximation_inittemper(), approximation_temper(), capacity_matrix(), Cmat, mechbclc::comp_sum(), compute_req_val(), compute_req_valt(), conductivity_matrix(), probdesct::convergcontrolt, Cp, timecontr::dtmin, timecontr::endtime(), probdesct::err, probdesc::errnr, f, probdescc::fcsolv, fillm(), fillv(), fluxest, fullnewtonc, fullnewtont, lhsrhs::give_lhs(), lhsrhst::give_lhs(), lhsrhs::give_rhs(), lhsrhst::give_rhs(), lhsrhst::give_tdlhs(), gmatrix::gmxv(), Gtm, Gtt, probdesct::hdbcont, probdesc::hdbcont, hdbcontr::hdbtype, incr_internal_forces(), timecontr::initialtimeincr(), transmat::initmaterialmodels(), mechmat::initmaterialmodels(), internal_fluxes(), internal_forces(), timecontr::isitimptime(), Kmat, lrhst, lsm_quad(), Lsrs, Lsrst, Mb, mefel_right_hand_side(), Mespr, Mesprc, Mesprt, Mm, Mp, Ndofm, Ndoft, timecontr::newtime(), probdesct::nii, probdesc::niilnr, probdesc::nlman, nullv(), timecontr::oldtime(), probdesc::otherstate, Out, Outt, p, print_close(), print_closet(), print_flush(), print_flusht(), print_init(), print_initt(), print_step(), print_step_forced(), print_stept(), hdbcontr::restore_stat(), hdbcontr::save_stat(), gmatrix::scalgm(), Smat, solution_correction(), slesolv::solve_system(), solver_restore(), solver_save(), solvert_restore(), solvert_save(), ss(), probdesc::ssle, probdesct::ssle, timecontr::starttime(), stiffness_matrix(), nonlinman::stmat, probdesc::strainstate, probdesc::stressstate, mechbclc::sumcomp, timecontr::take_values(), tangent_stiff, probdesct::time, probdesc::time, probdescc::time, probdesc::timecon, probdescc::timecon, probdesct::timecont, Tm, Tp, probdesct::trestype, trfel_right_hand_side(), probdesct::trsolv, mechmat::updateipval(), transmat::updateipval(), yes, probdesc::zero, and probdesct::zero.

Referenced by newton_raphson_parcoupl().

void newton_raphson_parcoupl_nonlin_new ( long  lcid  ) 

function solves system of nonlinear algebraic equations by Newton-Raphson method

Parameters:
lcid - load case id

19.5.2004 - modified by TKr 31.10.2005

Definition at line 2348 of file pcsolver.cpp.

References timecontr::actualbacktimeincr(), gmatrix::addgm(), probdesct::alpha, approximation(), approximation_humid(), approximation_temper(), capacity_matrix(), Cmat, mechbclc::comp_sum(), compute_req_val(), conductivity_matrix(), Cp, timecontr::endtime(), probdesct::err, f, probdesc::filename, filename_decomposition(), fluxest, fullnewtont, gmatrix::gmxv(), Gtm, Gtt, timecontr::initialtimeincr(), mechmat::initmaterialmodels(), internal_fluxes(), internal_forces(), timecontr::isitimptime(), Kmat, lhsrhs::lhs, lhsrhst::lhs, lrhst, Lsrs, Lsrst, Mb, mefel_right_hand_side(), Mesprc, Mesprt, Mm, Mp, Ndofm, Ndoft, timecontr::newtime(), probdesct::nii, nullv(), Out, Outdm, outdriverm::outfn, Outt, p, probdesc::path, print_close(), print_closet(), print_flush(), print_flusht(), print_init(), print_initt(), print_step(), print_stept(), lhsrhs::rhs, lhsrhst::rhs, gmatrix::scalgm(), Smat, solution_correction(), slesolv::solve_system(), solver_save(), ss(), probdescc::ssle, timecontr::starttime(), stiffness_matrix(), probdesc::suffix, timecontr::take_values(), lhsrhst::tdlhs, probdesct::time, probdesc::time, probdescc::time, probdesc::timecon, probdescc::timecon, probdesct::timecont, Tp, probdesct::trestype, trfel_right_hand_side(), probdesct::trsolv, mechmat::updateipval(), and probdesct::zero.

Referenced by par_newton_raphson_parcoupl_lin().

void newton_raphson_parcoupl_nonlin_old ( long  lcid  ) 
void solve_pcouplprob (  ) 

partially coupled problem transport problem influences mechanical problem

function solves partially coupled thermo-hydro-mechanical time-dependent problem transport problems are linear or nonlinear mechanical analysis is always nonlinear

Definition at line 26 of file pcsolver.cpp.

References Cp, newton_raphson_parcoupl_common_dt(), newtonc, print_err(), and probdescc::tnlinsol.

Referenced by solve_metr_problem().

void vector_assemb ( double *  c,
double *  m,
double *  t 
)

Definition at line 684 of file fcsolver.cpp.

References gtopology::cngtopcorr, Gtm, Gtt, and gtopology::ndof.

void vector_decomp ( double *  c,
double *  m,
double *  t 
)

Definition at line 699 of file fcsolver.cpp.

References gtopology::cngtopcorr, Gtm, Gtt, and gtopology::ndof.


Generated by  doxygen 1.6.2