SIFEL/METR/SRC/fcsolver.cpp File Reference

#include <math.h>
#include <string.h>
#include "fcsolver.h"
#include "globalc.h"
#include "global.h"
#include "globalt.h"
#include "globmat.h"
#include "globmatt.h"
#include "globmatc.h"
#include "gmatrix.h"
#include "gtopology.h"
#include "elemswitch.h"
#include "dloadcase.h"
#include "transprint.h"
#include "mechprint.h"
#include "backupsol.h"
#include "backupsolt.h"

Go to the source code of this file.

Functions
void	newton_raphson_coupl (long lcid)
void	newton_raphson_coupl_new (long lcid)
void	nonlinear_solver_coupl (long lcid)
void	solve_fcouplprob ()
void	vector_assemb (double *c, double *m, double *t)
void	vector_decomp (double *c, double *m, double *t)

---

Function Documentation

void newton_raphson_coupl

(

long

lcid

)

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

Parameters:

lcid

- load case id

5.4.2003

Definition at line 390 of file fcsolver.cpp.

References timecontr::actualbacktimeincr(), gmatrix::addgm(), probdescc::alpha, approximationc(), gnode::cn, compute_req_valt(), Cp, D0mat, D1mat, timecontr::endtime(), probdescc::errnr, probdescc::fcsolv, first_order_matrix(), fluxesc, fullnewtonc, lhsrhsc::give_lhs(), lhsrhsc::give_rhs(), lhsrhsc::give_tdlhs(), gmatrix::gmxv(), gtopology::gnodes, Gtm, Gtt, timecontr::initialtimeincr(), internal_gforces(), lhsrhst::lhs, lhsrhs::lhs, lhsrhsc::lhsi, lrhsc, Lsrs, Lsrsc, Lsrst, Mesprc, metr_right_hand_side(), Mp, Ndofc, timecontr::newtime(), probdescc::niilnr, nullv(), Outc, Outt, p, print_flush(), print_flusht(), print_init(), print_initt(), print_step(), print_stept(), probdescc::restype, lhsrhst::rhs, lhsrhs::rhs, gmatrix::scalgm(), solution_correction(), slesolv::solve_system(), ss(), probdescc::ssle, timecontr::starttime(), timecontr::take_values(), lhsrhst::tdlhs, probdesct::time, probdesc::time, probdescc::time, probdesc::timecon, probdescc::timecon, probdesct::timecont, Tp, probdescc::zero, and zero_order_matrix().

void newton_raphson_coupl_new

(

long

lcid

)

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

Parameters:

lcid

- load case id

24/4/2007, revised by TKr 12/9/2008

Definition at line 80 of file fcsolver.cpp.

References timecontr::actualbacktimeincr(), gmatrix::addgm(), probdescc::alpha, approximationc(), Cmu, compute_req_val(), compute_req_valt(), copy_data(), gmatrix::copygm(), Cp, D0mat, D1mat, probdesc::eigstrains, timecontr::endtime(), probdescc::errnr, probdescc::fcsolv, first_order_matrix(), fluxesc, fullnewtonc, lhsrhsc::give_lhs(), lhsrhsc::give_rhs(), lhsrhsc::give_tdlhs(), gmatrix::gmxv(), Gtu, probdesct::hdbcont, probdesc::hdbcont, init_trfel_mefel(), timecontr::initialtimeincr(), transmat::initmaterialmodels(), mechmat::initmaterialmodels(), coupmatu::initmaterialmodels(), internal_forces(), internal_gforces(), lhsrhst::lhs, lhsrhs::lhs, lhsrhsc::lhsi, lrhsc, Lsrs, Lsrsc, Lsrst, Mespr, Mesprc, Mesprt, metr_right_hand_side(), Mm, Mp, Ndofc, timecontr::newtime(), probdescc::niilnr, nullv(), Outc, p, print_close(), print_closet(), print_flush(), print_flusht(), print_init(), print_initt(), print_step(), print_stept(), hdbcontr::restore_stat(), probdescc::restype, lhsrhst::rhs, lhsrhs::rhs, gmatrix::scalgm(), solution_correction(), slesolv::solve_system(), ss(), probdescc::ssle, timecontr::starttime(), timecontr::take_values(), lhsrhst::tdlhs, probdesct::time, probdesc::time, probdescc::time, probdesc::timecon, probdescc::timecon, probdesct::timecont, Tm, Tp, updateval(), probdescc::zero, and zero_order_matrix().

Referenced by nonlinear_solver_coupl().

void nonlinear_solver_coupl

(

long

lcid

)

Definition at line 55 of file fcsolver.cpp.

References Cp, newton_raphson_coupl_new(), newtonc, and probdescc::tnlinsol.

Referenced by solve_fcouplprob().

void solve_fcouplprob

(

)

function solves fully coupled thermo-hydro-mechanical problem

JK

Definition at line 25 of file fcsolver.cpp.

References lhsrhst::lhs, lhsrhsc::lhs, lhsrhs::lhs, Lsrs, Lsrsc, Lsrst, nonlinear_solver_coupl(), lhsrhst::rhs, lhsrhsc::rhs, and lhsrhs::rhs.

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.