Dynamic responses to sinusoidal excitations of beams with frictional joints

We consider the dynamic responses of a beam with a frictional joint. The frictional force at the joint is modeled using the Coulomb friction model. The frictional force at the joint makes the nature of the boundary conditions at the joint uncertain. Therefore, this problem represents a type of nonlinear problems where the boundary conditions are coupled to the solutions. Using numerical integration of the resulting differential equations obtained by combining the finite element method and the Lagrange equations, we study the steady-state solutions of the system to sinusoidal excitations. We explore the dependence of the system responses to various parameters including the frictional force, the forcing frequency and the forcing amplitude. A result of special interest is the existence of an optimum friction force if the frictional joint is used to control the system response amplitude. We also examine the ways that friction affects the resonance frequency of the structure. Experiments are carried out, which agree qualitatively with the numerical results.     

The fractional-order SIS epidemic model with variable population size

In this work, we deal with the fractional-order SIS epidemic model with constant recruitment rate, mass action incidence and variable population size. The stability of equilibrium points is studied. Numerical solutions of this model are given. Numerical simulations have been used to verify the theoretical analysis.     

Elasticity solution of multi-span beams with variable thickness under static loads

This paper studies the stress and displacement distributions of continuously varying thickness multi-span beams simply supported at two ends and under static loads. The intermediate supports of the beam may be elastic and/or rigid in one or two directions. On the basis of the two-dimensional plane elasticity theory, the general solution of stress function, which exactly satisfies the governing differential equations and the simply supported boundary conditions, is deduced. In the present analysis, the reaction forces of the intermediate supports are regarded as the unknown external forces acting on the lower surface of the beam under consideration. The unknown coefficients in the solutions are determined by using the Fourier sinusoidal series expansions to the boundary conditions on the upper and lower surfaces of the beam and using the linear relations between reaction forces and displacements of the beam at intermediate supports. The solution obtained is exact and excellent convergence has been confirmed. Comparing the numerical results obtained from the proposed method to those obtained from the Euler beam theory, the Timoshenko beam theory and those obtained from the commercial finite element software ANSYS, high accuracy of the present method is demonstrated.     

β-Expectation tolerance regions for a generalized multivariate model with normal error variables

In this paper a procedure of construction of β-expectation tolerance regions in the framework of the structural method of inference has been developed. The procedure has been applied to the generalized multivariate model and the β-expectation tolerance region for this case has been constructed assuming the normal distribution for the error variables of the model.     

Existence of global BV solutions to a model of reacting Euler fluid with variable thermodynamics parameters

We investigate the Cauchy problems of Euler equations coupled with chemical reaction equation. The specific heat and adiabatic parameters of compressible fluid are introduced as functions with respect to reactant mass fraction in this paper. By means of fractional step wave-front tracking scheme, we derive the existence of global BV solution, and verify that it is indeed an entropy weak solution. This system is also regarded as balance law with resonant characteristic fields. Hence some decay condition is imposed on the source term, in order to handle the damping effect.     

Contact problem for a thin elastic layer with variable thickness: Application to sensitivity analysis of articular contact mechanics

In the framework of the recently developed asymptotic models for tibio-femoral contact incorporating frictionless elliptical contact interaction between thin elastic, viscoelastic, or biphasic cartilage layers, we apply an asymptotic modeling approach for analytical evaluating the sensitivity of crucial parameters in joint contact mechanics due to small variations in the thicknesses of the contacting cartilage layers. The four term asymptotic expansion for the normal displacement at the contact surface is explicitly derived, which recovers the corresponding solution obtained previously for the 2D case in the compressible case. It was found that to minimize the influence of the cartilage thickness non-uniformity on the force–displacement relationship, the effective thicknesses of articular layers should be determined from a special optimization criterion.     

A coradiant based scalarization to characterize approximate solutions of vector optimization problems with variable ordering structures

This paper investigates some properties of approximate efficiency in variable ordering structures where the variable ordering structure is given by a special set valued map. We characterize $\epsilon$-minimal and $\epsilon$- nondominated elements as approximate solutions of a multiobjective optimization problem with a variable ordering structure and give necessary and sufficient conditions for these solutions, via scalarization.     

A generalized model of an irreversible thermal Brownian refrigerator and its performance

A generalized model of an irreversible thermal Brownian refrigerator, which consists of Brownian particles moving in a periodic sawtooth potential with external forces and contacting with the alternating hot and cold reservoirs along the space coordinate, is established in this paper. The heat flows driven by both potential and kinetic energies of the particles as well as the heat leakage between the hot and cold reservoirs are taken into account. The optimum performance of the generalized model is analyzed using the theory and method of finite time thermodynamics. The analytical expressions for cooling load, coefficient of performance (COP) and power input of the Brownian refrigerator are derived. It is shown by numerical examples that due to the heat leakage between the heat reservoirs and heat flow via the change of kinetic energy of the particles, the Brownian refrigerator is always irreversible and the COP can never attain the Carnot COP. The influences of the heat leakage, the external force, barrier height of the potential, asymmetry of the sawtooth potential and temperature ratio of the heat reservoirs on the performance of the Brownian refrigerator are also investigated in detail. The effective regions of external force and barrier height of the potential in which the Brownian motor can operate as a refrigerator are determined. It is found that the performance of the Brownian refrigerator depends strictly on the design parameters. If these design parameters are properly chosen, the Brownian refrigerator can be controlled to operate in the optimal regimes. The results obtained herein about the general Brownian refrigerator model include those obtained in many previous literatures.     

Turing-Hopf bifurcation in the reaction-diffusion system with delay and application to a diffusive predator-prey model

The interactions of diffusion-driven Turing instability and delay-induced Hopf bifurcation always give rise to rich spatiotemporal dynamics. In this paper, we first derive the algorithm for the normal forms associated with the Turing-Hopf bifurcation in the reaction-diffusion system with delay, which can be used to investigate the spatiotemporal dynamical classification near the Turing-Hopf bifurcation point in the parameter plane. Then, we consider a diffusive predator-prey model with weak Allee effect and delay. Through investigating the dynamics of the corresponding normal form of Turing-Hopf bifurcation induced by diffusion and delay, the spatiotemporal dynamics near this bifurcation point can be divided into six categories. Especially, stable spatially homogeneous/inhomogeneous periodic solutions and steady states, coexistence of two stable spatially inhomogeneous periodic solutions, coexistence of two stable spaially inhomogeneous steady states and the transition from one kind of spatiotemporal patterns to another are found.     

The dynamics of an impulsive delay SI model with variable coefficients

An impulsive delayed SI model with variable coefficients and a nonlinear incidence is formulated and analyzed. By introducing three thresholds, we obtain sufficient conditions for eradication and permanence of the disease, respectively. It is shown that the conditions depend on time delay for both the global attractivity of the positive infection-free periodic solution and permanence of the model. Furthermore, our results indicate that the disease will disappear if the ratio of the maximum to minimum of the pulse vaccination rate is lager than some value. The main feature of this paper is that we introduce multi-delays and variable coefficients into the SI model, and exhibit a new method which is applied to investigate this model. Numerical results show that the system we considered has complex dynamics including periodic and quasi-periodic oscillations.     

Analysis of multivariate skew normal models with incomplete data

We establish computationally flexible methods and algorithms for the analysis of multivariate skew normal models when missing values occur in the data. To facilitate the computation and simplify the theoretic derivation, two auxiliary permutation matrices are incorporated into the model for the determination of observed and missing components of each observation. Under missing at random mechanisms, we formulate an analytically simple ECM algorithm for calculating parameter estimation and retrieving each missing value with a single-valued imputation. Gibbs sampling is used to perform a Bayesian inference on model parameters and to create multiple imputations for missing values. The proposed methodologies are illustrated through a real data set and comparisons are made with those obtained from fitting the normal counterparts.     

Global asymptotic stability of a generalized SIRS epidemic model with transfer from infectious to susceptible

In this paper, we propose a generalized SIRS epidemic model with varying total population size caused by the death rate due to the disease and transfer from infectious to susceptible, where the incidence rate employed includs a wide range of monotonic and concave incidence rates. Applying the geometric approach developed by Smith, Li and Muldowey, we prove that the endemic equilibrium is globally asymptotically stable provided that the rate of loss of immuity $\delta$ is in a critical interval $[\eta,\bar\delta)$ when the basic reproduction number $R_0$ is greater than unity.     

Codimension-two bifurcation analysis of the continuous stirred tank reactor model with delay

The aim of this paper is to research the dynamical behaviors of the continuous stirred tank reactor (CSTR) model with delay. Firstly, we discuss the situation that its related characteristic equation has a simple zero root and a pair of purely imaginary roots. Secondly, the center manifold method and the normal form method are used to reduce the equation of CSTR model. Finally, some characteristics about the CSTR model can be obtained. We analyze three different topological structure and give entire bifurcation diagrams and phase portraits, which are innovative phenomenon. At the end, we obtain the stable and unstable periodic solutions by numerical simulation.     

The application of the auxiliary equation technique to a generalized mKdV equation with variable coefficients

The auxiliary differential equation technique is employed to investigate a generalized mKdV equation with variable coefficients. The Jacobi elliptic function wave-like solutions of the equation are expressed under several circumstances. The degenerated soliton-like and trigonometric function solutions are discussed in detail as the modulus of the Jacobi elliptic wave-like solutions tends to 1 and 0, respectively.     

一类具有总人口变化含时滞的SIS流行病模型分析

§ 1　ＩｎｔｒｏｄｕｃｔｉｏｎＩｎａｎｅｐｉｄｅｍｉｏｌｏｇｉｃａｌｍｏｄｅｌ,ｉｆｔｈｅｒｅａｒｅｍａｎｙｄｉｓｅａｓｅ ｒｅｌａｔｅｄｄｅａｔｈｓｏｒｔｈｅｂｉｒｔｈｓａｒｅｎｏｔｂａｌａｎｃｅｄｂｙｔｈｅｄｅａｔｈｓ,ｔｈｅｎｔｈｅｐｏｐｕｌａｔｉｏｎｓｉｚｅｍｕｓｔｂｅａｓｓｕｍｅｄｔｏｂｅａｆｕｎｃｔｉｏｎｏｆｔｉｍｅ[1] .ＩｎＳＩＳｅｐｉｄｅｍｉｏｌｏｇｉｃｍｏｄｅｌｓ,ｓｕｓｃｅｐｔｉｂｌｅｓｂｅｃｏｍｅｉｎｆｅｃｔｅｄａｆｔｅｒｓｕｆｆｉｃｉｅｎｔｃｏｎｔａｃｔｗｉｔｈａｎｉｎｆｅｃ…     

A computational model of a bonus-malus system with a variable trend

The paper considers a model of a bonus-malus system with a variable trend in the case of a three-parameter gamma distribution. The parameters of the distribution and the value of the trend are selected based on statistical data. The indexes of the bonus-malus system are obtained.     

Bifurcation analysis of a generalized Gause model with prey harvesting and a generalized Holling response function of type III

In this paper we study a generalized Gause model with prey harvesting and a generalized Holling response function of type III: . The goal of our study is to give the bifurcation diagram of the model. For this we need to study saddle-node bifurcations, Hopf bifurcation of codimension 1 and 2, heteroclinic bifurcation, and nilpotent saddle bifurcation of codimension 2 and 3. The nilpotent saddle of codimension 3 is the organizing center for the bifurcation diagram. The Hopf bifurcation is studied by means of a generalized Liénard system, and for b=0 we discuss the potential integrability of the system. The nilpotent point of multiplicity 3 occurs with an invariant line and can have a codimension up to 4. But because it occurs with an invariant line, the effective highest codimension is 3. We develop normal forms (in which the invariant line is preserved) for studying of the nilpotent saddle bifurcation. For b=0, the reversibility of the nilpotent saddle is discussed. We study the type of the heteroclinic loop and its cyclicity. The phase portraits of the bifurcations diagram (partially conjectured via the results obtained) allow us to give a biological interpretation of the behavior of the two species.     

Numerical study of a diffusive epidemic model of influenza with variable transmission coefficient

A diffusive epidemic model for H1N1 influenza is formulated with a view to gain basic understanding of the virus behavior. All newborns are assumed to be susceptible. Mortality rate for infective individuals in the population is assumed to be greater than natural mortality rate. Latent, infectious and immune periods are assumed to be constants throughout this study. The numerical solutions of this model are carried out under three different initial populations distribution. In order to investigate the effect of the disease transmission coefficient on the spread of disease, β is taken to be constant as well as a function of seasonally varying time t and a function of spatial variable x  . The threshold quantity (_R0)$({\mathfrak{R}}_0)$ that governs the disease dynamics is derived. Numerical simulation shows that the system supports the existence of sustained and damped oscillations depending on initial populations distribution, the disease transmission rate and diffusion.