barelc Class Reference

#include <barelc.h>

Public Member Functions
	barelc (void)
void	eleminit (long eid)
void	giveloccoord (vector &x, vector &y, vector &lx)
void	intpointgrad (long eid)
void	intpointval (long eid)
void	lower_cap_coup_matrix (long eid, long ri, long ci, matrix &vm)
void	lower_cond_coup_matrix (long eid, long ri, long ci, matrix &vm)
void	lower_internal_fluxes (long ri, long ci, long eid, vector &ifl)
void	mainip_strains (long eid, long ri, long ci, vector &x, vector &s, vector &r)
void	res_lower_cap_coup_matrix (long eid, matrix &vm)
void	res_lower_cond_coup_matrix (long eid, matrix &vm)
void	res_lower_internal_fluxes (long eid, vector &elemif)
void	res_mainip_strains (long eid)
void	res_upper_cap_coup_matrix (long eid, matrix &vm)
void	res_upper_cond_coup_matrix (long eid, matrix &vm)
void	res_upper_cond_coup_vector (vector &f, long eid)
void	res_upper_internal_forces (long eid, vector &ifor)
void	res_volume_rhs_vector (vector &f, long eid, long lcid)
void	transf_matrix (ivector &nodes, matrix &tmat)
void	upper_cap_coup_matrix (long eid, long ri, long ci, matrix &vm)
void	upper_cond_coup_matrix (long eid, long ri, long ci, matrix &vm)
void	upper_cond_coup_vector (vector &tvm, vector &nodval, long eid, long ri, long ci)
void	upper_internal_forces (long ri, long ci, long eid, vector &ifl)
void	volume_rhs_vector (long lcid, long eid, long ri, long ci, vector &vrhs)
	~barelc (void)
Public Attributes
long *	dofe
	number of DOFs for particular medium
long **	intordvlm
	orders of integration of lower coupling matrices
long **	intordvum
	orders of integration of upper coupling matrices
long	mnb
	number of blocks in strain vector
long	mndofe
	number of degrees of freedom of mechanical part
long *	mordering
	array containing ordering of mechanical unknowns
long **	nipl
	numbers of integration points of lower coupling matrices
long **	nipu
	numbers of integration points of upper coupling matrices
long	nnemp
	number of nodes on one element in mechanical problem
long	nnetp
	number of nodes on one element in transport problem
long	ntm
	number of transported media
long	tndofe
	number of degrees of freedom of transport part
long	tnipl
long	tnipu
	total number of integration point
long	tnmcomp
	number of components of strain vector
long **	tordering
	array containing ordering of transport unknowns

Detailed Description

class barelc defines onedimensional element for mechanical-transport coupling

Definition at line 11 of file barelc.h.

Constructor & Destructor Documentation

barelc ( void )

Definition at line 17 of file barelc.cpp.

References Bar2d, Barq2d, Cp, dofe, intordvlm, intordvum, probdescc::lbb, Lbt, lin_lin, mnb, mndofe, mordering, barelq2d::nb, barel2d::nb, quadbart::ndofe, barelq2d::ndofe, linbart::ndofe, barel2d::ndofe, nipl, nipu, quadbart::nne, barelq2d::nne, linbart::nne, barel2d::nne, nnemp, nnetp, probdesct::ntm, ntm, onemedium, Qbt, quad_lin, quad_quad, probdescc::savemode, threemediacoup, probdesct::tmatt, barelq2d::tncomp, barel2d::tncomp, tndofe, tnipl, tnipu, tnmcomp, tordering, Tp, and twomediacoup.

~barelc ( void )

Definition at line 296 of file barelc.cpp.

References dofe, intordvlm, intordvum, mnb, mordering, nipl, nipu, ntm, and tordering.

Member Function Documentation

void eleminit ( long eid )

Definition at line 324 of file barelc.cpp.

References Ct, couptop::elements, intordvlm, elementc::intordvlm, intordvum, elementc::intordvum, mnb, nipl, elementc::nipl, nipu, elementc::nipu, and ntm.

void giveloccoord	(	vector &	x,
		vector &	y,
		vector &	lx
	)

Definition at line 352 of file barelc.cpp.

References vector::n, and sqr.

Referenced by lower_cap_coup_matrix(), lower_cond_coup_matrix(), lower_internal_fluxes(), res_mainip_strains(), upper_cap_coup_matrix(), upper_cond_coup_matrix(), upper_cond_coup_vector(), and upper_internal_forces().

void intpointgrad ( long eid )

function computes gradients in integration points from nodal values

Parameters:

eid

- element id

TKr, 22.3.2004

Definition at line 497 of file barelc.cpp.

References vector::a, ivector::a, allocv(), Cml, Cmu, Cp, Ct, destrv(), dofe, couptop::elements, elemvalues(), gauss_points(), transtop::give_code_numbers(), couptop::give_node_coord2d(), grad, quadbart::grad_matrix(), linbart::grad_matrix(), intordvlm, intordvum, elementc::ippl, elementc::ippu, probdescc::lbb, Lbt, lin_lin, mnb, mxv(), nnetp, probdesct::ntm, Qbt, quad_lin, quad_quad, probdescc::savemode, coupmatl::storegrad_cml(), coupmatu::storegrad_cmu(), tndofe, tordering, Tp, and Tt.

Referenced by intpointgradientsc().

void intpointval ( long eid )

function computes values at integration points from nodal values

Parameters:

eid

- element id

TKr, 22.3.3004

Definition at line 405 of file barelc.cpp.

References vector::a, ivector::a, allocv(), quadbart::approx(), linbart::approx(), intpointsc::av, Cml, Cmu, Cp, Ct, destrv(), dofe, couptop::elements, elemvalues(), gauss_points(), transtop::give_code_numbers(), intordvlm, intordvum, coupmatl::ip, coupmatu::ip, elementc::ippl, elementc::ippu, probdescc::lbb, Lbt, lin_lin, mnb, nnetp, probdesct::ntm, ntm, Qbt, quad_lin, quad_quad, probdescc::savemode, tndofe, tordering, Tp, and Tt.

Referenced by intpointvaluesc().

void lower_cap_coup_matrix	(	long	eid,
		long	ri,
		long	ci,
		matrix &	vm
	)

function computes lower coupling capacity matrix of 1D problems

Parameters:

	eid	- number of element
	ri,ci	- row and column indices of the computed block in the resulting matrix
	vm	- coupling matrix

JK, 10.4.2003, corrected by TKr 25.4.2007

Definition at line 919 of file barelc.cpp.

References Bar2d, Barq2d, bdbjac(), quadbart::bf_matrix(), linbart::bf_matrix(), Cml, Cp, Ct, barelq2d::dirvect(), couptop::elements, fillm(), gauss_points(), barelq2d::geom_matrix(), barel2d::geom_matrix(), mechcrsec::give_area(), mechtop::give_node_coord2d(), giveloccoord(), intordvlm, elementc::ippl, probdescc::lbb, Lbt, lin_lin, coupmatl::matcap(), Mc, mndofe, Mt, nnemp, nnetp, Qbt, quad_lin, quad_quad, probdescc::savemode, and tnmcomp.

Referenced by res_lower_cap_coup_matrix().

void lower_cond_coup_matrix	(	long	eid,
		long	ri,
		long	ci,
		matrix &	vm
	)

function computes lower coupling conductivity matrix of 1D problems

Parameters:

	lcid	- number of load case
	eid	- number of element
	ri,ci	- row and column indices of the computed block in the resulting matrix
	vm	- coupling matrix

JK, 10.4.2003, corrected by TKr 25.4.2007

Definition at line 789 of file barelc.cpp.

References Bar2d, Barq2d, bdbjac(), quadbart::bf_matrix(), linbart::bf_matrix(), Cml, Cp, Ct, barelq2d::dirvect(), couptop::elements, fillm(), gauss_points(), barelq2d::geom_matrix(), barel2d::geom_matrix(), mechcrsec::give_area(), mechtop::give_node_coord2d(), giveloccoord(), intordvlm, elementc::ippl, probdescc::lbb, Lbt, lin_lin, coupmatl::matcond(), Mc, mndofe, Mt, nnemp, nnetp, Qbt, quad_lin, quad_quad, probdescc::savemode, and tnmcomp.

Referenced by res_lower_cond_coup_matrix().

void lower_internal_fluxes	(	long	ri,
		long	ci,
		long	eid,
		vector &	ifl
	)

function computes internal fluxes

Parameters:

	ri,ci	- row and column indices
	eid	- element id
	ifl	- vector of internal fluxes

JK, 17.7.2005, corrected by TKr 25.4.2007

Definition at line 1347 of file barelc.cpp.

References addv(), linbart::approx(), Cml, cmulv(), coupmatl::computenlfluxes(), Cp, Ct, couptop::elements, fillv(), gauss_points(), transcrsec::give_area(), transtop::give_elemnodes(), transtop::give_node_coord2d(), coupmatl::givefluxes_cml(), giveloccoord(), quadbart::grad_matrix(), linbart::grad_matrix(), intordvlm, elementc::ippl, probdescc::lbb, Lbt, lin_lin, mndofe, mtxv(), nnetp, nodes, Qbt, quad_lin, quad_quad, probdescc::savemode, Tc, tnmcomp, and Tt.

Referenced by res_lower_internal_fluxes().

void mainip_strains	(	long	eid,
		long	ri,
		long	ci,
		vector &	x,
		vector &	s,
		vector &	r
	)

function computes strains in integration points of element

Parameters:

	lcid	- load case id
	eid	- element id
	ri	- row index
	ci	- column index

TKr, 22.3.2004

Definition at line 596 of file barelc.cpp.

References vector::a, allocm(), allocv(), Bar2d, Barq2d, Cml, Cmu, Cp, Ct, destrm(), destrv(), couptop::elements, gauss_points(), barelq2d::geom_matrix(), barel2d::geom_matrix(), intordvlm, intordvum, elementc::ippl, elementc::ippu, probdescc::lbb, lin_lin, mnb, mndofe, mxv(), quad_lin, quad_quad, probdescc::savemode, coupmatl::storestrain_cml(), and coupmatu::storestrain_cmu().

Referenced by res_mainip_strains().

void res_lower_cap_coup_matrix	(	long	eid,
		matrix &	vm
	)

function assembles lower coupling capacity matrices into one element matrix

Parameters:

	eid	-element id
	vm	- element lower coupling capacity matrix

JK, 17.7.2005

Definition at line 1089 of file barelc.cpp.

References quadbart::codnum(), linbart::codnum(), Cp, fillm(), probdescc::lbb, Lbt, lin_lin, lower_cap_coup_matrix(), mat_localize(), mnb, mndofe, mordering, nnetp, ntm, Qbt, and quad_lin.

Referenced by lower_cap_coupl_mat().

void res_lower_cond_coup_matrix	(	long	eid,
		matrix &	vm
	)

function assembles lower coupling conductivity matrices into one element matrix

Parameters:

	eid	- element id
	vm	- element lower coupling conductivity matrix

JK, 17.7.2005

Definition at line 1019 of file barelc.cpp.

References quadbart::codnum(), linbart::codnum(), Cp, fillm(), probdescc::lbb, Lbt, lin_lin, lower_cond_coup_matrix(), mat_localize(), mnb, mndofe, mordering, nnetp, ntm, Qbt, and quad_lin.

Referenced by lower_cond_coupl_mat().

void res_lower_internal_fluxes	(	long	eid,
		vector &	ifl
	)

function computes internal fluxes caused by mechanical effects

Parameters:

	eid	- element id
	ifl	- vector of internal fluxes

JK, 17.7.2005

Definition at line 1418 of file barelc.cpp.

References vector::a, allocv(), quadbart::codnum(), destrv(), fillv(), locglob(), lower_internal_fluxes(), ntm, Qbt, and tnmcomp.

Referenced by internal_coup_fluxes().

void res_mainip_strains ( long eid )

function computes resulting strains in integration points of element

Parameters:

eid

- element id

TKr, 22.3.2004

Definition at line 677 of file barelc.cpp.

References allocm(), allocv(), Barq2d, copyv(), Cp, Ct, destrm(), destrv(), barelq2d::dirvect(), eldispl(), mechtop::give_code_numbers(), couptop::give_elemnodes(), couptop::give_node_coord2d(), giveloccoord(), probdescc::lbb, lgvectortransf(), mechtop::locsystems(), mainip_strains(), mndofe, Mt, nnemp, nodes, probdesct::ntm, quad_lin, quad_quad, Tp, and transf_matrix().

Referenced by intpointstrainsc().

void res_upper_cap_coup_matrix	(	long	eid,
		matrix &	vm
	)

function assembles upper coupling capacity matrices into one element matrix

Parameters:

	eid	-element id
	vm	- element upper coupling capacity matrix

JK, 17.7.2005

Definition at line 1054 of file barelc.cpp.

References quadbart::codnum(), linbart::codnum(), Cp, fillm(), probdescc::lbb, Lbt, lin_lin, mat_localize(), mnb, mndofe, mordering, nnetp, ntm, Qbt, quad_lin, and upper_cap_coup_matrix().

Referenced by upper_cap_coupl_mat().

void res_upper_cond_coup_matrix	(	long	eid,
		matrix &	vm
	)

function assembles upper coupling conductivity matrices into one element matrix

Parameters:

	eid	- element id
	vm	- element upper coupling conductivity matrix

JK, 17.7.2005

Definition at line 984 of file barelc.cpp.

References quadbart::codnum(), linbart::codnum(), Cp, fillm(), probdescc::lbb, Lbt, lin_lin, mat_localize(), mnb, mndofe, mordering, nnetp, ntm, Qbt, quad_lin, and upper_cond_coup_matrix().

Referenced by upper_cond_coupl_mat().

void res_upper_cond_coup_vector	(	vector &	f,
		long	eid
	)

function computes

Parameters:

	f	- vector of right-hand side
	eid	- element id

JK, 15.7.2005

Definition at line 1200 of file barelc.cpp.

References vector::a, allocv(), quadbart::codnum(), linbart::codnum(), Cp, destrv(), dofe, fillv(), transtop::give_code_numbers(), globloc(), initialvalues(), probdescc::lbb, Lbt, lin_lin, ntm, Qbt, quad_lin, tndofe, Tt, and upper_cond_coup_vector().

Referenced by assemble_coup().

void res_upper_internal_forces	(	long	eid,
		vector &	ifor
	)

function computes internal forces caused by transport effects

Parameters:

	eid	- element id
	ifor	- vector of internal forces

JK, 17.7.2005

Definition at line 1324 of file barelc.cpp.

References addv(), fillv(), mnb, mndofe, and upper_internal_forces().

Referenced by internal_coup_forces().

void res_volume_rhs_vector	(	vector &	f,
		long	eid,
		long	lcid
	)

function assembles resulting element volume right-hand side

Parameters:

	eid	- element id
	lcid	- load case id
	f	- resulting volume right-hand side vector of one element

TKr, 09/03/2011

Definition at line 1525 of file barelc.cpp.

References vector::a, allocv(), quadbart::codnum(), linbart::codnum(), Cp, destrv(), probdescc::lbb, Lbt, lin_lin, locglob(), mnb, mndofe, Qbt, quad_lin, and volume_rhs_vector().

Referenced by volume_rhs_vectorc().

void transf_matrix	(	ivector &	nodes,
		matrix &	tmat
	)

transformation matrix x_g = T x_l

Definition at line 377 of file barelc.cpp.

References node::e1, node::e2, fillm(), matrix::m, Mt, ivector::n, mechtop::nodes, and node::transf.

Referenced by res_mainip_strains().

void upper_cap_coup_matrix	(	long	eid,
		long	ri,
		long	ci,
		matrix &	vm
	)

function computes upper coupling capacity matrix of 1D problems

Parameters:

	eid	- number of element
	ri,ci	- row and column indices of the computed block in the resulting matrix
	vm	- coupling matrix

JK, 9.1.2003, corrected by TKr 25.4.2007

Definition at line 854 of file barelc.cpp.

References Bar2d, Barq2d, bdbjac(), quadbart::bf_matrix(), linbart::bf_matrix(), Cmu, Cp, Ct, barelq2d::dirvect(), couptop::elements, fillm(), gauss_points(), barelq2d::geom_matrix(), barel2d::geom_matrix(), mechcrsec::give_area(), mechtop::give_node_coord2d(), giveloccoord(), intordvum, elementc::ippu, probdescc::lbb, Lbt, lin_lin, coupmatu::matcap(), Mc, mndofe, Mt, nnemp, nnetp, Qbt, quad_lin, quad_quad, probdescc::savemode, and tnmcomp.

Referenced by res_upper_cap_coup_matrix().

void upper_cond_coup_matrix	(	long	eid,
		long	ri,
		long	ci,
		matrix &	vm
	)

function computes upper coupling conductivity matrix of 1D problems

Parameters:

	lcid	- number of load case
	eid	- number of element
	ri,ci	- row and column indices of the computed block in the resulting matrix
	vm	- coupling matrix

JK, 10.4.2003, corrected by TKr 25.4.2007

Definition at line 723 of file barelc.cpp.

References Bar2d, Barq2d, bdbjac(), quadbart::bf_matrix(), linbart::bf_matrix(), Cmu, Cp, Ct, barelq2d::dirvect(), couptop::elements, fillm(), gauss_points(), barelq2d::geom_matrix(), barel2d::geom_matrix(), mechcrsec::give_area(), mechtop::give_node_coord2d(), giveloccoord(), intordvum, elementc::ippu, probdescc::lbb, Lbt, lin_lin, coupmatu::matcond(), Mc, mndofe, Mt, nnemp, nnetp, Qbt, quad_lin, quad_quad, probdescc::savemode, and tnmcomp.

Referenced by res_upper_cond_coup_matrix().

void upper_cond_coup_vector	(	vector &	tvm,
		vector &	nodval,
		long	eid,
		long	ri,
		long	ci
	)

function computes coupling vector of 1D problems

Parameters:

	vm	- coupling vector
	nodval	- vector of nodal values
	eid	- number of element
	ri,ci	- row and column indices of the computed block in the resulting matrix