Skip to main content

Displacement Calculation

For the displacement calculation the polygons resulting from these cuts is used

  • 1D – the element cut of the plane is expanded to the enclosing box
  • 2D – the element cut line is expanded to the top and bottom surfaces
  • 3D – the polygon on the cut plane

The displacement and rotations are based on a displacement plane described for each component jj by

uj=aj+bjy+cjzu_{j} = a_{j} + b_{j}y + c_{j}z

A least-squares fit across the nn points gives a set of equations (one for each of x,y,zx,y,z directions)

[nyiziyiyi2yiziziyizizi2]{ajbjcj}={ui,jui,jyiui,jzi}\begin{bmatrix} n & \sum_{}^{}y_{i} & \sum_{}^{}z_{i} \\ \sum_{}^{}y_{i} & \sum_{}^{}{y_{i}}^{2} & \sum_{}^{}{y_{i}z_{i}} \\ \sum_{}^{}z_{i} & \sum_{}^{}{y_{i}z_{i}} & \sum_{}^{}{z_{i}}^{2} \\ \end{bmatrix}\begin{Bmatrix} a_{j} \\ b_{j} \\ c_{j} \\ \end{Bmatrix} = \begin{Bmatrix} \sum_{}^{}u_{i,j} \\ \sum_{}^{}u_{i,j}y_{i} \\ \sum_{}^{}u_{i,j}z_{i} \\ \end{Bmatrix}

Then the displacement are

ua,x=axua,y=ayua,z=az\begin{aligned}{u}_{a,x} = a_{x} \\ {u}_{a,y} = a_{y} \\ {u}_{a,z} = a_{z}\end{aligned}

And rotations are about yy and zz come from the curvature terms in the displacement equation in the xx direction

θa,y=cx\theta_{a,y} = c_{x}
θa,z=bx\theta_{a,z} = b_{x}

The twisting rotation is calculated from

θa,x=1ni=1n((ui,zu~z)yi(ui,yu~y)zi(yi2+zi2))\theta_{a,x} = \frac{1}{n}\sum_{i = 1}^{n}\left( \frac{\left( u_{i,z} - {\widetilde{u}}_{z} \right)y_{i} - \left( u_{i,y} - {\widetilde{u}}_{y} \right)z_{i}}{\left( {y_{i}}^{2} + {z_{i}}^{2} \right)} \right)