private Matrix LocalStiffnessMatrixFunction(double ksi, double eta)
{
Matrix derivativeMatrix = GetLocalDerivativeMatrix(elementCurrent, ksi, eta);
JacobianRectangular J = new JacobianRectangular();
J.Element = elementCurrent;
return derivativeMatrix * ConstMatrix * Matrix.Transpose(derivativeMatrix) * J.GetJacobianDeterminant(ksi, eta);
}
protected Matrix LocalNonlinearMatrixFunction(double ksi, double eta)
{
Matrix derivativeMatrix = GetLocalDerivativeMatrix(elementCurrent, ksi, eta);
Matrix derivativeUMatrix = GetLocalUDerivativeMatrix(previousU, elementCurrent, ksi, eta);
JacobianRectangular J = new JacobianRectangular();
J.Element = elementCurrent;
return derivativeMatrix * derivativeUMatrix * Matrix.Transpose(derivativeMatrix) * J.GetJacobianDeterminant(ksi, eta);
}
private Matrix LocalMassMatrixFunction(double ksi, double eta)
{
Matrix baseFunctionsMatrix = GetLocalBaseFunctionsMatrix(ksi, eta);
JacobianRectangular J = new JacobianRectangular();
J.Element = elementCurrent;
return baseFunctionsMatrix * Matrix.Transpose(baseFunctionsMatrix) * J.GetJacobianDeterminant(ksi, eta);
}
private Matrix LocalVectorFunction(double ksi, double eta)
{
Matrix derivativeMatrix = GetLocalDerivativeMatrix(elementCurrent, ksi, eta);
JacobianRectangular J = new JacobianRectangular();
J.Element = elementCurrent;
Matrix VariableVectorOnElement = GetVariableVectorOnElement(elementCurrent, ksi, eta);
return derivativeMatrix * VariableVectorOnElement * J.GetJacobianDeterminant(ksi, eta);
}
