Stability of Composite Plates with a Dense System of Ribs in Two Directions

TL;DR

This paper introduces a fast analytical method for analyzing the stability of composite plates with dense ribs, avoiding resource-heavy simulations.

Contribution

A novel analytical approach using tolerance averaging simplifies stability analysis of ribbed composite plates with discontinuous coefficients.

Findings

The analytical method produces results matching FEM calculations but with much lower computational cost.

The transformed equations resemble classical plate equations, enabling straightforward analytical solutions.

Material and geometrical properties significantly influence critical forces under shear and compression.

Abstract

What are the main findings?
A simple analytical method of stability analysis of a thin plate with a dense system of ribs is derived and verified.Issues of uniaxial compression and shear in-plane loadings are investigated.All presented results are verified with FEM-based calculations, proving the correctness of the derived solution.

A simple analytical method of stability analysis of a thin plate with a dense system of ribs is derived and verified.

Issues of uniaxial compression and shear in-plane loadings are investigated.

All presented results are verified with FEM-based calculations, proving the correctness of the derived solution.

What are the implications of the main findings?
The derived analytical approach allows swift and reliable stability analysis of composites, contrary to FEM-based investigations, where the modelling of microstructure requires a lot of computational resources and long computation times.The performed investigations of the influence of material and geometrical properties on the shear stability of the plate allows us to predict critical forces without the necessity of performing time-consuming calculations.

The derived analytical approach allows swift and reliable stability analysis of composites, contrary to FEM-based investigations, where the modelling of microstructure requires a lot of computational resources and long computation times.

The performed investigations of the influence of material and geometrical properties on the shear stability of the plate allows us to predict critical forces without the necessity of performing time-consuming calculations.

This paper presents an easy-to-use analytical method for stability analysis of composite plates with dense bidirectional microstructure. The main characteristic feature of such a defined composite is that due to its periodic nature the obtainable governing partial differential equations are characterised by discontinuous, strongly oscillating coefficients. Such cases bring many difficulties during derivation of their solution. In order to simplify calculations, the initial governing equations are transformed with the use of the tolerance averaging technique, so a system of partial differential equations with constant coefficients is obtained. The most important finding of the presented work is that the form of the mentioned equations is similar to the classic equations, which describe the stability issue of the thin homogeneous plate. Consequently, the analytical solution to such issues is easily obtainable. Moreover, when compared to, for example, finite element method (FEM) analysis, it requires substantially less computation resources, which can be perceived as its superior feature. Therefore, the proposed method is convenient for engineering applications. In this paper, a comparative analysis of the results obtained from the proposed analytical models with the results obtained from the FEM has been carried out. The impact of materials and dimensions of microstructure on the values of critical normal and shear forces has also been analysed.

Figures

Full text

1. Introduction

1.1. Research Significance

Non-homogeneous plates have many applications in modern engineering issues, such as aircraft construction, space vehicles and even in traditional civil engineering. The most common ones are sandwich panels. A review of computational methods for such panels can be found in Hu et al. [1]. Another type of composite plates are plates with internal microstructure, which will be analysed in this paper.

There are many papers in the literature considering similar structures, such as stiffened, fibre-reinforced or functionally graded plates. The work of Deng et al. [2] was dedicated to the buckling analysis of unidirectionally stiffened panels. Teter and Kolakowski [3] and Piculin and Može [4] examined the stability of unidirectionally reinforced panels with thin ribs. The distances between the ribs in these works are large enough to effectively use the finite element method for analysis. The stability analysis of fibre-reinforced plates can be found in the works of Taheri and Ebrahimi [5], where the Mori–Tanaka method was used in the homogenization procedure and in the work [6] where a cantilevered multilayer plate of Marco fibre and carbon fibre composites was analysed. Additional issues connected with, e.g., delamination and manufacturing of carbon-fibre-reinforced polymers can be found in [7,8,9]. The literature concerning plates made of functionally graded materials (FGM) is full of various investigations of interesting cases, such as nonlinear vibrations [10,11], thermal buckling [12], dynamic stability [13], dynamics of FGM shells [14] or flexural behaviour of sandwich plates with FGM core [15].

Possible methods to analyse such plates are, among others, the homogenisation technique or finite element method. The latter is a labour-intensive process due to the necessity of detailing the element mesh to the microstructure of the composite. Tolerance averaging technique, presented in this work, is a special modelling procedure similar to homogenisation. It gives the opportunity to replace strongly oscillating discontinuous coefficients of differential equations with constant ones. Consequently, an easy-to-use method for the evaluation of critical forces of the composite is obtainable. Its superior feature is that it requires substantially less computation resources, than, for example, finite element method (FEM) analysis. Thus, the proposed method is convenient for engineering applications.

The aim of the presented work is to obtain and apply a mathematical macroscopic model describing stability of the thin plate with a dense system of ribs in two directions based on the tolerance averaging technique. General rules of modelling within this technique can be found in the literature, cf. Woźniak et al. [16,17]. Applications of this method for the analysis of various dynamic and thermomechanical problems can be found in multiple works. In the paper presented by Matysiak and Nagórko [18], propagation of harmonic waves in a periodically laminated composite is analysed. The analysis of dynamic behaviour of honeycomb-based composite solids can be performed according to the work of Wierzbicki and Woźniak [19]. In the paper [20], modelling of the dynamics of thin plates made of functionally graded material (FGM) is considered. The obtained averaged equations of motion are verified by comparing the results from the tolerance model with those from the finite element method. The dynamics of three-layer plates (sandwich plates) is analysed in the work of Marczak and Jędrysiak [21], while Domagalski et al. [22] analyse the vibrations of a periodic Timoshenko beam with nonlinearity of the von Karman type. Jędrysiak and Ostrowski [23] describe natural vibrations of a plate made of material of uncertain properties. The description of thermomechanical issues may be found in [24]. Ostrowski and Michalak [25] study an FGM conductor, where the heat conduction in a microheterogeneous cylinder is verified using the finite element method.

The application of the tolerance averaging technique to model stability issues of various periodic composites can be found in a series of articles. The stability of lightly wrinkled plates is analysed by Michalak [26]. Stability of thin periodically stiffened cylindrical shells is investigated in the work [27]. Marczak and Jędrysiak [28] study the stability of periodic beams interacting with elastic foundation, while the vibration and stability of viscoelastic periodic beams based on a periodic foundation with damping is considered in [29]. What is more, the application of the tolerance averaging technique to obtain a macroscopic model describing the stability of plates made of FGM is shown in [30] and is used in the paper [31], where the stability issue of a functionally variable annular plate cooperating with an elastic microheterogenous subsoil is described.

It should be noted that in the abovementioned works the cell dimension has to be much larger than the thickness of the structure. The fundamental distinctive feature of the presented work is that the considered plate is reinforced by a bidirectional rib system, which microstructure length parameter l has a dimension similar to the plate thickness h ( $[eqn]$ ).

1.2. Purpose and Scope of Work

In this work, composite plates with a dense system of ribs are analysed, cf. Figure 1. The main characteristic feature of this contribution is that its periodicity cell dimension of the microstructure is comparable to the thickness of the plate. The main purpose of this work is to obtain and apply a macroscopic mathematical model describing the stability of the plate under normal and shear loadings. A tolerance-averaging technique is used to build the model, which does not contain equations with discontinuous and strongly oscillating coefficients. Eventually, the results obtained from the tolerance-averaging technique are compared with the ones derived from finite element method (FEM). The analysis of the impact of materials and dimensions on critical shear forces is described as well.

2. Methods of Calculation

2.1. Description of the Considered Composite

Let us analyse a thin plate in a Cartesian coordinate system $[eqn]$ with characteristic dimensions $[eqn]$ and $[eqn]$ along the $[eqn]$ and $[eqn]$ axis, respectively, cf. Figure 1. The considered object is a thin composite with constant thickness h and a dense bidirectional system of ribs. Throughout the article, ribs should be understood as those parts of the thin plate which are made of material having superior properties, when compared to the rest of the composite (the so-called matrix).

Let us assume that it is possible to distinguish two families of thin homogeneous ribs with widths $[eqn]$ and $[eqn]$ , respectively, which axes intersect at right a angle, cf. Figure 1. The areas between the ribs are filled with homogeneous matrix material. Additionally, let us state that it is possible to distinguish m ribs along the $[eqn]$ axis and $[eqn]$ ribs along the $[eqn]$ axis in the whole structure, where $[eqn]$ Those ribs have a constant width within a specific direction and are distributed uniformly throughout the whole plate. The microstructure parameter $[eqn]$ (l1 and l2 are periodicity cell dimensions) is assumed to be sufficiently small when compared to the dimensions of the plate midplane, but at the same time it is comparable to the plate thickness h ( $[eqn]$ ). A macroscopic model of dynamics of such a type of plate can be found in the work of Michalak and Rabenda [32].

Let us introduce the concept of a basic periodicity cell; $[eqn]$ as the smallest repeatable part of the plate, where $[eqn]$ and $[eqn]$ are the cell’s dimensions. Moreover, it is assumed that the whole plate is made of isotropic materials with a constant value of Poisson’s ratio. Setting x ≡ (x1, x2) and $[eqn]$ , the undeformed plate occupies the region $[eqn]$ , where $[eqn]$ is the plate’s midplane and $[eqn]$ is the constant plate’s thickness.

2.2. Formulation of the Preliminary Governing Equations

Let us start with the direct description of the considered plate in the framework of the well-known second order non-linear theory of thin plates. The displacement field of the arbitrary point of the plate is assumed in the form

$$[eqn]$$

where $[eqn]$ , $[eqn]$ , are displacements along the $[eqn]$ -axis directions and $[eqn]$ , $[eqn]$ are in-plane displacements caused by in-plane loadings.

Denoting by $[eqn]$ the metric tensor, the Ricci tensor by $[eqn]$ , setting $[eqn]$ and the gradient operators $[eqn]$ , where $[eqn]$ and $[eqn]$ are first- and second-order gradients, a system of governing equations of linear approximated theory for thin plates (cf. [33]) is obtained:

• (i)strain-displacement relations:

$$[eqn]$$

(ii)strain energy averaged over the plate thickness:

$$[eqn]$$

where $[eqn]$ is the tensile stiffness and $[eqn]$ with $[eqn]$ is the bending stiffness.

(iii)work of external forces:

$$[eqn]$$

where $[eqn]$ , $[eqn]$ , are external loadings applied to the midplane of the considered composite along the $[eqn]$ -axis, respectively.

In order to derive governing equations of the considered plate we shall define the stationary action functional:

$$[eqn]$$

where Lagrangean $[eqn]$ . From stationary action principle ( $[eqn]$ ) we obtain governing equations:

$$[eqn]$$

where generalised forces are defined as:

$$[eqn]$$

It can be observed that this direct description leads to a system of partial differential equations with discontinuous and highly oscillating coefficients, which are difficult to use in multiple issues of modern engineering. The above equations are used as a preliminary formulation for the tolerance averaging modelling procedure.

2.3. Basic of the Tolerance Modelling

In order to derive averaged equations for the composite plate under consideration we applied the tolerance averaging approach [16,17]. This technique is based on the concept of indiscernibility (cf. [34,35]) and tolerance relations. Let us mention here some basic concepts of this approach, such as an averaging operator, a tolerance periodic function, a slowly varying function and a highly oscillating function.

The fundamental concept of the modelling technique is the averaging operation:

$$[eqn]$$

We shall refer to Formula (8) as an averaging of arbitrary integrable function $[eqn]$ for every $[eqn]$ .

The main concept of the tolerance averaging approach is that values of functions defined on a region $[eqn]$ can be determined only with a certain accuracy $[eqn]$ . Let $[eqn]$ stand for an arbitrary positive number and let $[eqn]$ be a linear normed space. Tolerance relation $[eqn]$ for a certain $[eqn]$ is defined by:

$$[eqn]$$

where $[eqn]$ is called the tolerance parameter. Such a relation can be treated as the already mentioned indiscernibility relation between $[eqn]$ and $[eqn]$ .

Periodic approximation. Let $[eqn]$ be the Sobolev space for fixed $[eqn]$ . Function $[eqn]$ , $[eqn]$ , $[eqn]$ will be referred to as the periodic approximation of $[eqn]$ in $[eqn]$ (where $[eqn]$ - $[eqn]$ -th gradient in $[eqn]$ ). For $[eqn]$ we define $[eqn]$ , $[eqn]$ .

Tolerance periodic function. Function $[eqn]$ is called the tolerance periodic function (with respect to cell $[eqn]$ and tolerance parameter $[eqn]$ ), $[eqn]$ , if for $[eqn]$ the following conditions are met:

$$[eqn]$$

Slowly varying function. Function $[eqn]$ will be called the slowly varying function (with respect to the cell $[eqn]$ and tolerance parameter $[eqn]$ ), and denoted by $[eqn]$ , if for $[eqn]$ the following conditions are met:

$$[eqn]$$

It can be observed that periodic approximation $[eqn]$ of $[eqn]$ in $[eqn]$ is a constant function for every $[eqn]$ . In other words, if $[eqn]$ , then:

$$[eqn]$$

Highly oscillating function. Function $[eqn]$ is called the highly oscillating function (with respect to the cell $[eqn]$ and tolerance parameter $[eqn]$ ), and denoted by $[eqn]$ , if for $[eqn]$ the following conditions are held:

$$[eqn]$$

which leads to the following conclusion:

$$[eqn]$$

Let us by $[eqn]$ denote a highly oscillating function, $[eqn]$ , defined on $[eqn]$ , continuous together with gradient $[eqn]$ . Its second derivative, $[eqn]$ is piecewise continuous and bounded. Function $[eqn]$ is called the fluctuation shape function of the second kind, if it depends on $[eqn]$ as a parameter and satisfies conditions:

$$[eqn]$$

The set of all fluctuation shape functions of the second kind is denoted by $[eqn]$ .

2.4. Tolerance Averaging Technique in the Modelling of Microperiodic Plate

In this section the averaged equations of thin plate are presented. Since the derivation of such equations require a specific use of the tolerance averaging technique, those investigations can be treated as an original content of this work.

The starting point for the derivation of the averaged governing equations for the plate under consideration are Equations (1)–(7). The modelling technique is based on the tolerance averaging approximation and on the micro-macro decomposition of the displacement fields given by:

$$[eqn]$$

for $[eqn]$ , $[eqn]$ and $[eqn]$ . The basic tolerance modelling assumption states that macro-displacements $[eqn]$ and amplitudes of displacements fluctuation $[eqn]$ are slowly varying functions together with all partial derivatives. Functions $[eqn]$ are the basic unknowns. Functions $[eqn]$ are known fluctuation shape functions, which are dependent on microstructure length parameter $[eqn]$ . These functions describe the shape of displacements fluctuations in cell $[eqn]$ caused by the microstructure of the plate.

Let $[eqn]$ be a periodic approximation of $[eqn]$ in cell $[eqn]$ , respectively. Due to the fact that $[eqn]$ are tolerance-periodic functions, it can be stated that the periodic approximation of $[eqn]$ and their derivatives in $[eqn]$ , $[eqn]$ can be defined as:

$$[eqn]$$

for every $[eqn]$ and almost every $[eqn]$ .

By setting $[eqn]$ and $[eqn]$ into Lagrangian $[eqn]$ we can assume that $[eqn]$ . Hence the periodic approximation of $[eqn]$ in every $[eqn]$ we denote by $[eqn]$ . In order to derive the governing equations we shall define tolerance-averaged Lagrangian $[eqn]$ :

$$[eqn]$$

Substituting the right-hand sides of Equation (15) into (3), with use of the tolerance averaging technique, we arrive at the following formula for the strain energy averaged over the cell $[eqn]$ :

$$[eqn]$$

The averaged energy of the external loadings over the cell $[eqn]$ is defined as:

$$[eqn]$$

From the principle of stationary action of the averaged Lagrangean $[eqn]$ we derive equations, which describe

(i)plane stress state:

$$[eqn]$$

where averaged normal forces are equal to

$$[eqn]$$

(ii)bending state:

$$[eqn]$$ $$[eqn]$$

where

$$[eqn]$$

We can observe that from Equation (22b), direct representation of oscillation amplitudes $[eqn]$ can be obtained. Let $[eqn]$ stand for a linear transformation operator such that $[eqn]$ . Thus amplitudes of oscillation are defined as

$$[eqn]$$

Eventually, by denoting

$$[eqn]$$

we obtain a governing equation for stability of the plate:

$$[eqn]$$

The above equation has an identical form as the stability equation for thin orthotropic plate. Hence, it is solvable with well-known methods of structural mechanics, which is very useful in practice. Let us once again emphasise that, contrary to classic description, the coefficients in Equation (25) are functional but smooth coefficients.

3. Calculation Example and Results

3.1. Common Assumptions

In order to verify the correctness of the proposed tolerance model some benchmark analysis should be made. Let us consider a rectangular composite plate with constant width of the ribs and simply supported on all edges.

In the calculation example, the stability of a square plate is investigated under either uniaxial compression or shear stress. In all further calculations it is assumed that the matrix is made of concrete having Young’s modulus $[eqn]$ and Poisson’s ratio $[eqn]$ , while the ribs are made of steel for which $[eqn]$ and $[eqn]$ . Additionally, let us state that the considered plate is a biperiodic square plate with characteristic dimensions $[eqn]$ , thickness $[eqn]$ and size of the cell $[eqn]$ . The considered structure is analysed with the use of both presented approaches: tolerance averaging technique and FEM.

3.2. Initial Remarks on the Used Tolerance Model

During the whole tolerance modelling procedure, there are several assumptions which have to be made. The important one is the form of fluctuation shape functions. Our further analysis is restricted to the case where it is decided to express the microstructural behaviour of the composite with two shape functions. Both of them are assumed as a certain product, formulated as follows:

$$[eqn]$$

where $[eqn]$ .

(a)Fluctuation shape function $[eqn]$ (cf. Figure 2) is defined as

$$[eqn]$$ $$[eqn]$$

where $[eqn]$ .

(b)Fluctuation shape function $[eqn]$ (cf. Figure 3) is defined as

$$[eqn]$$ $$[eqn]$$

where $[eqn]$ .

In the case of composite plate with constant width of the ribs the stability Equation (25) of the tolerance model transforms into

$$[eqn]$$

where for square plate and $[eqn]$ and $[eqn]$ we have

$$[eqn]$$

In the above formulas it is assumed that Poisson’s ratio $[eqn]$ and $[eqn]$ , $[eqn]$ as the stiffness of the matrix and ribs, respectively. In such form, the derived tolerance model of the composite is ready to be used in further calculations.

3.3. Finite Element Modelling of the Composite

In order to assess the quality of results derived within the averaged model of a thin plate with a dense system of ribs, the FEM formulation of the same structure is performed, cf. Figure 4. The whole composite is created as a two-dimensional shell plate with a thickness of 0.10 m. Several different geometries are investigated, which differ from each other with the width of ribs $[eqn]$ and $[eqn]$ . Obviously, different sets of material properties (elastic moduli and Poisson’s ratios) are applied to both ribs and matrix in every calculation case.

In all calculations it is assumed that the considered plates are rectangular and simply supported on all for edges, which within FEM modelling should be understood as no perpendicular displacements on the whole boundary of the composite. Such a case is ideal for the assessment of quality of the averaged model, as its solution is easily obtainable using mathematical methods.

Let us make several remarks concerning the quality of FEM analysis. As has been stated, the whole structure is modelled with shell elements (S4R and S3) with thickness equal to 0.10 m. Multiple mesh sizes were considered in order to provide reliable results with reasonable computing times. Eventually, all investigated structures were modelled with the use of approximately 100,000–180,000 elements, depending on the ribs’ dimensions. Within such approach the convergence of results is assured.

3.4. Calculation Example I—Stability Under Uniaxial Compression

In this section the stability of biperiodic rectangular plate compressed in one direction $[eqn]$ is investigated, cf. Figure 5.

In such a case, the stability Equation (31) is reduced to the form

$$[eqn]$$

where $[eqn]$ . The above issue is solved in a way which is similar to the known method for a simply supported rectangular homogeneous plate. Let use the following form of solution:

$$[eqn]$$

where $[eqn]$ , $[eqn]$ and $[eqn]$ are coefficients of the sine series corresponding to certain buckling shapes defined by parameters $[eqn]$ , $[eqn]$ . In the next step, after substituting (34) into Equation (33) we obtain

$$[eqn]$$

where $[eqn]$ . The critical force for the m-th and n-th buckling mode is equal to

$$[eqn]$$

It is clear that for $[eqn]$ we obtain the lowest value of the critical force, which is related to the first mode of buckling.

The results of the tolerance model (TM), together with its benchmark results of FEM for the following widths of ribs: $[eqn]$ , are shown in Table 1. We can also notice that the values of critical forces for $[eqn]$ are identical to those for a classic homogeneous plate made of concrete with the appropriate value of the elastic modulus.

The results presented in Table 1 require special comment. It can be noticed that the relative error increases with the width of the ribs. It is caused by the growing influence of microheterogeneities in the calculations of critical force. The presented modelling procedure of TM clearly indicates that within this solution the existence of ribs with certain width is averaged on the whole periodicity cell. Such an approach allows us to simplify the initial governing equations so that they are easily solvable.

On the other hand, these transformations should be always perceived as a kind of approximation, the consequences of which become clear when compared to the explicit modelling of the rib in FEM, cf. Figure 6. In reality, it can be observed that the ribs (made of materials having superior properties when compared to the matrix) carry most of the in-plane loadings. Obviously, the averaging operation in the case of such disproportions must lead to noticeable discrepancies, which become larger with higher widths of ribs. Consequently, the wider the ribs, the lesser precision of the averaged model, which is noticeable in the results presented in Table 1.

Nevertheless, it should be emphasised that for the considered biperiodic plate with a microstructure, the differences in the results between TM and the FEM model are smaller than 10% when $[eqn]$ .

3.5. Calculation Example II—Shear Stability

In the current calculations, shear stability of a square plate is analysed. The critical value of the uniform shear load on a simply supported rectangular plate (Figure 7) is determined. Material properties and geometry of the composite are assumed as in Section 3.1. Considering only the strain energy of bending $[eqn]$ and the work of the shear load $[eqn]$ , the critical value of the load is determined from the condition

$$[eqn]$$

The strain energy of the considered microstructural plate is equal to

$$[eqn]$$

During the deformation of the plate, the work of the shear load can be expressed as

$$[eqn]$$

The solution to the stated issue is assumed in the form

$$[eqn]$$

where $[eqn]$ , $[eqn]$ . Substituting (40) into Formulas (39) and (38), respectively, we obtain

$$[eqn]$$ $$[eqn]$$

where $[eqn]$ and $[eqn]$ must be odd numbers. For the even values of $[eqn]$ or $[eqn]$ the work of the shear loading $[eqn]$ , and consequently no buckling modes are obtainable. From the energy equilibrium condition (37) the expression for the critical forces for the m-th and n-th buckling mode is obtained:

$$[eqn]$$

It can be noticed that it is no longer true that the lowest values of the critical force are obtained for the lowest values of parameters m, p, n, q. In such a case, the smallest critical value of the load is obtained by levelling the derivatives of each factor $[eqn]$ to zero. Subsequently, a system of linear equations for the unknowns $[eqn]$ are received. By limiting the values of the indicators $[eqn]$ and $[eqn]$ , this system of equations has the form

$$[eqn]$$

where $[eqn]$ and

$$[eqn]$$

Eventually, a system of eigenvalues of the system of Equation (44) are evaluated, which corresponds to the critical value of the shear force:

$$[eqn]$$

The tolerance model (TM) is verified for the following widths of ribs: $[eqn]$ . Results of both TM and FEM analysis are shown in Table 2.

Results presented in Table 2 have similar characteristic features as those presented in Table 1. It is clear that the precision of the TM is decreasing in the case of structures with higher disproportions in material properties, which in this calculation are represented by composites with wider ribs. Similarly, as in the case of uniaxial compression, the differences in the results of critical shear forces between TM and the FEM models are smaller than 10% when $[eqn]$ .

3.6. Calculation Example III—Influence of Material and Geometrical Properties on Shear Stability of the Plate

In this section, the influence of certain geometry and material properties on the shear stability is investigated. Let us remember that all calculations are performed for a simply supported square plate with dimensions $[eqn]$ and thickness $[eqn]$ . Additionally, all considered cases are characterised by ribs made of steel $[eqn]$ .

Firstly, the influence of the coefficient $[eqn]$ (d—rib’s width, l—cell’s dimension) on the critical value of shear forces is evaluated. The range of the assumed ratio is assumed as $[eqn]$ . The calculations are made for two different matrix materials: typical concrete with Young’s modulus of $[eqn]$ and light concrete with $[eqn]$ . Figure 8 shows the critical values of shear forces versus ratio $[eqn]$ .

The nature of the graphs in Figure 8 is almost linear (slightly exponential). It can be seen that the ratio $[eqn]$ increases as the parameter $[eqn]$ increases. On the other hand, the difference between the values $[eqn]$ is almost constant. It can be proved that the values of the critical forces for $[eqn]$ derived in TM are identical to those for a homogeneous plate made of concrete with the appropriate value of the $[eqn]$ modulus.

The influence of the matrix material on the values of the critical forces is shown in Figure 9. The calculations were performed for three values of the proportion of the rib material in the cell: $[eqn]$ . The matrix material varies from $[eqn]$ to $[eqn]$ .

In Figure 9, the linear dependence of the critical shear force value on the matrix elasticity modulus is easily noticeable. This relation is different for different $[eqn]$ . Table 3 shows explicit relations between the modulus of elasticity of the matrix, parameter $[eqn]$ and critical shear force values in considered cases.

Figure 10 presents the ratio of the critical force of a composite plate to the critical force of a homogeneous concrete plate $[eqn]$ depending on the share of the rib material in the cell ( $[eqn]$ ). These relations are evaluated for three ratios of the matrix stiffness to the rib stiffness $[eqn]$ , which are obtained by altering elastic modulus of the matrix ( $[eqn]$ ). Functions: $[eqn]$ for $[eqn]$ , $[eqn]$ for $[eqn]$ and $[eqn]$ for $[eqn]$ are drawn in Figure 10.

It can be seen that for a smaller proportion of the matrix stiffness in the cell, an increase in the width of the ribs causes a faster increase in the value $[eqn]$ . Such an observation is compliant with our engineering intuition, as the introduction of strong ribs into a relatively weak plate (low $[eqn]$ ratio) should result in a rapid increase in the critical force of the considered composite. For higher ratios of $[eqn]$ this reinforcement is not so meaningful, as the stronger matrix itself is capable of withstanding higher in-plane loadings without buckling.

4. Conclusions and Future Research

Let us present several of the most important conclusions which can be drawn from the presented work.

The problem of stability in the thin plate with a dense system of microstructures is described by partial differential equations (PDE) with highly oscillating and discontinuous coefficients. Therefore, the tolerance averaging technique is applied to obtain averaged PDEs with functional but smooth coefficients.The obtained averaged PDEs have the same structure as the equations which describe the stability of orthotropic homogeneous plates. Hence, the solution of specific boundary problems of stability of certain composite plates can be obtained using simple mathematical methods. Such convenience in the modelling of complicated composites is superior to other known methods.The verification of the averaged equations of the tolerance model was carried out using the finite element method. The verification process of the model equations passed satisfactorily. In the case of plates with the share of rib material in the cell $[eqn]$ , the differences between the values of the critical forces obtained from both methods did not exceed 15%. Much more similar results of TM and FEM are obtained for lesser values of $[eqn]$ . Consequently, it can be stated that the presented approach can be applied in multiple engineering issues.For plates which are simply supported on all edges, the relation between the shear critical force and the modulus of elasticity of the matrix $[eqn]$ is linear. This makes it possible to easily determine the value of the critical force for different values of $[eqn]$ without the need to solve differential equations. Consequently, the presented solution allows us to save a great deal of effort and computational resources.The influence of the rib width $[eqn]$ on the critical shear force is particularly advantageous in plates with a low $[eqn]$ rate. In such a case, an introduction of a dense system of ribs allows us to significantly reinforce the whole composite. This is particularly important when the considered plates are expected to have good thermal insulation properties, which usually are reserved for materials with lesser density and poor elastic modulus $[eqn]$ .

The presented analysis can be significantly expanded by multiple issues: the stability analysis of microperiodic plates situated on microheterogeneous foundations, the analysis of FGM plates with a dense system of ribs or plates having different shapes (e.g., circular). All of those issues can also be investigated with the use of experiments, which should be considered as the ultimate verification of the presented analysis.