Nonlinear damping in a micromechanical oscillator

Nonlinear elastic effects play an important role in the dynamics of microelectromechanical systems (MEMS). A Duffing oscillator is widely used as an archetypical model of mechanical resonators with nonlinear elastic behavior. In contrast, nonlinear dissipation effects in micromechanical oscillators are often overlooked. In this work, we consider a doubly clamped micromechanical beam oscillator, which exhibits nonlinearity in both elastic and dissipative properties. The dynamics of the oscillator is measured in both frequency and time domains and compared to theoretical predictions based on a Duffing-like model with nonlinear dissipation. We especially focus on the behavior of the system near bifurcation points. The results show that nonlinear dissipation can have a significant impact on the dynamics of micromechanical systems. To account for the results, we have developed a continuous model of a geometrically nonlinear beam-string with a linear Voigt–Kelvin viscoelastic constitutive law, which shows a relation between linear and nonlinear damping. However, the experimental results suggest that this model alone cannot fully account for all the experimentally observed nonlinear dissipation, and that additional nonlinear dissipative processes exist in our devices.     

Dissipation-induced instabilities and symmetry

The paradox of destabilization of a conservative or non-conservative system by small dissipation,or Ziegler’s paradox(1952),has stimulated a growing interest in the sensitivity of reversible and Hamiltonian systems with respect to dissipative perturbations.Since the last decade it has been widely accepted that dissipation-induced instabilities are closely related to singularities arising on the stability boundary,associated with Whitney’s umbrella.The first explanation of Ziegler’s paradox was given(much earlier)by Oene Bottema in 1956.The aspects of the mechanics and geometry of dissipation-induced instabilities with an application to rotor dynamics are discussed.     

土的本构方程与热力学

介绍一种基于热力学理论建立土力学本构方程的一般性理论框架. 这一方法利用两个势函数即自由能势函数和耗散势函数(或屈服函数)以及固定的过程和框架, 建立土的本构方程. 简要介绍了建立热力学本构方程中所用到的热力学内变量理论, 利用Legendre变换建立了热力学势函数之间以及各耗散函数与屈服函数之间的关系;利用自由能势函数和耗散势函数(或屈服函数)建立土的本构方程及其具体步骤. 最后讨论了土力学本构方程研究的意义以及它和应用之间的关系.      

Non-linear viscoelastodynamic equations of three-dimensional rotating structures in finite displacement and finite element discretization

This paper deals with the non-linear viscoelastodynamics of three-dimensional rotating structure undergoing finite displacement. In addition, the non-linear dynamics is studied with respect to geometrical and mechanical perturbations. On part of the boundary of the structure, a rigid body displacement field is applied which moves the structure in a rotation motion. A time-dependent Dirichlet condition is applied to another part of the boundary. For instance, this corresponds to the cycle step of a helicopter rotor blade. A surface force field is applied to the third part of the boundary and depends on the time history of the structural displacement field. For example, this might corresponds to general unsteady aerodynamics forces applied to the structure. The objective of this paper is to model the non-linear dynamic behavior of such a rotating viscoelastic structure undergoing finite displacements, and to allow small geometrical and mechanical (mass, constitutive equations) perturbations analysis to be performed. The model is constructed by the introduction of a reference configuration which is deduced from the non-linear steady boundary value problem. A constitutive equation deduced from the Coleman and Noll theory concerning the viscoelasticity in finite displacement is used. Thereafter, the weak formulation of the boundary value problem is constructed and discretized using the finite element method. In order to simplify the mathematical study of the equations, multilinear forms are introduced in the algebraic calculation and their mathematical properties are presented.     

Relationship between plastic and intrinsic dissipations

In this paper, the relationship between the plastic and intrinsic dissipations is addressed within the normality structure of [Rice, J.R., 1971. Inelastic constitutive relations for solids: an integral variable theory and its application to metal plasticity. J. Mech. Phys. Solids 19, 433–455; Rice, J.R., 1975. Continuum mechanics and thermodynamics of plasticity in relation to microscale deformation mechanisms. In: Argon, A.S. (Ed.), Constitutive Equations in Plasticity. MIT Press, Cambridge, MA, pp. 23–79.] It is shown that the plastic dissipation is generally not equal to the intrinsic dissipation. Within the normality structure, the microscale and macroscale thermodynamic fluxes and forces are related by the conditions of energy and dissipation equivalence. If the plastic dissipation is required to be equal to the intrinsic dissipation, J₂ potential and the Levy–Mises equation are recovered from the condition of dissipation equivalence for incompressible plastic flows.     

Nonlinear longitudinal forced vibration of a hysteretic bar: An analytical solution

This work investigates the nonlinear longitudinal forced vibrations of a bar with hysteretic dynamics typical of rocks and man-made geomaterials. The material properties are described by a constitutive relation that includes energy dissipation effects consistently with hysteretic dynamics. The harmonic balance method and a perturbation technique that uses the material's modulus defect as perturbation parameter are employed to solve the equation of motion. The spatial dependence of the modulus defect, which is disregarded in earlier solutions of this problem, is duly accounted for. According to this model, the resonance frequency shift is proportional to the amplitude of the excitation, while the inverse of the quality factor increases as the square root of the latter. Further, the spatial variation of the modulus defect is shown to affect the modulation of the fundamental component along the bar. The nonlinear spectral components are of odd order, with amplitude proportional to the maximum value of the modulus defect and to the square of the excitation's amplitude, and decreases non-monotonically with increasing harmonic order. The motivation for this work is twofold. Analytical models may improve our understanding of the dynamics of hysteretic materials in general, and of the mutual interaction of material defects in particular. Secondly, this work establishes a bench-mark result on a mathematically simple hysteretic system against which alternative mathematical approaches, possibly concerning material with complex constitutive relations, could be tested.     

Newtonian Limit of Maxwell Fluid Flows

In this paper, we revise Maxwell’s constitutive relation and formulate a system of first-order partial differential equations with two parameters for compressible viscoelastic fluid flows. The system is shown to possess a nice conservation–dissipation (relaxation) structure and therefore is symmetrizable hyperbolic. Moreover, for smooth flows we rigorously verify that the revised Maxwell’s constitutive relations are compatible with Newton’s law of viscosity.     

基于Riesz势空间分数阶算子的非局部粘弹性力学元件

将幂函数引入Eringen非局部线粘弹性本构,导出Riesz势形式的应力-应变关系。利用该关系,构造非局部弹簧和非局部阻尼器两类元件;利用元件的串联和并联,建立非局部Kelvin和非局部Maxwell粘弹性模型,推导模型的松弛模量和蠕变柔量。进一步,给出非局部粘弹性模型在生物组织超声波耗散建模中的应用。     

Crystallization under external pressure

The influence of external pressure on the dynamics of crystallization is examined by considering a solid–liquid two-phase system in a cylinder closed by a piston. The dynamic equations are derived using three methods, namely, Rational Thermodynamics (Liu’s procedure), the Matrix Model, and the general equation for the nonequilibrium reversible-irreversible coupling (GENERIC) formalism. The constitutive relation for the multiphase pressure on the piston is identical for all three methods, whereas some aspects of the result for the phase change dynamics differ. The rational thermodynamics treatment constrains the phase change dynamics of only those structural variables that enter into the dissipation inequality, whereas the other two formalisms make statements about the phase change of all structural variables. Nevertheless, all three methods show, first, that the phase change happens instantaneously at constant volume and, second, where morphological detail can be built into the model without violating thermodynamic principles. It is discussed how an appropriate choice of the morphological variables allows one to incorporate impingement of crystals and depletion zone effects, as well as to distinguish crystal shapes.     

Transition waves in bistable structures. II. Analytical solution: wave speed and energy dissipation

We consider dynamics of chains of rigid masses connected by links described by irreversible, piecewise linear constitutive relation: the force-elongation diagram consists of two stable branches with a jump discontinuity at the transition point. The transition from one stable state to the other propagates along the chain and excites a complex system of waves. In the first part of the paper (Cherkaev et al., 2004, Transition waves in bistable structures. I. Delocalization of damage), the branches could be separated by a gap where the tensile force is zero, the transition wave was treated as a wave of partial damage. Here we assume that there is no zero-force gap between the branches. This allows us to obtain steady-state analytical solutions for a general piecewise linear trimeric diagram with parallel and nonparallel branches and an arbitrary jump at the transition. We derive necessary conditions for the existence of the transition waves and compute the speed of the wave. We also determine the energy of dissipation which can be significantly increased in a structure characterized by a nonlinear discontinuous constitutive relation. The considered chain model reveals some phenomena typical for waves of failure or crushing in constructions and materials under collision, waves in a structure specially designed as a dynamic energy absorber and waves of phase transitions in artificial and natural passive and active systems.     

Maximum norm wellposedness of nonlinear kinematic hardening models

The Chaboche model represents a widely popular rate independent elastoplastic stress-strain law of nonlinear kinematic hardening type. We consider it as a special case of a certain system of ordinary differential equations coupled to a hysteresis nonlinearity which expresses the principle of maximal dissipation of energy. For this system, we prove wellposedness with respect to the maximum norm, that is, continuous dependence of the stress on the strain (or vice versa) with respect to perturbations of small amplitude. Received: May 15, 1997; accepted: September 7, 1997     

Dissipation in a dilute suspension of spheres in a second-order fluid

We examine the effective medium properties of a dilute suspension of spheres in a second-order fluid under linear shear. Since the second-order fluid is the first step toward the general viscoelastic fluid, the results obtained may provide a qualitative feel for the problem in which the suspending fluid obeys a more complicated (and realistic) constitutive relation.The dissipation in the medium is calculated by determining the rate of working by surface forces; this is compared to the dissipation in a homogeneous fluid to give the effective properties. The results show that the term linear in volume fraction increases the corresponding rheological coefficient, just as in the Newtonian case. It is to be noted that the second-order dissipation is zero for simple shear and other weak flows, whereas for strong flows the small correction may increase or decrease the overall dissipation.     

Thermoelasticity without energy dissipation

This paper deals with thermoelastic material behavior without energy dissipation; it deals with both nonlinear and linear theories, although emphasis is placed on the latter. In particular, the linearized theory of thermoelasticity discussed possesses the following properties: (a) the heat flow, in contrast to that in classical thermoelasticity characterized by the Fourier law, does not involve energy dissipation; (b) a constitutive equation for an entropy flux vector is determined by the same potential function which also determines the stress; and (c) it permits the transmission of heat as thermal waves at finite speed. Also, a general uniqueness theorem is proved which is appropriate for linear thermoelasticity without energy dissipation.     

Linear Dynamics of Perturbations in Flows with Constant Horizontal Shear

The dynamics of perturbations in shallow water and incompressible stratified fluid flows with constant horizontal shear are described using the nonmodal analysis. It is shown that the shear flow perturbations can be divided into two classes on the basis of the potential vorticity: rapidly oscillating wave perturbations with zero potential vorticity and slow vortex perturbations with nonzero potential vorticity. In the cases of weak and strong shear the main features of the dynamics of wave and vortex perturbations are studied analytically (using the WKBG method) and numerically. It is shown that for large times the wave perturbation energy increases linearly, i.e., the shear flow is algebraically unstable due to the growth of rapid wave perturbations. This instability can be of importance in processes of turbulence development and surface and internal wave generation.     

On the Iwan models for lap-type bolted joints

The paper presents mathematical modeling of the non-linear constitutive relation for bolted joints in the framework of the Kragelsky-Demkin theory of rough contact. It is shown that this approach, which maintains the tribology-related features of bolted joint interfaces, leads to a singular Iwan distribution density. In particular, we show that the Iwan density is expressed in terms of the height distribution density of the surface asperities, whereas its singular exponent is determined by the shape exponent of the surface asperities. Following this, constitutive relations for lap joints and the corresponding backbone (force-deflection) curves are obtained. Finally, Masing's hypothesis is applied and Goodman's relation for energy dissipation is recovered in order to describe the effects of cyclic loading. The two cases of a rough surface in contact with a flat surface and of two contacting rough surfaces are treated separately.     

On kinetics of hysteresis

An initiation criterion and a kinetic relation describing hysteresis in shape memory alloys consistent with non-equilibrium thermodynamics are proposed. The initiation of phase transition begins always at zero driving force, defined as the negative derivative of free energy with respect to the volume fraction of high strain phase. The kinetic relation is obtained from the dissipation potential which is proportional to the magnitude of the volume fraction rate times the newly formed volume fraction. The proposed theory turns out to be rate-independent, but history dependent, and can describe all features of hysteresis loops observed in experiments.      

On a general concept for the analysis of crack growth and material damage

For finite strain dynamics a variational model of crack evolution is formulated within the generalized oriented continuum methodology. In this approach position- and direction-dependent deformation and strain measures are used to describe the (macro)motion of the body with defects, which may evolve relative to the moving body. The inelastic behaviour of continua with evolving defects is represented by phenomenological equations including the transversal crack evolution. A strain-induced crack propagation criterion is defined by the difference between the strain energy release rate and the rate of the surface energy of the crack. A possible nucleation of microcracks in terms of the average drag coefficient of the crack configuration is proposed. Based on the crack growth criterion presented in this paper, the kinking of cracks is investigated using variational concepts. A constitutive damage model of Kachanov's type accounting for the crack density is derived in terms of the free energy functional and a dissipation potential.     

On multiscale significance of Rice’s normality structure

The normality structure proposed by [Rice, J.R., 1971. Inelastic constitutive relations for solids: an integral variable theory and its application to metal plasticity. J. Mech. Phys. Solids 19, 433–455.] provides a minimal framework of multiscale thermodynamics. As shown in this paper, Rice’s multiscale thermodynamic formalism is exactly consistent with Ziegler’s essential notion [Ziegler, H., 1977. An Introduction to Thermomechanics, North-Holland, Amsterdam.] that the entire constitutive response is determined by the knowledge of two scalar potential functions: an energy function and a dissipation function. In Rice’s multiscale thermodynamic formulation, the variational equation relating macroscale and microscale thermodynamic fluxes and forces plays a central role and ensures the equality between microscale and macroscale dissipation rate. The variational equation can be further reformulated into a principle of maximum equivalent dissipation. Based on the variation equation, the transformation from microscale to macroscale is characterized by two linear transformations with the same corresponding matrix.     

Coherent vorticity dynamics and dissipation in a utility-scale wind turbine wake with uniform inflow

The vorticity dynamics and its relationship to dissipation in the wake of a utility-scale wind turbine are investigated through large-eddy simulation. The vorticity dynamics is assessed through the enstrophy, which is related to the turbulent dissipation. The averaged enstrophy and turbulent dissipation are shown to be quantitatively similar in the wake. Using temporal phase averaging, the vorticity fluctuations are decomposed into coherent and random fluctuations with respect to the frequency of the tip vortices. The enstrophy in the tip vortices is dominated by coherent fluctuations, while the coherent fluctuations of root vortices are immediately saturated by the random vorticity fluctuations of the unstable hub vortex. The coherent strain rate has significant differences compared to the coherent enstrophy within one diameter downwind of blade tip, but the random enstrophy and strain rate are relatively similar. Differences in coherent enstrophy and strain rate decrease further from the rotor.     

Dynamics of a gravity stonewall

Treated is the dynamics of a gravity stonewall. The wall is excited by a transient damped periodic oscillation simulating an earthquake. The model adopts a stick-slip friction constitutive law. Sensitivity of energy dissipation to parameters such as number of blocks, friction coefficient, sticktion and slipping stiffness and excitation amplitude and frequency is determined. A 2-D model of the monolithic wall is also analyzed to compare displacement and shear stress of the two constructions.