trlinaxisym Class Reference

#include <trlinaxisym.h>

Public Member Functions
double	approx (vector &areacoord, vector &nodval)
double	approx_nat (double xi, double eta, vector &nodval)
void	bf_matrix (matrix &n, vector &areacoord)
void	boundary_flux (vector &v, long lcid, long eid, long leid, long ri, long ci)
void	capacity_matrix (long eid, long ri, long ci, matrix &cm)
void	codnum (long *cn, long ri)
void	conductivity_matrix (long lcid, long eid, long ri, long ci, matrix &km)
void	convection_vector (vector &v, long lcid, long eid, long leid, long ri, long ci)
void	edge_integral (long edg, vector &x, vector &y, long intord, vector &gp, vector &w, vector &coef, matrix &km)
void	edgenodeval (long edg, vector &nodval, vector &list)
void	give_approx_fun (double &f, vector &areacoord, long i)
void	grad_matrix (matrix &gm, vector &b, vector &c)
void	internal_fluxes (long lcid, long eid, vector &ifl)
void	intpointgrad (long eid)
void	intpointother (long eid)
void	intpointval (long eid)
void	ipcoordblock (long eid, long ri, long ci, double **coord)
void	nod_others (long lcid, long eid, long ri, long ci)
void	quantity_source_vector (vector &sv, vector &nodval, long eid, long ri, long ci)
void	res_boundary_flux (vector &f, long lcid, long eid, long leid)
void	res_capacity_matrix (long eid, matrix &cm)
void	res_conductivity_matrix (long eid, long lcid, matrix &km)
void	res_convection_vector (vector &f, long lcid, long eid, long leid)
void	res_internal_fluxes (long eid, vector &elemif)
void	res_quantity_source_vector (vector &sv, vector &nodval, long lcid, long eid)
void	res_transmission_vector (vector &f, long lcid, long eid, long leid)
double	total_integral (long eid, vector &nodval)
void	transf_coeff (long edg, vector &coeff, vector &list, long eid, long ri, long ci, long ipp, bocontypet *bc)
void	transf_flux (long edg, vector &coeff, vector &list, vector &trc, vector &trr, long eid, long ri, long ci, long ipp, bocontypet *bc)
void	transf_val (long edg, vector &nodval, vector &list, vector &trc, vector &trr, long eid, long ri, long ci, long ipp, bocontypet *bc)
void	transmission_matrix (long lcid, long eid, long ri, long ci, matrix &km)
void	transmission_vector (vector &v, long lcid, long eid, long leid, long ri, long ci)
	trlinaxisym (void)
	~trlinaxisym (void)
Public Attributes
long **	dofe
	number of DOFs on the element
long	intordb
	order of integration on boundaries
long **	intordcm
	order of integration of capacity matrix
long **	intordkm
	order of integration of conductivity matrix
long	napfun
	number of approximated functions on the element
long	ncomp
	number of components of stress and strain tensors
long	ndofe
	total number of DOFs on the element
long	ned
	number of edges
long **	nip
	number of integration points for conductivity matrix
long	nne
	number of nodes on one element
long	nned
	number of nodes on one edge
long	ntm
	number of transported matter
long **	ordering

Detailed Description

class trlinaxisym defines triangular element with linear approximation functions for transport problems for axisymmetric solution

Definition at line 11 of file trlinaxisym.h.

Constructor & Destructor Documentation

trlinaxisym ( void )

Definition at line 12 of file trlinaxisym.cpp.

References dofe, intordb, intordcm, intordkm, napfun, ncomp, ndofe, ned, nip, nne, nned, nomedium, probdesct::ntm, ntm, onemedium, ordering, print_err(), probdesct::savemode, threemediacoup, probdesct::tmatt, Tp, and twomediacoup.

~trlinaxisym ( void )

Definition at line 105 of file trlinaxisym.cpp.

References dofe, intordcm, intordkm, nip, ntm, and ordering.

Member Function Documentation

double approx	(	vector &	areacoord,
		vector &	nodval
	)

function approximates function defined by nodal values

Parameters:

	areacoord	- vector containing area coordinates
	nodval	- nodal values

JK, 25.9.2001

Definition at line 163 of file trlinaxisym.cpp.

References f, and scprd().

Referenced by capacity_matrix(), conductivity_matrix(), edge_integral(), internal_fluxes(), intpointother(), intpointval(), quantity_source_vector(), and total_integral().

double approx_nat	(	double	xi,
		double	eta,
		vector &	nodval
	)

function approximates function defined by nodal values

Parameters:

	xi,eta	- natural coordinates
	nodval	- nodal values

Definition at line 177 of file trlinaxisym.cpp.

References f, and scprd().

Referenced by ipcoordblock().

void bf_matrix	(	matrix &	n,
		vector &	areacoord
	)

function returns matrix of base function

Parameters:

	n	- array of approximation functions
	areacoord	- area coordinates

JK, 5.11.2001

Definition at line 421 of file trlinaxisym.cpp.

References fillm().

Referenced by capacity_matrix(), conductivity_matrix(), edge_integral(), and quantity_source_vector().

void boundary_flux	(	vector &	v,
		long	lcid,
		long	eid,
		long	leid,
		long	ri,
		long	ci
	)

function computes contributions to the boundary flux from transmission vector

{Gamma_3} N^T c_{tr} N dGamma * nodal_external_value

Parameters:

	v	- vector of boundary fluxes
	lcid	- load case id
	eid	- element id
	leid	- loaded element id
	ri,ci	- row and column indices

TKr, 28.2.2004

Definition at line 1220 of file trlinaxisym.cpp.

References addv(), edge_integral(), edgenodeval(), transtop::elements, loadcaset::elemload, fillm(), fillv(), gauss_points(), loadelt::give_bc(), transtop::give_elemnodes(), transtop::give_node_coord2d(), loadelt::give_nodval(), loadelt::give_trc(), loadelt::give_trr(), intordkm, elementt::ipp, transbclc::lc, mxv(), ned, nne, nned, nodes, nodval(), probdesct::savemode, Tb, Tp, transf_flux(), and Tt.

Referenced by res_boundary_flux().

void capacity_matrix	(	long	eid,
		long	ri,
		long	ci,
		matrix &	cm
	)

function computes capacity matrix for one medium

Parameters:

	eid	- element id
	ri,ci	- row and column indices
	cm	- capacity matrix

JK, 5.11.2001

Definition at line 453 of file trlinaxisym.cpp.

References matrix::a, approx(), bf_matrix(), transmat::capcoeff(), dofe, transtop::elements, fillm(), gauss_points_tr(), transcrsec::give_density(), transtop::give_elemnodes(), transtop::give_node_coord2d(), intordcm, intordkm, elementt::ipp, matrix::m, matrix::n, nne, nnj(), nodes, probdesct::savemode, Tc, Tm, Tp, and Tt.

Referenced by res_capacity_matrix().

void codnum	(	long *	cn,
		long	ri
	)

Definition at line 123 of file trlinaxisym.cpp.

References nne, and ordering.

Referenced by res_boundary_flux(), res_capacity_matrix(), res_conductivity_matrix(), res_convection_vector(), res_internal_fluxes(), res_quantity_source_vector(), and res_transmission_vector().

void conductivity_matrix	(	long	lcid,
		long	eid,
		long	ri,
		long	ci,
		matrix &	km
	)

function computes conductivity matrix of one medium

Parameters:

	lcid	- load case id
	eid	- element id
	km	- conductivity matrix

JK, 5.11.2001

Definition at line 528 of file trlinaxisym.cpp.

References matrix::a, approx(), bdbj(), bdbjac(), bf_matrix(), dofe, transtop::elements, fillm(), transtop::give_elemnodes(), transtop::give_node_coord2d(), grad_matrix(), elementt::ipp, transmat::matcond(), transmat::matcond2(), ncomp, nne, nodes, plsb(), plsc(), reallocm(), probdesct::savemode, Tm, Tp, elementt::transi, transmission_matrix(), and Tt.

Referenced by res_conductivity_matrix().

void convection_vector	(	vector &	v,
		long	lcid,
		long	eid,
		long	leid,
		long	ri,
		long	ci
	)

function computes contribution to the convection vector

{Gamma_2} N^T N dGamma * nodal_flux_values

Parameters:

	v	- array of nodal fluxes
	lcid	- load case id
	eid	- element id
	leid	- id of loaded element
	ri,ci	- row and column indices

JK, 19.8.2004

Definition at line 1015 of file trlinaxisym.cpp.

References addv(), edge_integral(), edgenodeval(), transtop::elements, loadcaset::elemload, fillm(), fillv(), gauss_points(), loadelt::give_bc(), transtop::give_elemnodes(), transtop::give_node_coord2d(), loadelt::give_nodval(), intordkm, elementt::ipp, transbclc::lc, mxv(), ned, nne, nned, nodes, nodval(), probdesct::savemode, Tb, Tp, and Tt.

Referenced by res_convection_vector().

void edge_integral	(	long	edg,
		vector &	x,
		vector &	y,
		long	intord,
		vector &	gp,
		vector &	w,
		vector &	coef,
		matrix &	km
	)

function integrates N^T c N over edges

Parameters:

	edg	- number of edge
	x,y	- arrays of nodal coordinates
	intord	- order of numerical integration
	gp	- array of coordinates of integration points
	w	- array of weights
	coef	- array of nodal values of coefficient
	km	- resulting matrix

JK

Definition at line 1290 of file trlinaxisym.cpp.

References matrix::a, approx(), bf_matrix(), intordb, jac1d_2d(), matrix::m, matrix::n, nne, and nnj().

Referenced by boundary_flux(), convection_vector(), transmission_matrix(), and transmission_vector().

void edgenodeval	(	long	edg,
		vector &	nodval,
		vector &	list
	)

function picks up nodal values on required edges

Parameters:

	edg	- number of required edge
	nodval	- array of nodal values
	list	- array of nodal values defined on all edges

JK, 19.8.2004

Definition at line 1472 of file trlinaxisym.cpp.

References ivector::a, fillv(), lintriangle_edgnod(), and nned.

Referenced by boundary_flux(), convection_vector(), and transmission_vector().

void give_approx_fun	(	double &	f,
		vector &	areacoord,
		long	i
	)

function returns one approximation function evaluated in required point

Parameters:

	f	- value of approximation function
	areacoord	- area coordinates
	i	- number of approximation function

JK, 1.2.2003

Definition at line 439 of file trlinaxisym.cpp.

void grad_matrix	(	matrix &	gm,
		vector &	b,
		vector &	c
	)

function assembles gradient of matrix of base functions

Parameters:

	gm	- gradient matrix
	b,c	- array containing node coordinates

JK, 25.9.2001

Definition at line 506 of file trlinaxisym.cpp.

References fillm(), and nne.

Referenced by conductivity_matrix(), internal_fluxes(), and intpointgrad().

void internal_fluxes	(	long	lcid,
		long	eid,
		vector &	ifl
	)

function computes internal fluxes of 1D problems for one transported matter finite element with linear approximation functions

Parameters:

	lcid	- number of load case
	eid	- number of element
	ifl	- vector of internal fluxes

JK, 31.3.2002

Definition at line 640 of file trlinaxisym.cpp.

References addv(), approx(), cmulv(), transmat::computenlfluxes(), dofe, transtop::elements, transtop::give_elemnodes(), transtop::give_node_coord2d(), transmat::givefluxes(), grad_matrix(), elementt::ipp, mtxv(), ncomp, nne, nodes, plsb(), plsc(), probdesct::savemode, Tm, Tp, and Tt.

Referenced by res_internal_fluxes().

void intpointgrad ( long eid )

function computes values in integration points from nodal values

eid - element id

JK, 25.9.2001

Definition at line 261 of file trlinaxisym.cpp.

References vector::a, ivector::a, dofe, transtop::elements, elemvalues(), gauss_points_tr(), transtop::give_code_numbers(), transtop::give_node_coord2d(), grad, grad_matrix(), intordcm, intordkm, elementt::ipp, mxv(), ncomp, ndofe, nne, probdesct::ntm, ordering, plsb(), plsc(), reallocv(), probdesct::savemode, transmat::storegrad(), Tm, Tp, and Tt.

Referenced by intpointgradients().

void intpointother ( long eid )

function approximates nodal values of array other to integration points

Parameters:

eid

- element id

JK, 17.9.2005

Definition at line 344 of file trlinaxisym.cpp.

References vector::a, approx(), transtop::elements, gauss_points_tr(), transtop::give_elemnodes(), intordcm, intordkm, transmat::ip, elementt::ipp, nodet::ncompother, nne, nodalotherval(), transtop::nodes, nodes, probdesct::ntm, intpointst::other, reallocv(), probdesct::savemode, Tm, Tp, and Tt.

Referenced by intpointothers().

void intpointval ( long eid )

function computes values in integration points from nodal values

eid - element id

JK, 25.9.2001

Definition at line 195 of file trlinaxisym.cpp.

References vector::a, ivector::a, approx(), intpointst::av, dofe, transtop::elements, elemvalues(), gauss_points_tr(), transtop::give_code_numbers(), intordcm, intordkm, transmat::ip, elementt::ipp, ndofe, nne, probdesct::ntm, ordering, reallocv(), probdesct::savemode, Tm, Tp, and Tt.

Referenced by intpointvalues().

void ipcoordblock	(	long	eid,
		long	ri,
		long	ci,
		double **	coord
	)

function returns coordinates of integration points

Parameters:

	eid	- element id
	ri	- row index
	ci	- column index
	coord	- array of coordinates

19.1.2002

Definition at line 141 of file trlinaxisym.cpp.

References vector::a, approx_nat(), gauss_points_tr(), transtop::give_node_coord2d(), intordkm, nne, and Tt.

void nod_others	(	long	lcid,
		long	eid,
		long	ri,
		long	ci
	)

function computes other values in nodes of element

Parameters:

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

10.5.2002

Definition at line 958 of file trlinaxisym.cpp.

References vector::a, ivector::a, transtop::elements, elemvalues(), transtop::give_code_numbers(), transtop::give_elemnodes(), transmat::givecompother(), transmat::givencompother(), h(), elementt::ipp, ndofe, nne, nod, ntm, other, reallocv(), probdesct::savemode, Tm, Tp, and Tt.

void quantity_source_vector	(	vector &	sv,
		vector &	nodval,
		long	eid,
		long	ri,
		long	ci
	)

function computes source vector of one matter on one element

{Omega} N^T N d Omega . s

Parameters:

	sv	- source vector of one matter
	lcid	- load case id
	eid	- element id
	ri,ci	- row and column indices of the block (ri must be equal to ci)

JK, 4.10.2001

Definition at line 597 of file trlinaxisym.cpp.

References matrix::a, addv(), approx(), bf_matrix(), dofe, fillm(), gauss_points_tr(), transtop::give_node_coord2d(), intordcm, mxv(), nne, nnj(), and Tt.

Referenced by res_quantity_source_vector().

void res_boundary_flux	(	vector &	f,
		long	lcid,
		long	eid,
		long	leid
	)

function assembles resulting element boundary flux vector

Parameters:

	f	- resulting boundary flux vector of one element
	lcid	- load case id
	eid	- element id
	leid	- element id

TKr, 28.2.2004

Definition at line 925 of file trlinaxisym.cpp.

References vector::a, boundary_flux(), codnum(), dofe, transtop::elements, locglob(), ntm, reallocv(), elementt::transi, and Tt.

Referenced by elem_transmission_flux().

void res_capacity_matrix	(	long	eid,
		matrix &	cm
	)

function assembles resulting element capacity matrix

Parameters:

	eid	- element id
	cm	- resulting capacity matrix of one element

JK, 6.1.2002

Definition at line 719 of file trlinaxisym.cpp.

References capacity_matrix(), codnum(), probdesct::diagcap, dofe, mat_localize(), ntm, reallocm(), and Tp.

Referenced by capacmat(), and res_internal_fluxes().

void res_conductivity_matrix	(	long	eid,
		long	lcid,
		matrix &	km
	)

function assembles resulting element conductivity matrix

Parameters:

	eid	- element id
	lcid	- load case id
	km	- resulting conductivity matrix of one element

JK, 6.1.2002

Definition at line 692 of file trlinaxisym.cpp.

References codnum(), conductivity_matrix(), dofe, mat_localize(), ntm, and reallocm().

Referenced by conductmat().

void res_convection_vector	(	vector &	f,
		long	lcid,
		long	eid,
		long	leid
	)

function assembles resulting element convection vector

Parameters:

	f	- resulting convection vector of one element
	lcid	- load case id
	eid	- element id
	leid	- element id

JK, 6.1.2002

Definition at line 765 of file trlinaxisym.cpp.

References vector::a, cmulv(), codnum(), convection_vector(), dofe, locglob(), and reallocv().

Referenced by elem_neumann_vector().

void res_internal_fluxes	(	long	eid,
		vector &	elemif
	)

function assembles resulting element internal fluxes vector

Parameters:

	eid	- element id
	elemif	- resulting internal fluxes vector of one element

JK, 6.1.2002

Definition at line 845 of file trlinaxisym.cpp.

References vector::a, codnum(), dofe, gtopology::give_code_numbers(), Gtt, internal_fluxes(), locglob(), mxv(), ndofe, nodalderivatives(), ntm, reallocm(), reallocv(), res_capacity_matrix(), and subv().

Referenced by internal_fluxes().

void res_quantity_source_vector	(	vector &	sv,
		vector &	nodval,
		long	lcid,
		long	eid
	)

function assembles resulting element source vector

Parameters:

	sv	- resulting source vector of one element
	lcid	- load case id
	eid	- element id

JK, 6.1.2002

Definition at line 824 of file trlinaxisym.cpp.

References vector::a, codnum(), dofe, locglob(), quantity_source_vector(), and reallocv().

Referenced by source_vector().

void res_transmission_vector	(	vector &	f,
		long	lcid,
		long	eid,
		long	leid
	)

function assembles resulting element transmission vector

Parameters:

	f	- resulting transmission vector of one element
	lcid	- load case id
	eid	- element id
	leid	- element id

JK, 6.1.2002

Definition at line 793 of file trlinaxisym.cpp.

References vector::a, codnum(), dofe, transtop::elements, locglob(), ntm, reallocv(), elementt::transi, transmission_vector(), and Tt.

Referenced by elem_newton_vector().

double total_integral	(	long	eid,
		vector &	nodval
	)

function computes element quantity integral

Parameters:

	eid	- element id
	nodval	- vector of quantity nodal values

Return values:

f

- element quantity integral

TKr, 30.1.2004

Definition at line 881 of file trlinaxisym.cpp.

References approx(), f, gauss_points_tr(), transtop::give_elemnodes(), transtop::give_node_coord2d(), nne, nodes, and Tt.

void transf_coeff	(	long	edg,
		vector &	coeff,
		vector &	list,
		long	eid,
		long	ri,
		long	ci,
		long	ipp,
		bocontypet *	bc
	)

Definition at line 1386 of file trlinaxisym.cpp.

References ivector::a, fillv(), transtop::give_elemnodes(), lintriangle_edgnod(), nne, nned, nodes, Tm, transmat::transmission_transcoeff(), and Tt.

Referenced by transmission_matrix().

void transf_flux	(	long	edg,
		vector &	coeff,
		vector &	list,
		vector &	trc,
		vector &	trr,
		long	eid,
		long	ri,
		long	ci,
		long	ipp,
		bocontypet *	bc
	)

Definition at line 1352 of file trlinaxisym.cpp.

References ivector::a, fillv(), transtop::give_elemnodes(), lintriangle_edgnod(), nne, nned, nodes, Tm, transmat::transmission_flux(), and Tt.

Referenced by boundary_flux().

void transf_val	(	long	edg,
		vector &	nodval,
		vector &	list,
		vector &	trc,
		vector &	trr,
		long	eid,
		long	ri,
		long	ci,
		long	ipp,
		bocontypet *	bc
	)

Parameters:

	edg	- number of required edge
	nodval	- array of transformed nodal values
	list	- array of nodal values defined on all edges
	trc	-
	trr	-
	ri,ci	- row and column indices
	ipp	- integration point number
	bc	- array defining boundary conditions