# Radioss Coordinate System

- Global system ($X$, $Y$, and $Z$)
- Natural system (isoparametric frame) ($\xi ,\eta ,\zeta $)
- Local element coordinate system ($x$, $y$, and $z$)

($X$, $Y$, and $Z$) - Global Cartesian fixed system:

($\xi ,\eta ,\zeta $) - Natural system (non-normalized coordinate system).

- $\xi $ is from middle point of Line 14 to middle point of Line 23.
- $\eta $ is from middle point of Line 12 to middle point of Line 34.
- Plane ($\xi ,\eta $) is in the middle surface of shell element and $\zeta $ is normal of the middle surface.

- $z$ is normal of middle surface.
- ($x$ and $y$) are in the middle surface
- $x$ and $y$ are positioned so that they have same angle between $x$ and $\xi $, $y$ and $\eta $

The origin of ($\xi ,\eta ,\zeta $) and ($x$, $y$, and $z$) are the same as it is at the intersection point of middle point line.

## 3-node Shell Element

($X$, $Y$, and $Z$) - Global Cartesian fixed system

($\xi ,\eta ,\zeta $) - Natural system (non-normalized coordinate system).

- $\xi $ is from Node 1 to Node 2.
- $\eta $ is from Node 1 to Node 3.
- Plane ($\xi ,\eta $) is in the middle surface of shell element and $\zeta $ is normal of the middle surface.

- $z$ is normal of middle surface.
- $x$ is from Node 1 to Node 2.
- $y$ is orthogonal to $x$ and ($x$ and $y$) are in the middle surface.

## Solids and Thick Shells (hexa)

- Global system ($X$, $Y$, and $Z$)
- Natural system ($r$, $s$, and $t$)
- Local element coordinate system ($x$, $y$, and $z$)
- Material system

($X$, $Y$, and $Z$) - Global Cartesian fixed system

- $r$ is from the center of surface (1, 2, 6, and 5) to center of surface (4, 3, 7, and 8)
- $s$ is from the center of surface (1, 2, 3, and 4) to center of surface (5, 6, 7, and 8)
- $t$ is from the center of surface (1, 4, 8, and 5) to center of surface (2, 3, 7, and 6)

- $r$ is also from middle point of Line 1' and 2' to middle point of Line 3' and 4'.
- $t$ is also from middle point of Line 1' and 2' to middle point of Line 3' and 4'
- $n$ is normal of middle surface (1', 2', 3', and 4')

($x$, $y$, and $z$) - Local coordinate system (orthogonal, normalized elemental coordinate system).

Local coordinate system in middle surface (1', 2', 3', and 4') is the same as the local coordinate system in middle surface (1, 2, 3, and 4) for shell element. r in solid is the same as $\xi $ in shell element.

## Tetra Elements

- $r$ is from node 4 to node 1
- $s$ is from node 4 to node 2
- $t$ is from node 4 to node 3

## Material System

(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|

V_{X} |
V_{Y} |
V_{Z} |
$$\varphi $$ |

Use vector $V$ and angle $\varphi $ material direction 1 ($m1$) can be defined (along normal n project vector $V$ to middle surface and get vector $V\text{'}$. Rotate angle $$\varphi $$ of vector $V\text{'}$ then get material direction $m1$. Material direction $m1$ is normally the fiber direction. For composite, a different ply could be defined with one vector $V$ and different $$\varphi $$.

- $n$ is normal of shell middle surface

`I`

_{orth}can determine the relative orientation of the material system.

`I`_{orth}=0 (default): The orthotropic direction follows the local co-rotational reference. The angle between $x$ and $m1$ is constant during the simulation. Internal force is computed in local frame and then rotated to the global system. This formulation is more accurate, if a large rotation occurs.`I`_{orth}=1: The orthotropic direction is attached to the local isoparametric frame. The angle between $\xi $ and $m1$ is updated during the simulation. It is updated in a way that projection of vector $m1$ to $\xi $ and $\eta $ is always constant during the simulation. Pure shear could not well descripted with this method, but traction could well be described. So, this method usually is used to defined fiber direction in airbag.

`I`

_{orth}), like shell elements. In /PROP/TYPE6 (SOL_ORTH) use the option

`I`

_{P}to determine the reference plane.

`I`_{P}=0- Use
`Skew_ID` `I`_{P}=1- Plane ($r$,$s$) + angle $$\varphi $$
`I`_{P}=2- Plane ($s$,$t$) + angle $$\varphi $$
`I`_{P}=3- Plane ($t$,$r$) + angle $$\varphi $$
`I`_{P}=11- Plane ($r$,$s$) + orthogonal projection of references vector $$V$$ on plane ($r$,$s$)
`I`_{P}=12- Plane ($s$,$t$) + orthogonal projection of references vector $$V$$ on plane ($s$,$t$)
`I`_{P}=13- Plane ($t$,$r$) + orthogonal projection of references vector $$V$$ on plane ($t$,$r$)

- Definition is the same for any
`I`_{solid}and`I`_{frame}parameters - In the simplest case, material directions $m1$, $m2$ and
$m3$ directly with skew (
`I`_{P}=0) are recommended - For
`I`_{P}> 0 the isoparametric, non-orthogonal system $r$, $s$, and $t$, is used to determine material directions.- First material axis $m1$ is determined according to
`I`_{P}. - For example, for
`I`_{P}=1 - The first material axis $m1$ and $m2$ is orthogonal and rotated by angle $$\varphi $$ in the ($r\text{'}$ and $s\text{'}$) plane.
- The third material axis $m3$ is normal of $m1$ and $m2$ plane (vector product of $m1$ and $m2$).
- The ($r\text{'}$, $s\text{'}$, and $t\text{'}$) system is orthogonal and it is generated from non-orthogonal isoparametric system ($r$, $s$, and $t$).

- First material axis $m1$ is determined according to

Depending on `I`_{solid} and `I`_{frame} parameters, three definitions of systems are used in
Radioss for hexa elements (8-noded bricks) using
/PROP/TYPE6 (SOl_ORTH).

### Global System Definition

- Definition 1: Solids,
`I`_{solid}=1, 2, 17 +`I`_{frame}=0, 1 (default)Global system is used, no element system (non-co-rotational formulation) available.

### Element System Definition

- Definition 2: Solids,
`I`_{solid}=1, 2, 17 +`I`_{frame}=2Element system (with

`I`_{frame}=2 co-rotational formulation) is used. - Definition 3: Solids,
`I`_{solid}=14 or 24`I`_{frame}parameter has no effect. Element system is used, and co-rotational formulation defined already.

`I`

_{orth}) keeps the same orientation with respect to the local (co-rotating) frame and is therefore also co-rotating.