Vertical Stand

1 Description

Figure 1: CAD Geometry

The vertical stand (see xx--CrossReference--dft--fig:cad--xx) represents a structural part of a machine tool. On one of its surfaces a pair of guide rails is located. Caused by a machining process a tool slide is moving on these rails. The machining process produces a certain amount of heat which is transported through the structure into the vertical stand. This heat source is considered to be a temperature input at the guide rails. The induced temperature field, denoted by $T$ is modeled by the heat equation

$c_p\rho\frac{\partial{T}}{\partial{t}}=\Delta T=0$

with the boundary conditions

$\lambda\frac{\partial T}{\partial n}=f \qquad\qquad\qquad$ on $\Gamma_{slide}$ (surface where the tool slide is moving on the guide rails),

describing the heat transfer between the tool slide and the vertical stand. The heat transfer to the ambiance is given by the fixed Robin-type boundary condition

$\lambda\frac{\partial T}{\partial n}=\kappa_i(T-T_i^{ext})$ on $\Gamma_{surf}$ (remaining boundaries).

The motion driven temperature input and the associated change in the temperature field $T$ lead to deformations $u$ within the stand structure. Further, it is assumed that no external forces are induced to the system such that the deformation is purely driven by the change of temperature. Since the mechanical behavior of the machine stand is assumed to be much faster than the propagation of the thermal field, it is sufficient to consider the stationary linear elasticity equations

$\begin{align} -\operatorname{div}(\sigma(u)) &= 0&\text{ on }\Omega,\\ \varepsilon(u) &= {\mathbf{C}}^{-1}:\sigma(u)+\beta(T-T_{ref})I_d&\text{ on }\Omega,\\ {\mathbf{C}}^{-1}\sigma(u) &=\frac{1+\nu}{E_u}\sigma(u)-\frac{\nu}{E_u}\text{tr}(\sigma(u))I_d&\text{ on }\Omega,\\ \varepsilon(u) &= \frac{1}{2}(\nabla u+\nabla u^T)&\text{on }\Omega. \end{align}$

Geometrical dimensions:

Stand: Width ($x$ direction): $519mm$, Height ($y$ direction): $2\,010mm$, Depth ($z$ direction): $480mm$

Guide rails: $y\in [519, 2\,004] mm$

Slide: Width: $430 mm$, Height: $500mm$, Depth: $490 mm$

1.1 Discretized Model

The solid model has been generated and meshed in ANSYS. For the spatial discretization the finite element method with linear Lagrange elements has been used and is implemented in FEniCS. The resulting system of ordinary differential equations (ODE), representing the thermal behavior of the stand, reads

$\begin{align} E \frac{\partial}{\partial t} T(t) &= A(t)T(t) + B(t)z(t), \\ T(0) &= T_0, \end{align}$

with $t>0$ and a system dimension of $n=16\,626$ degrees of freedom and $m=6$ inputs. Note that $A(.)\in\mathbb{R}^{n\times n}$ and $B(.)\in\mathbb{R}^{n\times m}$ are time-dependent matrix-valued functions. That is, the underlying model is represented by a linear time-Varying (LTV) state-space system. More precisely, here the time dependence originates from the change of the boundary condition on $\Gamma_{slide}$ due to the motion of the tool slide. The system input is, according to the boundary conditions, given by

$z_i=\begin{cases} q, i=1,\\ \kappa_i T_i^{ext}, i=2,\dots,6 \end{cases}.$

The discretized stationary elasticity model becomes

$Mu(t)=KT(t).$

For the observation of the displacements in single points/regions of interest an output equation of the form

$y(t) = C u(t)$

is given.

Exploiting the one-sided coupling of the temperature and deformation fields, and reorganizing the elasticity equation in the form $u(t)=M^{-1}KT(t)$, the heat equation and the elasticity model can easily be combined via the output equation. Finally, the thermo-elastic control system is of the form

$\begin{align} E \frac{\partial}{\partial t} T(t) &= A(t)T(t) + B(t)z(t), \\ y(t)& = \tilde{C}T(t),\\ T(0) &= T_0, \end{align}$

where the modified output matrix $\tilde{C}=CM^{-1}KT(t)$ includes the entire elasticity information.

The motion of the tool slide and the associated variation of the affected input boundary are modeled by two different system representations. The following specific model representations have been developed and investigated in , .

1.2 Switched Linear System

Figure 2: Schematic segmentation

For the model description as a switched linear system, the guide rails of the machine stand are modeled as fifteen equally distributed horizontal segments with a height of $99mm$ (see a schematic depiction in xx--CrossReference--dft--fig:segm--xx). Any of these segments is assumed to be completely covered by the tool slide if its midpoint (in y-direction) lies within the height of the slide. On the other hand, each segment whose midpoint is not covered is treated as not in contact and therefore the slide always covers five to six segments at each time. Still, the covering of six segments does not have a significant effect on the behavior of the temperature and displacement fields. Due to that and in order to keep the number of subsystems small, this scenario will be neglected. Than, in fact eleven distinct, discrete boundary condition configurations for the stand model that are prescribed by the geometrical dimensions of the segmentation and the tool slide are defined. These distinguishable setups define the subsystems

$\begin{align} E\dot{T}&=A_{\alpha}T+B_{\alpha}z,\\ y&=\tilde{C}T, \end{align}$

of the switched linear system , where $\alpha$ is a piecewise constant function of time, which takes its value from the index set $\mathcal{J}=\{1,\dots,11\}$. To be more precise, the switching signal $\alpha$ implicitly maps the slide position to the number of the currently active subsystem.

1.3 Linear Parameter-Varying System

For the parametric model description, the FE nodes located at the guide rails are clustered with respect to their y-coordinates. This results in $233$ distinct layers in y-direction. According to these layers, the matrices $A(t)=A(\mu(t)), B(t)=B(\mu(t))$ are defined in a parameter-affine representation of the form

$\begin{align} A(\mu) = A_0+f_1(\mu)A_1+...+f_{m_A}(\mu)A_{m_A},\\ B(\mu) = B_0+g_1(\mu)B_1+...+g_{m_B}(\mu)B_{m_B} \end{align}$

with the scalar functions $f_i,g_j\in\{0,1\}, i=1,...,m_A, j=1,...,m_B$ selecting the active layers, covered by the tool slide. The matrix $A_0$ represents the discretization of the Laplacian $\Delta$, as well as the discrete portions of $\kappa_i T$ from the Robin-type boundaries, associated to the exchange with the ambience given by the inputs $z_i,~i=2,...,6$. The remaining summands $A_j,~j=1,...,m_A$ denote the discretization of $\kappa_i T$ from the moving Robin-type boundaries.

Note that in general the number of summands of these representations need not to be equal. Still, according to the number of layers, for this example, it holds $m_A=m_B=233$. For more details on parametric models, see e.g., and the references therein.

Then, the final linear parameter-varying (LPV) reformulation of the above LTV system reads

$\begin{align} E\dot{T}&=A(\mu)T+B(\mu)z,\\ y&=\tilde{C}T. \end{align}$

2 Acknowledgement & Origin

The base model was developed , in the Collaborative Research Centre Transregio 96 Thermo-Energetic Design of Machine Tools funded by the Deutsche Forschungsgemeinschaft .

3 Data

3.1 Switched System Data

The data file VertStand_SLS.tar.gz contains the matrices

$\begin{align} E\in\mathbb{R}^{n\times n}, \tilde{C}\in\mathbb{R}^{q\times n}, A_\alpha\in\mathbb{R}^{n\times n}, B_\alpha\in\mathbb{R}^{n\times m}, \alpha=1,\dots,11. \end{align}$

defining the subsystems of the switched linear system. The matrices $A_\alpha$ and $B_\alpha$ are numbered according to the slide position in descending order (1 - uppermost slide position / 11 - lowest slide position).

System dimensions:

$n=16\,626, m=6, q=27$

3.2 Parametric System Data

The data file VertStand_PAR.tar.gz contains the matrices

$E,A_j\in\mathbb{R}^{n\times n},j=1,...,234, B_1\in\mathbb{R}^{n\times 233}, B_2\in\mathbb{R}^{n\times 5}$ and $\tilde{C}\in\mathbb{R}^{q\times n}, q=27, n=16\,626,$

as well as a file ycoord_layers.mtx containing the y-coordinates of the layers located on the guide rails.

Here $B_1$ contains all columns corresponding to the different layers and $B_2$ correlates to the boundaries where the external temperatures act on.

In order to set up the parameter dependent matrices $A(\mu),B(\mu)$ the active matrices $A_i$ and columns $B_1(:,i)$ associated to the covered layers have to be identified via the current position of the tool slide and the y-coordinates of the different layers given in ycoord_layers.mtx. Then, $B(\mu)$ has to be set up in the form $B(\mu)=[\sum_{i\in id_{active}}\!\!\!\!\!B_1(:,i),B_2]$. Similar to the switched systems approach, the layers are numbered in descending order.

4 Citation

To cite this benchmark, use the following references:

• For the benchmark itself and its data:

The MORwiki Community, Vertical Stand. MORwiki - Model Order Reduction Wiki, 2018. http://modelreduction.org/index.php/Vertical_Stand

@MISC{morwiki_vertstand,
  author =       {{The MORwiki Community}},
  title =        {Vertical Stand},
  howpublished = {{MORwiki} -- Model Order Reduction Wiki},
  url =          {http://modelreduction.org/index.php/Vertical_Stand},
  year =         2014
}

• For the background on the benchmark:

@Article{morLanSB14,
  author =       {Lang, N. and Saak, J. and Benner, P.},
  title =        {Model Order Reduction for Systems with Moving Loads},
  journal =      {at-Automatisierungstechnik},
  volume =       62,
  number =       7,
  pages =        {512--522},
  year =         2014,
  doi =          {10.1515/auto-2014-1095}
} 

6 Contact

Jens Saak