Asymptotic and spectral results on the initial boundary‐value problem for the coupled bending‐torsion vibration model (which is important in such areas of engineering sciences as bridge and tall building designs, aerospace and oil pipes modeling, etc.) are presented. The model is given by a system of two hyperbolic partial differential equations equipped with a three‐parameter family of non‐self‐adjoint (linear feedback type) boundary conditions modeling the actions of self‐straining actuators. The system is rewritten in the form of the first‐order evolution equation in a Hilbert space of a four‐component Cauchy data. It is shown that the dynamics generator is a matrix differential operator with compact resolvent, whose discrete spectrum splits asymptotically into two disjoint subsets called the α‐branch and the β‐branch, respectively. Precise spectral asymptotics for the eigenvalues from each branch as the number of an eigenvalue tends to ∞ have been derived. It is also shown that the leading asymptotical term of the α‐branch eigenvalue depends only on the
 NSFPAR ID:
 10420587
 Publisher / Repository:
 WileyBlackwell
 Date Published:
 Journal Name:
 Studies in Applied Mathematics
 Volume:
 150
 Issue:
 4
 ISSN:
 00222526
 Page Range / eLocation ID:
 p. 9961025
 Format(s):
 Medium: X
 Sponsoring Org:
 National Science Foundation
More Like this

Mathematical analysis of the well known model of a piezoelectric energy harvester is presented. The harvester is designed as a cantilever Timoshenko beam with piezoelectric layers attached to its top and bottom faces. Thin, perfectly conductive electrodes are covering the top and bottom faces of the piezoelectric layers. These electrodes are connected to a resistive load. The model is governed by a system of three partial differential equations. The first two of them are the equations of the Timoshenko beam model and the third one represents Kirchhoff’s law for the electric circuit. All equations are coupled due to the piezoelectric effect. We represent the system as a single operator evolution equation in the Hilbert state space of the system. The dynamics generator of this evolution equation is a nonselfadjoint matrix differential operator with compact resolvent. The paper has two main results. Both results are explicit asymptotic formulas for eigenvalues of this operator, i.e., the modal analysis for the electrically loaded system is performed. The first set of the asymptotic formulas has remainder terms of the order O ( 1 n ) , where n is the number of an eigenvalue. These formulas are derived for the model with variable physical parameters. The second set of the asymptotic formulas has remainder terms of the order O ( 1 n 2 ) , and is derived for a less general model with constant parameters. This second set of formulas contains extra term depending on the electromechanical parameters of the model. It is shown that the spectrum asymptotically splits into two disjoint subsets, which we call the α branch eigenvalues and the θ branch eigenvalues. These eigenvalues being multiplied by “i” produce the set of the vibrational modes of the system. The α branch vibrational modes are asymptotically located on certain vertical line in the left half of the complex plane and the θ branch is asymptotically close to the imaginary axis. By having such spectral and asymptotic results, one can derive the asymptotic representation for the mode shapes and for voltage output. Asymptotics of vibrational modes and mode shapes is instrumental in the analysis of control problems for the harvester.more » « less

null (Ed.)The distribution of natural frequencies of the Euler–Bernoulli beam resting on elastic foundation and subject to an axial force in the presence of several damping mechanisms is investigated. The damping mechanisms are: ( i ) an external or viscous damping with damping coefficient ( − a 0 ( x )), ( ii ) a damping proportional to the bending rate with the damping coefficient a 1 ( x ). The beam is clamped at the left end and equipped with a fourparameter (α, β, κ 1 , κ 2 ) linear boundary feedback law at the right end. The 2 × 2 boundary feedback matrix relates the control input (a vector of velocity and its spacial derivative at the right end) to the output (a vector of shear and moment at the right end). The initial boundary value problem describing the dynamics of the beam has been reduced to the first order in time evolution equation in the state Hilbert space of the system. The dynamics generator has a purely discrete spectrum (the vibrational modes). Explicit asymptotic formula for the eigenvalues as the number of an eigenvalue tends to infinity have been obtained. It is shown that the boundary control parameters and the distributed damping play different roles in the asymptotical formulas for the eigenvalues of the dynamics generator. Namely, the damping coefficient a 1 and the boundary controls κ 1 and κ 2 enter the leading asymptotical term explicitly, while damping coefficient a 0 appears in the lower order terms.more » « less

Classical statistical mechanics has long relied on assumptions such as the equipartition theorem to understand the behavior of the complicated systems of many particles. The successes of this approach are well known, but there are also many wellknown issues with classical theories. For some of these, the introduction of quantum mechanics is necessary, e.g., the ultraviolet catastrophe. However, more recently, the validity of assumptions such as the equipartition of energy in classical systems was called into question. For instance, a detailed analysis of a simplified model for blackbody radiation was apparently able to deduce the Stefan–Boltzmann law using purely classical statistical mechanics. This novel approach involved a careful analysis of a “metastable” state which greatly delays the approach to equilibrium. In this paper, we perform a broad analysis of such a metastable state in the classical Fermi–Pasta–Ulam–Tsingou (FPUT) models. We treat both the αFPUT and βFPUT models, exploring both quantitative and qualitative behavior. After introducing the models, we validate our methodology by reproducing the wellknown FPUT recurrences in both models and confirming earlier results on how the strength of the recurrences depends on a single system parameter. We establish that the metastable state in the FPUT models can be defined by using a single degreeoffreedom measure—the spectral entropy (η)—and show that this measure has the power to quantify the distance from equipartition. For the αFPUT model, a comparison to the integrable Toda lattice allows us to define rather clearly the lifetime of the metastable state for the standard initial conditions. We next devise a method to measure the lifetime of the metastable state tm in the αFPUT model that reduces the sensitivity to the exact initial conditions. Our procedure involves averaging over random initial phases in the plane of initial conditions, the P1Q1 plane. Applying this procedure gives us a powerlaw scaling for tm, with the important result that the power laws for different system sizes collapse down to the same exponent as Eα2→0. We examine the energy spectrum E(k) over time in the αFPUT model and again compare the results to those of the Toda model. This analysis tentatively supports a method for an irreversible energy dissipation process suggested by Onorato et al.: fourwave and sixwave resonances as described by the “wave turbulence” theory. We next apply a similar approach to the βFPUT model. Here, we explore in particular the different behavior for the two different signs of β. Finally, we describe a procedure for calculating tm in the βFPUT model, a very different task than for the αFPUT model, because the βFPUT model is not a truncation of an integrable nonlinear model.more » « less

null (Ed.)Abstract Painful herniated discs are treated surgically by removing extruded nucleus pulposus (NP) material (nucleotomy). NP removal through enzymatic digestion is also commonly performed to initiate degenerative changes to study potential biological repair strategies. Experimental and computational studies have shown a decrease in disc stiffness with nucleotomy under single loading modalities, such as compressiononly or bendingonly loading. However, studies that apply more physiologically relevant loading conditions, such as compression in combination with bending or torsion, have shown contradicting results. We used a previously validated bone–disc–bone finite element model (Control) to create a Nucleotomy model to evaluate the effect of dual loading conditions (compression with torsion or bending) on intradiscal deformations. While disc joint stiffness decreased with nucleotomy under single loading conditions, as commonly reported in the literature, dual loading resulted in an increase in bending stiffness. More specifically, dual loading resulted in a 40% increase in bending stiffness under flexion and extension and a 25% increase in stiffness under lateral bending. The increase in bending stiffness was due to an increase and shift in compressive stress, where peak stresses migrated from the NP–annulus interface to the outer annulus. In contrast, the decrease in torsional stiffness was due to greater fiber reorientation during compression. In general, large radial strains were observed with nucleotomy, suggesting an increased risk for delamination or degenerative remodeling. In conclusion, the effect of nucleotomy on disc mechanics depends on the type and complexity of applied loads.more » « less

Abstract Quartz crystals with zircon inclusions were synthesized using a pistoncylinder apparatus to experimentally evaluate the use of inclusions in “soft” host minerals for elastic thermobarometry. Synthesized zircon inclusion strains and, therefore, pressures (
P _{inc}) were measured using Raman spectroscopy and then compared with the expected inclusion strains and pressures calculated from elastic models. Measured inclusion strains and inclusion pressures are systematically more tensile than the expected values and, thus, recalculated entrapment pressures are overestimated. These discrepancies are not caused by analytical biases or assumptions in the elastic models and strain calculations. Analysis shows that inclusion strain discrepancies progressively decrease with decreasing experimental temperature in the αquartz field. This behavior is consistent with inelastic deformation of the host–inclusion pairs induced by the development of large differential stresses during experimental cooling. Therefore, inclusion strains are more reliable for inclusions trapped at lower temperature conditions in the αquartz field where there is less inelastic deformation of the host–inclusion systems. On the other hand, entrapment isomekes of zircon inclusions entrapped in the βquartz stability field plot along the α–β quartz phase boundary, suggesting that the inclusion strains were mechanically reset at the phase boundary during experimental cooling and decompression. Therefore, inclusions contained in soft host minerals can be used for elastic thermobarometry and inclusions contained in βquartz may provide constraints on theP –T at which the host–inclusion system crossed the phase boundary during exhumation.