Abstract
Kirchhoff–Love plate theory is widely used in structural engineering. In this paper, efficient and accurate numerical algorithms are developed to solve a generalized Kirchhoff–Love plate model subject to three common physical boundary conditions: (i) clamped; (ii) simply supported; and (iii) free. The generalization stems from the inclusion of additional physics to the classical Kirchhoff–Love model that accounts for bending only. We solve the model equation by discretizing the spatial derivatives using second-order finite-difference schemes, and then advancing the semi-discrete problem in time with either an explicit predictor–corrector or an implicit Newmark-Beta time-stepping algorithm. Stability analysis is conducted for the schemes, and the results are used to determine stable time steps in practice. A series of carefully chosen test problems are solved to demonstrate the properties and applications of our numerical approaches. The numerical results confirm the stability and 2nd-order accuracy of the algorithms and are also comparable with experiments for similar thin plates. As an application, we illustrate a strategy to identify the natural frequencies of a plate using our numerical methods in conjunction with a fast Fourier transformation power spectrum analysis of the computed data. Then we take advantage of one of the computed natural frequencies to simulate the interesting physical phenomena known as resonance and beat for a generalized Kirchhoff–Love plate.
Similar content being viewed by others
References
Niiranen J, Kiendl J, Niemi AH, Reali A (2017) Isogeometric analysis for sixth-order boundary value problems of gradient-elastic Kirchhoff plates. Comput Method Appl Mech Eng 316:328–348
Niiranen J, Niemi AH (2017) Variational formulations and general boundary conditions for sixth-order boundary value problems of gradient-elastic Kirchhoff plates. Eur J Mech Solids 61:164–179
Flügge W (2013) Viscoelasticity. Springer, Berlin
Vlasov VZ (1966) Beams, plates and shells on elastic foundation. Israel Program for Scientific Translation
Koiter WT, Simmonds JG (1973) Foundations of shell theory. In: Becker E, Mikhailov GK (eds) Theoretical and applied mechanics. Springer, Berlin, pp 150–176
Leissa A (1969) Vibration of plates. Tech Rep NASA-SP-160, NASA
Love AEH (1888) The small free vibrations and deformation of a thin elastic shell. Philos Trans R Soc A 179:491–546
Reissner E (1976) On the theory of transverse bending of elastic plates. Int J Solids Struct 12:545–554
Timoshenko S, Woinowsky-Krieger S (1959) Theory of plates and shells, 2nd edn. McGraw-Hill, New York
Mindlin RD (1951) Influence of rotatory inertia and shear on flexural motions of isotropic, elastic plates. J Appl Mech 18(1):31–38
Koiter WT (1960) A consistent first approximation in the general theory of thin elastic shells. In: Proceedings of the IUTAM symposium on the theory of thin elastic shells, pp. 12–33. North–Holland, Amsterdam
Canic S, Tambaca J, Guidoboni G, Mikelic A, Hartley CJ, Rosenstrauch D (2006) Modeling viscoelastic behavior of arterial walls and their interaction with pulsatile blood flow. SIAM J Appl Math 67:164–193
Arnold DN, Madureira AL, Zhang S (2002) On the range of applicability of the Reissner–Mindlin and Kirchhoff–Love plate bending models. J Elast 67:171–185
Banks JW, Henshaw WD, Schwendeman DW (2014) An analysis of a new stable partitioned algorithm for FSI problems. Part II: incompressible flow and structural shells. J Comput Phys 268:399–416
Li L, Henshaw WD, Banks JW, Schwendeman DW, Main GA (2016) A stable partitioned FSI algorithm for incompressible flow and deforming beams. J Comput Phys 312:272–306
Szilard R (2004) Theories and applications of plate analysis: classical numerical and engineering methods. Wiley, Chichester
Bilbao S (2008) A family of conservative finite difference schemes for the dynamical von Karman plate equations. Numer Methods Partial Differ Equ 24(1):193–216
Ji H, Li L (2019) Numerical methods for thermally stressed shallow shell equations. J Comput Appl Math 362:626–652
Batoz JL (1982) An explicit formulation for an efficient triangular plate-bending element. Int J Numer Meth Eng 18:1077–1089
Bécache E, Derveaux G, Joly P (2005) An efficient numerical method for the resolution of the Kirchhoff–Love dynamic plate equation. Numer Methods Partial Differ Equ 21(2):323–348
Bischoff M, Bletzinger KU, Wall W, Ramm E (2004) Models and finite elements for thin-walled structures 2
da Veiga LB, Niiranen J, Stenberg R (2007) A family of \({C}^0\) finite elements for Kirchhoff plates I: error analysis. SIAM J Numer Anal 45:2047–2071
da Veiga LB, Niiranen J, Stenberg R (2008) A family of \({C}^0\) finite elements for Kirchhoff plates II: numerical results. Comput Method Appl Mech Eng 197(21):1850–1864
Huang J, Huang X, Xu Y (2011) Convergence of an adaptive mixed finite element method for Kirchhoff plate bending problems. SIAM J Numer Anal 49:574–607
Ibrahimbegović A (1993) Quadrilateral finite elements for analysis of thick and thin plates. Comput Method Appl Mech Eng 110:195–209
Oñate E, Zárate F (2000) Rotation-free triangular plate and shell elements. Int J Numer Meth Eng 47:557–603
Perotti L, Bompadre A, Ortiz M (2013) Automatically inf-sup compliant diamond-mixed finite elements for Kirchhoff plates. Int J Numer Meth Eng 96:405–424
Frangi A, Guiggiani M (1999) Boundary element analysis of Kirchhoff plates with direct evaluation of hypersingular integrals. Int J Numer Meth Eng
Benson D, Bazilevs Y, Hsu M, Hughes T (2010) Isogeometric shell analysis: The Reissner–Mindlin shell. Comput Method Appl Mech Eng 199(5):276–289
da Veiga LB, Hughes TJR, Kiendl J, Lovadina C, Niiranen J, Reali A, Speleers H (2015) A locking-free model for Reissner–Mindlin plates: analysis and isogeometric implementation via NURBS and triangular NURPS. Math Models Methods Appl Sci 25(08):1519–1551
Kiendl J, Bletzinger KU, Linhard J, Wüchner R (2009) Isogeometric shell analysis with Kirchhoff–Love elements. Comput Method Appl Mech Eng 198(49):3902–3914
Kiendl J, Hsu MC, Wu MC, Reali A (2015) Isogeometric Kirchhoff–Love shell formulations for general hyperelastic materials. Comput Method Appl Mech Eng 291:280–303
Zou Z, Hughes T, Scott M, Sauer R, Savitha E (2021) Galerkin formulations of isogeometric shell analysis: alleviating locking with Greville quadratures and higher-order elements. Comput Method Appl Mech Eng 380:113757
Zou Z, Scott M, Miao D, Bischoff M, Oesterle B, Dornisch W (2020) An isogeometric Reissner–Mindlin shell element based on Brézier dual basis functions: overcoming locking and improved coarse mesh accuracy. Comput Method Appl Mech Eng 370:113283
Henshaw WD (2006) A high-order accurate parallel solver for Maxwell’s equations on overlapping grids. SIAM J Sci Comput 28(5):1730–1765
Lambert JD (1973) Computational methods in ordinary differential equations. Cambridge University Press, New York
Papkova IV, Awrejcewicz NJ, Krysko AV, Krylova EY, Krysko VA (2020) Theory of flexible mesh-type shallow Kirchhoff-love structures based on the modified couple stress theory, vol 134, chap. 17. Springer, Cham, pp 331–344
Newmark NM (1959) A method of computation for structural dynamics. Proc Am Soc Civil Eng 85:67–74
Li M, Guan X, Mao S (2014) New error estimates of the Morley element for the plate bending problems. J Comput Appl Math 263:405–416
Ming W, Xu J (2006) The Morley element for fourth order elliptic equations in any dimensions. Numer Math 103:155–169
Morley LSD (1968) The triangular equilibrium element in the solution of plate bending problems. Aero Q 19:149–169
Ciarlet PG (2002) The finite element method for elliptic problems. SIAM, Philadelphia
Bukač M, Čanič S, Glowinski R, Tambača J, Quaini A (2013) Fluid-structure interaction in blood flow capturing non-zero longitudinal structure displacement. J Comput Phys 235:515–541
Bukač M, Čanič S, Muha B (2015) A partitioned scheme for fluid-composite structure interaction problems. J Comput Phys 281:493–517
Bukač M, Čanič S, Tambača J, Wang Y (2019) Fluid-structure interaction between pulsatile blood flow and a curved stented coronary artery on a beating heart: a four stent computational study. Comput Method Appl Mech Eng 350:679–700
Tuan PH, Wen CP, Chiang PY, Yu YT, Liang HC, Huang KF, Chen YF (2015) Exploring the resonant vibration of thin plates: reconstruction of Chladni patterns and determination of resonant wave numbers. J Acoust Soc Am 137:2113–2123
Chugh A (2007) Natural vibration characteristics of gravity structures. Int J Numer Anal Meth Geomech 31:607–648
Chesshire GS, Henshaw WD (1990) Composite overlapping meshes for the solution of partial differential equations. J Comput Phys 90(1):1–64
Acknowledgements
This research was supported by the Louisiana Board of Regents Support Fund under contract No. LEQSF(2018-21)-RD-A-23. L. Li is grateful to Professor W.D. Henshaw of Rensselaer Polytechnic Institute (RPI) for helpful conversations. Portions of this research were conducted with high performance computational resources provided by the Louisiana Optical Network Infrastructure (http://www.loni.org).
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendix A: Nodal line patterns for the eigenvalue problem
Appendix A: Nodal line patterns for the eigenvalue problem
We show the results of the eigenvalue problem (31) here. Nodal lines of the first 25 eigenmodes (with multiplicity) for the square plate with clamped edges and the annular plate with simply supported boundaries are shown in Figs. 18 and 19, respectively. The eigenmodes plotted for each degenerated pair are arbitrary so they can be asymmetric.
Rights and permissions
About this article
Cite this article
Nguyen, D.T.A., Li, L. & Ji, H. Stable and accurate numerical methods for generalized Kirchhoff–Love plates. J Eng Math 130, 6 (2021). https://doi.org/10.1007/s10665-021-10163-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s10665-021-10163-x