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 solutio...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.展开更多
Ocean energy has progressively gained considerable interest due to its sufficient potential to meet the world’s energy demand,and the blade is the core component in electricity generation from the ocean current.Howev...Ocean energy has progressively gained considerable interest due to its sufficient potential to meet the world’s energy demand,and the blade is the core component in electricity generation from the ocean current.However,the widened hydraulic excitation frequency may satisfy the blade resonance due to the time variation in the velocity and angle of attack of the ocean current,even resulting in blade fatigue and destructively interfering with grid stability.A key parameter that determines the resonance amplitude of the blade is the hydrodynamic damping ratio(HDR).However,HDR is difficult to obtain due to the complex fluid-structure interaction(FSI).Therefore,a literature review was conducted on the hydrodynamic damping characteristics of blade-like structures.The experimental and simulation methods used to identify and obtain the HDR quantitatively were described,placing emphasis on the experimental processes and simulation setups.Moreover,the accuracy and efficiency of different simulation methods were compared,and the modal work approach was recommended.The effects of key typical parameters,including flow velocity,angle of attack,gap,rotational speed,and cavitation,on the HDR were then summarized,and the suggestions on operating conditions were presented from the perspective of increasing the HDR.Subsequently,considering multiple flow parameters,several theoretical derivations and semi-empirical prediction formulas for HDR were introduced,and the accuracy and application were discussed.Based on the shortcomings of the existing research,the direction of future research was finally determined.The current work offers a clear understanding of the HDR of blade-like structures,which could improve the evaluation accuracy of flow-induced vibration in the design stage.展开更多
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 a...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.展开更多
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 ...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.展开更多
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 m...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.展开更多
The first order differential matrix equations of the host shell and constrained layer for a sandwich rotational shell are derived based on the thin shell theory.Employing the layer wise principle and first order shear...The first order differential matrix equations of the host shell and constrained layer for a sandwich rotational shell are derived based on the thin shell theory.Employing the layer wise principle and first order shear deformation theory, only considering the shearing deformation of the viscoelastic layer, the integrated first order differential matrix equation of a passive constrained layer damping rotational shell is established by combining with the normal equilibrium equation of the viscoelastic layer.A highly precise transfer matrix method is developed by extended homogeneous capacity precision integration technology.The numerical results show that present method is accurate and effective.展开更多
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 con...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.展开更多
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）, th...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.展开更多
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 burni...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.展开更多
The outstanding comprehensive mechanical properties of newly developed hybrid lattice structures make them useful in engineering applications for bearing multiple mechanical loads.Additive-manufacturing technologies m...The outstanding comprehensive mechanical properties of newly developed hybrid lattice structures make them useful in engineering applications for bearing multiple mechanical loads.Additive-manufacturing technologies make it possible to fabricate these highly spatially programmable structures and greatly enhance the freedom in their design.However,traditional analytical methods do not sufficiently reflect the actual vibration-damping mechanism of lattice structures and are limited by their high computational cost.In this study,a hybrid lattice structure consisting of various cells was designed based on quasi-static and vibration experiments.Subsequently,a novel parametric design method based on a data-driven approach was developed for hybrid lattices with engineered properties.The response surface method was adopted to define the sensitive optimization target.A prediction model for the lattice geometric parameters and vibration properties was established using a backpropagation neural network.Then,it was integrated into the genetic algorithm to create the optimal hybrid lattice with varying geometric features and the required wide-band vibration-damping characteristics.Validation experiments were conducted,demonstrating that the optimized hybrid lattice can achieve the target properties.In addition,the data-driven parametric design method can reduce computation time and be widely applied to complex structural designs when analytical and empirical solutions are unavailable.展开更多
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 iterativ...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 .展开更多
The present study investigates the wavespace of Highly Contrasted Structures(HCS)and Highly Dissipative Structures(HDS)by wave-based models.The Asymptotic Homogenization Method(AHM),exploits the asymptotic Zig-Zag mod...The present study investigates the wavespace of Highly Contrasted Structures(HCS)and Highly Dissipative Structures(HDS)by wave-based models.The Asymptotic Homogenization Method(AHM),exploits the asymptotic Zig-Zag model and homogenization technique to compute the bending wavenumbers via a 6th-order equation.The General Laminate Model(GLM)employs Mindlin’s displacement field to establish displacement-constraint relationships and resolves a quadratic Eigenvalue Problem(EVP)of the dispersion relation.The Wave Finite Element(WFE)scheme formulates the Nonlinear Eigenvalue Problem(NEP)for waves in varying directions and tracks complex wavenumbers using Weighted Wave Assurance Criteria(WWAC).Two approaches are introduced to estimate the Damping Loss Factor(DLF)of HDS,with the average DLF calculated by the modal density at various angles where non-homogeneity is present.Evaluation of robustness and accuracy is made by comparing the wavenumbers and DLF obtained from AHM and GLM with WFE.WFE is finally extended to a sandwich metastructure with a non-homogeneous core,and the Power Input Method(PIM)with Finite Element Method(FEM)data is employed to assess the average DLF,demonstrating an enhanced DLF compared to layered configurations with the same material portion,indicating increased energy dissipation due to the bending-shear coupling effects.展开更多
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 fri...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.展开更多
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 rel...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.展开更多
<span style="font-family:Verdana;">In this paper, for the initial and boundary value problem of beams with</span> <span style="font-family:Verdana;">structural damping, by introdu...<span style="font-family:Verdana;">In this paper, for the initial and boundary value problem of beams with</span> <span style="font-family:Verdana;">structural damping, by introducing intermediate variables, the original </span><span style="font-family:Verdana;">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</span><span style="font-family:Verdana;"><span style="font-family:Verdana;"><span style="font-family:Verdana;">.</span></span></span><span><span><span style="font-family:Verdana;"> Numerical examples are provided to confirm the theoretical results. In the end, we test the value of <em>δ</em></span><span style="font-family:Verdana;"> to observe its influence on the model.</span></span></span>展开更多
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 variabl...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.展开更多
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 coeffic...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.展开更多
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 propose...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.展开更多
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 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%.展开更多
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 forc...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.展开更多
文摘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.
基金Supported by the National Natural Science Foundation of China(Nos.52222904 and 52309117)China Postdoctoral Science Foundation(Nos.2022TQ0168 and 2023M731895).
文摘Ocean energy has progressively gained considerable interest due to its sufficient potential to meet the world’s energy demand,and the blade is the core component in electricity generation from the ocean current.However,the widened hydraulic excitation frequency may satisfy the blade resonance due to the time variation in the velocity and angle of attack of the ocean current,even resulting in blade fatigue and destructively interfering with grid stability.A key parameter that determines the resonance amplitude of the blade is the hydrodynamic damping ratio(HDR).However,HDR is difficult to obtain due to the complex fluid-structure interaction(FSI).Therefore,a literature review was conducted on the hydrodynamic damping characteristics of blade-like structures.The experimental and simulation methods used to identify and obtain the HDR quantitatively were described,placing emphasis on the experimental processes and simulation setups.Moreover,the accuracy and efficiency of different simulation methods were compared,and the modal work approach was recommended.The effects of key typical parameters,including flow velocity,angle of attack,gap,rotational speed,and cavitation,on the HDR were then summarized,and the suggestions on operating conditions were presented from the perspective of increasing the HDR.Subsequently,considering multiple flow parameters,several theoretical derivations and semi-empirical prediction formulas for HDR were introduced,and the accuracy and application were discussed.Based on the shortcomings of the existing research,the direction of future research was finally determined.The current work offers a clear understanding of the HDR of blade-like structures,which could improve the evaluation accuracy of flow-induced vibration in the design stage.
基金the National Natural Science Foundation of China (Grant No.50579035)
文摘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.
基金National Natural Science Foundation of ChinaUnder Grant No. 50578047, 50338020 China Ministry ofEducation (Program for New Century Excellent Talents inUniversity) China Ministry of Science and Technology UnderGrant No.2003AA602150
文摘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.
基金National Natural Science Foundation under Grant No. 51179093National Basic Research Program of China under Grant No. 2011CB013602Program for New Century Excellent Talents in University under Grant No.NCET-10-0531
文摘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.
基金supported by the National Natural Science Foundation of China (No.10662003)Educational Commission of Guangxi Province of China (No.200807MS109)
文摘The first order differential matrix equations of the host shell and constrained layer for a sandwich rotational shell are derived based on the thin shell theory.Employing the layer wise principle and first order shear deformation theory, only considering the shearing deformation of the viscoelastic layer, the integrated first order differential matrix equation of a passive constrained layer damping rotational shell is established by combining with the normal equilibrium equation of the viscoelastic layer.A highly precise transfer matrix method is developed by extended homogeneous capacity precision integration technology.The numerical results show that present method is accurate and effective.
基金Supported by National Key Basic Research Program of China(973 Program,Grant No.2014CB046405)State Key Laboratory of Fluid Power and Mechatronic Systems(Zhejiang University)Open Fund Project(Grant No.GZKF-201502)Hebei Military and Civilian Industry Development Funds Projects of China(Grant No.2015B060)
文摘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.
基金supported by a project from the Youth Science Foundation of the National Natural Science Foundation of China (11104089)
文摘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.
文摘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.
基金supported by National Natural Science Foundation of China(Grant No.52375380)National Key R&D Program of China(Grant No.2022YFB3402200)the Key Project of National Natural Science Foundation of China(Grant No.12032018).
文摘The outstanding comprehensive mechanical properties of newly developed hybrid lattice structures make them useful in engineering applications for bearing multiple mechanical loads.Additive-manufacturing technologies make it possible to fabricate these highly spatially programmable structures and greatly enhance the freedom in their design.However,traditional analytical methods do not sufficiently reflect the actual vibration-damping mechanism of lattice structures and are limited by their high computational cost.In this study,a hybrid lattice structure consisting of various cells was designed based on quasi-static and vibration experiments.Subsequently,a novel parametric design method based on a data-driven approach was developed for hybrid lattices with engineered properties.The response surface method was adopted to define the sensitive optimization target.A prediction model for the lattice geometric parameters and vibration properties was established using a backpropagation neural network.Then,it was integrated into the genetic algorithm to create the optimal hybrid lattice with varying geometric features and the required wide-band vibration-damping characteristics.Validation experiments were conducted,demonstrating that the optimized hybrid lattice can achieve the target properties.In addition,the data-driven parametric design method can reduce computation time and be widely applied to complex structural designs when analytical and empirical solutions are unavailable.
文摘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 .
基金supported by the Natural Sciences and Engineering Research Council of Canada-Discovery Grant(individual)Program(No.NSEC-DG#355433-2009)funded by the LabEx CeLyA(Centre Lyonnais d’Acoustique,No.ANR-10-LABX-0060)of Universite?de Lyon。
文摘The present study investigates the wavespace of Highly Contrasted Structures(HCS)and Highly Dissipative Structures(HDS)by wave-based models.The Asymptotic Homogenization Method(AHM),exploits the asymptotic Zig-Zag model and homogenization technique to compute the bending wavenumbers via a 6th-order equation.The General Laminate Model(GLM)employs Mindlin’s displacement field to establish displacement-constraint relationships and resolves a quadratic Eigenvalue Problem(EVP)of the dispersion relation.The Wave Finite Element(WFE)scheme formulates the Nonlinear Eigenvalue Problem(NEP)for waves in varying directions and tracks complex wavenumbers using Weighted Wave Assurance Criteria(WWAC).Two approaches are introduced to estimate the Damping Loss Factor(DLF)of HDS,with the average DLF calculated by the modal density at various angles where non-homogeneity is present.Evaluation of robustness and accuracy is made by comparing the wavenumbers and DLF obtained from AHM and GLM with WFE.WFE is finally extended to a sandwich metastructure with a non-homogeneous core,and the Power Input Method(PIM)with Finite Element Method(FEM)data is employed to assess the average DLF,demonstrating an enhanced DLF compared to layered configurations with the same material portion,indicating increased energy dissipation due to the bending-shear coupling effects.
基金the National Natural Science Foundation of China(50275030)the"111" Project (B07018)
文摘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.
文摘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.
文摘<span style="font-family:Verdana;">In this paper, for the initial and boundary value problem of beams with</span> <span style="font-family:Verdana;">structural damping, by introducing intermediate variables, the original </span><span style="font-family:Verdana;">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</span><span style="font-family:Verdana;"><span style="font-family:Verdana;"><span style="font-family:Verdana;">.</span></span></span><span><span><span style="font-family:Verdana;"> Numerical examples are provided to confirm the theoretical results. In the end, we test the value of <em>δ</em></span><span style="font-family:Verdana;"> to observe its influence on the model.</span></span></span>
基金supported by the China Important National Science and Technology Specific Projects (No2011ZX05024-001-02)
文摘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.
基金National Natural Science Foundation of China under Grant No.51078032the Visiting Scholar Foundation of China Scholarship Councilthe Center for Infrastructure Engineering Studies at Missouri University of Science and Technology
文摘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.
文摘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.
基金National Natural Science Foundation of China under Grant No.51468005 and 51368008Guangxi Natural Science Foundation under Grant No.2014GXNSFAA118315the Innovative Research Team Program of Guangxi University of Science and Technology(2015)
文摘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%.
文摘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.
