Designing Actuation Concepts for Adaptive Slabs with Integrated Fluidic Actuators Using Inﬂuence Matrices

: Previous work has shown that ﬂoor slabs make up most of the material mass of building structures and are typically made of reinforced concrete. Considering the associated resource consumption and greenhouse gas emissions, new approaches are needed in order to reduce the built environment’s impact on the ongoing climate crisis. Various studies have demonstrated that adaptive building structures offer a potential solution for reducing material resource consumption and associated emissions. Adaptive structures have the ability to improve load-bearing performance by speciﬁcally reacting to external loads. This work applies the concept of adaptive structures to reinforced concrete slabs through the integration of ﬂuidic actuators into the cross-section. The optimal integration of actuators in reinforced concrete slabs is a challenging interdisciplinary design problem that involves many parameters. In this work, actuation inﬂuence matrices are extended to slabs and used as an analysis and evaluation tool for deriving actuation concepts for adaptive slabs with integrated ﬂuidic actuators. To deﬁne requirements for the actuator concept, a new procedure for the selection of actuation modes, actuator placement and the computation of actuation forces is developed. This method can also be employed to compute the required number of active elements for a given load case. The new method is highlighted in a case study of a 2 m × 2 m ﬂoor.


Previous Work
The construction industry is responsible for around 60% of global resource consumption, 50% of global waste generation, approximately 40% to 50% of global emissions of greenhouse gases and more than 35% of global energy consumption [1,2]. Identifying potential savings and thus reducing resource consumption, energy expenditure and emissions is essential for the sustainable development of our built environment.
Previous analyses based on [3][4][5] have shown that floor slabs in building structures contribute significantly to the total mass. Reinforced concrete slabs are one of the most widespread components due to the low cost and ease of construction [6]. Reinforced concrete construction involves high process-related greenhouse gas emissions, which is a cause of concern for the ongoing climate crisis. In addition, concrete structures require considerable amounts of sand, which is becoming increasingly scarce [7,8]. The load transfer in floor slabs takes place primarily through bending. Usually, this means that serviceability limit states, mostly deflection limits, are the governing criteria for sizing [9]. For example, it may be necessary to increase the thickness of a slab in order to limit bending-induced deformations and thus prevent damage to adjacent structural and non-structural elements. Optimisation strategies for passive reinforced concrete slabs quickly lead to a shift in the

New Contribution
The generation of forces and strokes with the use of auxiliary energy can be ach through a variety of functional principles. Divided into sub-functions, an actuator co of an energy converter, a conductor and supplier [38,39]. The actuator concept is u stood as the ideal combination of these three sub-functions.
In two-way slabs, bending moments from external loads do not just occur in direction. The actuation concept of integrated fluidic actuators presented for uniaxia transfer must therefore be extended to multi-axial load transfer. For this purpose pairs in several spatial directions above the neutral plane are considered. Pushing face pairs apart through actuation induces moments in the x, y and xy directions. T tuation concept is influenced by different design parameters such as the control obj the actuation modes, the adaptation level, the actuator placement and the computat the required actuation forces.
The design of an adaptive structure is a multi-disciplinary and iterative proc which various factors have to be taken into account. In this work, actuation concep not developed as usual by employing stock actuators. The following method can be to create a list of requirements for the actuator concept. Conversely, the method a lows the requirements of the actuator concept to be taken into account. An optimal tive structure, in terms of material and emission efficiency, can be achieved using su integrated approach whereby the actuation and actuator concept are harmonised each other.
This work builds on and extends the methods developed in [22] for the integ of fluidic actuators in reinforced concrete structures and [34] concerning the use of tion influence matrices for the design of the actuation concepts respectively of the vidual design parameters. New contributions offered by this work are: 1. The elaboration of a new actuation concept based on fluidic actuators that c integrated into floor slabs to control the response under loading. 2. Extension of actuation influence matrices to two-way slabs. 3. Computation of actuation forces and design spaces for actuator placement bas influence matrices. 4. A new methodology that makes it possible to co-design actuation and actuato cepts for adaptive floor slabs. The opposing steel sheets (pair of faces) are pushed apart by hydraulic pressure. Numerical and experimental investigations of integrated fluidic actuators in reinforced concrete beams have shown potential for the manipulation of deformations [18,22,36,37].

New Contribution
The generation of forces and strokes with the use of auxiliary energy can be achieved through a variety of functional principles. Divided into sub-functions, an actuator consists of an energy converter, a conductor and supplier [38,39]. The actuator concept is understood as the ideal combination of these three sub-functions.
In two-way slabs, bending moments from external loads do not just occur in the x-direction. The actuation concept of integrated fluidic actuators presented for uniaxial load transfer must therefore be extended to multi-axial load transfer. For this purpose, face pairs in several spatial directions above the neutral plane are considered. Pushing these face pairs apart through actuation induces moments in the x, y and xy directions. The actuation concept is influenced by different design parameters such as the control objective, the actuation modes, the adaptation level, the actuator placement and the computation of the required actuation forces.
The design of an adaptive structure is a multi-disciplinary and iterative process in which various factors have to be taken into account. In this work, actuation concepts are not developed as usual by employing stock actuators. The following method can be used to create a list of requirements for the actuator concept. Conversely, the method also allows the requirements of the actuator concept to be taken into account. An optimal adaptive structure, in terms of material and emission efficiency, can be achieved using such an integrated approach whereby the actuation and actuator concept are harmonised with each other.
This work builds on and extends the methods developed in [22] for the integration of fluidic actuators in reinforced concrete structures and [34] concerning the use of actuation influence matrices for the design of the actuation concepts respectively of the individual design parameters. New contributions offered by this work are: 1.
The elaboration of a new actuation concept based on fluidic actuators that can be integrated into floor slabs to control the response under loading.

2.
Extension of actuation influence matrices to two-way slabs.

3.
Computation of actuation forces and design spaces for actuator placement based on influence matrices.

4.
A new methodology that makes it possible to co-design actuation and actuator concepts for adaptive floor slabs.

Influence Matrices
The model adopted for plate structures is based on Reissner-Mindlin Plate Theory with bilinear shape functions on a quadrilateral element with 4 corner nodes and 4 Gauss nodes as shown in Figure 2 and [40][41][42].

Influence Matrices
The model adopted for plate structures is based on Reissner-Mindlin Plate Theory with bilinear shape functions on a quadrilateral element with 4 corner nodes and 4 Gauss nodes as shown in Figure 2 and [40][41][42]. Therefore, the element stiffness matrix has 12 by 12 entries. The order of the entries of the element deformation ( ) and force ( ) vectors are as follows: Assuming a slab with DOF degrees of freedom and i # active elements, the equilibrium and compatibility conditions of an adaptive slab can be written as (cf. [34]): Here and in the following, the superscript # indicates the total number of the respective counting variable. The only exception is the designation DOF, as in this work no counting variable is used for it. K ∊ ℝ DOF x DOF is the stiffness matrix, d ∊ ℝ DOF is the deformation vector, and f ∊ ℝ DOF is the force vector which is separated into forces for the passive state fpas and forces of the active state fact. Assuming small strains and displacements, a distinction can be made between three system states [11,43]: the "passive" state, which is the response of the structure to conventional loads (including self-weight and dead load), the "active" system state, which is the response of the structure caused by the actuation forces, and the "adaptive" state, which is taken as the superposition of passive and active states. Using integrated fluidic actuators, forces are applied through a pair of faces located eccentrically to the cross-section neutral plane. A certain thickness of concrete is assumed to be covering the actuator, in order to achieve a homogeneous stress distribution and corrosion protection for the actuator. Figure 3 shows a single floor slab element of planar dimensions , , , . Indicated in red is the (pressurised) load application area with a height of ℎ . It is assumed that an actuation pushes the opposing faces pairwise apart, resulting in a pair of moments, similarly to what is shown in [22]. At this stage, the actuator concept is not fully defined since some of the parameters are an output of the process. Therefore, only moment pairs and no pressurised faces will be considered in the following. Therefore, the element stiffness matrix has 12 by 12 entries. The order of the entries of the element deformation d (e) and force f (e) vectors are as follows: Assuming a slab with DOF degrees of freedom and i # active elements, the equilibrium and compatibility conditions of an adaptive slab can be written as (cf. [34]): Here and in the following, the superscript # indicates the total number of the respective counting variable. The only exception is the designation DOF, as in this work no counting variable is used for it. K ∈ R DOF × DOF is the stiffness matrix, d ∈ R DOF is the deformation vector, and f ∈ R DOF is the force vector which is separated into forces for the passive state f pas and forces of the active state f act . Assuming small strains and displacements, a distinction can be made between three system states [11,43]: the "passive" state, which is the response of the structure to conventional loads (including self-weight and dead load), the "active" system state, which is the response of the structure caused by the actuation forces, and the "adaptive" state, which is taken as the superposition of passive and active states.
Using integrated fluidic actuators, forces are applied through a pair of faces located eccentrically to the cross-section neutral plane. A certain thickness of concrete is assumed to be covering the actuator, in order to achieve a homogeneous stress distribution and corrosion protection for the actuator. Figure 3 shows a single floor slab element of planar dimensions l a,x , l a,y . Indicated in red is the (pressurised) load application area A a with a height of h a . It is assumed that an actuation pushes the opposing faces pairwise apart, resulting in a pair of moments, similarly to what is shown in [22]. At this stage, the actuator concept is not fully defined since some of the parameters are an output of the process. Therefore, only moment pairs and no pressurised faces will be considered in the following. As long as the force pairs applied on the red faces ( Figure 3) are only short-circuited along the actuator itself, the distance between the pair of moments can be varied in order to counteract the bending moment caused by the external load as needed. In this case, short-circuiting means that the resulting reaction forces cancel each other out through a connecting structure without leading to forces affecting the surrounding region of concrete. Three actuation modes are considered as shown in Figure 4, whereby the biaxial mode is a superposition of the two uniaxial modes. The labels 45° and 135° refer to the axis of the resulting moments with reference to the coordinate system shown in Figure 4. The force vector for the active state depends on the actuation mode employed ( Figure 4). The force vector in the active state can be written as follows (cf. [20]): The vector u contains the actuation input (see also Section 2.3). The matrix A ∊ ℝ # is denoted as the actuator allocation matrix. It collates the forces produced by the unit actuation of each individual element through a particular actuation mode column-wise: For example, the column vector °| corresponds to the uniaxial actuation mode 135° for element 1:  As long as the force pairs applied on the red faces ( Figure 3) are only short-circuited along the actuator itself, the distance between the pair of moments can be varied in order to counteract the bending moment caused by the external load as needed. In this case, short-circuiting means that the resulting reaction forces cancel each other out through a connecting structure without leading to forces affecting the surrounding region of concrete. Three actuation modes are considered as shown in Figure 4, whereby the biaxial mode is a superposition of the two uniaxial modes. The labels 45 • and 135 • refer to the axis of the resulting moments with reference to the coordinate system shown in Figure 4. The force vector for the active state depends on the actuation mode employed ( Figure 4).  As long as the force pairs applied on the red faces ( Figure 3) are only short-circuited along the actuator itself, the distance between the pair of moments can be varied in order to counteract the bending moment caused by the external load as needed. In this case, short-circuiting means that the resulting reaction forces cancel each other out through a connecting structure without leading to forces affecting the surrounding region of concrete. Three actuation modes are considered as shown in Figure 4, whereby the biaxial mode is a superposition of the two uniaxial modes. The labels 45° and 135° refer to the axis of the resulting moments with reference to the coordinate system shown in Figure 4. The force vector for the active state depends on the actuation mode employed ( Figure 4). The force vector in the active state can be written as follows (cf. [20]): The vector u contains the actuation input (see also Section 2.3). The matrix A ∊ ℝ # is denoted as the actuator allocation matrix. It collates the forces produced by the unit actuation of each individual element through a particular actuation mode column-wise: For example, the column vector °| corresponds to the uniaxial actuation mode 135° for element 1:  The force vector in the active state can be written as follows (cf. [20]): The vector u contains the actuation input (see also Section 2.3). The matrix A ∈ R DOF × i # is denoted as the actuator allocation matrix. It collates the forces produced by the unit actuation of each individual element through a particular actuation mode column-wise: For example, the column vector a 135 • |1 corresponds to the uniaxial actuation mode 135 • for element 1: The force vector for the passive state can be generated in the same way, except that the entries of the A matrix must correspond to those of an out-of-plane load. In keeping with the concept of influence surfaces [33], only the active state is considered and f pas is set to zero. Furthermore, only a unit actuation input is applied for all active elements. At this stage, the actuation input vector u can be formulated as the identity matrix of dimensions R i # × i # , which can therefore be omitted [34]. Replacing Equation (3) in (4) yields the equation to compute the actuation influence matrix for displacements: The column vectors e d|i ∈ R DOF quantify the influence on nodal displacements (global coordinates) caused by unit actuation of the associated active element with the actuation mode that is being considered. To observe only the translational displacements w, the influence matrix of displacements E d ∈ R DOF × i # can be reduced to E d(w) ∈ R n # × i # by simply selecting every third global degree of freedom (cf. Equation (1)), n # denotes the total number of element corner nodes. Each column vector e d(w)|i ∈ R n # quantifies the influence of a single active element on the nodal translational displacements w of the whole structure.
To obtain the actuation influence matrices for stress resultants, the strain-displacement matrices B b ∈ R 3 × 12 and B s ∈ R 3 × 12 as well as the stress-strain material matrices C b ∈ R 3 × 3 and C s ∈ R 2 × 2 are required for bending and shear respectively [40,42]. The formulation given in Equations (9)-(12) is for a single element following [40]. The coordinates must be transformed to the isoparametric coordinates of the gauss element according to [40,42].
A gauss node stress resultants vector can be calculated for the moments according to Equation (13) and for the shear forces according to Equation (14).
By iterating through the 4 gauss nodes while using the isoparametric coordinates (current gauss node coordinates) ξ and η in the strain-displacement matrices of Equations (13) and (14), the stress resultant components (m x ; m y ; m xy ; ν x ; ν y ) for the gaussian elements are obtained in turn. An example for m (e ) x and ν (e ) x is given in Equations (15) and (16) [40,42].
After calculating the stress resultants for each gauss element in turn, the values at the element corner nodes are obtained using bilinear extrapolation [42]. To avoid discontinuities between adjacent elements, the nodal stresses are averaged. For the averaging, the same weight is assigned to all elements that meet at a node [42]. Averaging reduces each stress resultant vector to one value per node, i.e., m x|i , m y|i , m xy|i , ν x|i , ν x|i ∈ R n # , where n # denotes the total number of element corner nodes and i the respective active element.
To obtain actuation influence matrices for bending moments E M ∈ R n # × i # or shear forces E V ∈ R n # × i # , the procedure in Equations (13)-(16) (including averaging) must be performed for each column vector e d|i of the influence matrix E d . The stress resultant vectors generated are then collated column-wise into matrices as expressed in Equations (17) and (18).
To gain an understanding of the load-bearing behaviour of two-way slabs, it is useful to analyse the principal moments. With the calculated influence matrices E M x , E M y and E M xy , the influence matrices E m 1 and E m 2 for the two principal moment directions can be calculated. The influence matrices E m 1 and E m 2 for the two principal moments are obtained by successively collating the result of Equation (19) [44] applied to each column vector of

Determining Suitable Actuation Modes
The optimal actuation mode depends on the control objective, e.g., reduction of deformations or homogenisation of bending moments. The selection of suitable actuation modes can be carried out on the basis of the computed influences. Summing the rows of the corresponding influence column vectors for each actuation mode k as indicated in Equations (20) e act|k = e act|k,1 e act|k,2 · · · e act|k,i · · · e act|k,i # −1 e act|k,i # For each actuation mode k, the row vectors e act|k,i ∈ R i # are collated vertically into a matrix E act ∈ R k # × i # . Each row contains the summed influence of actuation mode k while each column entry is the summed influence of an active element for the considered control objective at all element corner nodes, cf. Equation (22).
Using E act , the actuation mode with the highest influence can be determined for each active element. The maximum value per column indicates which actuation mode has the highest influence on the selected control objective at the respective active element (column index). In the following, the index of the actuation mode that has maximum influence at the respective active element i is denoted kmax.
The optimal actuation mode varies for the individual elements. Therefore, a new influence matrix E d,comb is assembled from the influence column vectors corresponding to the elements and actuation modes with the highest summed influence, using E act . For further processing (Section 2.3) and to visualise the selected actuation modes, it is also necessary to generate a new actuator allocation matrix of the combined actuation modes. For this purpose, the columns (active elements) of E d and A corresponding to the actuation mode kmax on the selected control objective, are collated into the combined matrices E d,comb and A comb respectively, as expressed in Equations (23) and (24).

Actuation Load, Adaption Level and Adaptive State
As described in [34], influence matrices can also be used to calculate the actuation forces required for a chosen adaptive nominal state y ada . According to [33], Maxwell's theorem of reciprocal displacements also applies to elastic slabs. For example, consider the deformations in the passive state at a generic point P(x,y) caused by the external load and the deformations generated by a unit load (pair of moments) applied at a point A(x,y). The actuation moment pair required to completely reduce the effect of the external load can be obtained from the ratio between the deformations in the passive and active state at point P. The approach presented in [34] to determine the actuation forces for truss and frame structures can be extended to floor slabs: Often the deformations caused by the permanent load (self-weight + dead load) are compensated by other means (e.g., pre-cambering). In these cases, it is possible to isolate the effect of variable loads (e.g., live load) by excluding the permanent load from the external load in the passive state f pas .
Since influence matrices are typically not square and are therefore not invertible, the Moore-Penrose pseudoinverse, denoted with (·) + , is employed to compute the required actuation input u * . The pseudoinverse is a generalisation of the inverse matrix to singular and non-square matrices [45]. Solving Equation (25) for u * gives: The adaptive system state y * ada achievable with the actuation input u * is then obtained by substituting u * into Equation (25), as expressed in Equation (27). The same applies to the vector of actuation forces. Since the actuation modes considered in this work do not provide out-of-plane forces, there is no corresponding entry in A comb for the translational degree of freedom, and therefore only actuation moments are included in the vector of actuation forces f * act .

Actuator Placement
When all structural elements are actuated, the adaptive state can typically be reached with high accuracy. However, this requires a high number of actuators. To reduce the complexity of the actuation system, it is sensible to reduce the number of actuators as much as possible. The reduction of actuators (i.e., the actuator placement) can be done based on the influence of the active element towards the attainment of the target adaptive state (Equation (20)). An iterative process is carried out by removing each element in turn and evaluating the effect on the target adaptive state.
The elements with the lowest influences are removed first. The minimum number of active elements is selected based on the desired control objective. For example, if vertical displacements are to be compensated, an upper (x u ) and lower (x l ) limit for the translational displacement w of each node is defined. The removal continues for elements with increasing influence until all elements are removed. It is not aborted after the upper or lower limit is reached in order to leave open the possibility (in a further step, which is not dealt with in this work) of analysing and discussing the procedure (cf. Equation (29), Section 4).
Due to the local effect that characterises the actuation with integrated fluidic actuators [22], it is useful to specify a removal direction. For example, if the external load is distributed on the floor slab, low-influence elements should be removed from the outside inwards to the diagonals and towards the centre of the slab. In cases where the external load is only distributed over one of the quadrants, the elements of the loaded quadrant should be removed last.
Elements are removed by filling zero entries into the corresponding column in E d,comb . Each time an element is removed, the adaptive state y * ada,r (Equation (27)) is computed with the updated E d,comb and collated into an output matrix Y * ada ∈ R DOF × r # where r # denotes the total number of removed elements. Y * ada = y * ada,1 y * ada,2 · · · y * ada,r · · · y * ada,r # −1 y * ada,r # Due to the order of the degrees of freedom (Equation (1)), the matrix for the adaptive state for translational displacements Y * ada(w) ∈ R n # × i # can be determined by filtering every third column entry of Y * ada . The minimum number of active elements is equal to the total number of active elements minus the lowest column index from Y * ada(w) in which the target displacement bounds are exceeded. The reduced influence matrix E d,comb,r , which corresponds to the configuration with the minimum number of actuators, is used to compute the adaptive state and the required actuation moments using Equations (26)-(28).

Defining Pressure Levels
As described in Section 2.1, actuation moments are produced by applying pressure to pairs of faces, resulting in an actuation moment (cf. [22], Figure 3). Given f * act from Equation (28), the calculation of the pressure p * act is based on Equation (30) which was formulated and experimentally validated in [22] for beams. A a,x is the load application area in the x direction and z a is the inner lever arm (Figure 3).
The moment curve between a pair of actuation moments is assumed to be constant. As discussed in Section 2.1, this constant moment curve can be induced across multiple elements.
Adjacent elements that are actuated through approximately the same pressure are thought of as forming a pressure level. Consider the example of biaxial actuation shown in Figure 5. Since the actuation moments at adjacent nodes cancel each other out (Figure 5a), the combined elements are equivalent to a larger active element (cf. Figure 5b). This shifts the pairs of faces to the outer edges of the larger active element. This means that the number of active elements does not have to be the same as the actual number of actuators. The active elements are employed to define the active region in which the actuators will be placed. Furthermore, grouping active elements through pressure levels is useful for limiting the necessary number of valves of the hydraulic unit. The number of pressure levels can be chosen, for example, on the basis of the available hydraulic unit. The pressure levels correspond to an even distribution of the pressures over the maximum pressure range. The pressure range corresponds to the maximum required pressure or a selected pressure limit. A pressure limit can be chosen, for example, on the basis of requirements from the actuator concept as well as on the basis of material properties such as the concrete compressive strength. When reducing the number of active elements to obtain the actuator placement region, the remaining active elements are required to apply higher pressure to compensate for the missing elements. In some cases, this can lead to an extraordinary increase in the highest pressure levels. If a significantly high step occurs between individual pressures after sorting them, all pressures beyond the selected step are limited to the nearest smaller pressure level. The moment curve between a pair of actuation moments is assumed to be cons As discussed in Section 2.1, this constant moment curve can be induced across mul elements.
Adjacent elements that are actuated through approximately the same pressur thought of as forming a pressure level. Consider the example of biaxial actuation sh in Figure 5. Since the actuation moments at adjacent nodes cancel each other out (Fi 5a), the combined elements are equivalent to a larger active element (cf. Figure 5b). shifts the pairs of faces to the outer edges of the larger active element. This means tha number of active elements does not have to be the same as the actual number of actua The active elements are employed to define the active region in which the actuators be placed. Furthermore, grouping active elements through pressure levels is usefu limiting the necessary number of valves of the hydraulic unit. The number of pres levels can be chosen, for example, on the basis of the available hydraulic unit. The pres levels correspond to an even distribution of the pressures over the maximum pres range. The pressure range corresponds to the maximum required pressure or a sele pressure limit. A pressure limit can be chosen, for example, on the basis of requirem from the actuator concept as well as on the basis of material properties such as the con compressive strength. When reducing the number of active elements to obtain the a tor placement region, the remaining active elements are required to apply higher pres to compensate for the missing elements. In some cases, this can lead to an extraordi increase in the highest pressure levels. If a significantly high step occurs between ind ual pressures after sorting them, all pressures beyond the selected step are limited t nearest smaller pressure level.

Results
In the following, the application of the method formulated in Section 2 is applie a case study of a simply supported slab. A summary of this procedure can be foun Appendix A. Dimensions and material properties, indicated in Table 1, have been sele on the basis of planned experimental testing and available resources. The control obje is the reduction of displacements. The nodal translational displacements are constra along all edges. The floor slab is discretised in 400 elements as shown in Figure 6.

Results
In the following, the application of the method formulated in Section 2 is applied to a case study of a simply supported slab. A summary of this procedure can be found in Appendix A. Dimensions and material properties, indicated in Table 1, have been selected on the basis of planned experimental testing and available resources. The control objective is the reduction of displacements. The nodal translational displacements are constrained along all edges. The floor slab is discretised in 400 elements as shown in Figure 6. The element size was chosen based on the mesh convergence study shown in Figure 7. A comparison of the numerically determined translational displacements and moments in xand y-direction at slab midspan with the analytically calculated values according to Czerny [46] yields a deviation of around 10%. Due to the relatively small values (cf. Section 3.4), this is considered to be sufficiently accurate for the presentation of the method. The simulations have been performed in Python 3.7.8. element size was chosen based on the mesh convergence study shown in Figure 7. A parison of the numerically determined translational displacements and moments in x y-direction at slab midspan with the analytically calculated values according to Cz [46] yields a deviation of around 10%. Due to the relatively small values (cf. Section this is considered to be sufficiently accurate for the presentation of the method. The ulations have been performed in Python 3.7.8.

Influence Matrices
It is useful to compare the individual plots of the columns of the influence matrices for each actuation mode in the case of the chosen control objective, the columns of the influence matrix ( ) . To compute the resulting displacement response in the passive state, a similar influence matrix can be built by applying a unitary out-of-plane load at every element corner node of a single element. Figure 8 shows an example of three plotted columns for the passive (a-c) and active state (d-f). In this example three application points are chosen: the centre of the slab, the centre of a diagonal and the principal axes related to one quadrant.

Influence Matrices
It is useful to compare the individual plots of the columns of the influence matrices for each actuation mode in the case of the chosen control objective, the columns of the influence matrix E d(w) . To compute the resulting displacement response in the passive state, a similar influence matrix can be built by applying a unitary out-of-plane load at every element corner node of a single element. Figure 8 shows an example of three plotted columns for the passive (a-c) and active state (d-f). In this example three application points are chosen: the centre of the slab, the centre of a diagonal and the principal axes related to one quadrant. In the passive system state, the highest influences are indicated in dark blue (negative deflection) and, in the active state, they are indicated in bright yellow (positive deflection). The response is more pronounced in the immediate vicinity of the load application points or position of the active elements (marked with pink dots and arrows) and then decreases rapidly. The influence extends over a smaller region when load application points (or the position of the active element) are considered closer to the edges. In the passive system state, the highest influences are indicated in dark blue (negative deflection) and, in the active state, they are indicated in bright yellow (positive deflection). The response is more pronounced in the immediate vicinity of the load application points or position of the active elements (marked with pink dots and arrows) and then decreases rapidly. The influence extends over a smaller region when load application points (or the position of the active element) are considered closer to the edges.  Figure 9 shows the plot of the summed influences on translational displacements w under biaxial actuation (a), 45 • (b) and 135 • (c) uniaxial actuation (cf. Section 2.2) for all elements. In red areas, the influence on the control objective is the largest, in green areas it is negligible and in purple areas it is negative, i.e., actuation causes downward displacements.

Combining Actuation Modes
CivilEng 2022, 3 822 Figure 9. Summed influences over all elements on the displacements w of a biaxial actuation (a), 45° (b) as well as 135° (c) uniaxial actuation. Figure 9 shows that the biaxial actuation mode is most efficient in the centre of the slab but has almost no effect in the corners. The uniaxial actuation modes affect the deformations along the diagonals and mostly towards the corners, according to their alignment with the principal moment m1. It is important to align uniaxial actuation correctly with respect to the principal moment m1. The principal moments change sign at the corners, so actuation has a negative effect on the displacements w. Figure 10 shows the actuation modes selected after applying Equations (20)- (24).

Actuation Load and Adaptation Level
To obtain measurable displacements under the considered dimensions and material properties (Table 1), a distributed load of 15 kN/m² is applied. The active moments required to compensate displacements w as well as rotations and are computed with Equations (25)- (27) for a distribution of actuation modes according to Figure 10. The required actuation pressures are indicated in Figure 11. If all the coloured elements in  Figure 9 shows that the biaxial actuation mode is most efficient in the centre of the slab but has almost no effect in the corners. The uniaxial actuation modes affect the deformations along the diagonals and mostly towards the corners, according to their alignment with the principal moment m 1 . It is important to align uniaxial actuation correctly with respect to the principal moment m 1 . The principal moments change sign at the corners, so actuation has a negative effect on the displacements w. Figure 10 shows the actuation modes selected after applying Equations (20)- (24).  Figure 9 shows that the biaxial actuation mode is most efficient in the centre of the slab but has almost no effect in the corners. The uniaxial actuation modes affect the deformations along the diagonals and mostly towards the corners, according to their alignment with the principal moment m1. It is important to align uniaxial actuation correctly with respect to the principal moment m1. The principal moments change sign at the corners, so actuation has a negative effect on the displacements w. Figure 10 shows the actuation modes selected after applying Equations (20)- (24).

Actuation Load and Adaptation Level
To obtain measurable displacements under the considered dimensions and material properties (Table 1), a distributed load of 15 kN/m² is applied. The active moments required to compensate displacements w as well as rotations and are computed with Equations (25)-(27) for a distribution of actuation modes according to Figure 10. The required actuation pressures are indicated in Figure 11. If all the coloured elements in

Actuation Load and Adaptation Level
To obtain measurable displacements under the considered dimensions and material properties (Table 1), a distributed load of 15 kN/m 2 is applied. The active moments required to compensate displacements w as well as rotations ϕ x and ϕ y are computed with Equations (25)- (27) for a distribution of actuation modes according to Figure 10. The required actuation pressures are indicated in Figure 11. If all the coloured elements in Figure 10 are active, the maximum required pressure is approximately 76 bar. The required pressure is the highest in the proximity of the sign change of the principal moment m 1 . The pressure is computed with Equation (30). The chosen values A a,x = 0.003 m 2 (0.1 m × 0.03 m) and z a = 0.015 m correspond to the maximum size per element to obtain a concrete cover of 2 cm.

Number of Necessary Active Elements
If a permissible range is defined for the displacements w, not all active elements are required. For the case study under consideration, the application of the procedure described in Section 2.4 leads to a required number of active elements of 232 (out of a total of 400). Figure 12 shows a zoomed-out view of the matrix ( ) * ∊ ℝ # # . The individual entries that lie within the permissible deformation range of = 1 × 10 −5 to =−1 × 10 −5 are highlighted in red. The representation is structured in such a way that one element is removed per column from left to right. From top to bottom, the node number is listed in ascending order. The set permissible range can no longer be complied with when an active element is removed from the area where the principal moment m1 is close to zero.

Number of Necessary Active Elements
If a permissible range is defined for the displacements w, not all active elements are required. For the case study under consideration, the application of the procedure described in Section 2.4 leads to a required number of active elements of 232 (out of a total of 400). Figure 12 shows a zoomed-out view of the matrix Y * ada(w) ∈ R n # × i # . The individual entries that lie within the permissible deformation range of x u = 1 × 10 −5 to x l =−1 × 10 −5 are highlighted in red. The representation is structured in such a way that one element is removed per column from left to right. From top to bottom, the node number is listed in ascending order. The set permissible range can no longer be complied with when an active element is removed from the area where the principal moment m 1 is close to zero.
The maximum value of the displacement in the passive state w pas is 3.4 × 10 −4 m. When all elements are active, the maximum displacement in the adaptive state w ada is 0.032 × 10 −4 m (a reduction of 99.9%). After removing 168 active elements ( Figure 13) the maximum displacement in the adaptive state is only 0.074 × 10 −4 m (a reduction of 97.8%), due to the strict chosen displacement bounds.

Defining Pressure Levels
A reduction in the number of active elements can result in a large increase in the required actuation pressure. Figure 14a shows the required actuation pressures to satisfy the displacement bounds set in Section 3.4. Compared to the case in which all elements are active (Figure 11), the maximum pressure increases from 76 bar to 164 bar. To further reduce the number of actuators, elements with a similar pressure are combined through pressure levels (Section 2.5). In addition, since the pressure only increases significantly for some of the elements, a pressure limit is introduced. Six pressure levels are defined, since the hydraulic unit that will be used for testing has six valves. For the general pressure levels, this results in steps of roughly 10 bar. The pressure level of the few active elements with extreme values is manually limited to the highest pressure level. The actuation pressure plot for this case is shown in Figure 14b. The maximum displacement for this variation is 0.42 × 10 −4 m (w pas = 3.4 × 10 −4 m, reduction of 87.6%). are active (Figure 11), the maximum pressure increases from 76 bar to 164 bar. To further reduce the number of actuators, elements with a similar pressure are combined through pressure levels (Section 2.5). In addition, since the pressure only increases significantly for some of the elements, a pressure limit is introduced. Six pressure levels are defined, since the hydraulic unit that will be used for testing has six valves. For the general pressure levels, this results in steps of roughly 10 bar. The pressure level of the few active elements with extreme values is manually limited to the highest pressure level. The actuation pressure plot for this case is shown in Figure 14b. The maximum displacement for this variation is 0.42 × 10 −4 m (wpas = 3.4 × 10 −4 m, reduction of 87,6%). Figure 14. (a) Actuation pressures after reduction of the active elements. (b) Actuation pressures using six pressure levels, including pressure limit.

Discussion
The concept of integrated fluidic actuation has already been validated on a largescale beam in [18]. In this study, the underlying approach of counteracting moments caused by external loads with opposing active moments has been extended to floor slabs. The directionality of the load transfer in two-way slabs increases the complexity. Influence matrices are employed to determine the location of the active elements and the required actuation forces to satisfy a chosen control objective.
The simplification of the 2-dimensional plate theory does not allow an evaluation of the stress in thickness direction since the stress curve from pressurised eccentric faces does not correspond to the linear stress curve of the plate theory. However, the knowledge gained regarding the orientation and intensity of the actuation moments, as well as the possibility of identifying actuator design spaces (pattern from pressure levels) with reduced computation time, justifies this approach. Stresses in the thickness direction can be evaluated separately, since individual variations still have to be checked during the development of the actuator concept.
In general, practical implementation will have to be tested to appreciate the feasibility of the method proposed in this paper. This includes the necessary proof that the distance between the pairs of moments can be increased by an appropriate actuator structure.

Discussion
The concept of integrated fluidic actuation has already been validated on a large-scale beam in [18]. In this study, the underlying approach of counteracting moments caused by external loads with opposing active moments has been extended to floor slabs. The directionality of the load transfer in two-way slabs increases the complexity. Influence matrices are employed to determine the location of the active elements and the required actuation forces to satisfy a chosen control objective.
The simplification of the 2-dimensional plate theory does not allow an evaluation of the stress in thickness direction since the stress curve from pressurised eccentric faces does not correspond to the linear stress curve of the plate theory. However, the knowledge gained regarding the orientation and intensity of the actuation moments, as well as the possibility of identifying actuator design spaces (pattern from pressure levels) with reduced computation time, justifies this approach. Stresses in the thickness direction can be evaluated separately, since individual variations still have to be checked during the development of the actuator concept.
In general, practical implementation will have to be tested to appreciate the feasibility of the method proposed in this paper. This includes the necessary proof that the distance between the pairs of moments can be increased by an appropriate actuator structure. The exact implementation (1 actuator = 1 active element or 1 pressure level field) would mean that there is no concrete between the individual pressurised faces. In addition, the actuator requires a construction space (enclosure, connecting elements, etc.) where no force can be applied on the concrete. This issue can be solved by allowing gaps between the active elements (or pressure fields). The applied pressure is distributed in plane, resulting in actuation moments that are not just directly at the pressurised faces. The distances between individual actuators can be selected based on the decay lengths of the applied in-plane stresses (cf. [22]). These adjustments lead to a smaller load application area A a . In practical implementation, the maximum pressure required for displacement compensation will therefore be higher, with expected values of about 100 bar. In the case study presented in this work, the use of pressure barriers has a larger impact on the deviation from the desired displacement bounds than the reduction of active elements. This means that if a bigger area (with more combined active elements) is actuated, smaller inaccuracies can also be compensated for.
The effect of the actuation modes on actuator placement and determination of the minimum number of active elements for different control objectives will be the subject of future work. A review of Y * ada(w) shows that there are still some active elements that can be removed without affecting compliance with the deformation criteria. Further work will also analyse and improve the removal process.
Although slabs are typically dimensioned based on deformation limits and very good results for deformation reduction at relatively low pressures have been demonstrated in this study, it remains to be proven that the basic goal behind the adaptive slab, i.e., savings in mass and resources, can be achieved. Especially since no statement can be made about the mass of the actuators at this point in the design phase.

Conclusions
This study presents an extension of the use of influence matrices from trusses and beams to two-way slabs with integrated fluidic actuators. The actuation influence matrices are employed to place the actuators as well as to determine preliminary requirements for the actuator concept. It has been shown that the characteristic value of the summed influences is a useful indicator to choose suitable actuation modes. Through the determination of pressure levels, optimal position and forces of the actuators can be obtained for a chosen control objective including the distance between the pair of faces (design space) of the fluidic actuators. Through the load case independent analysis, it was shown that biaxial actuation in the centre of the slab and uniaxial actuation along the diagonals is advantageous. Furthermore, the paper illustrates that actuation in the corner areas is not necessary for displacement compensation of a uniformly distributed load. This highlights that adaptive slabs require specifically designed actuators instead of industry standard ones to reduce the number of actuators and to achieve a good performance. Moreover, it was presented that with a gradation of just six pressure levels, displacements from a uniformly distributed load can be reduced by more than 85%. Future work will test the feasibility of the presented method in this paper and general practical applicability of adaptive slabs through experimental testing.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Conflicts of Interest:
The authors declare no conflict of interest.

Glossary
Sub-and Superscripts: # indicates the total number of the respective counting variable * indicates that the vector or matrix refers to the required target state DOF degrees of freedom (e) single element, evaluation at the element corner nodes (e') single "gauss" element, evaluation at the gauss nodes i counting variable for the active elements i = 1, 2, . . . , i # k counting variable for the actuation modes kmax index of the actuation mode that has maximum influence at the respective active element n counting variable for the element corner nodes n = 1, 2, . . . , n # r counting variable for the number of removed elements r = 1, 2, . . . , r # Vectors: a i ∈R DOF column vector of the actuator allocation matrix. Including the actuation mode, e.g., corresponds to the uniaxial actuation mode, e.g., a 135 • |1 for element 1 d (e) . ∈ R 12 element deformation vector d ∈ R DOF global deformation vector e act|k ∈R i # row vector of the summed rows of the corresponding influence column vectors of each element and actuation mode on the control objective e d|i ∈ R DOF column vector of the actuation influence matrix for displacements e d(w)|i ∈R n # column vector of the actuation influence matrix for translational displacements w e m 1,2 ∈R n # column vector of actuation influence matrices for the two principal moment directions f ∈ R DOF global force vector f (e) ∈ R 12 element force vector f * act ∈ R DOF force vector of the required actuation moments f * act ∈ R DOF global force vector for the active state f pas ∈ R DOF global force vector for the passive state m (e') ∈ R 4 gauss node stress resultant vector for the bending moments p * act ∈ R DOF actuation pressure vector u * ∈ R DOF required actuation input vector, for a target state u ∈ R DOF actuation input vector ν (e') ∈ R 4 gauss node stress resultant vector for the shear forces y act ∈ R DOF output vector for the active nominal state y * ada ∈ R DOF output vector for the adaptive system state achievable with the actuation input u * ∈ R DOF output vector for the adaptive system state achievable with the actuation input u * and the current number of active elements y ada ∈ R DOF output vector for the adaptive nominal state y pas ∈ R DOF output vector for the passive nominal state Matrices: A ∈R DOF x i # actuator allocation matrix A comb ∈R DOF x i # actuator allocation matrix of the combined actuation modes B b ∈ R 3 x 12 strain-displacement matrices for bending B s ∈ R 3 x 12 strain-displacement matrices for shear C b ∈ R 3 x 3 stress-strain material matrices for bending C s ∈ R 2 x 2 stress-strain material matrices for shear E act ∈R k # x i # actuation influence matrix for the summed influences of each element and actuation mode on the control objective E d ∈R DOF x i # actuation influence matrix for displacements E d,comb ∈R DOF x i # actuation influence matrix for displacements of the combined actuation modes with the highest summed influence E d,comb,r ∈R DOF x i # actuation influence matrix for displacements of the combined actuation modes with the highest summed influence and the minimum number of active elements to stay within the target state bounds E d(w) ∈ R n # x i # actuation influence matrix for translational displacements w E M ∈ R n # x i # actuation influence matrix for bending moments E m 1 and E m 2 ∈R n # x i # actuation influence matrices for the two principal moment directions E V ∈ R n # x i # actuation influence matrix for shear forces K ∈ R DOF x DOF global stiffness matrix k counting variable for the actuation modes Y * ada ∈R DOF x i # output matrix for the adaptive system state achievable with the actuation input u * and the current number of active elements Y * ada(w) ∈R n # x i # output matrix for the translational displacements in the adaptive system state achievable with the actuation input u * and the current number of active elements Symbols: (·) + Moore-Penrose pseudoinverse ϕ x and ϕ y rotational degrees of freedom A a load application area e act|k,i summed influence of a single element on the control objective for the respective actuation mode e act|n,i influence of a single element on the control objective at the respective element corner node h a height of the load application area l a,x and l a,y edge length of a single slab element in x and y direction. Applies also when using pressure levels to form a single larger active element. l x and l y span width of the slab in x and y direction w translational displacement x u upper limit for the translational displacement w of each node x l lower limit for the translational displacement w of each node z a inner lever

Appendix A
A summary of the procedure described in Section 2.1-2.5 is given in Figure A1 as a flow chart. First of all, the influence matrices for each actuation mode need to be calculated. Summing the rows of the influence column vectors yields a characteristic value to determine suitable actuation modes. After assigning a load case and setting an adaptive target state, the adaptive system state can be calculated. By removing active elements with the lowest influence in turn the minimum number of active elements can be determined. Assigning pressure levels and setting a pressure limit leads to possible distances for the actuation faces.
flow chart. First of all, the influence matrices for each actuation mode need to be calculated. Summing the rows of the influence column vectors yields a characteristic value to determine suitable actuation modes. After assigning a load case and setting an adaptive target state, the adaptive system state can be calculated. By removing active elements with the lowest influence in turn the minimum number of active elements can be determined. Assigning pressure levels and setting a pressure limit leads to possible distances for the actuation faces. Figure A1. Summary of the procedure described in Sections 2.1-2.5. Figure A1. Summary of the procedure described in Section 2.1, Section 2.2, Section 2.3, Section 2.4, Section 2.5.