public override Vector IntegralArgument_Fe(Vector Xi)//override Fe
{
int NNPE = ElementNodes.Length;
Vector N;
Vector Temp1 = new Vector(NNPE);
Vector Temp2 = new Vector(NNPE);
double DsDXi = Calculate_NDLine_DsDXi(Xi, ElementNodes, out N);
double t_x = Interpolator.Interpolate_Variable(ForceX.Values, N);
double t_y = Interpolator.Interpolate_Variable(ForceY.Values, N);
Temp1 = (t_x * DsDXi) * N;
Temp2 = (t_y * DsDXi) * N;
Vector BigFe = new Vector(NNPE * 2);
for (int i = 0; i < NNPE; i++)
{
BigFe.Values[i * 2] = Temp1.Values[i];
BigFe.Values[i * 2 + 1] = Temp2.Values[i];
}
return(BigFe);
}
public Vector Calculate_U(Vector Xi)
{
Vector[] Ue = Get_ElementNodal_U(Element_ND.U);
return(Interpolator.Interpolate_Variable(Xi, Ue));
}
public double Calculate_U(Vector Xi)
{
double[] Ue = Get_ElementNodal_U();
return(Interpolator.Interpolate_Variable(Xi, Ue));
}