Commit 754bac5d by lburger

Adapatation of VolUniDirShell and VolCylShell jacobian transformations.

parent 6c9dbcb4
Pipeline #1031 passed with stage
in 44 seconds
......@@ -483,7 +483,7 @@ double Transformation(int Dim, int Type, struct Element * Element, MATRIX3x3 * J
"Valid parameters: Radius1 Radius2 CenterX CenterY CenterZ Power 1/Infinity");
}
}
else if(Type == JACOBIAN_CYL_SHELL){
else if(Type == JACOBIAN_VOL_CYL_SHELL){
if(Element->JacobianCase->NbrParameters >= 3)
Axis = (int)Element->JacobianCase->Para[2];
if(Element->JacobianCase->NbrParameters >= 4)
......@@ -571,7 +571,7 @@ double Transformation(int Dim, int Type, struct Element * Element, MATRIX3x3 * J
break;
}
}
else if(Type == JACOBIAN_CYL_SHELL){
else if(Type == JACOBIAN_VOL_CYL_SHELL){
switch(Axis) {
case 1: R = sqrt(SQU(Y-Cy)+SQU(Z-Cz)); YR = (Y-Cy)/R; ZR = (Z-Cz)/R; dRdy = (Y-Cy)/R; dRdz = (Z-Cz)/R; break;
case 2: R = sqrt(SQU(X-Cx)+SQU(Z-Cz)); XR = (X-Cx)/R; ZR = (Z-Cz)/R; dRdx = (X-Cx)/R; dRdz = (Z-Cz)/R; break;
......@@ -1084,7 +1084,7 @@ double JacobianVolCylShell3D(struct Element * Element, MATRIX3x3 * Jac)
double DetJac1, DetJac2;
DetJac1 = JacobianVol3D(Element, &Jac1);
DetJac2 = Transformation(_3D, JACOBIAN_CYL_SHELL, Element, &Jac2);
DetJac2 = Transformation(_3D, JACOBIAN_VOL_CYL_SHELL, Element, &Jac2);
Get_ProductMatrix(_3D, &Jac1, &Jac2, Jac);
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment