artisynth.core.femmodels

Class IntegrationPoint3d

java.lang.Object

artisynth.core.femmodels.IntegrationPoint3d

Direct Known Subclasses:

MFreeIntegrationPoint3d

public class IntegrationPoint3d
extends java.lang.Object

This class stores readonly and transient information for each integration point within a 3D element. To avoid excessive memory requirements, this data will generally be allocated statically within the each element class. The readonly data presents no problem in this regard, while the transient data will be valid only while stiffness update computations are being performed for the element.

Data that needs to be stored and retained per-element is declared in the IntegrationData3d class.

Constructor Summary

Constructors

Constructor and Description
IntegrationPoint3d(int nnodes)
IntegrationPoint3d(int nnodes, int npvals)
IntegrationPoint3d(int nnodes, int npvals, double s0, double s1, double s2, double w)

Method Summary

Modifier and Type	Method and Description
void	computeCoordsForRender(float[] coords, FemNode3d[] nodes)
double	computeGradient(Matrix3d F, FemNode3d[] nodes, Matrix3d invJ0) Computes the current deformation gradient at this integration point.
void	computeGradientForRender(Matrix3d Fmat, FemNode3d[] nodes, Matrix3d invJ0)
double	computeInverseJacobian(Matrix3d invJ, FemNode3d[] nodes) Computes the current inverse Jacobian at this integration point.
double	computeInverseRestJacobian(Matrix3d invJ0, FemNode3d[] nodes)
void	computeJacobian(Matrix3d J, FemNode3d[] nodes) Computes the current Jacobian at this integration point.
void	computeJacobian(Matrix3d J, FemNode3d[] nodes, FemElement.ElementClass type)
void	computeJacobian(Matrix3d J, FemNode3d[] nodes, FemElement.ElementClass type, FemNode3d.CoordType ctype)
double	computeJacobianDeterminant(FemNode3d[] nodes) Computes the determinant of the current Jacobian at this integration point.
void	computePosition(Point3d pos, FemElement3dBase elem)
void	computePosition(Point3d pos, FemNode3d[] nodes)
void	computeRestJacobian(Matrix3d J0, FemNode3d[] nodes)
void	computeRestPosition(Point3d pos, FemElement3dBase elem)
void	computeRestPosition(Point3d pos, FemNode3d[] nodes)
Vector3d[]	computeShapeGradient(Matrix3d invJ) Computes and returns the gradient dN/dx of the shape functions, given an inverse Jacobian.
void	computeShapeGradient(Matrix3d invJ, Vector3d[] out) Computes and returns the gradient dN/dx of the shape functions, given an inverse Jacobian.
static IntegrationPoint3d	create(FemElement3dBase elem, double s0, double s1, double s2, double w) Create an integration point for a given element at a specific set of natural coordinates.
Vector3d	getCoords()
FemElement.ElementClass	getElementClass()
Vector3d[]	getGNs()
int	getNumber() Returns the number of this integration point.
VectorNd	getPressureWeights()
Vector3d[]	getShapeGradient() Returns the gradient dN/dx of the shape functions.
VectorNd	getShapeWeights()
double	getWeight()
void	setCoords(double s0, double s1, double s2)
void	setNumber(int num) Used internally by the system to set the number for this integration point.
void	setPressureWeights(VectorNd vals)
void	setShapeGrad(int i, Vector3d dNds)
void	setShapeWeights(VectorNd vals)
void	setWeight(double w)
Vector3d[]	updateShapeGradient(Matrix3d invJ) Updates and returns the gradient dN/dx of the shape functions, given an inverse Jacobian.

Methods inherited from class java.lang.Object

equals, getClass, hashCode, notify, notifyAll, toString, wait, wait, wait

Constructor Detail

IntegrationPoint3d

public IntegrationPoint3d(int nnodes,
                          int npvals,
                          double s0,
                          double s1,
                          double s2,
                          double w)

IntegrationPoint3d

public IntegrationPoint3d(int nnodes)

IntegrationPoint3d

public IntegrationPoint3d(int nnodes,
                          int npvals)

Method Detail

getElementClass

public FemElement.ElementClass getElementClass()

getNumber

public int getNumber()

Returns the number of this integration point. This will be in the range 0 to numi-1, where numi is the number of integration points for the associated elements.

setNumber

public void setNumber(int num)

Used internally by the system to set the number for this integration point.

create

public static IntegrationPoint3d create(FemElement3dBase elem,
                                        double s0,
                                        double s1,
                                        double s2,
                                        double w)

Create an integration point for a given element at a specific set of natural coordinates.

Parameters:

elem - element to create the integration point for

s0 - first coordinate value

s1 - second coordinate value

s2 - third coordinate value

w - weight

Returns:

new integration point

setCoords

public void setCoords(double s0,
                      double s1,
                      double s2)

getCoords

public Vector3d getCoords()

getWeight

public double getWeight()

setWeight

public void setWeight(double w)

setShapeWeights

public void setShapeWeights(VectorNd vals)

getShapeWeights

public VectorNd getShapeWeights()

setShapeGrad

public void setShapeGrad(int i,
                         Vector3d dNds)

setPressureWeights

public void setPressureWeights(VectorNd vals)

getPressureWeights

public VectorNd getPressureWeights()

computeJacobian

public void computeJacobian(Matrix3d J,
                            FemNode3d[] nodes)

Computes the current Jacobian at this integration point.

Parameters:

J - returns the Jacobian

nodes - FEM nodes, used to obtain the element node positions

computeJacobian

public void computeJacobian(Matrix3d J,
                            FemNode3d[] nodes,
                            FemElement.ElementClass type)

computeJacobian

public void computeJacobian(Matrix3d J,
                            FemNode3d[] nodes,
                            FemElement.ElementClass type,
                            FemNode3d.CoordType ctype)

computeRestJacobian

public void computeRestJacobian(Matrix3d J0,
                                FemNode3d[] nodes)

computeInverseRestJacobian

public double computeInverseRestJacobian(Matrix3d invJ0,
                                         FemNode3d[] nodes)

computeJacobianDeterminant

public double computeJacobianDeterminant(FemNode3d[] nodes)

Computes the determinant of the current Jacobian at this integration point.

Parameters:

nodes - FEM nodes, used to obtain the element node positions

Returns:

determinant of the Jacobian

computeGradient

public double computeGradient(Matrix3d F,
                              FemNode3d[] nodes,
                              Matrix3d invJ0)

Computes the current deformation gradient at this integration point.

Parameters:

F - returns the deformation gradient

nodes - FEM nodes, used to obtain the element node positions

invJ0 - inverse rest position Jacobian

Returns:

determinant of F

computeInverseJacobian

public double computeInverseJacobian(Matrix3d invJ,
                                     FemNode3d[] nodes)

Computes the current inverse Jacobian at this integration point.

Parameters:

invJ - returns the inverse Jacobian

nodes - FEM nodes, used to obtain the element node positions

Returns:

determinant of the Jacobian

computeGradientForRender

public void computeGradientForRender(Matrix3d Fmat,
                                     FemNode3d[] nodes,
                                     Matrix3d invJ0)

computePosition

public void computePosition(Point3d pos,
                            FemNode3d[] nodes)

computePosition

public void computePosition(Point3d pos,
                            FemElement3dBase elem)

computeRestPosition

public void computeRestPosition(Point3d pos,
                                FemElement3dBase elem)

computeRestPosition

public void computeRestPosition(Point3d pos,
                                FemNode3d[] nodes)

computeCoordsForRender

public void computeCoordsForRender(float[] coords,
                                   FemNode3d[] nodes)

computeShapeGradient

public Vector3d[] computeShapeGradient(Matrix3d invJ)

Computes and returns the gradient dN/dx of the shape functions, given an inverse Jacobian. The gradient is supplied as an array of 3-vectors, one for each shape function, because this is more convenient for computation.

Parameters:

invJ - inverse Jacobian

Returns:

shape function gradient (can be modifed).

computeShapeGradient

public void computeShapeGradient(Matrix3d invJ,
                                 Vector3d[] out)

Computes and returns the gradient dN/dx of the shape functions, given an inverse Jacobian. The gradient is supplied as an array of 3-vectors, one for each shape function, because this is more convenient for computation.

Parameters:

invJ - inverse Jacobian

out - shape function gradient to populate

updateShapeGradient

public Vector3d[] updateShapeGradient(Matrix3d invJ)

Updates and returns the gradient dN/dx of the shape functions, given an inverse Jacobian. The gradient is supplied as an array of 3-vectors, one for each shape function, because this is more convenient for computation.

Parameters:

invJ - inverse Jacobian

Returns:

shape function gradient (must not be modifed).

getShapeGradient

public Vector3d[] getShapeGradient()

Returns the gradient dN/dx of the shape functions. This is supplied as an array of 3-vectors, one for each shape function, because this is more convenient for computation.

Returns:

shape function gradient (must not be modifed).

getGNs

public Vector3d[] getGNs()