Analysis of nonlinear steady state vibration of a multi-degree-of-freedom system using component mode synthesis method

In this paper, an analytical approach for nonlinear forced vibration of a multi-degree-of-freedom system is proposed using the component mode synthesis method. The whole system is divided into some components and a nonlinear modal equation of each component is derived using the free-interface vibration modes. The modal equations of all components and the conjunction conditions are solved simultaneously, and then the modal responses of components are derived. Finally, the dynamic responses of the whole system can be obtained. The degrees of freedom of modal equations can be reduced when the lower vibration modes are only adopted in each component. As a numerical example, a nine-degree-of-freedom system is considered, in which all spring have cubic type nonlinearity. As a result, it is shown that when there are no rigid modes in components, the compliance by the proposed method agrees very well with the exact one even if the lower vibration modes of components are only adopted. The other hand, in the case with rigid modes in components, the compliance has a little error compared with the exact result. It is recognized that the method proposed is very effective in the case without rigid modes in components for the actual application.     

A combined finite element-stiffness equation transfer method for steady state vibration response analysis of structures

An extended finite element transfer matrix method, in combination with stiffness equation transfer, is applied to dynamic response analysis of the structures under periodic excitations. In the present method, the transfer of state vectors from left to right in a combined finite element-transfer matrix (FE-TM) method is changed into the transfer of general stiffness equations of every section from left to right. This method has the advantages of reducing the order of standard transfer equation systems, and minimizing the propagation of round-off errors occurring in recursive multiplication of transfer and point matrices. Furthermore, the drawback that in the ordinary FE-TM method, the number of degrees of freedom on the left boundary be the same on the right boundary, is now avoided. A FESET program based on this method using microcomputers is developed. Finally, numerical examples are presented to demonstrate the accuracy as well as the potential of the proposed method for steady state vibration response analysis of structures.     

Discontinuous Galerkin method coupled with an integral representation for solving the three-dimensional time-harmonic Maxwell equations

In this paper, we present a mathematical and numerical studies of the three-dimensional time-harmonic Maxwell equations. The problem is solved by a discontinuous Galerkin DG method coupled with an integral representation. This study was completed by some numerical tests to justify the effectiveness of the proposed approach. The numerical simulation was done by an iterative solver implemented in FORTRAN.     

Equilibrium state and non-equilibrium steady state in an isolated human system

The principle of increasing entropy （PIE） is commonly considered as a universal physical law tbr natural systems. It also means that a non-equilibrium steady state （NESS） must not appear in any isolated natural systems. Here we experimentally investigate an isolated human social system with a clustering effect. We report that the PIE cannot always hold, and that NESSs can come to appear. Our study highlights the role of human adaptability in the PIE, and makes it possible to study human social systems by using some laws originating from traditional physics.     

A stabilized MLPG method for steady state incompressible fluid flow simulation

In this paper, the meshless local Petrov–Galerkin (MLPG) method is extended to solve the incompressible fluid flow problems. The streamline upwind Petrov–Galerkin (SUPG) method is applied to overcome oscillations in convection-dominated problems, and the pressure-stabilizing Petrov–Galerkin (PSPG) method is applied to satisfy the so-called Babuška–Brezzi condition. The same stabilization parameter τ(τ_SUPG = τ_PSPG) is used in the present method. The circle domain of support, linear basis, and fourth-order spline weight function are applied to compute the shape function, and Bubnov–Galerkin method is applied to discretize the PDEs. The lid-driven cavity flow, backward facing step flow and natural convection in the square cavity are applied to validate the accuracy and feasibility of the present method. The results show that the stability of the present method is very good and convergent solutions can be obtained at high Reynolds number. The results of the present method are in good agreement with the classical results. It also seems that the present method (which is a truly meshless) is very promising in dealing with the convection- dominated problems.     

Topology optimization using the improved element-free Galerkin method for elasticity

The improved element-free Galerkin(IEFG) method of elasticity is used to solve the topology optimization problems.In this method, the improved moving least-squares approximation is used to form the shape function. In a topology optimization process, the entire structure volume is considered as the constraint. From the solid isotropic microstructures with penalization, we select relative node density as a design variable. Then we choose the minimization of compliance to be an objective function, and compute its sensitivity with the adjoint method. The IEFG method in this paper can overcome the disadvantages of the singular matrices that sometimes appear in conventional element-free Galerkin(EFG) method. The central processing unit(CPU) time of each example is given to show that the IEFG method is more efficient than the EFG method under the same precision, and the advantage that the IEFG method does not form singular matrices is also shown.     

Free vibration and buckling analysis of clamped skew sandwich plates by the Galerkin method

Galerkin's variational method has been used in the past by several investigators [1–3] to solve bending problems of clamped skew plates. In this paper the suitability of the Galerkin method for solution of problems of buckling under the action of in-plane forces and of free vibration of skew plates is studied. The method is first applied to investigate the problems for clamped rectangular sandwich plates. After the validity of the method has been established, the method is then extended to analyze similar problems for clamped skew sandwich plates. The governing differential equations for the skew sandwich plates are obtained by transforming the corresponding differential equations in Cartesian coordinates into skew co-ordinates. The parameters considered herein for the buckling and free vibration behaviour of the skew sandwich plates are the aspect ratio of the plate, Poisson's ratio, skew angle and various shearing stiffnesses of the core. Simplicity and quick convergence is the advantage of the method in comparison with other much more laborious numerical methods requiring extensive computer facilities.     

Analytical solution of Everett integral using Lorentzian Preisach function approximation

In this paper, we take in to account the Lorentzian function, already proposed as analytical function approximation in scalar Preisach-type modeling of hysteresis of soft magnetic materials. Preliminarily, we point out the properties of the Lorentzian function and the physical and mathematical meaning of its parameters. Successively, we show how the use of the Lorentzian function approximation allows to solve in complete analytical way the Everett's integral. In particular, we present in the paper the analytical expression in closed form of the first magnetization curve, of the symmetric and non-symmetric minor loops, and of the first- and second-order reversal curves. In addition, we show the use of the complete analytical formulas of the symmetric magnetic loops above-mentioned, applied to a simple identification procedure of the Lorentzian function parameters, by the knowledge of the measured major loop. Finally, in order to show the practical use of the analytical expression found, some computation examples and comparisons with experimental data are shown.     

Normal modes of a continuous system with quadratic and cubic non-linearities

A power series solution is presented for the free vibrations of simply supported beams resting on elastic foundation having quadratic and cubic non-linearities. The time-dependence is assumed harmonic and the problem is posed as a non-linear eigenvalue problem. The spatial variable is transformed into an independent variable that satisfies the boundary conditions. This permits a power series expansion of the beam motion in terms of the new variable. A recurrence relation is obtained from the governing equation and used in conjunction with the Rayleigh energy principle to compute the natural frequencies. The results show that, for a first order approximation, only the lower frequencies and first mode shape are significantly affected by the cubic non-linearity.     

Galerkin method for steady-state response of nonlinear forced vibration of axially moving beams at supercritical speeds

The present paper investigates the steady-state periodic response of an axially moving viscoelastic beam in the supercritical speed range. The straight equilibrium configuration bifurcates in multiple equilibrium positions in the supercritical regime. It is assumed that the excitation of the forced vibration is spatially uniform and temporally harmonic. Under the quasi-static stretch assumption, a nonlinear integro-partial-differential equation governs the transverse motion of the axially moving beam. The equation is cast in the standard form of continuous gyroscopic systems via introducing a coordinate transform for non-trivial equilibrium configuration. For a beam constituted by the Kelvin model, the primary resonance is analyzed via the Galerkin method under the simply supported boundary conditions. Based on the Galerkin truncation, the finite difference schemes are developed to verify the results via the method of multiple scales. Numerical simulations demonstrate that the steady-state periodic responses exist in the transverse vibration and a resonance with a softening-type behavior occurs if the external load frequency approaches the linear natural frequency in the supercritical regime. The effects of the viscoelastic damping, external excitation amplitude, and nonlinearity on the steady-state response amplitude for the first mode are illustrated.     

伽辽金有限元法在超导磁悬浮受力分析中的应用

超导块材和永磁铁相互作用形成悬浮系统。讨论了对超导块材进行网格划分的方法,利用矢量电位法建立了Nd-FeB永磁铁和超导体相互作用的电磁场数学方程。利用伽辽金法形成有限元线性方程组。求解过程中采用Kim模型作为迭代条件,解出超导体电流密度分布,并进行受力分析。     

电各向异性色散介质电磁散射的三维递推卷积-时域有限差分方法分析

根据递推卷积原理，将磁化等离子体的频域介电系数过渡到时域，通过引入时域复数极化率张量和时域复数电位移矢量，得到了磁化等离子体的三维时域有限差分方法迭代式. 为了验证该方法，用它计算了非磁等离子体球的后向雷达散射截面，与移位算子法结果符合很好. 应用该方法计算和分析了磁化等离子体球的电磁波散射，发现其后向散射时域波形明显出现了交叉极化分量. 关键词： 递推卷积 磁化等离子体 电磁散射 时域有限差分方法     

Dynamics and non-equilibrium steady state in a system of coupled harmonic oscillators

A system of two coupled oscillators, each of them coupled to an independent reservoir, is analysed. The analytical solution of the non-rotating wave master equation is obtained in the high-temperature and weak coupling limits. No thermal entanglement is found in the high-temperature limit. In the weak coupling limit the system converges to an entangled non-equilibrium steady state. A critical temperature for the appearance of quantum correlations is found.     

Calculation of polarization using a density functional method with localized charge

A density functional method, which represents the total charge density as a sum of self-consistently determined localized densities, is described. While this approach is generally less accurate than conventional band-structure methods, it offers a relatively simple interpretation of polarization and related properties. The method is illustrated with results for NaCl, MgO, and AlP.     

Using the Galerkin method to construct an analytical solution for the integral equation of a certain diffraction problem

Grodno State University. Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Radiofizika, Vol. 31, No. 11, pp. 1355–1361, November, 1988.     

Nonequilibrium steady state transport of collective-qubit system in strong coupling regime

We investigate the steady state photon transport in a nonequilibrium collective-qubit model. By adopting the noninteracting blip approximation, which is applicable in the strong photon–qubit coupling regime, we describe the essential contribution of indirect qubit–qubit interaction to the population distribution, mediated by the photonic baths. The linear relations of both the optimal flux and noise power with the qubits system size are obtained. Moreover, the inversed power-law style for the finite-size scaling of the optimal photon–qubit coupling strength is exhibited, which is proposed to be universal.     

A general expression of the steady state composition for multicomponent system during ion sputtering

A general expression of the steady state composition during sputtering has been derived based on the preferential sputtering model. The resulting formula has been applied to calculate the steady state concentration in a AgAuPd ternary system and the calculated values are in agreement, at least qualitatively, with the experimental data.     

A method for estimating time-dependent acoustic cross-sections of bubbles and bubble clouds prior to the steady state

Models for the acoustic cross-sections of gas bubbles undergoing steady-state pulsation in liquid have existed for some time. This article presents a theoretical scheme for estimating the cross-sections of single bubbles, and bubble clouds, from the start of insonation onward. In this period the presence of transients can significantly alter the cross-section from the steady-state value. The model combines numerical solutions of the Herring-Keller model with appropriate damping values to calculate the extinction cross-section of a bubble as a function of time in response to a continuous harmonic sound field (it is also shown how the model can be adapted to estimate the time-dependent scatter cross-section). The model is then extended to determine the extinction cross-section area of multiple bubbles of varying population distributions assuming no bubble-bubble interactions. The results have shown that the time taken to reach steady state is dependent on the closeness of the bubble to resonance, and on the driving pressure amplitude. In the response of the population as a whole, the time to reach steady state tends to decrease with increasing values of the driving pressure amplitude; and with the increasing values of the ratio of the numbers of bubbles having radii much larger than resonance to the number of resonant bubbles. The implications of these findings for the use of acoustic pulses are explored.     

Irreversible thermodynamics of the steady state with uniform disequilibrium

It is demonstrated that familiar equations of contemporary irreversible thermodynamics, which are currently believed to require a foundation in statistical microscopic theory, can be derived by purely macroscopic reasoning, provided the concept of uniform disequilibrium is accepted. The consequences may be important in relieving the customary limitation of irreversible thermodynamics to near equilibrium conditions.