quadrilatc Class Reference

#include <quadrilatc.h>

Public Member Functions
void	dmatblockcol (long ri, long ci, matrix &d, matrix &dd)
void	dmatblockrow (long ri, long ci, matrix &d, matrix &dd)
void	eleminit (long eid)
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	mainip_strains (long eid, long ri, long ci, vector &x, vector &y, vector &r)
	quadrilatc (void)
void	res_lower_cap_coup_matrix (long eid, matrix &vm)
void	res_lower_cond_coup_matrix (long eid, matrix &vm)
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	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)
	~quadrilatc (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 *	mncomp
	array of numbers of components in strain vectors
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 quadrilatc defines quadrilateral element for problems of mechanical-transport coupling

Definition at line 13 of file quadrilatc.h.

Constructor & Destructor Documentation

quadrilatc ( void )

Definition at line 14 of file quadrilatc.cpp.

References Cp, dofe, intordvlm, intordvum, probdescc::lbb, lin_lin, Lqt, mnb, mncomp, mndofe, mordering, planeelemqq::nb, planeelemlq::nb, planeelemqq::ncomp, planeelemlq::ncomp, quadquadrilatt::ndofe, planeelemqq::ndofe, quadlineart::ndofe, planeelemlq::ndofe, nipl, nipu, quadquadrilatt::nne, planeelemqq::nne, quadlineart::nne, planeelemlq::nne, nnemp, nnetp, probdesct::ntm, ntm, onemedium, Pelq, Peqq, Qqt, quad_lin, quad_quad, probdescc::savemode, threemediacoup, probdesct::tmatt, planeelemqq::tncomp, planeelemlq::tncomp, tndofe, tnipl, tnipu, tnmcomp, tordering, Tp, and twomediacoup.

~quadrilatc ( void )

Definition at line 278 of file quadrilatc.cpp.

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

Member Function Documentation

void dmatblockcol	(	long	ri,
		long	ci,
		matrix &	d,
		matrix &	dd
	)

function assembles blocks of stiffness matrix of material

Parameters:

	ri	- row index
	ci	- column index
	d	- stiffness matrix of material
	dd	- required block of stiffness matrix of material

Definition at line 334 of file quadrilatc.cpp.

References fillm(), and mnb.

Referenced by upper_cap_coup_matrix(), upper_cond_coup_matrix(), and upper_cond_coup_vector().

void dmatblockrow	(	long	ri,
		long	ci,
		matrix &	d,
		matrix &	dd
	)

function assembles blocks of stiffness matrix of material

Parameters:

	ri	- row index
	ci	- column index
	d	- stiffness matrix of material
	dd	- required block of stiffness matrix of material

Definition at line 362 of file quadrilatc.cpp.

References fillm(), and mnb.

Referenced by lower_cap_coup_matrix(), and lower_cond_coup_matrix().

void eleminit ( long eid )

Definition at line 297 of file quadrilatc.cpp.

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

void intpointgrad ( long eid )

function computes gradients in integration points from nodal values

Parameters:

eid

- element id

TKr, 01/10/2012

Definition at line 489 of file quadrilatc.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, quadquadrilatt::grad_matrix(), quadlineart::grad_matrix(), intordvlm, intordvum, elementc::ippl, elementc::ippu, probdescc::lbb, lin_lin, Lqt, mnb, mxv(), nnetp, probdesct::ntm, Qqt, quad_lin, quad_quad, probdescc::savemode, coupmatl::storegrad_cml(), coupmatu::storegrad_cmu(), tndofe, tordering, Tp, and Tt.

void intpointval ( long eid )

function computes values at integration points from nodal values

Parameters:

eid

- element id

TKr, 01/10/2012

Definition at line 393 of file quadrilatc.cpp.

References vector::a, ivector::a, allocv(), quadquadrilatt::approx(), quadlineart::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, lin_lin, Lqt, mnb, nnetp, probdesct::ntm, ntm, Qqt, quad_lin, quad_quad, probdescc::savemode, tndofe, tordering, Tp, and Tt.

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

function computes lower coupling capacity matrix of 2D 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

Definition at line 951 of file quadrilatc.cpp.

References quadquadrilatt::approx(), quadlineart::approx(), bdbjac(), quadquadrilatt::bf_matrix(), quadlineart::bf_matrix(), Cml, Cp, Ct, dmatblockrow(), dofe, couptop::elements, fillm(), gauss_points(), planeelemqq::geom_matrix(), planeelemlq::geom_matrix_block(), transtop::give_elemnodes(), mechtop::give_node_coord2d(), transcrsec::give_thickness(), intordvlm, elementc::ippl, probdescc::lbb, lin_lin, Lqt, coupmatl::matcap(), mncomp, mndofe, Mt, nnemp, nnetp, nodes, Pelq, Peqq, Qqt, quad_lin, quad_quad, probdescc::savemode, Tc, tnmcomp, and Tt.

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 2D 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

Definition at line 803 of file quadrilatc.cpp.

References quadquadrilatt::approx(), quadlineart::approx(), bdbjac(), quadquadrilatt::bf_matrix(), quadlineart::bf_matrix(), Cml, Cp, Ct, dmatblockrow(), dofe, couptop::elements, fillm(), gauss_points(), planeelemqq::geom_matrix(), planeelemlq::geom_matrix_block(), transtop::give_elemnodes(), mechtop::give_node_coord2d(), transcrsec::give_thickness(), intordvlm, elementc::ippl, probdescc::lbb, lin_lin, Lqt, coupmatl::matcond(), mncomp, mndofe, Mt, nnemp, nnetp, nodes, Pelq, Peqq, Qqt, quad_lin, quad_quad, probdescc::savemode, Tc, tnmcomp, and Tt.

Referenced by res_lower_cond_coup_matrix().

void mainip_strains	(	long	eid,
		long	ri,
		long	ci,
		vector &	x,
		vector &	y,
		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, 01/10/2012

Definition at line 593 of file quadrilatc.cpp.

References vector::a, allocm(), allocv(), Cml, Cmu, Cp, Ct, destrm(), destrv(), couptop::elements, gauss_points(), planeelemqq::geom_matrix(), planeelemlq::geom_matrix(), planeelemlq::geom_matrix_block(), intordvlm, intordvum, elementc::ippl, elementc::ippu, probdescc::lbb, lin_lin, mnb, mncomp, mndofe, mxv(), Pelq, Peqq, 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 1135 of file quadrilatc.cpp.

References allocm(), quadquadrilatt::codnum(), quadlineart::codnum(), Cp, destrm(), dofe, fillm(), probdescc::lbb, lin_lin, lower_cap_coup_matrix(), Lqt, mat_localize(), mnb, mndofe, mordering, nnetp, ntm, Qqt, 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 1064 of file quadrilatc.cpp.

References allocm(), quadquadrilatt::codnum(), quadlineart::codnum(), Cp, destrm(), dofe, fillm(), probdescc::lbb, lin_lin, lower_cond_coup_matrix(), Lqt, mat_localize(), mnb, mndofe, mordering, nnetp, ntm, Qqt, and quad_lin.

Referenced by lower_cond_coupl_mat().

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 680 of file quadrilatc.cpp.

References allocm(), allocv(), copyv(), Cp, Ct, destrm(), destrv(), eldispl(), mechtop::give_code_numbers(), couptop::give_elemnodes(), couptop::give_node_coord2d(), probdescc::lbb, lgvectortransf(), lin_lin, mechtop::locsystems(), mainip_strains(), mndofe, Mt, nnemp, nodes, probdesct::ntm, Pelq, Peqq, quad_lin, quad_quad, Tp, planeelemqq::transf_matrix(), and planeelemlq::transf_matrix().

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 1099 of file quadrilatc.cpp.

References allocm(), quadquadrilatt::codnum(), quadlineart::codnum(), Cp, destrm(), dofe, fillm(), probdescc::lbb, lin_lin, Lqt, mat_localize(), mnb, mndofe, mordering, nnetp, ntm, Qqt, 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 1026 of file quadrilatc.cpp.

References allocm(), quadquadrilatt::codnum(), quadlineart::codnum(), Cp, destrm(), dofe, fillm(), probdescc::lbb, lin_lin, Lqt, mat_localize(), mnb, mndofe, mordering, ntm, Qqt, quad_lin, quad_quad, 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	-
	eid	- element id

JK, 19.7.2005

Definition at line 1257 of file quadrilatc.cpp.

References vector::a, allocv(), quadquadrilatt::codnum(), quadlineart::codnum(), Cp, destrv(), dofe, fillv(), transtop::give_code_numbers(), globloc(), initialvalues(), probdescc::lbb, lin_lin, Lqt, mnb, ntm, Qqt, quad_lin, tndofe, Tt, and upper_cond_coup_vector().

Referenced by assemble_coup().

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

function computes upper coupling capacity matrix of 2D 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

Definition at line 876 of file quadrilatc.cpp.

References quadquadrilatt::approx(), quadlineart::approx(), bdbjac(), quadquadrilatt::bf_matrix(), quadlineart::bf_matrix(), Cmu, Cp, Ct, dmatblockcol(), dofe, couptop::elements, fillm(), gauss_points(), planeelemqq::geom_matrix(), planeelemlq::geom_matrix_block(), transtop::give_elemnodes(), mechtop::give_node_coord2d(), transcrsec::give_thickness(), intordvum, elementc::ippu, probdescc::lbb, lin_lin, Lqt, coupmatu::matcap(), mncomp, mndofe, Mt, nnemp, nnetp, nodes, Pelq, Peqq, Qqt, quad_lin, quad_quad, probdescc::savemode, Tc, tnmcomp, and Tt.

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 2D 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

Definition at line 729 of file quadrilatc.cpp.

References quadquadrilatt::approx(), quadlineart::approx(), bdbjac(), quadquadrilatt::bf_matrix(), quadlineart::bf_matrix(), Cmu, Cp, Ct, dmatblockcol(), dofe, couptop::elements, fillm(), gauss_points(), planeelemqq::geom_matrix(), planeelemlq::geom_matrix_block(), transtop::give_elemnodes(), mechtop::give_node_coord2d(), transcrsec::give_thickness(), intordvum, elementc::ippu, probdescc::lbb, lin_lin, Lqt, coupmatu::matcond(), mncomp, mndofe, Mt, nnemp, nnetp, nodes, Pelq, Peqq, Qqt, quad_lin, quad_quad, probdescc::savemode, Tc, tnmcomp, and Tt.

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 2D 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