Moving-Particle Semi-Implicit Method for Vortex Patterns and Roll Damping of 2D Ship Sections

A meshless method, Moving-Particle Semi-hnplicit Method （MPS） is presented in this paper to simulate the rolling of different 2D ship sections. Sections S. S. 0.5, S.S. 5.0 and S. S. 7.0 of series 60 with CB = 0.6 are chosen for the simulation. It shows that the result of MPS is very close to results of experiments or mesh-numerical simulations. In the simulation of MPS, vortices are found periodically in bilges of ship sections. In section S. S. 5.0 and section S. S. 7.0, which are close to the middle ship, two little vortices are found at different bilges of the section, in section S. S. 0.5, which is close to the bow, only one big vortex is found at the bottom of the section, these vortices patterns are consistent with the theory of Ikeda. The distribution of shear stress and pressure on the rolling hull of ship section is calculated. When vortices are in bilges of the section, the sign clmnge of pressure can be found, but in section S. S. 0.5, there is no sign change of pressure because only one vortex in the bottom of the section. With shear stress distribution, it can be found the shear stress in bilges is bigger than that at other part of the ship section. As the free surface is considered, the shear stress of both sides near the free surface is close to zero and even sign changed.

Stability of average acceleration method for structures with nonlinear damping

The energy approach is used to theoretically verify that the average acceleration method （AAM）, which is unconditionally stable for linear dynamic systems, is also unconditionally stable for structures with typical nonlinear damping, including the special case of velocity power type damping with a bilinear restoring force model. Based on the energy approach, the stability of the AAM is proven for SDOF structures using the mathematical features of the velocity power function and for MDOF structures by applying the virtual displacement theorem. Finally, numerical examples are given to demonstrate the accuracy of the theoretical analysis.

Accuracy of the half-power bandwidth method with a third-order correction for estimating damping in multi-DOF systems

A third-order correction was recently suggested to improve the accuracy of the half-power bandwidth method in estimating the damping of single DOF systems.This paper analyzes the accuracy of the half-power bandwidth method with the third-order correction in damping estimation for multi-DOF linear systems.Damping ratios in a two-DOF linear system are estimated using its displacement and acceleration frequency response curves,respectively.A wide range of important parameters that characterize the shape of these response curves are taken into account.Results show that the third-order correction may greatly improve the accuracy of the half-power bandwidth method in estimating damping in a two-DOF system.In spite of this,the half-power bandwidth method may significantly overestimate the damping ratios of two-DOF systems in some cases.

A time domain induced polarization relaxation time spectrum inversion method based on a damping factor and residual correction

Relaxation time spectra （RTS） derived from time domain induced polarization data （TDIP） are helpful to assess oil reservoir pore structures. However, due to the sensitivity to the signal-to-noise ratio （SNR）, the inversion accuracy of the traditional singular value decomposition （SVD） inversion method reduces with a decrease of SNR. In order to enhance the inversion accuracy and improve robustness of the inversion method to the SNR, an improved inversion method, based on damping factor and spectrum component residual correction, is proposed in this study. The numerical inversion results show that the oscillation of the RTS derived from the SVD method increased with a decrease of SNR, which makes it impossible to get accurate inversion components. However, the SNR has little influence on inversion components of the improved method, and the RTS has high inversion accuracy and robustness. Moreover, RTS derived from core sample data is basically in accord with the pore-size distribution curve, and the RTS derived from the actual induced polarization logging data is smooth and continuous, which indicates that the improved method is practicable.

Resonance and Bifurcation of Fractional Nonlinear Systems with Power Damping Term for Robot Grinding

A fractional nonlinear system with power damping term is introduced to study the forced vibration system in order to solve the resonance and bifurcation problems between grinding wheel and steel bar during robot grinding.The robot,grinding wheel and steel bar are reduced to a spring-damping second-order system model.The implicit function equations of vibration amplitude of the dynamic system with coulomb friction damping,linear damping,square damping and cubic damping are obtained by average method.The stability of the system is analyzed and explained,and the stability condition of the system is proposed.Then,the amplitude-frequency characteristic curves of the system under different fractional differential orders,nonlinear stiffness parameters,fractional differential term coefficients and external excitation amplitude are analyzed.It is shown that the fractional differential term in the dynamic system is the damping characteristic.Then the influence of four kinds of damping on the vibration amplitude of the system under the same parameter is investigated and it is proved that the cubic damping suppresses the vibration of the system to the maximum extent.Finally,based on the idea that the equilibrium point of the system is the constant part of the Fourier series expansion term,the bifurcation behavior caused by the change of damping parameters in linear damping,square damping and cubic damping systems with different values of fractional differential order is investigated.

Force Control Compensation Method with Variable Load Stiffness and Damping of the Hydraulic Drive Unit Force Control System

Each joint of hydraulic drive quadruped robot is driven by the hydraulic drive unit（HDU）, and the contacting between the robot foot end and the ground is complex and variable, which increases the difficulty of force control inevitably. In the recent years, although many scholars researched some control methods such as disturbance rejection control, parameter self-adaptive control, impedance control and so on, to improve the force control performance of HDU, the robustness of the force control still needs improving. Therefore, how to simulate the complex and variable load characteristics of the environment structure and how to ensure HDU having excellent force control performance with the complex and variable load characteristics are key issues to be solved in this paper. The force control system mathematic model of HDU is established by the mechanism modeling method, and the theoretical models of a novel force control compensation method and a load characteristics simulation method under different environment structures are derived, considering the dynamic characteristics of the load stiffness and the load damping under different environment structures. Then, simulation effects of the variable load stiffness and load damping under the step and sinusoidal load force are analyzed experimentally on the HDU force control performance test platform, which provides the foundation for the force control compensation experiment research. In addition, the optimized PID control parameters are designed to make the HDU have better force control performance with suitable load stiffness and load damping, under which the force control compensation method is introduced, and the robustness of the force control system with several constant load characteristics and the variable load characteristics respectively are comparatively analyzed by experiment. The research results indicate that if the load characteristics are known, the force control compensation method presented in this paper has positive compensation effects on the load characteristics variation, i.e., this method decreases the effects of the load characteristics variation on the force control performance and enhances the force control system robustness with the constant PID parameters, thereby, the online PID parameters tuning control method which is complex needs not be adopted. All the above research provides theoretical and experimental foundation for the force control method of the quadruped robot joints with high robustness.

Numerical Research on The Nozzle Damping Effect by A Wave Attenuation Method

Nozzle damping is one of the most important factors in the suppression of combustion instability in solid rocket motors.For an engineering solid rocket motor that experiences combustion instability at the end of burning,a wave attenuation method is proposed to assess the nozzle damping characteristics numerically.In this method,a periodic pressure oscillation signal which frequency equals to the first acoustic mode is superimposed on a steady flow at the head end of the chamber.When the pressure oscillation is turned off,the decay rate of the pressure can be used to determine the nozzle attenuation constant.The damping characteristics of three other nozzle geometries are numerically studied with this method under the same operating condition.The results show that the convex nozzle provides more damping than the conical nozzle which in turn provides more damping than the concave nozzle.All the three nozzles have better damping effect than that of basic nozzle geometry.At last,the phase difference in the chamber is analyzed,and the numerical pressure distribution satisfies well with theoretical distribution.

A CLASS OF r-POINT (r+1)ST-ORDER A-STABLE ONE-BLOCK METHODS WITH DAMPING AT THE INFINITE POINT

This paper presents a class of r-point (r+1)st-order A-stable one-block methods with damping at the infinite point (DIAOB r,r+1 ). Under the conditions of the same order, A-stability, operation count (at each iterative step) and storage space are the same as the methods in , the methods in the paper improve the stability in a neighborhood at the infinite point. And, by using the OOPI method , it possesses much faster rate of convergence for solving systems of nonlinear equations produced by the DIAOB r,r+1 .

A Damping Characteristics Calculation Method of Metal Dry Friction Isolators

The dry friction ring-type vibration isolator is considered as an isotropic continuous medium. A method of dry friction hysteresis loop calculation is proposed based on friction force analysis of contact beam. The friction force is modeled as an equivalent distributed moment to use the finite element method （FEM） to calculate the dry friction vibration isolator hysteresis loop, so the damping characteristics can be obtained. A comparison of the hysteresis loop calculation results and the experimental results shows the average relative error is 2.7 %, it proves the calculation method is feasible.

Sound Radiation Analysis of Damping Structure Based on Response Surface Method

The methods of modifying dimension and shape of structure, or covering damping material are effective to reduce structure-borne noise, while these methods are based on the knowledge of qualitative and quantitative relationship between sound radiation and design parameters. In order to decrease the complexity of the problem, response surface method(RSM) was utilized to analyze and optimize the vibro-acoustic properties of the damping structure. A simple case was illustrated to demonstrate the capabilities of the developed procedure. A structure-born noise problem was approximated by a series of polynomials using RSM. Three main performances were considered, i.e. sound radiation power, first order modal frequency and total mass. Consequently, the response surface model not only gives the direction of design modification, it also leads to an optimal design of complex systems.

Mixed Finite Volume Element Method for Vibration Equations of Beam with Structural Damping

In this paper, for the initial and boundary value problem of beams with structural damping, by introducing intermediate variables, the original fourth-order problem is transformed into second-order partial differential equations, and the mixed finite volume element scheme is constructed, and the existence, uniqueness and convergence of the scheme are analyzed. Numerical examples are provided to confirm the theoretical results. In the end, we test the value of δ to observe its influence on the model.

THE AUTO-ADJUSTABLE DAMPING METHOD FORSOLVING NONLINEAR EQUATIONS

The general approach for solving the nonlinear equations is linearizing the equations and forming various iterative procedures, then executing the numerical simulation. For the strongly nonlinear problems, the solution obtained in the iterative process is always difficult, even divergent due to the numerical instability. It can not fulfill the engineering requirements. Newton's method and its variants can not settle this problem. As a result, the application of numerical simulation for the strongly nonlinear problems is limited. An auto-adjustable damping method has been presented in this paper. This is a further improvement of Newton's method with damping factor. A set of vector of damping factor is introduced. This set of vector can be adjusted continuously during the iterative process in accordance with the judgement and adjustment. An effective convergence coefficient and quichening coefficient are employed to relax the restricted requirements for the initial values and to shorten the iterative process. Then, the numerical stability will be ensured for the solution of complicated strongly nonlinear equations. Using this method, some complicated strongly nonlinear heat transfer problems in airplanes and aeroengines have been numerically simulated successfully. It can be used for the numerical simulation of strongly nonlinear problems in engineering such as nonlinear hydrodynamics and aerodynamics, heat transfer and structural dynamic response etc.

Variable damping constraint tomography and its application in VSP Data

First arrival travel time tomography has achieved wide application. However, tomographic resolution is insufficient because geometry constraints cause rays to be unevenly distributed in the velocity model. The variable damping constraint method adopts uneven priori information to match uneven data distribution which can lessen the correlation between velocity correction values and ray coverage density. In this paper, we combine the variable damping constraint with a smoothness constraint which is added into the regularization equations in velocity inversion to avoid instability caused by only using the variable damping constraint method. The alpha-trimmed-mean filter is used to smooth and denoise intermediate results in the velocity inversion process. We use the LSQR algorithm to enhance the convergence rate and suppress error propagation in solving linear equations. In this paper, we apply the proposed tomographic method to perform velocity inversion using VSP data. The application in recovery test of the checkerboard model and velocity inversion of real VSP data show that the variable damping constraint method can improve tomographic quality because it can solve the effects of uneven ray coverage. In addition, the examples show that the tomographic result near geophones is much more reliable than other areas in the velocity model.

Suboptimal Rayleigh damping coefficients in seismic analysis of viscously-damped structures

An optimization method for the consistent evaluation of two Rayleigh damping coefficients is proposed. By minimizing an objective function such as an error term of the peak displacement of a structure, the two coefficients can be determined with response spectral analysis. The optimization method degenerates into the conventional method used in current practices when only two modes of vibration are included in the objective function. Therefore, the proposed method with all significant modes included for simplicity in practical applications results in suboptimal damping coefficients. The effects of both spatial distribution and frequency content of excitations as well as structural dynamic characteristics on the evaluation of Rayleigh damping coefficients were investigated with a five-story building structure. Two application examples with a 62-story high-rise building and a 840 m long cable-stayed bridge under ten earthquake excitations demonstrated the accuracy and effectiveness of the proposed method to account for all of the above effects.

Implementation of configuration dependent stiffness proportional damping for the dynamics of rigid multi-block systems

The distinct element method(DEM)has been used successfully for the dynamic analysis of rigid block sys- tems.One of many difficulties associated with DEM is modeling of damping.In this paper,new procedures are proposed for the damping modeling and its numerical implementation in distinct element analysis of rigid muhi-block systems.The stiff- ness proportional damping is constructed for the prescribed damping ratio,based on the non-zero fundamental frequency ef- fective during the time interval while the boundary conditions remain essentially constant.At this time interval,the funda- mental frequency can be estimated without complete eigenvalue analysis.The damping coefficients will vary while the damp- ing ratio remains the same throughout the entire analysis.A new numerical procedure is developed to prevent unnecessary energy loss that can occur during the separation phases.These procedures were implemented in the development of the dis- tinet element method for the dynamic analyses of piled multi-block systems.The analysis results |or the single-block and two-block systems were in a good agreement with the analytic predictions.Applications to the seismic analyses of piled four- block systems revealed that the new procedures can make a significant difference and may lead to much-improved results.

Anharmonic 1D actuator model including electrostatic and Casimir forces with fractional damping perturbed by an external force

We modeled a one-dimensional actuator including the Casimir and electrostatic forces perturbed by an external force with fractional damping. The movable electrode was assumed to oscillate by an anharmonic elastic force originated from Murrell- Mottram or Lippincott potential. The nonlinear equations have been solved via the Adomian decomposition method. The behavior of the displacement of the electrode from equilibrium position, its velocity and acceleration were described versus time. Also, the changes of the displacement have been investigated according to the frequency of the external force and the voltage of the electrostatic force. The convergence of the Adomian method and the effect of the orders of expansion on the displacement versus time, frequency, and voltage were discussed. The pull-in parameter was obtained and compared with the other models in the literature. This parameter was described versus the equilibrium position and anharmonicity constant.

Equivalent damping of SDOF structure with Maxwell damper

To predict the maximum earthquake response of an SDOF structure with a Maxwell fluid damper or supplemental brace-viscous damper system using the seismic design response spectrum technique,a new approach is presented to determine the first-and second-order equivalent viscous damping and stiffness,the peak responses,and the damper force of the above structure.Based on the fact that the dynamic characteristics of a general linear viscoelastically damped structure are fully determined by its free vibration properties and the relaxation time constants of a Maxwell fluid damper and supplemental brace-viscous damper system in engineering practice are all small,the method of improved multiple time scales and the equivalent criterion in which all free vibration properties are the same are used to obtain the first-and second-order equivalent viscous damping and stiffness of the above structure in closed form.The accuracy of the proposed method is higher and significantly better than that of the modal strain energy method.Furthermore,in the parametric range of the requirements of the Chinese ＂Code for Seismic Design of Buildings＂,the error of the proposed second-order equivalent system for the abovementioned engineering structure is not more than 0.5%.

Dynamic analysis of a flexible rotor supported by ball bearings with damping rings based on FEM and lumped mass theory

A dynamic model of a flexible rotor supported by ball bearings with rubber damping rings was proposed by combining the finite element and the mass-centralized method.In the proposed model,the rotor was built with the Timoshenko beam element,while the supports and bearing outer rings were modelled by the mass-centralized method.Meanwhile,the influences of the rotor’s gravity,unbalanced force and nonlinear bearing force were considered.The governing equations were solved by precise integration and the Runge-Kutta hybrid numerical algorithm.To verify the correctness of the modelling method,theoretical and experimental analysis is carried out by a rotor-bearing test platform,where the error rate between the theoretical and experimental studies is less than 10%.Besides that,the influence of the rubber damping ring on the dynamic properties of the rotor-bearing coupling system is also analyzed.The conclusions obtained are in agreement with the real-world deployment.On this basis,the bifurcation and chaos behaviors of the coupling system were carried out with rotational speed and rubber damping ring’s stiffness.The results reveal that as rotational speed increases,the system enters into chaos by routes of crisis,quasi-periodic and intermittent bifurcation.However,the paths of crisis,quasi-periodic bifurcation,and Hopf bifurcation to chaos were detected under the parameter of rubber damping ring’s stiffness.Additionally,the bearing gap affects the rotor system’s dynamic characteristics.Moreover,the excessive bearing gap will make the system’s periodic motion change into chaos,and the rubber damping ring’s stiffness has a substantial impact on the system motion.

Numerical analysis of added mass and damping of floating production,storage and offloading system

An integral equation approach is utilized to in- vestigate the added mass and damping of floating produc- tion, storage and offloading system （FPSO system）. Finite water depth Green function and higher-order boundary ele- ment method are used to solve integral equation. Numeri- cal results about added mass and damping are presented for odd and even mode motions of FPSO. The results show ro- bust convergence in high frequency range and can be used in wave load analysis for FPSO designing and operation.

Nonlinear performance of classical damping

The performance of a classical damping matrix, constructed either from the use of initial structural properties or current structural properties, in the step-by-step solution of a nonlinear multiple degree of freedom （MDOF） system is analytically evaluated. The analytical results are confirmed by numerical examples. Consequently, some conclusions are drawn from these analytical results that might be considered as rough guidelines for practical applications. It is found that a classical damping matrix constructed from initial structural properties is adequate for practical applications, since it has approximately the same damping effect as obtained by current structural properties and is more efficient in terms of computing.