Mathematical Problems in Engineering
Volume 2018 (2018), Article ID 3989250, 9 pages
https://doi.org/10.1155/2018/3989250

Research Article

Reissner Plates with Plastic Behavior: Probability of Failure

Ernesto Pineda-León,¹ Alejandro Rodríguez-Castellanos,2,3 Dante Tolentino,¹ José Manuel Rosales-Juárez,³ Ivan Felix-González,² and Supriyono⁴

¹Escuela Superior de Ingeniería y Arquitectura, SEPI-ESIA, Instituto Politécnico Nacional, Unidad Profesional Adolfo López Mateos, s/n, Ciudad de México, Mexico
²Instituto Mexicano del Petróleo, Eje Central Lázaro Cárdenas 152, 07730 Gustavo A. Madero, Ciudad de México, Mexico
³Escuela Superior de Ingeniería Mecánica y Eléctrica, SEPI-ESIME, Instituto Politécnico Nacional, Unidad Profesional Adolfo López Mateos, s/n, Ciudad de México, Mexico
⁴Department of Mechanical Engineering, Universitas Muhammadiyah Surakarta, Surakarta, Indonesia

Correspondence should be addressed to Alejandro Rodríguez-Castellanos; arcastel@imp.mx

Received 9 February 2018; Revised 11 April 2018; Accepted 16 April 2018; Published 24 May 2018

Academic Editor: Anna M. Gil-Lafuente

Copyright © 2018 Ernesto Pineda-León et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The current paper shows the application of the boundary element method for the analysis of plates under shear stress causing plasticity. In this case, the shear deformation of a plate is considered by means of Reissner’s theory. The probability of failure of a Reissner’s plate due to a proposed index plastic behavior is calculated taking into account the uncertainty in mechanical and geometrical properties. The problem is developed in three dimensions. The classic plasticity’s theory is applied and a formulation for initial stresses that lead to the boundary integral equations due to plasticity is also used. For the plasticity calculation, the von Misses criterion is used. To solve the nonlinear equations, an incremental method is employed. The results show a relatively small failure probability () for the ranges of loads between . However, for values between , the probability of failure increases significantly. Consequently, for , the plate failure is imminent. The results are compared to those that were found in the literature and the agreement is good.

1. Introduction

All the materials have an elastic limit state; when the load is above this limit state, the material gets into the plastic deformation. This plastic strain depends not only on the final state of stresses but also on the load history. The plastic strain is nonlinear and is irreversible; hence, the physical relations that define the process are more complicated than those that define the elastic behavior.

Plates are commonly used in civil and mechanical engineering. They are used to fabricate vehicles, airplanes, equipment, slabs, and so forth. The main references with information about the elastic theory with strain due to shear on plates are due to Reissner [1], Mindlin [2], and Timoshenko and Goodier [3]. On the other hand, plasticity theory has been reported by Mendelson [4], Hill [5], and Prager [6].

Additionally, the study of plates has been diversified and more realistic conditions have been studied for its analysis and design. Thermal effects, geometric nonlinearity, explosive loads, plasticity, creep, and probabilistic behavior among others are the kind of nonlinear problems that have been considered for plate analysis.

Recently, several formulations based on Reissner-Mindlin plate theory have been developed for the study of foundation plates considering thermal effects [7]. These authors also reported analytical solutions for plates subjected to bending and free vibration. Additional studies on the behavior of foundation plates based on the Reissner-Mindlin theory can be found in [8], Wen and Aliabadi [9], and Jankovski and Skarzauskas [10]. In this last reference, the boundary element method was formulated in the Laplace domain to solve the problem. Also, Nosier and Fallah [11] studied the effect of mechanical and thermal loads on the axisymmetric and asymmetric behavior of circular plates, where material properties, boundary conditions, plate thickness, and other factors were studied and discussed.

On the other hand, Reissner-Mindlin theory has been applied to the study of nonlinear problems of large deformations on plates. With this idea, Xiao-Yan as well as Mao-kuang and Xiuxi [13] developed an interesting application of the plate boundary element method with geometric nonlinearity in isotropic plates, including the effects of shear deformation. Moreover, recently, Benatta et al. [14] studied plates subjected to pressure loads, where geometric nonlinearity was incorporated into the strain-displacement equations based on the Von-Karman theory. For this case, the results suggest that the precision obtained by this formulation is comparable with high-order theories.

Moreover, Reissner-Mindlin theory has been successfully applied to problems dealing with cracks, composites, and effects of concentrated and uniform loads on plates [15–17], including the study of oceanic plates [18]. In the latter case, the authors numerically solved the flexion equations associated with the Reissner-Mindlin theory and the application was focused on the study of the bending of the Nazca oceanic plate. This fact reveals the power of this theory and how it can be applied to highly complex cases.

Also, there are different factors that modify the response of a structural element; such factors are associated with its mechanical and geometrical properties; that is, there is a probability that the manufacturer will provide a structural element with both geometric and mechanical characteristics different from those specified, which generates uncertainty, since the structural element is designed based on manufacturer's specifications. The loads to which the structural element will be subjected are variable in magnitude, occurrence, direction, and so on. This leads us to consider as stochastic variables all the parameters that intervene in the response of the structural element.

The application of BEM to probabilistic analysis has been successfully applied to solve elastostatic problems [19–22], crack problems [23, 24], and elastic-linear problems [25–27]. Moreover, the BEM method was used for solving optimization problems considering the cost as a parameter of design [28]. Unfortunately, none of the mentioned approaches using BEM have been used to obtain the probability of failure based on a plastic behavior index that considers the relationship between demand and capacity of the structural element. In the present work, the BEM method is used in Reissner plates with plastic behavior and the probability of failure is obtained considering the mechanical and geometrical uncertainties.

In Section 2, the boundary element method for Reissner plates is formulated. The validation of the formulation is included in Section 3. In Section 4, a failure probability approach is given to the analysis of Reissner plates and finally Section 5 includes an application example.

2. Formulation of the Plastic Problem

Our formulation is based on Supriyono and Aliabadi [29] and is closely related to Dirgantara and Aliabadi [30]. Furthermore, our research intends to extend the formulation of Supriyono and Aliabadi [29] including the probability of failure of Reissner plates with plastic behavior.

Assuming the reference system of Figure 1(a) and the definition of bending and shear displacements as (where ), displacements in the plane as () and their associated tractions can be seen in Figures 1(b) and 1(c); for this case, the strain-displacement relations can be written as

where is transverse shear strain; is bending deformation; and is strain in the middle plane of the element.

Figure 1: (a) Definition of the reference system. (b) Displacements and tractions due to bending and shear. (c) Displacements and tractions in the plane.

On one hand, if we consider a differential element of a plate with dimensions subject to a uniform load (see Figure 2), the equilibrium equations for moments and forces, ignoring higher-order terms, can be expressed as

where results from bending stresses, results from shear stresses, and results from membrane stresses.

(a)

(b)

Figure 2: Load’s equilibrium on an element.

The notation used in Figure 2 means that is used to denote + ; the notation represents + .

On the other hand, Reissner [1] obtained the stress-strain relationships for flexural plates, which are

Based on Hooke’s law for plane stresses, the resulting stress-strain relationships for membrane behavior can be expressed as follows:

In terms of displacements, (7)–(9) are written as

where tension stiffness flexural stiffness, and shear stiffness. The Kronecker delta is represented by .

2.1. Integral Boundary Formulation for Reissner Plates

According to the equations of equilibrium (4)–(6), it is possible to express the following integral equation:

where and are weight functions and is the domain where integration takes place. Equation (13) is the integral representation associated with the governing equation for the results of bending and shear stress, while (14) is associated with the governing equation for the resultant membrane stresses.

2.2. Focus on Initial Stresses

In this section, we introduce the basic equation of elastoplastic problems in continuous media. Plasticity is a phenomenon independent of time and admits increases in load as a variation from one state to another. For example, the stress for the load step () can be represented as

In general terms, the nonlinearity of the material for flexural plates is considered by the incorporation of plastic deformations due to flexion and membrane [12]. Thus, strain ratios can be generically defined as

where the superscripts “” and “” represent the elastic and plastic parts, respectively.

On the other hand, the results of stresses can be established as

Equations (19), (20), and (21) represent the resulting stresses due to the momentum and membrane moment, respectively. The stress-displacement relations can be represented by the phenomenon of plasticity as

2.3. Integral Equations of Elastoplastic Strain

Here, we formulate the integral displacement equations for Reissner plates in materials with nonlinear behavior (i.e., plastic behavior). We begin by considering the equations of equilibrium (4)–(6) to have the following relations:

where and are the weight functions and is the domain where the integration takes place. The above equations (25) and (26) together with Betti’s reciprocal theorem are used to obtain the integral equations of displacement.

2.4. Integral Representation of Rotation and Deflection

The integral representing the ratio of rotations () and the deflection () can be obtained by integrating (25) by parts and applying Green’s identity to obtain

where and are the field points on the boundary and the domain, respectively, and , , and are the fundamental solutions of displacements, tractions, and strains, respectively.

2.5. Integral Representations of In-Plane Displacements

The integral representations related to the governing equations for the resulting membrane stresses can be obtained by integrating (26) by parts, and applying Green’s identity gives

where , , and are the fundamental solutions of displacements, tractions, and strains, respectively.

2.6. Integral Boundary Equations of Displacement

In order to obtain solutions for points on the boundary, it is necessary to consider the limit in such a way that the boundary equation of the displacements can be written as

where represents the principal value of Cauchy, are the source and field points, respectively, on the boundary, and and are the jump terms. The values of and are equal to and , respectively, when is located on a smooth boundary.

2.7. Integral Equations of Elastoplastic Stresses

The results of elastoplastic stresses on the points of the domain can be evaluated by differential equations with respect to the source point . Finally, the integral equation for elastoplastic stresses can be written as

For the result of stress due to moments,

resulting from shear forces, and

For the result of membrane stresses where the kernels are a linear combination of the first derivative of and , the kernels , are a linear combination of the first derivative of and . The kernels , are a linear combination of the first derivative of and .

We used BIE for in-plane stress/strain because it is considered that the plastic zone is not only because of bending but also because of in-plane load. On the other hand, the relation between the boundary integral equations for plate bending and plane strain is included in (31), (32), and (33) of the paper, where the last integral of each equations involves the nonlinear strain. Equation (31) represents stresses due to moments that cause plate bending, (32) denotes the stresses resulting from shear forces, and (33) characterizes the membrane stresses.

The above equations represent the integral boundary formulation for Reissner plates with large deformation and plastic behavior.

3. Verification of the Method

In order to verify the application and precision of the method presented in this paper, the problem solved by Karam and Telles [12] has been selected. This consists of a square plate (simply supported plate at its boundary) with dimensions such that the ratio , where is the plate thickness and is the plate width; the plate is subjected to a uniform load in the vertical direction. The plate properties are  GPa, , and  MPa. For these problems, the material is considered with elasto-perfectly plastic behavior (no hardening is considered). In Figure 3 we use the nondimensional parameters and , where and . The variable represents the uniform load applied to the plate. Figure 3 shows the results obtained by the present formulation and those obtained by Karam and Telles [12]. Good agreement between both formulations is observed. This figure shows the dimensionless deflection calculated at the center of the plate for several values of load, which reaches a maximum value of 24.97. The results clearly illustrate a linear behavior and subsequently the plastic behavior of the plate.

Figure 3: Deflection in the center of a square plate with uniform load. The symbols represent the results of Karam and Telles [12] and the continuous line represents those obtained by the present formulation.

4. Probability of Failure in a Structural Element

A structural element or system is designed to work in service conditions (elastic range) below the load for which it was designed (i.e., own weight, dead load, live load, earthquake, wind, and waves) for a certain interval of time. If an extraordinary load occurs before or after said time interval, there is a high probability that the structural element incurs in the inelastic range and the failure could occur. According to the above, it is important to consider such loads within the analysis in order to estimate the probability of failure of the structural element. One of the advantages of estimating the probability of failure is that it allows the redesign of the structural element so that it implicitly considers all possible magnitudes of loads, where a structural element would be safety. The disadvantage is that such element would be robust and costly.

The probability of failure of a structural element can be obtained by evaluating both the load solicitations (demand) and the ultimate capacity. The relationship between demand and capacity expressed in terms of a damage indicator (i.e., stress and deformation) is called damage index. The damage index has the advantage of normalizing the damage indicator with values between 0 and 1 and it permits quantifying the probability of no damage, the probability of damage, and the probability of failure. For the purpose of establishing a damage index, several authors have proposed damage indexes for structures subjected to earthquakes [31–36]. With the aim to estimate the probability of failure in Reissner’s plates, a plastic behavior index that relates the relationship between the requested dimensionless displacement and the ultimate displacement is proposed as follows:

where is the requested dimensionless displacement; is the dimensionless displacement, where the plastic behavior of the structural element begins; is the ultimate capacity value of the structural element. According to the above, values of indicate that the structural element has not entered the plastic range and values of indicate the failure of the structural element.

Based on the plastic behavior index, , it is possible to quantify empirically the probability of failure by a sample of Monte Carlo simulation as the sum of failures, , between the number of cases, ; that is, . Therefore, the mean of the natural logarithm, , and the standard deviation of the natural logarithm, , are obtained as follows:

If failures of the structural element follow a lognormal distribution [37], they can be fitted theoretically as follows:

where is the normal distribution function; is the standard normal probability distribution function.

5. Methodology

The probability of failure of the plate due to loads is studied by means of a Monte Carlo simulation. The methodology is described in Figure 4.

Figure 4: Flowchart for obtaining the probability of failure.

6. Application Example

In this section, the failure probability of a plate is evaluated considering the uncertainties in its mechanical (yielding, ultimate stress, and modulus of elasticity) and geometrical (thickness, length, and width) properties. The nominal properties of the plate are the following: 0.1 m long by 0.1 m wide, with the thickness being 0.05 m. The yield stress is 248.108 MPa, the ultimate stress is 400.111 MPa, and the modulus of elasticity is 200055.66 MPa. Table 1 shows the statistical parameters reported by Mansour et al. [38] and Galambos and Ravindra [39].

Table 1: Statistical parameters used in the present study.

Variable
Mean value
Coefficient of variation (CV)
Type of distribution

1.05
0.044
Lognormal

0.988
0.046
Lognormal

0.987
0.076
Lognormal

1.3
0.124
Normal

1.05
0.075
Normal

The plastic behavior index and the probability of failure are obtained considering the procedure described in the previous section. 55 cases were simulated considering the geometric and mechanical uncertainties. It was verified that a larger number of samples did not change the statistical results. Figure 5 shows the plastic behavior index (see (34)) considering the simulated cases. It is noted that the variability of the curves of the plastic behavior index is due to the fact of considering the uncertainties in mechanical and geometric properties. According to the above, different values can be observed, where , that is, different failure scenarios of the structural element. For , values are between 1.26 and 2.38, which correspond to 0.023 and 0.039 m, respectively. This range is similar to the displacement failure values reported by Rudrapatna et al. [40], where the authors studied a plate with similar geometric and mechanical characteristics using the finite element method and validated the results by experimentally testing the plate up to failure.

Figure 5: Plastic behavior index for the 20 considered cases.

Figure 6 shows the probability of empirical and theoretical failure. It is noticed that the points of the empirical failure probability plot follow a lognormal distribution, so the assumption made in (37) is verified. Figure 6 displays that the probability of failure for is null. However, for values of , failure probability increases as increases: and , respectively. This range of failure probability can occur when the uncertainty in mechanical and geometric properties is considered. The failure event is imminent for , since empirical and theoretical probabilities of failure are the same: .

Figure 6: Probability of empirical and theoretical failure.

7. Conclusions

In the present work, the application of the boundary element method is used to solve the nonlinear problem of plasticity in Reissner’s plates. We use the classical theory of plasticity and a formulation of initial stresses, which allow the formulation of integral boundary equations due to plasticity. For the calculation of the plastic zone, the von Mises criterion is used. To verify the validity and accuracy of our formulation, the results were compared with those obtained by Karam and Telles [12]. The plastic behavior index proposed in this work to quantify the probability of failure has been validated against published data.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors are grateful for the support granted by CONACYT under the Project CB-2014/240413 and they acknowledge Instituto Politécnico Nacional for the support.

References

E. Reissner, “On a Variational Theorem in Elasticity,” Journal of Mathematics and Physics, vol. 29, pp. 90–95, 1950.

R. D. Mindlin, “Influence of Rotatory Inertia and Shear on Flexural Motions of Isotropic, Elastic Plates,” Journal of Applied Mechanics, vol. 18, pp. 31–38, 1951.

S. Timoshenko and J. N. Goodier, Theory of Elasticity, McGraw-Hill International Editions, Springer, Berlin, Germany, 3rd edition, 1970.

A. Mendelson, Plasticity: Theory and Application, Macmillan, New york, NY, USA, 1968.

R. Hill, The Mathematical Theory of Plasticity, Clarendon Press, Oxford, UK, 1950.

W. Prager, An Introduction to Plasticity, Addison-Wesley, Amsterdam, Netherlands, 1959.

H.-S. Shen, J. Yang, and L. Zhang, “Dynamic response of Reissner-Mindlin plates under thermomechanical loading and resting on elastic foundations,” Journal of Sound and Vibration, vol. 232, no. 2, pp. 309–329, 2000.

P. H. Wen, “The fundamental solution of Mindlin plates resting on an elastic foundation in the Laplace domain and its applications,” International Journal of Solids and Structures, vol. 45, no. 3-4, pp. 1032–1050, 2008.

P. H. Wen and M. H. Aliabadi, “Boundary element formulations for Mindlin plate on an elastic foundation with dynamic load,” Engineering Analysis with Boundary Elements, vol. 33, no. 10, pp. 1161–1170, 2009.

V. Jankovski and V. Skaržauskas, “The physically nonlinear analysis of circular plate on deformable foundation,” The Baltic Journal of Road and Bridge Engineering, vol. 6, no. 1, pp. 59–66, 2011.

A. Nosier and F. Fallah, “Reformulation of Mindlin-Reissner governing equations of functionally graded circular plates,” Acta Mechanica, vol. 198, no. 3-4, pp. 209–233, 2008.

V. J. Karam and J. C. F. Telles, “Nonlinear material analysis of Reissners plates, Plate Bending Analysis with Boundary Element, Computational Mechanics Publications,” Nonlinear material analysis of Reissners plates, Plate Bending Analysis with Boundary Element, Computational Mechanics Publications, 1998.

X.-Y. Lei, M.-K. Huang, and X. Wang, “Geometrically nonlinear analysis of a Reissner type plate by the boundary element method,” Computers & Structures, vol. 37, no. 6, pp. 911–916, 1990.

M. A. Benatta, A. Kaci, A. Tounsi, M. S. Houari, K. Bakhti, and E. A. Adda Bedia, “Nonlinear bending analysis of functionally graded plates under pressure loads using a four variable refined plate theory,” International Journal of Computational Methods, vol. 11, no. 4, Article ID 1350062, 1350062, 16 pages, 2014.

F. G. Yuan and S. Yang, “Asymptotic crack-tip fields in an anisotropic plate subjected to bending, twisting moments and transverse shear loads,” Composites Science and Technology, vol. 60, no. 12-13, pp. 2489–2502, 2000.

J. G. Simmonds, “Pointwise errors in the classical and in Reissner's linear theory of plates, especially for concentrated loads,” Journal of Elasticity: The Physical and Mathematical Science of Solids, vol. 23, no. 2-3, pp. 219–232, 1990.

W. Yu and D. H. Hodges, “A simple thermopiezoelastic model for smart composite plates with accurate stress recovery,” Smart Materials and Structures, vol. 13, no. 4, pp. 926–938, 2004.

P. Manŕiquez, E. Contreras-Reyes, and A. Osses, “Lithospheric 3-D flexure modelling of the oceanic plate seaward of the trench using variable elastic thickness,” Geophysical Journal International, vol. 196, no. 2, pp. 681–693, 2013.

M. Ettouney, H. Benaroya, and J. Wright, “Boundary element methods in probabilistic structural analysis (PBEM),” Applied Mathematical Modelling, vol. 13, no. 7, pp. 432–441, 1989.

I. Kaljević and S. Saigal, “Stochastic boundary elements in elastostatics,” Computer Methods Applied Mechanics and Engineering, vol. 109, no. 3-4, pp. 259–280, 1993.

M. Kaminski, “Stochastic second-order BEM perturbation formulation,” Engineering Analysis with Boundary Elements, vol. 23, no. 2, pp. 123–129, 1999.

M. Kamiński, “Application of the generalized perturbation-based stochastic boundary element method to the elastostatics,” Engineering Analysis with Boundary Elements, vol. 31, no. 6, pp. 514–527, 2007.

X. Huang and M. H. Aliabadi, “Probabilistic fracture mechanics by the boundary element method,” International Journal of Fracture, vol. 171, no. 1, pp. 51–64, 2011.

G. Chen, S. Rahman, and Y. H. Park, “Shape sensitivity and reliability analyses of linear-elastic cracked structures,” International Journal of Fracture, vol. 112, no. 3, pp. 223–246, 2001.

C. Su and C. Zheng, “Probabilistic fracture mechanics analysis of linear-elastic cracked structures by spline fictitious boundary element method,” Engineering Analysis with Boundary Elements, vol. 36, no. 12, pp. 1828–1837, 2012.

C. Su, S. Zhao, and H. Ma, “Reliability analysis of plane elasticity problems by stochastic spline fictitious boundary element method,” Engineering Analysis with Boundary Elements, vol. 36, no. 2, pp. 118–124, 2012.

C. Su and J. Xu, “Reliability analysis of Reissner plate bending problems by stochastic spline fictitious boundary element method,” Engineering Analysis with Boundary Elements, vol. 51, pp. 37–43, 2015.

D. W. Kim and B. M. Kwak, “Reliability-based shape optimization of two-dimensional elastic problems using BEM,” Computers & Structures, vol. 60, no. 5, pp. 743–750, 1996.

Supriyono and M. H. Aliabadi, “Boundary element method for shear deformable plates with combined geometric and material nonlinearities,” Engineering Analysis with Boundary Elements, vol. 30, no. 1, pp. 31–42, 2006.

T. Dirgantara and M. H. Aliabadi, “Elastoplastic boundary element method for shear deformable shells,” Engineering Structures, vol. 45, pp. 62–67, 2012.

H. Banon and D. Veneziano, “Seismic safety of reinforced concrete members and structures,” Earthquake Engineering & Structural Dynamics, vol. 10, no. 2, pp. 179–193, 1982.

M. S. L. Roufaiel and C. Meyer, “Analytical modeling of hysteretic behavior of R/C frames,” Journal of Structural Engineering (United States), vol. 113, no. 3, pp. 429–444, 1987.

E. Cosenza, G. Manfredi, and R. Ramasco, “The use of damage functionals in earthquake engineering: A comparison between different methods,” Earthquake Engineering & Structural Dynamics, vol. 22, no. 10, pp. 855–868, 1993.

J. M. Bracci, A. M. Reinhorn, and J. B. Mander, “Deterministic model for seismic damage evaluation of reinforced concrete structures,” Tech. Rep., State University of New York at Buffalo, 1989.

Y.-J. Park and A. H.-S. Ang, “Mechanistic seismic damage model for reinforced concrete,” Journal of Structural Engineering, vol. 111, no. 4, pp. 722–739, 1985.

G. S. Kamaris, G. D. Hatzigeorgiou, and D. E. Beskos, “A new damage index for plane steel frames exhibiting strength and stiffness degradation under seismic motion,” Engineering Structures, vol. 46, pp. 727–736, 2013.

E. Rosenblueth and L. Esteva, “Reliability basis for some Mexican codes,” in in Probabilistic Design of Reinforced Concrete Buildings, ACI Publication SP-31, American Concrete Institute, p. 41, MI, Farmington Hills, 1972.

A. E. Mansour, H. Y. Jan, C. I. Zigelman, Y. N. Chen, and S. J. Harding, “Implementation of reliability methods to marine structures,” Transactions Society of Naval Architects and Marine Engineers, vol. 92, pp. 353–382, 1984.

T. V. Galambos and T. V. Ranvindra, Properties of steel for use in LRFD Journal of the Structural Division104, pp. 1459–1468, 1978.

N. S. Rudrapatna, R. Vaziri, and M. D. Olson, “Deformation and failure of blast-loaded square plates,” International Journal of Impact Engineering, vol. 22, no. 4, pp. 449–467, 1999.