CN-122020902-A - Edge plate damping device optimization method based on Monte Carlo and Latin hypercube

CN122020902ACN 122020902 ACN122020902 ACN 122020902ACN-122020902-A

Abstract

The invention provides a flange plate damping device optimization method based on Monte Carlo simulation and Latin hypercube design, which is characterized in that a certain order of blade mode is selected to carry out damper structure optimization design based on blade mode vibration mode, a blade harmonic response frequency range is determined, a dual mode synthesis method is used for processing a blade group substructure, an integral structure matrix is combined through a substructure reduction matrix, a parameter matrix and a two-dimensional contact surface friction model are obtained, a dynamic equation set with dry friction nonlinearity is constructed, a dynamic equation set with dry friction nonlinearity is obtained, a amplitude frequency characteristic curve is solved through a harmonic balance method and an arc length extension method, monte Carlo simulation and Latin hypercube design is compiled, parameter optimization analysis aiming at the best damping effect of a damping device is completed, and an optimal parameter set is obtained to further optimize an original damping structure. The model is processed by the structure-based method, so that the vibration response of the blade group with dry friction damping under simple harmonic excitation can be predicted.

Inventors

ZHOU YUHANG
ZHANG DEJIA
SUN JINGGUO
QI WENKAI
WANG ZHENKUN
JI CHEN
SUN YANBO

Assignees

中国船舶集团有限公司第七〇三研究所

Dates

Publication Date: 20260512
Application Date: 20260202

Claims (7)

1. A method for optimizing a marginal plate damping device based on Monte Carlo and Latin hypercube is characterized by comprising the steps of selecting a certain order vibration mode of a blade to perform damper optimization design, obtaining analysis object parameters and node information, establishing a dynamic equation set with dry friction nonlinearity, solving a blade amplitude-frequency characteristic curve, optimizing analysis of multiple parameters of a damping structure and optimizing an original damping structure, and specifically comprises the following steps: s1, selecting a certain order of blade mode shape based on the blade mode shape to perform damper structure optimization design, and determining a blade resonance response frequency range; s2, processing the sub-structure of the blade group based on a dual mode synthesis method, and reducing a matrix combination overall structure matrix through the sub-structure; S3, constructing a dynamic equation set with dry friction nonlinearity based on the parameter matrix obtained in the S2 and the two-dimensional contact surface friction model; S4, solving an amplitude-frequency characteristic curve through a harmonic balance method and an arc length extension method based on the dynamic equation set with the dry friction nonlinearity obtained in the step S3; S5, based on Monte Carlo simulation and Latin hypercube design, writing a multi-parameter Monte Carlo simulation analysis program, and completing parameter optimization analysis aiming at the optimal damping effect of the damping device; S6, further optimizing the original damping structure based on the optimal parameter set obtained in the S4.
2. The optimization method of the flange plate damping device based on the hypercube of Monte Carlo and Latin as set forth in claim 1, wherein in the optimization design of the damper by selecting a certain order of vibration mode of the blade in S1, the blade modal analysis is performed by establishing a single-blade finite element model.
3. The optimization method of the flange plate damping device based on the hypercube of Monte Carlo and Latin according to claim 1, wherein the parameter matrix comprises the mass and rigidity matrix of the blade and the damper in the parameter matrix and the node information of the analysis object, a finite element model of two blades and the damper is established through finite element analysis software, the parameter matrix and the node information of the substructure are derived and recombined, and the finite element analysis software adopts ANSYS APDL.
4. The optimization method of the flange plate damping device based on the hypercube of Monte Carlo and Latin according to claim 1, wherein in the step S3, a dynamics equation set with dry friction nonlinearity is established, wherein the dynamics equation set is a differential equation set, and the differential equation set is mapped to a nonlinear algebraic equation set under a frequency domain through Fourier Galerkin mapping.
5. The optimization method of the flange plate damping device based on the hypercube of Monte Carlo and Latin according to claim 1, wherein in the S4 solving of the amplitude-frequency characteristic curve of the blade, the solution under the initial frequency is solved by a harmonic balance method, wherein the nonlinear force and jacobian matrix solving method is to establish a functional relation between the normal force and tangential force on the contact surface and the relative displacement through a single-point two-dimensional friction model, then discretizing the contact surface to divide a plurality of contact units, obtaining the total contact force and the contact stiffness by superposing the contact force and the contact stiffness of each contact unit in actual calculation, and solving the solution under the initial frequency through a Newton iteration method. And finally, solving an amplitude-frequency characteristic curve by an arc length extension method.
6. The optimization method of the flange plate damping device based on the Monte Carlo and Latin hypercube is characterized in that in the multi-parameter optimization analysis, a turbine blade set vibration response calculation program based on a model reduction method is firstly compiled through S2-S4, input parameters of the turbine blade set vibration response calculation program comprise friction factors, modal damping ratios, left and right flange plate angles, contact coefficients, damping block design quality, installation radius, rotor rotating speed, contact pair number, harmonic balance method harmonic retention number and periodic sampling point number, and then a multi-parameter Monte Carlo simulation analysis program is compiled by combining Monte Carlo simulation and Latin hypercube design. The multiparameter Monte Carlo simulation analysis program is written in MATLAB language.
7. The optimization method of the flange plate damping device based on the Monte Carlo and Latin hypercube of claim 1, wherein in the multi-parameter optimization analysis, the design parameters of the flange plate damping device are friction factors, left and right flange plate angle angles, contact stiffness and damping fin design quality.

Description

Edge plate damping device optimization method based on Monte Carlo and Latin hypercube Technical Field The invention belongs to the technical field of gas turbines, and particularly relates to an optimization method of a flange plate damping device based on Monte Carlo and Latin hypercube. Background The problem of fatigue fracture caused by blade vibration is one of the core problems of modern marine gas turbine blade strength designs. In order to improve the reliability and availability of the gas turbine and the development requirement of the blades with wide temperature range and long service life, the dynamic characteristic analysis and vibration reduction design of the gas turbine blades are required. The dry friction damping has the advantages of good vibration damping performance, simple structure, high reliability and the like, and has wide application in the related fields. The turbine rotor blade is taken as one of key parts and is arranged in high-temperature high-pressure airflow, the vibration problem is particularly remarkable, and the vibration is more likely to cause structural fatigue and fracture, so that the vibration suppression of the turbine blade is particularly urgent. The flange plate damping device is taken as an important part widely applied to the vibration reduction of turbine blades, and high requirements are put on the materials and the shapes of the damping device. Disclosure of Invention The invention aims to provide an optimization method of a flange plate damping device based on Monte Carlo and Latin hypercube. A method for optimizing a marginal plate damping device based on Monte Carlo and Latin hypercube includes the steps of selecting a certain order vibration mode of a blade to conduct damper optimization design, obtaining analysis object parameters and node information, establishing a dynamic equation set with dry friction nonlinearity, solving a blade amplitude-frequency characteristic curve, optimizing analysis of multiple parameters of a damping structure and optimizing an original damping structure, and specifically comprises the following steps: s1, selecting a certain order of blade mode shape based on the blade mode shape to perform damper structure optimization design, and determining a blade resonance response frequency range; s2, processing the sub-structure of the blade group based on a dual mode synthesis method, and reducing a matrix combination overall structure matrix through the sub-structure; S3, constructing a dynamic equation set with dry friction nonlinearity based on the parameter matrix obtained in the S2 and the two-dimensional contact surface friction model; S4, solving an amplitude-frequency characteristic curve through a harmonic balance method and an arc length extension method based on the dynamic equation set with the dry friction nonlinearity obtained in the step S3; S5, based on Monte Carlo simulation and Latin hypercube design, writing a multi-parameter Monte Carlo simulation analysis program, and completing parameter optimization analysis aiming at the optimal damping effect of the damping device; S6, further optimizing the original damping structure based on the optimal parameter set obtained in the S4. Further, in the step S1, a certain-order vibration mode of the blade is selected for damper optimization design, and blade mode analysis is performed by establishing a single-blade finite element model. Further, in the step of obtaining the parameter matrix and the node information of the analysis object, the parameter matrix comprises the mass and the rigidity matrix of the blade and the damper, the finite element models of the two blades and the damper are established through finite element analysis software, the parameter matrix and the node information of the substructure are derived and recombined, and the finite element analysis software adopts ANSYS APDL. Further, in the step S3, a dynamics equation set with dry friction nonlinearity is established, wherein the dynamics equation set is a differential equation set, and the differential equation set is mapped to a nonlinear algebraic equation set under a frequency domain through fourier galy mapping. Further, in the step S4 of solving the amplitude-frequency characteristic curve of the blade, a solution under an initial frequency is solved through a harmonic balance method, wherein the nonlinear force and jacobian matrix solving method is to establish a functional relation between a normal force and a tangential force on a contact surface and relative displacement through a single-point two-dimensional friction model, then discretizing the contact surface to divide a plurality of contact units, in actual calculation, the contact force and the contact stiffness of each contact unit are overlapped to obtain an overall contact force and contact stiffness, and the solution under the initial frequency is solved through a Newton iteration method. And finally, solving an amplitude-frequ