[go: up one dir, main page]

CN119087806B - A parameter optimization method for nonlinear shock dampers - Google Patents

A parameter optimization method for nonlinear shock dampers

Info

Publication number
CN119087806B
CN119087806B CN202411200384.1A CN202411200384A CN119087806B CN 119087806 B CN119087806 B CN 119087806B CN 202411200384 A CN202411200384 A CN 202411200384A CN 119087806 B CN119087806 B CN 119087806B
Authority
CN
China
Prior art keywords
coefficient
zero
order term
damping force
nonlinear
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202411200384.1A
Other languages
Chinese (zh)
Other versions
CN119087806A (en
Inventor
吴九林
张炜
邓卓俊
车吉星
姜伟
陈学东
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huazhong University of Science and Technology
Original Assignee
Huazhong University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Huazhong University of Science and Technology filed Critical Huazhong University of Science and Technology
Priority to CN202411200384.1A priority Critical patent/CN119087806B/en
Publication of CN119087806A publication Critical patent/CN119087806A/en
Application granted granted Critical
Publication of CN119087806B publication Critical patent/CN119087806B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Vibration Prevention Devices (AREA)

Abstract

本申请属于振动抗冲控制技术领域,具体公开了一种非线性抗冲阻尼器的参数优化方法,包括:确定非线性抗冲阻尼器的控制目标、基本参数以及载荷和载体受到冲击后的初速度;将控制目标、基本参数以及初速度根据非线性动力学方程获取速度衰减为零时的最大位移,将最大位移结合抗冲行程的理想值构建非线性阻尼系数设计不等式,并确定非线性阻尼系数的零次项系数和一次项系数;基于抗冲行程和阻尼力过渡系数对零次项系数和一次项系数进行优化,直至抗冲行程和阻尼力过渡系数均达到理想值,得到优化后的零次项系数和一次项系数。通过本申请,能够提高非线性抗冲阻尼器的参数优化效率。

The present application belongs to the field of vibration anti-shock control technology and specifically discloses a parameter optimization method for a nonlinear anti-shock damper, including: determining the control target, basic parameters, and initial velocity of the load and carrier after the nonlinear anti-shock damper is subjected to an impact; obtaining the maximum displacement when the velocity decays to zero according to the nonlinear dynamic equation based on the control target, basic parameters, and initial velocity; constructing a nonlinear damping coefficient design inequality based on the maximum displacement and the ideal value of the anti-shock stroke, and determining the zero-order term coefficient and the first-order term coefficient of the nonlinear damping coefficient; optimizing the zero-order term coefficient and the first-order term coefficient based on the anti-shock stroke and the damping force transition coefficient until both the anti-shock stroke and the damping force transition coefficient reach the ideal value, thereby obtaining the optimized zero-order term coefficient and the first-order term coefficient. Through the present application, the parameter optimization efficiency of the nonlinear anti-shock damper can be improved.

Description

Parameter optimization method of nonlinear impact damper
Technical Field
The application belongs to the technical field of vibration impact control, and particularly relates to a parameter optimization method of a nonlinear impact damper.
Background
In the service and precision manufacturing of moving carrier equipment, the impact phenomenon is ubiquitous and a key factor for limiting the service performance and precision manufacturing precision of the equipment. In order to effectively absorb and dissipate the energy of external impact, reduce the adverse effect thereof on the system, various types of impact devices have been developed and applied to practical engineering, and play a very important role in ensuring the performance of the system and improving the safety.
At present, in engineering, a constant design is often adopted for damping coefficients of the damper for impact resistance, but the constant damping coefficients cannot exert the energy consumption effect of the damper to the greatest extent, and the variable damping coefficient type impact resistance damper parameter design has the problems of low universality, high complexity and the like. That is, the current design method may not be suitable for various application scenarios, and the parameter design and optimization process is complicated, and complex parameters and conditions need to be processed, resulting in lower overall efficiency of parameter optimization of the nonlinear impact damper.
Therefore, how to solve the defect of low parameter optimization efficiency caused by low universality and high complexity of the nonlinear impact damper design method in the prior art is a technical problem to be solved currently.
Disclosure of Invention
Aiming at the defects of the prior art, the application aims to provide a parameter optimization method of a nonlinear impact damper, which aims to solve the problem of low parameter optimization efficiency of the existing impact damper.
In a first aspect, the present application provides a method for optimizing parameters of a nonlinear impact damper, comprising:
Determining a control target of the nonlinear impact damper, basic parameters and an initial speed of the load and the carrier after the load and the carrier are impacted, wherein the control target comprises an impact stroke and a damping force transition coefficient, the basic parameters comprise rigidity and load mass, and the initial speed is determined according to an impact signal acquired in the field;
Obtaining maximum displacement when the speed attenuation is zero according to the control target, the basic parameters and the initial speed and the nonlinear dynamics equation, constructing a nonlinear damping coefficient design inequality by combining the maximum displacement with an ideal value of the impact stroke, and determining a zero order coefficient and a first order coefficient of the nonlinear damping coefficient;
And optimizing the zero order coefficient and the primary coefficient based on the impact stroke and the damping force transition coefficient until the impact stroke and the damping force transition coefficient reach ideal values, and obtaining the optimized zero order coefficient and primary coefficient.
The method comprises the steps of obtaining a control target, a basic parameter and an initial speed, systematically and standardizing the design of the impact damper to improve the universality of the design, substituting the control target, the basic parameter and the initial speed into a nonlinear dynamic equation, calculating the maximum displacement, constructing a design inequality, reducing the complexity of an optimization process and the uncertainty of parameter design through an optimization design process, optimizing a zero order coefficient and a first order coefficient until an ideal value is achieved, improving the adaptability and the optimization effect of the parameter design and optimization through an iterative optimization process, simultaneously reducing the difficulty caused by complex design steps, optimizing the parameter through convenient optimization steps, ensuring clear and operable overall design flow, and further improving the parameter design and the optimization efficiency of the impact damper.
Optionally, the optimizing method of the zero-order term coefficient and the first-order term coefficient includes:
Substituting the updated zero-order term coefficient and the updated primary term coefficient into a constructed discrete dynamics equation to obtain displacement response and velocity response, and obtaining a damping force change curve according to the displacement response and the velocity response;
And judging the control target by utilizing the displacement response and damping force change curve, continuously optimizing the nonlinear damping coefficient until the impact stroke and the damping force transition coefficient reach ideal values, and determining the zero-order coefficient and the first-order coefficient after optimization.
Optionally, determining the control target by using the displacement response and the damping force variation curve and continuously optimizing the nonlinear damping coefficient until the impact stroke and the damping force transition coefficient reach ideal values, and determining the optimized zero-order term coefficient and the optimized first-order term coefficient includes:
According to the design inequality, an initial zero order term coefficient and a primary term coefficient are determined, the primary term coefficient is fixed, the zero order term coefficient is continuously increased, the zero order term coefficient and the primary term coefficient are input into the discrete dynamics equation, a displacement response curve is obtained, and until the maximum value of the displacement response reaches the ideal value of the impact stroke;
Determining a fixed relation between a zero order term coefficient and a primary term coefficient, and under the condition that the fixed relation is kept unchanged, continuously increasing the primary term coefficient and continuously reducing the zero order term coefficient to obtain a speed response curve;
And obtaining a change curve of the damping force along with the displacement by utilizing the displacement response curve and the speed response curve until the damping force transition coefficient reaches an ideal value, and obtaining the optimized zero-order term coefficient and the optimized primary term coefficient.
Optionally, the expression of the nonlinear damping coefficient is:
C=C0+C1x
Where x is the relative displacement of the carrier and load, c is the nonlinear damping coefficient, c 0 is the zero order term coefficient, and c 1 is the first order term coefficient.
Optionally, the damping force transition coefficient is determined according to a ratio of an initial damping force to a maximum damping force at an initial position;
The damping force is determined by the following formula:
wherein F c is the damping force, Is the relative speed of movement between the carrier and the load;
when x=0, the damping force is the initial damping force.
Optionally, the obtaining process of the maximum displacement specifically includes:
acquiring a relative velocity motion equation according to a nonlinear power equation, wherein the relative velocity motion equation is shown in the following formula:
Wherein, the For the relative movement speed at the moment t after being impacted, k is the rigidity, m is the load mass, and v 0 is the initial speed;
When the relative movement speed is zero, the maximum displacement is determined as shown in the following formula:
where t 0 is the time required for the speed decay to zero.
Optionally, the nonlinear dynamical equation is:
Wherein, the Is the relative motion acceleration between the carrier and the load, k is the stiffness, and m is the load mass.
Optionally, the discrete kinetic equation is represented by the formula:
c1nx2(t)+(m+kn2+c0n)x(t)-(2m+C0n)x(t-n)-C1nx(t)x(t-n)+mx(t-2n)=0
where n is the interval of discrete time, x (t) is the relative displacement at that moment, x (t-n) is the relative displacement from before n at that moment, and x (t-2 n) is the relative displacement from before 2n at that moment.
In a third aspect, the application provides an electronic device comprising at least one memory for storing a program, and at least one processor for executing the program stored in the memory, the processor being adapted to perform the method described in the first aspect or any one of the possible implementations of the first aspect when the program stored in the memory is executed.
In a fourth aspect, the present application provides a computer readable storage medium storing a computer program which, when run on a processor, causes the processor to perform the method described in the first aspect or any one of the possible implementations of the first aspect.
In a fifth aspect, the application provides a computer program product which, when run on a processor, causes the processor to perform the method described in the first aspect or any one of the possible implementations of the first aspect.
It will be appreciated that the advantages of the second to fifth aspects may be found in the relevant description of the first aspect, and are not described here again.
In general, the above technical solutions conceived by the present application have the following beneficial effects compared with the prior art:
(1) The application makes the design of the impact damper systemize and standardize by obtaining control targets (such as impact stroke and damping force transition coefficient), basic parameters (rigidity and load quality) and initial speed to promote the universality of the design, calculates the maximum displacement by substituting the control targets, the basic parameters and the initial speed into a nonlinear dynamics equation, constructs a design inequality, reduces the complexity of the optimization process and the uncertainty of the parameter design by optimizing the design process, optimizes the zero order coefficient and the first order coefficient until reaching ideal values, improves the adaptability and the optimization effect of the parameter design and optimization by iterative optimization process, reduces the difficulty brought by complex design steps, optimizes the parameters by convenient optimization steps, has clear and operable overall design flow, and further promotes the parameter design and optimization efficiency of the nonlinear impact damper.
(2) The application can ensure that the impact damper can effectively control the displacement in the impact process in practical application by setting the ideal value of the impact stroke and optimizing the nonlinear damping coefficient, thereby protecting the load and the carrier from overlarge impact. Through the design inequality and the optimization process, the variation characteristic of the damping force can be adjusted, so that the damping force has better transition performance in the impact process.
(3) The nonlinear impact damper parameter optimization method provided by the application has high universality, and can be used for core parameter designs of various types of impact dampers such as output control strategy designs of electric impact dampers and fluid viscosity control rate designs of magnetorheological impact dampers.
Drawings
FIG. 1 is a schematic flow chart of a method for optimizing parameters of a nonlinear impact damper according to an embodiment of the present application;
FIG. 2 is a schematic diagram of a dynamics model of a nonlinear impact damper of an embodiment of the present application;
FIG. 3 is a second flow chart of a method for optimizing parameters of a nonlinear impact damper according to an embodiment of the present application;
FIG. 4 is a graph of the relative displacement response of a nonlinear impact damper provided by an embodiment of the present application after parameter preselection;
FIG. 5 is a graph of the relative velocity response of a nonlinear impact damper provided by an embodiment of the present application after parameter optimization;
FIG. 6 is a graph showing damping force versus displacement after the optimal design of nonlinear impact damper parameters according to an embodiment of the present application;
FIG. 7 is a graph of damping coefficient after the optimal design of nonlinear impact damper parameters provided by the embodiment of the application;
Fig. 8 is a schematic structural diagram of an electronic device according to an embodiment of the present application.
Detailed Description
The present application will be described in further detail with reference to the drawings and examples, in order to make the objects, technical solutions and advantages of the present application more apparent. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the application.
The term "and/or" is used herein to describe an association relationship of associated objects, and means that there may be three relationships, for example, a and/or B, and that there may be three cases where a exists alone, while a and B exist together, and B exists alone. The symbol "/" herein indicates that the associated object is or is a relationship, e.g., A/B indicates A or B.
The terms "first" and "second" and the like in the description and in the claims are used for distinguishing between different objects and not for describing a particular sequential order of objects. For example, the first response message and the second response message, etc. are used to distinguish between different response messages, and are not used to describe a particular order of response messages.
In embodiments of the application, words such as "exemplary" or "such as" are used to mean serving as an example, instance, or illustration. Any embodiment or design described herein as "exemplary" or "e.g." in an embodiment should not be taken as preferred or advantageous over other embodiments or designs. Rather, the use of words such as "exemplary" or "such as" is intended to present related concepts in a concrete fashion.
In the description of the embodiments of the present application, unless otherwise specified, the meaning of "plurality" means two or more, for example, a plurality of processing units means two or more processing units and the like, and a plurality of elements means two or more elements and the like.
Next, the technical scheme provided in the embodiment of the present application is described.
Referring to fig. 1, the present application provides a parameter optimization method of a nonlinear impact damper, comprising:
S101, determining a control target, basic parameters and an initial speed of the nonlinear impact damper after the load and the carrier are impacted, wherein the control target comprises an impact stroke and a damping force transition coefficient, the basic parameters comprise rigidity and load mass, and the initial speed is determined according to an impact signal acquired in the field;
S102, acquiring a maximum displacement when the speed attenuation is zero according to a control target, basic parameters and an initial speed and a nonlinear dynamics equation, constructing a nonlinear damping coefficient design inequality by combining the maximum displacement with an ideal value of an impact stroke, and determining a zero-order term coefficient and a primary term coefficient of the nonlinear damping coefficient;
S103, optimizing the zero order coefficient and the primary coefficient based on the impact stroke and the damping force transition coefficient until the impact stroke and the damping force transition coefficient reach ideal values, and obtaining the optimized zero order coefficient and primary coefficient.
Specifically, it should be noted that the nonlinear characteristic of the impact damper refers to the variation of the damping coefficient with the displacement, which results in nonlinear relationship between the damping force and the velocity.
Referring to fig. 2, control targets, system base parameters and initial velocity of the nonlinear impact damper are determined, two of the main control targets, namely impact stroke and damping force transition coefficient. Basic parameters of a vibration damping impact system include the stiffness (k) and the load mass (m) of the system.
The impact stroke refers to the maximum displacement allowed by the system when impacted, while the damping force transition coefficient describes the smoothness of the damping force as a function of displacement. The damping force transition coefficient η is a ratio of the damping force F c0 at the initial position to the maximum damping force F cmax, and a smaller ratio means that the impact process is flatter.
Stiffness determines how resistant the system is to displacement changes, while load mass affects the system's response to impact. And determining the initial speed of the load and the carrier after being impacted according to the impact signals acquired in the field.
In S102 described above, a nonlinear damping coefficient design inequality is constructed. Firstly, combining the ideal values of the maximum displacement and the impact stroke, and constructing a design inequality of a nonlinear damping coefficient. The zero order coefficient c 0 and the first order coefficient c 1 of the nonlinear damping coefficient are determined, and these coefficients will influence the curve of the damping force as a function of displacement. The relation between the nonlinear damping coefficient c and the zero-order term coefficient c 0 and the first-order term coefficient c 1 is as follows:
C=C0+C1x
Where x is the relative displacement of the carrier and load, c is the nonlinear damping coefficient, c 0 is the zero order term coefficient, and c 1 is the first order term coefficient.
Finally, the damping coefficient is optimized through the step S103. Based on the impact stroke and damping force transition coefficients, the zero order coefficient and the first order coefficient are optimized, and the coefficients are continuously adjusted through an iteration process in the embodiment so as to meet the design requirements. And through the optimization process, the impact stroke and the damping force transition coefficient reach ideal values. The resulting optimized zero order coefficient and one order coefficient will be used to design a nonlinear impact damper to ensure that the performance of the system upon impact meets a predetermined control objective.
Further, the optimization method of the zero-order term coefficient and the first-order term coefficient comprises the following steps:
Substituting the updated zero-order term coefficient and the updated primary term coefficient into a constructed discrete dynamics equation to obtain displacement response and velocity response, and obtaining a damping force change curve according to the displacement response and the velocity response;
And judging the control target by utilizing the displacement response and damping force change curve, continuously optimizing the nonlinear damping coefficient until the impact stroke and the damping force transition coefficient reach ideal values, and determining the zero-order coefficient and the first-order coefficient after optimization.
Further, judging the control target by using the displacement response and damping force change curve and continuously optimizing the nonlinear damping coefficient until the impact stroke and the damping force transition coefficient reach ideal values, and determining the optimized zero-order term coefficient and the optimized primary term coefficient, including:
According to the design inequality, an initial zero order term coefficient and a primary term coefficient are determined, the primary term coefficient is fixed, the zero order term coefficient is continuously increased, the zero order term coefficient and the primary term coefficient are input into the discrete dynamics equation, a displacement response curve is obtained, and until the maximum value of the displacement response reaches the ideal value of the impact stroke;
Determining a fixed relation between a zero order term coefficient and a primary term coefficient, and under the condition that the fixed relation is kept unchanged, continuously increasing the primary term coefficient and continuously reducing the zero order term coefficient to obtain a speed response curve;
And obtaining a change curve of the damping force along with the displacement by utilizing the displacement response curve and the speed response curve until the damping force transition coefficient reaches an ideal value, and obtaining the optimized zero-order term coefficient and the optimized primary term coefficient.
The specific process is as follows:
Firstly, a nonlinear damping coefficient is preliminarily determined according to a design inequality, a first order coefficient is fixed, and a zero order coefficient is gradually increased. And inputting the first order term coefficient and the adjusted zero order term coefficient into a discrete dynamics equation for numerical simulation. And obtaining displacement response curves under different zero-order term coefficients through simulation. And continuously adjusting and optimizing the zero-order term coefficient until the maximum value of the displacement response reaches the ideal value of the impact stroke.
And secondly, determining a fixed relation between the zero-order term coefficient and the first-order term coefficient, wherein the relation is as follows:
A=c0+0.5c1xi
Where A is a fixed value and x i is the ideal value for the impact stroke.
While keeping the above fixed relationship unchanged, the first order term coefficients are increased and the zero order term coefficients are reduced accordingly. Using the adjusted coefficients, they are again input into the discrete kinetic equation, resulting in a velocity response and a displacement response. The damping force is calculated from the velocity response curve and the displacement response curve.
The damping force is determined by the following formula:
wherein F c is the damping force, Is the relative speed of movement between the carrier and the load;
when x=0, the damping force is the initial damping force.
And drawing a curve of the damping force changing along with displacement according to the calculated damping force, and optimizing the damping coefficient until the ideal damping force transition coefficient is reached. And evaluating a damping force curve under the current coefficient, and checking whether an ideal value of the damping force transition coefficient is met. If the damping force transition coefficient does not reach the ideal value, further iterative optimization of the value of the zero order coefficient and the first order coefficient is required. And finally obtaining the optimal zero-order term coefficient and the primary term coefficient which meet the ideal values of the impact stroke and the damping force transition coefficient through an iteration process.
According to the embodiment of the application, through accurately adjusting the nonlinear damping coefficient, the performance of the impact damper can reach the expected impact stroke and damping force transition coefficient when the impact damper is impacted, through an iterative optimization process, the optimal damping coefficient can be found, so that the performance of the system is optimized when the impact damper is impacted, the impact resistance and the response speed of the system are improved, the zero-order coefficient and the first-order coefficient are optimized until the ideal value is reached, through an iterative optimization process, the adaptability and the effect of the design are improved, meanwhile, the difficulty caused by complex design steps is reduced, the parameter optimization is carried out through convenient optimization steps, the whole design flow is clear and operable, and the parameter design and the optimization efficiency of the impact damper are improved.
Optionally, the obtaining process of the maximum displacement specifically includes:
acquiring a relative velocity motion equation according to a nonlinear power equation, wherein the relative velocity motion equation is shown in the following formula:
Wherein, the For the relative movement speed at the moment t after being impacted, k is the rigidity, m is the load mass, and v 0 is the initial speed;
When the relative movement speed is zero, the maximum displacement is determined as shown in the following formula:
where t 0 is the time required for the speed decay to zero.
Optionally, the nonlinear dynamical equation is:
Wherein, the Is the relative motion acceleration between the carrier and the load, k is the stiffness, and m is the load mass.
Further, the maximum displacement is combined with the ideal value of the impact stroke to construct a nonlinear damping coefficient design inequality, as shown in the following formula:
Where x i is the ideal value for the impact stroke, t i is one-fourth of the system reciprocation period, and t 0 is the time required for the velocity decay to zero.
It should be noted that, under the condition that the system is free from damping after being impacted, the system can reciprocate with the inherent period, at the moment of a quarter inherent period t i, the relative displacement is at a peak value, while the existence of nonlinear damping can accelerate the consumption of impact energy, reduce the maximum value x 0 of the relative displacement, and meanwhile, the required time t 0 is smaller than t i.
Optionally, the discrete kinetic equation is represented by the formula:
c1nx2(t)+(m+kn2+c0n)x(t)-(2m+c0n)x(t-n)-c1nx(t)x(t-n)+mx(t-2n)=0
where n is the interval of discrete time, x (t) is the relative displacement at that moment, x (t-n) is the relative displacement from before n at that moment, and x (t-2 n) is the relative displacement from before 2n at that moment.
Referring to fig. 3, a complete flowchart of a parameter optimization method of a nonlinear impact damper according to an embodiment of the present application includes the steps of:
S1, defining ideal values x i and eta i of control targets (an impact stroke x 0 and a damping force transition coefficient eta);
Identifying the load mass m and the rigidity k of the vibration reduction and impact resistance system;
estimating the initial speed v 0 of the load/carrier after being impacted according to the impact signals acquired in the field;
s2, calculating a displacement expression when the speed decay is zero according to a nonlinear dynamics equation;
Constructing a nonlinear damping coefficient design inequality by combining an ideal value x i of a first control target;
Preliminarily selecting nonlinear damping coefficients c 0 and c 1;
S3, c 0 is increased while c 1 is fixed, and the displacement response curve is solved by substituting the built discrete dynamics equation;
determine whether the maximum displacement is near the ideal value of the impact stroke?
Increasing c 1 and decreasing c 0 on the basis of ensuring the sum of c 0 and 0.5c 1xi to be unchanged;
substituting into a discrete dynamics equation to solve the speed and displacement response, so as to obtain a curve of the damping force changing along with the displacement;
Judging whether the damping force transition coefficient is close to the ideal value?
Finally determining a nonlinear damping coefficient;
and finally, combining the output function of the specific type of damper to design the core parameters of the nonlinear impact damper.
Next, the present application will be described in detail with a typical impact vibration reduction system in which the system load mass m=100 kg, the stiffness k=3948n/m, the initial velocity v 0 =0.098m/s after impact, the first control target ideal value x i =5 mm, the second control target ideal value η i =0.3, and the quarter natural period t i =0.25 s.
As shown in fig. 4, a response graph of the relative displacement after the pre-selection of the nonlinear impact damper parameters provided in any of the foregoing embodiments is obtained by referring to the above steps to obtain the design inequality of 0.0025c 1+c0 >1466.5, and after pre-selection of c 0=1096.5、c1 =200000, it can be seen from the response graph of the relative displacement that the first control target x 0 is close to the ideal value x i.
As shown in fig. 5, in the graph of the relative velocity response after the optimization design of the nonlinear impact damper parameters provided in any of the foregoing embodiments, after the parameter optimization is performed with reference to step S3, it is finally determined that c 0=346.5、c1 =500000, and as seen from the relative velocity response curve, the velocity decay of 0.13S after the impact is zero, and then the nonlinear impact damper is slowly moved to the initial equilibrium position by means of the spring restoring force.
As shown in fig. 6, after the damping force versus displacement curve after the nonlinear impact damper parameter optimization design provided in any of the foregoing embodiments is optimized with reference to step S3, the damping force is increased and then decreased in the impact stage, and the ratio of the initial damping force to the maximum damping force is 0.29, which is close to the ideal value of the second control target, so that a gentle impact can be achieved.
As shown in fig. 7, a damping coefficient curve chart after the optimal design of the nonlinear impact damper parameters provided in any one of the foregoing embodiments is shown, and based on the parameter curve, by combining a displacement sensor and a speed sensor, core parameter designs of various types of impact dampers such as force control strategy designs of electric impact dampers and fluid viscosity control rate designs of magnetorheological impact dampers can be realized.
Referring to fig. 8, based on the method in the above embodiment, an electronic device is provided in an embodiment of the present application, which may include a processor 810, a communication interface (Communications Interface) 820, a memory 830, and a communication bus 840, where the processor 810, the communication interface 820, and the memory 830 complete communication with each other through the communication bus 840. The processor 810 may invoke logic instructions in the memory 830 to perform the methods of the embodiments described above.
Further, the logic instructions in the memory 830 described above may be implemented in the form of software functional units and may be stored in a computer-readable storage medium when sold or used as a stand-alone product. Based on this understanding, the technical solution of the present application may be embodied essentially or in a part contributing to the prior art or in a part of the technical solution, in the form of a software product stored in a storage medium, comprising several instructions for causing a computer device (which may be a personal computer, a server, a network device, etc.) to perform all or part of the steps of the method according to the embodiments of the present application.
Based on the method in the above embodiment, the embodiment of the present application provides a computer-readable storage medium storing a computer program, which when executed on a processor, causes the processor to perform the method in the above embodiment.
Based on the method in the above embodiments, an embodiment of the present application provides a computer program product, which when run on a processor causes the processor to perform the method in the above embodiments.
It is to be appreciated that the processor in embodiments of the present application may be a central processing unit (central processing unit, CPU), other general purpose processor, digital signal processor (DIGITAL SIGNAL processor, DSP), application Specific Integrated Circuit (ASIC), field programmable gate array (field programmable GATE ARRAY, FPGA) or other programmable logic device, transistor logic device, hardware components, or any combination thereof. The general purpose processor may be a microprocessor, but in the alternative, it may be any conventional processor.
The steps of the method in the embodiment of the present application may be implemented by hardware, or may be implemented by executing software instructions by a processor. The software instructions may be comprised of corresponding software modules that may be stored in random access memory (random access memory, RAM), flash memory, read-only memory (ROM), programmable ROM (PROM), erasable programmable ROM (erasable PROM, EPROM), electrically Erasable Programmable ROM (EEPROM), registers, hard disk, removable disk, CD-ROM, or any other form of storage medium known in the art. An exemplary storage medium is coupled to the processor such the processor can read information from, and write information to, the storage medium. In the alternative, the storage medium may be integral to the processor. The processor and the storage medium may reside in an ASIC.
In the above embodiments, it may be implemented in whole or in part by software, hardware, firmware, or any combination thereof. When implemented in software, may be implemented in whole or in part in the form of a computer program product. The computer program product includes one or more computer instructions. When loaded and executed on a computer, produces a flow or function in accordance with embodiments of the present application, in whole or in part. The computer may be a general purpose computer, a special purpose computer, a computer network, or other programmable apparatus. The computer instructions may be stored in or transmitted across a computer-readable storage medium. The computer instructions may be transmitted from one website, computer, server, or data center to another website, computer, server, or data center by a wired (e.g., coaxial cable, fiber optic, digital Subscriber Line (DSL)), or wireless (e.g., infrared, wireless, microwave, etc.). The computer readable storage medium may be any available medium that can be accessed by a computer or a data storage device such as a server, data center, etc. that contains an integration of one or more available media. The usable medium may be a magnetic medium (e.g., floppy disk, hard disk, tape), an optical medium (e.g., DVD), or a semiconductor medium (e.g., solid State Drive (SSD)), etc.
It will be appreciated that the various numerical numbers referred to in the embodiments of the present application are merely for ease of description and are not intended to limit the scope of the embodiments of the present application.
It will be readily appreciated by those skilled in the art that the foregoing description is merely a preferred embodiment of the application and is not intended to limit the application, but any modifications, equivalents, improvements or alternatives falling within the spirit and principles of the application are intended to be included within the scope of the application.

Claims (7)

1.一种非线性抗冲阻尼器的参数优化方法,其特征在于,包括:1. A parameter optimization method for a nonlinear impact damper, comprising: 确定非线性抗冲阻尼器的控制目标、基本参数以及载荷和载体受到冲击后的初速度;所述控制目标包括抗冲行程和阻尼力过渡系数,所述基本参数包括刚度和载荷质量,所述初速度是根据实地采集的冲击信号确定的;Determining the control objectives, basic parameters, and initial velocities of the load and carrier after impact of the nonlinear anti-impact damper; the control objectives include the anti-impact stroke and the damping force transition coefficient; the basic parameters include stiffness and load mass; and the initial velocity is determined based on the impact signal collected in the field; 将所述控制目标、基本参数以及初速度根据非线性动力学方程获取速度衰减为零时的最大位移,将所述最大位移结合所述抗冲行程的理想值构建非线性阻尼系数设计不等式,并确定非线性阻尼系数的零次项系数和一次项系数;The control target, basic parameters, and initial velocity are used in accordance with a nonlinear dynamic equation to obtain a maximum displacement when the velocity decays to zero, the maximum displacement is combined with an ideal value of the anti-impact stroke to construct a nonlinear damping coefficient design inequality, and the zero-order term coefficient and the first-order term coefficient of the nonlinear damping coefficient are determined; 基于所述抗冲行程和阻尼力过渡系数对所述零次项系数和一次项系数进行优化,直至所述抗冲行程和阻尼力过渡系数均达到理想值,得到优化后的零次项系数和一次项系数;Optimizing the zero-order term coefficient and the linear term coefficient based on the anti-impact stroke and the damping force transition coefficient until both the anti-impact stroke and the damping force transition coefficient reach ideal values, thereby obtaining optimized zero-order term coefficient and linear term coefficient; 所述零次项系数和一次项系数的优化方法包括:The optimization method of the zero-order term coefficient and the first-order term coefficient includes: 将更新后的零次项系数和一次项系数代入已构建的离散动力学方程中得到位移响应和速度响应,根据所述位移响应和速度响应得到阻尼力变化曲线;Substituting the updated zero-order term coefficient and the first-order term coefficient into the constructed discrete dynamic equation to obtain a displacement response and a velocity response, and obtaining a damping force variation curve according to the displacement response and the velocity response; 利用所述位移响应和阻尼力变化曲线对所述控制目标进行判断并不断优化所述非线性阻尼系数,直至所述抗冲行程和阻尼力过渡系数均达到理想值,确定优化后的零次项系数和一次项系数;The control target is judged by using the displacement response and damping force change curve and the nonlinear damping coefficient is continuously optimized until the anti-impact stroke and the damping force transition coefficient both reach ideal values, and the optimized zero-order term coefficient and first-order term coefficient are determined; 所述非线性阻尼系数的表达式为:The expression of the nonlinear damping coefficient is: 其中,是载体和载荷的相对位移,c是非线性阻尼系数,为零次项系数,为一次项系数;in, is the relative displacement between the carrier and the load, c is the nonlinear damping coefficient, is the coefficient of the zero-order term, is the coefficient of the first-order term; 所述阻尼力过渡系数是根据初始位置的初始阻尼力与最大阻尼力的比值确定的;The damping force transition coefficient is determined based on the ratio of the initial damping force at the initial position to the maximum damping force; 所述阻尼力通过以下公式确定:The damping force is determined by the following formula: 其中,为阻尼力,是载体和载荷之间的相对运动速度;in, is the damping force, is the relative speed between the carrier and the load; 时,阻尼力为初始阻尼力;when When , the damping force is the initial damping force; 不等式如以下公式所示:The inequality is expressed as follows: 其中,为一次项系数,为抗冲行程的理想值,为零次项系数,为系统往复运动周期的四分之一,为速度衰减为零所需要的时间,为载荷质量,为刚度。in, is the coefficient of the first-order term, is the ideal value of the impact stroke, is the coefficient of the zero-order term, is one quarter of the system's reciprocating motion period, is the time required for the velocity to decay to zero, is the load mass, For stiffness. 2.根据权利要求1所述的参数优化方法,其特征在于,利用所述位移响应和阻尼力变化曲线对所述控制目标进行判断并不断优化所述非线性阻尼系数,直至所述抗冲行程和阻尼力过渡系数均达到理想值,确定优化后的零次项系数和一次项系数,包括:2. The parameter optimization method according to claim 1, characterized in that the control target is judged by using the displacement response and damping force change curve and the nonlinear damping coefficient is continuously optimized until the anti-impact stroke and the damping force transition coefficient both reach ideal values, and the optimized zero-order term coefficient and first-order term coefficient are determined, comprising: 根据所述设计不等式,确定初始的零次项系数和一次项系数,将所述一次项系数固定,不断增大所述零次项系数,将所述零次项系数和一次项系数输入至所述离散动力学方程中,得到位移响应曲线,直至所述位移响应的最大值达到所述抗冲行程的理想值;According to the design inequality, determining an initial zero-order term coefficient and a linear term coefficient, fixing the linear term coefficient, continuously increasing the zero-order term coefficient, inputting the zero-order term coefficient and the linear term coefficient into the discrete dynamics equation to obtain a displacement response curve until the maximum value of the displacement response reaches the ideal value of the anti-impact stroke; 确定零次项系数和一次项系数的固定关系,在保持所述固定关系不变的情况下,不断增大所述一次项系数并不断减小所述零次项系数,获取速度响应曲线;Determining a fixed relationship between a zero-order term coefficient and a linear term coefficient, and while maintaining the fixed relationship, continuously increasing the linear term coefficient and continuously decreasing the zero-order term coefficient to obtain a speed response curve; 利用位移响应曲线和速度响应曲线得到阻尼力随位移的变化曲线,直至所述阻尼力过渡系数达到理想值,得到优化后的零次项系数和一次项系数。The displacement response curve and the velocity response curve are used to obtain a curve of the damping force varying with the displacement, until the damping force transition coefficient reaches an ideal value, thereby obtaining an optimized zero-order term coefficient and a first-order term coefficient. 3.根据权利要求1所述的参数优化方法,其特征在于,所述最大位移的获取过程具体包括:3. The parameter optimization method according to claim 1, wherein the process of obtaining the maximum displacement specifically comprises: 根据非线性动力方程获取相对速度运动方程;如以下公式所示:The relative velocity motion equation is obtained based on the nonlinear dynamic equation; as shown in the following formula: 其中,为受到冲击后时刻的相对运动速度,为刚度,为载荷质量,为初速度;in, After the impact The relative speed of movement at any moment, is the stiffness, is the load mass, is the initial velocity; 当相对运动速度为零时,确定所述最大位移,如以下公式所示:When the relative motion speed is zero, the maximum displacement is determined as shown in the following formula: 其中,为速度衰减为零所需要的时间。in, The time required for the velocity to decay to zero. 4.根据权利要求1或3所述的参数优化方法,其特征在于,所述非线性动力学方程为:4. The parameter optimization method according to claim 1 or 3, wherein the nonlinear dynamic equation is: 其中,是载体和载荷之间的相对运动加速度,为刚度,为载荷质量。in, is the relative acceleration between the carrier and the load, is the stiffness, is the load mass. 5.根据权利要求2所述的参数优化方法,其特征在于,所述离散动力学方程如下公式所示:5. The parameter optimization method according to claim 2, wherein the discrete dynamics equation is expressed as follows: 其中,为离散时间的间隔,是此刻的相对位移,是距离此刻前的相对位移,是距离此刻2n前的相对位移。in, is a discrete time interval, is the relative displacement at this moment, It's from this moment The relative displacement before is the relative displacement 2 n years ago from this moment. 6.一种电子设备,其特征在于,包括:6. An electronic device, comprising: 至少一个存储器,用于存储计算机程序;at least one memory for storing a computer program; 至少一个处理器,用于执行所述存储器存储的程序,当所述存储器存储的程序被执行时,所述处理器用于执行如权利要求1-5任一所述的方法。At least one processor is used to execute the program stored in the memory. When the program stored in the memory is executed, the processor is used to execute the method according to any one of claims 1 to 5. 7.一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,其特征在于,当所述计算机程序在处理器上运行时,使得所述处理器执行如权利要求1-5任一所述的方法。7. A computer-readable storage medium storing a computer program, wherein when the computer program runs on a processor, the processor is caused to execute the method according to any one of claims 1 to 5.
CN202411200384.1A 2024-08-29 2024-08-29 A parameter optimization method for nonlinear shock dampers Active CN119087806B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202411200384.1A CN119087806B (en) 2024-08-29 2024-08-29 A parameter optimization method for nonlinear shock dampers

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202411200384.1A CN119087806B (en) 2024-08-29 2024-08-29 A parameter optimization method for nonlinear shock dampers

Publications (2)

Publication Number Publication Date
CN119087806A CN119087806A (en) 2024-12-06
CN119087806B true CN119087806B (en) 2025-09-26

Family

ID=93691932

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202411200384.1A Active CN119087806B (en) 2024-08-29 2024-08-29 A parameter optimization method for nonlinear shock dampers

Country Status (1)

Country Link
CN (1) CN119087806B (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113297907A (en) * 2021-04-23 2021-08-24 东南大学 Nonlinear damping identification method based on data driving under pulse excitation
CN115879331A (en) * 2023-03-01 2023-03-31 山东科技大学 Kelvin model-based spring-damping vibration attenuation structure parameter optimization analysis method

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6212466B1 (en) * 2000-01-18 2001-04-03 Yamaha Hatsudoki Kabushiki Kaisha Optimization control method for shock absorber
CN102587530B (en) * 2012-02-17 2013-12-25 华中科技大学 Design method of viscous liquid damper for adjacent building structures
CN118484939A (en) * 2024-05-17 2024-08-13 中国铁道科学研究院集团有限公司铁道建筑研究所 A nonlinear optimization method for mechanical characteristics of vibration isolation elements

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113297907A (en) * 2021-04-23 2021-08-24 东南大学 Nonlinear damping identification method based on data driving under pulse excitation
CN115879331A (en) * 2023-03-01 2023-03-31 山东科技大学 Kelvin model-based spring-damping vibration attenuation structure parameter optimization analysis method

Also Published As

Publication number Publication date
CN119087806A (en) 2024-12-06

Similar Documents

Publication Publication Date Title
CN110121171B (en) Trust prediction method based on exponential smoothing method and gray model
CN110162739B (en) RFFKLMS Algorithm Weight Update Optimization Method Based on Variable Forgetting Factor
CN112462611B (en) Sliding friction modeling method for precision electromechanical system
JP7517448B2 (en) GENERATION METHOD, EVALUATION METHOD, GENERATION DEVICE, EVALUATION DEVICE, GENERATION PROGRAM, AND EVALUATION PROGRAM
CN112446487A (en) Method, device, system and storage medium for training and applying neural network model
WO2024088134A1 (en) Channel prediction method and apparatus
CN119087806B (en) A parameter optimization method for nonlinear shock dampers
CN116992806A (en) Automatic optimization method for large-scale analog integrated circuit based on self-attention mechanism
CN113031435A (en) Digital double-frequency extended state observer and disturbance observation method
CN119620961A (en) Cache management method, electronic device, storage medium and program product
CN119204222A (en) A large model reasoning acceleration method, device and medium
CN109491619A (en) Caching data processing method, device and system
CN110765360B (en) Text topic processing method and device, electronic equipment and computer storage medium
CN116680490A (en) Prefetching method and system for server cache
CN116866127A (en) Method, device, equipment and storage medium for determining weight coefficient of predistorter
CN113298256B (en) Adaptive curve learning method, device, computer equipment and storage medium
CN115618609A (en) Method and device for determining parameters of dynamic vibration absorber, electronic equipment and storage medium
CN114280939A (en) A model predictive control method, device, equipment and storage medium
CN118296765B (en) Motor topological structure optimization method and device based on prediction coefficient analysis
JP6588933B2 (en) Language model construction device, method and program
CN109726413B (en) Method and system for accelerating circuit optimization
CN120276682B (en) Memory management method, system and device for memory database and memory database server
CN111460732A (en) A Construction Method of Nonlinear Model of Planar Motor
CN113626720B (en) Recommendation method, device, electronic device and storage medium based on action pruning
CN118862592B (en) Method for constructing digital model of fan foundation based on temperature field and stress field analysis

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant