TY - JOUR
T1 - Spectral Analysis of a Mixed Method for Linear Elasticity
AU - ZHONG, Xiang
AU - QIU, Weifeng
PY - 2023
Y1 - 2023
N2 - The purpose of this paper is to analyze a mixed method for the linear elasticity eigenvalue problem, which approximates numerically the stress, displacement, and rotation, by piecewise (k + 1), k, and (k + 1) th degree polynomial functions (k ≥ 1), respectively. The numerical eigenfunction of stress is symmetric. By the discrete H1-stability of numerical displacement, we prove an O (hk+2) approximation to the L2-orthogonal projection of the eigenspace of exact displacement for the eigenvalue problem with a proper regularity assumption. Thus via postprocessing, we obtain a better approximation to the eigenspace of exact displacement for the eigenproblem than conventional methods. We also prove that numerical approximation to the eigenfunction of stress is locking free with respect to the Poisson ratio. We introduce a hybridization to reduce the mixed method to a condensed eigenproblem and prove an O (h2) initial approximation (independent of the inverse of the elasticity operator) of the eigenvalue for the nonlinear eigenproblem by using the discrete H1 -stability of numerical displacement, while only an O (h) approximation can be obtained if we use the traditional inf-sup condition. Finally, we report some numerical experiments. © 2023 The Author(s).
AB - The purpose of this paper is to analyze a mixed method for the linear elasticity eigenvalue problem, which approximates numerically the stress, displacement, and rotation, by piecewise (k + 1), k, and (k + 1) th degree polynomial functions (k ≥ 1), respectively. The numerical eigenfunction of stress is symmetric. By the discrete H1-stability of numerical displacement, we prove an O (hk+2) approximation to the L2-orthogonal projection of the eigenspace of exact displacement for the eigenvalue problem with a proper regularity assumption. Thus via postprocessing, we obtain a better approximation to the eigenspace of exact displacement for the eigenproblem than conventional methods. We also prove that numerical approximation to the eigenfunction of stress is locking free with respect to the Poisson ratio. We introduce a hybridization to reduce the mixed method to a condensed eigenproblem and prove an O (h2) initial approximation (independent of the inverse of the elasticity operator) of the eigenvalue for the nonlinear eigenproblem by using the discrete H1 -stability of numerical displacement, while only an O (h) approximation can be obtained if we use the traditional inf-sup condition. Finally, we report some numerical experiments. © 2023 The Author(s).
KW - linear elasticity
KW - eigenvalue problem
KW - error estimates
KW - mixed methods
UR - http://www.scopus.com/inward/record.url?scp=85169684602&partnerID=8YFLogxK
UR - https://www.scopus.com/record/pubmetrics.uri?eid=2-s2.0-85169684602&origin=recordpage
U2 - 10.1137/22M148611X
DO - 10.1137/22M148611X
M3 - RGC 21 - Publication in refereed journal
VL - 61
SP - 1885
EP - 1917
JO - SIAM Journal on Numerical Analysis
JF - SIAM Journal on Numerical Analysis
SN - 0036-1429
IS - 4
ER -