Ship body vibration suppression parameter determination method and device, computer equipment, readable storage medium and program product

Abstract

The application relates to a ship body vibration suppression parameter determination method and device, computer equipment, a computer readable storage medium and a computer program product. The method comprises the following steps: acquiring a coupling finite element model of an equipment, a base and a ship body structure; performing power flow analysis on the coupling finite element model to obtain a power flow analysis result; identifying a key transmission path based on the power flow analysis result; and determining path stiffness matching parameters of path blocking and integrated frequency parameters of the equipment, the base and the ship body based on the key transmission path. The method can improve the vibration suppression effect and facilitate modification of the parameters in the later stage.

Description

Technical Field

[0001] This application relates to the field of marine engineering vibration prevention and suppression technology, and in particular to a method, device, computer equipment, computer-readable storage medium and computer program product for determining hull vibration suppression parameters. Background Technology

[0002] Vibration control of marine engineering equipment is a key technology for ensuring platform comfort, equipment reliability, and structural safety. From a physical perspective, the core of vibration control lies in the effective management of the vibration energy transmission process. Vibration energy is generated from excitation sources (such as rotating or reciprocating machinery like diesel engines, compressors, and pump sets), passes through equipment support structures, elastic bases, and local hull structures, and is ultimately transmitted to the overall hull structure and radiated outwards, resulting in vibration and noise problems.

[0003] Traditional technologies mostly focus on the independent treatment of the vibration source or the receiving end, such as using equipment vibration isolators to reduce the vibration source output, or adding damping materials to the hull structure to suppress the response.

[0004] However, these methods fail to be proactively designed based on the physical nature of energy transfer, resulting in energy being transferred through unknown paths and poor vibration suppression. Summary of the Invention

[0005] Therefore, it is necessary to provide a method, apparatus, computer equipment, computer-readable storage medium, and computer program product for determining ship vibration suppression parameters that can improve the vibration suppression effect, in order to address the above-mentioned technical problems.

[0006] In a first aspect, this application provides a method for determining ship vibration suppression parameters, the method comprising:

[0007] Obtain a coupled finite element model of the equipment, base, and hull structure;

[0008] Power flow analysis was performed on the coupled finite element model to obtain the power flow analysis results;

[0009] Based on the power flow analysis results, key transmission paths are identified;

[0010] Based on the key transmission path, the path stiffness matching parameters for path blocking and the integrated frequency parameters of the equipment, base, and hull are determined.

[0011] In one embodiment, the power flow analysis results include the energy transfer rate of each candidate path; the step of performing power flow analysis based on the coupled finite element model to obtain power flow analysis results includes:

[0012] Obtain the discretized structural elements in the coupled finite element model, and calculate the power flow density vector field of each structural element;

[0013] Divergence calculations are performed on each of the power flux density vector fields to determine key nodes, and candidate paths are determined based on the key nodes.

[0014] Based on the power flux density vector fields of each candidate path, the energy flux is calculated along the candidate paths to obtain the energy transfer rate of each candidate path.

[0015] In one embodiment, identifying critical transmission paths based on power flow analysis results includes:

[0016] The energy transfer flux of each candidate path is determined based on the energy transfer rate of each candidate path.

[0017] The candidate paths whose energy transfer flux is greater than the flux threshold are identified as the critical transfer paths.

[0018] In one embodiment, determining the path stiffness matching parameters for path blocking based on the critical transmission path includes:

[0019] For each of the critical transmission paths, the base in the critical transmission path is divided into series impedance units;

[0020] Based on the impedance unit, determine the input and output terminals of the base;

[0021] A one-dimensional impedance transmission matrix model is established based on the vibration transmission between the input terminal and the output terminal.

[0022] The parameters of the base in the one-dimensional impedance transfer matrix model are adjusted until the objective function meets the preset conditions. The adjusted base parameters are then used as path stiffness matching parameters. The objective function is determined based on the input impedance and output impedance of the base at the device excitation frequency.

[0023] In one embodiment, adjusting the parameters of the base in the one-dimensional impedance transfer matrix model until the objective function satisfies a preset condition, and using the adjusted base parameters as path stiffness matching parameters, includes:

[0024] If the objective function does not meet the preset conditions, the structural parameters of the base in the one-dimensional impedance transfer matrix model are changed to obtain the updated structural parameters of the base; wherein, the updated structural parameters are used to increase the input impedance of the base and / or decrease the output impedance of the base.

[0025] Based on the updated structure parameters, determine the new objective function;

[0026] If the new objective function satisfies the preset conditions, the updated structural parameters are used as the path stiffness matching parameters.

[0027] In one embodiment, based on the critical transmission path, the integrated frequency planning parameters for the equipment, base, and hull are determined, including:

[0028] Modal frequencies are extracted from the coupled finite element model based on the key transmission path, wherein the modal frequencies include the device disturbance frequency, the first-order frequency of the base, and at least one first-order frequency of the base;

[0029] Obtain the ship's hull modal frequencies; wherein, the ship's hull modal frequencies include the ship's local modal frequencies and the ship's overall modal frequencies; the ship's local modal frequencies are extracted from the coupled finite element model; the ship's overall modal frequencies are determined by simplifying the ship as a continuous beam;

[0030] A first avoidance ratio is determined based on the device disturbance frequency and the first-order frequency of the base; wherein the first avoidance ratio is greater than or equal to a first threshold, which is used to prevent the device excitation from triggering resonance in the base;

[0031] A second avoidance ratio is determined based on at least one first-order frequency of the base and the modal frequency of the hull, wherein the second avoidance ratio is greater than or equal to a second threshold, used to prevent the base vibration from exciting overall hull resonance; the first threshold is less than the second threshold.

[0032] Secondly, this application also provides a device for determining ship vibration suppression parameters, the device comprising:

[0033] The acquisition module is used to acquire coupled finite element models of the equipment, base, and hull structure.

[0034] The analysis module is used to perform power flow analysis on the coupled finite element model and obtain the power flow analysis results;

[0035] The identification module is used to identify critical transmission paths based on power flow analysis results;

[0036] The determination module is used to determine the path stiffness matching parameters of the path blockage and the integrated frequency parameters of the equipment, base, and hull based on the key transmission path.

[0037] Thirdly, this application also provides a computer device, including a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to implement the steps of the above-described method.

[0038] Fourthly, this application also provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the above-described method.

[0039] Fifthly, this application also provides a computer program product, including a computer program that, when executed by a processor, implements the steps of the above-described method.

[0040] The aforementioned method, apparatus, computer equipment, computer-readable storage medium, and computer program product for determining ship vibration suppression parameters first obtain a coupled finite element model of the equipment, base, and ship structure; then, power flow analysis is performed on the coupled finite element model to obtain the power flow analysis results; secondly, based on the power flow analysis results, critical transmission paths are identified; finally, based on the critical transmission paths, path stiffness matching parameters for path blocking and integrated frequency parameters of the equipment, base, and ship are determined. On the one hand, through path stiffness matching parameters and integrated frequency parameters along the vibration energy transmission path, energy is reflected back to the equipment (vibration source) instead of heading towards the ship (receiving end); on the other hand, integrated frequency parameters prevent resonance between the equipment and the base / ship, preventing vibration energy from being amplified by the ship structure itself, thereby suppressing vibration. Attached Figure Description

[0041] To more clearly illustrate the technical solutions in the embodiments of this application or related technologies, the drawings used in the description of the embodiments of this application or related technologies will be briefly introduced below. Obviously, the drawings described below are only some embodiments of this application. For those skilled in the art, other related drawings can be obtained based on these drawings without creative effort.

[0042] Figure 1 This is a flowchart illustrating a method for determining hull vibration suppression parameters in one embodiment.

[0043] Figure 2 This is a schematic diagram of the power flow analysis based on a coupled finite element model in one embodiment;

[0044] Figure 3 This is a flowchart illustrating the process of determining path stiffness matching parameters for path blocking based on a critical transmission path in one embodiment.

[0045] Figure 4 This is a flowchart illustrating the process of adjusting the parameters of the base in a one-dimensional impedance transfer matrix model in one embodiment.

[0046] Figure 5 This is a flowchart illustrating the process of determining integrated frequency planning parameters for equipment, base, and hull based on a critical transmission path in one embodiment.

[0047] Figure 6 This is a structural block diagram of a device for determining ship vibration suppression parameters in one embodiment;

[0048] Figure 7 This is an internal structural diagram of a computer device in one embodiment. Detailed Implementation

[0049] To make the objectives, technical solutions, and advantages of this application clearer, the following detailed description is provided in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the scope of this application.

[0050] In one embodiment, such as Figure 1 As shown, a method for determining ship vibration suppression parameters is provided. This embodiment illustrates the application of this method to a terminal. It is understood that this method can also be applied to a server, and to a system including both a terminal and a server, and implemented through interaction between the terminal and the server. In this embodiment, the method includes the following steps S102 to S108. Wherein:

[0051] Step S102: Obtain the coupled finite element model of the equipment, base, and hull structure.

[0052] Among them, equipment refers to equipment that generates vibration, that is, rotating or reciprocating mechanical equipment on offshore platforms, such as diesel engines, compressors, pump sets, etc., which are the source of vibration energy. The excitation force (main frequency and harmonics) generated during their operation will be transmitted outward.

[0053] The base is an intermediate elastic support structure (such as an elastic base or vibration isolator) located between the equipment and the hull. It serves as both the mounting foundation for the vibration source and a critical link in the transmission path of vibration energy from the equipment to the hull.

[0054] The hull structure includes the local hull structure (the part connected to the base) and the overall hull structure. It is the final receiver of vibration. After the vibration energy is transmitted to the hull, it will cause structural vibration and radiate noise to the outside.

[0055] In practical applications, vibration energy is transmitted along the sequence of equipment-base-local hull structure-overall hull structure.

[0056] Optionally, a coupled finite element model is established at the terminal, consisting of the equipment (crude oil pump), the elastic base, and the local hull structure. This coupled finite element model accurately characterizes the equipment's mass characteristics, the base stiffness distribution, the dynamic characteristics of the local hull structure, and the contact stiffness characteristics of each connection interface. For example, the coupled finite element model includes: Pump mass characteristics: total weight 8.5t, center of mass height 0.95m; Base stiffness distribution: I-beam dimensions H200×200×8×12, elastic modulus 206GPa; Local hull structure: deck plate thickness 14mm, longitudinal girder spacing 600mm; Connection interfaces: stiffness of 4 vibration isolators (dynamic stiffness 2.2e6N / m), and bolt contact stiffness between the base and the deck (modeled as contact elements).

[0057] Step S104: Perform power flow analysis on the coupled finite element model to obtain the power flow analysis results.

[0058] Optionally, the terminal performs power flow analysis based on a coupled finite element model to calculate the transmission path distribution and energy flux of vibration energy in the structural unit, thereby obtaining the power flow analysis results.

[0059] Step S106: Based on the power flow analysis results, identify the critical transmission paths.

[0060] The asymmetry of the structure, the distribution of materials, and the differences in the location of connection points can lead to uneven energy distribution, resulting in more than one critical transmission path.

[0061] Optionally, based on the power flow analysis results such as the distribution of transmission paths and the magnitude of energy flux, the terminal can identify at least one critical transmission path and further locate the critical vibration transmission components in each critical transmission path, such as the main support component of the base, the connection interface of the equipment base, the connection interface of the base and the hull, and the high stress concentration components on the critical transmission path.

[0062] Furthermore, the terminal establishes a digital archive of the critical transmission path, including the topology of the critical transmission path, the design parameters of each component, the target impedance characteristics, and the energy flux benchmark value.

[0063] Step S108: Based on the critical transmission path, determine the path stiffness matching parameters for path blocking and the integrated frequency parameters of the equipment, base, and hull.

[0064] The path stiffness matching parameter is used to actively block the transmission of vibration energy, so that the energy is reflected back to the equipment side inside the base and cannot be transmitted to the hull. The integrated frequency parameter is used to ensure that the excitation frequency of the equipment is offset from the frequency of the base / hull, avoiding resonance and preventing the vibration energy from being amplified by the structure itself.

[0065] Optionally, the terminal determines the path stiffness matching parameters for path blocking based on parameters on the critical transmission path. This includes optimizing the cross-sectional dimensions, material distribution, and connection method of the base support component using a path impedance transmission model, and using the optimized cross-sectional dimensions, material distribution, and connection method of the base support component as the path stiffness matching parameters. This optimization can be achieved by introducing impedance abrupt change structures (such as variable cross-sections, material transition zones, flexible connections, etc.) into the critical vibration transmission component.

[0066] Optionally, the terminal determines the integrated frequency parameters of the equipment, base, and hull based on the parameters on the critical transmission path, including: the terminal jointly determines the integrated frequency parameters to prevent resonance of the equipment, base, and hull based on the parameters on the critical transmission path and the established coupled finite element model of the equipment, base, and hull structure.

[0067] Furthermore, based on path stiffness matching parameters and integrated frequency parameters, the terminal performs full-process quality control on the physical process of the critical transmission path, including traceability of path stiffness during the manufacturing process: after key processes in the manufacturing of critical components such as the base (e.g., welding completion, heat treatment completion, machining completion), on-site spot checks are conducted using the following methods:

[0068] Vibration testing hammer method: using a force hammer for excitation, measuring the acceleration response, calculating the frequency response function, and extracting the local dynamic stiffness of the component.

[0069] Static stiffness testing method: Apply a known static load, measure the displacement response, and calculate the static stiffness value.

[0070] Compare the measured stiffness values ​​with the design target values ​​to ensure that the deviations are within the allowable range (generally, the dynamic stiffness deviation should not exceed ±15%, and the static stiffness deviation should not exceed ±10%). Mark, rework, or scrap any components that exceed the tolerances, and trace the cause of the manufacturing deviations.

[0071] Establish a path stiffness traceability archive for the manufacturing stage, and record the measured data, inspection time, inspection personnel, and deviation handling records for each sampled component.

[0072] This also includes ensuring the integrity of the installation process path: A "Process Specification for Installation Interfaces" is formulated and implemented. This specification includes: contact surface treatment requirements: specifying the surface roughness and flatness requirements for the contact surfaces between the equipment base and the pedestal, and between the pedestal and the hull; fit rate inspection standards: using the coloring method or lead pressing method to inspect the contact surface fit rate, requiring the fit rate of the main load-bearing areas to be no less than 75%; connecting bolt preload control: using an intelligent torque wrench or hydraulic tensioner, applying preload in stages according to the specified preload sequence, and recording the actual preload value; installation quality inspection: after installation, using the hammer impact method or vibration exciter method to test the overall dynamic stiffness after installation to verify whether it meets the design requirements. This ensures that the interface contact stiffness after installation meets design expectations and does not introduce new vibration "short circuit" paths.

[0073] Furthermore, vibration sensors (accelerometers or velocity sensors) and / or strain gauges are deployed on key vibration transmission components to construct a path condition monitoring network. This monitoring includes: long-term continuous monitoring of the vibration acceleration, velocity, or dynamic stress response of key components; capturing changes in the path's dynamic characteristics during transient processes such as equipment start-up and shutdown, and operating condition switching; and periodically conducting transmission rate tests to evaluate the evolution of the overall path transmission characteristics.

[0074] The terminal establishes an evolution model of path transmission rate over time, operating conditions, and environmental conditions. This evolution model is based on the following factors: interface stiffness degradation caused by bolt preload relaxation; component stiffness changes caused by material fatigue; changes in cross-sectional properties caused by corrosion; and fretting wear of the contact surface caused by long-term vibration.

[0075] When monitoring data indicates that the transmission rate of a certain path has abnormally increased beyond the preset threshold, the terminal triggers a tiered warning: based on the warning and path characteristics, it automatically generates targeted maintenance instructions, such as the specific maintenance location.

[0076] In addition, a full lifecycle data management platform is established at the terminal to collect and structure-store path-related data from all stages of design, manufacturing, installation, and operation and maintenance. In-depth mining of operation and maintenance data is conducted: common degradation patterns on transmission paths are identified for different offshore platforms and equipment types; the correlation between degradation patterns and initial design parameters, manufacturing processes, and installation quality is analyzed; the impact weight of various factors on the long-term performance of the path is quantified; and the analysis results are fed back to the determination of path stiffness matching parameters and integrated frequency parameters of equipment, base, and hull, such as updating candidate paths, optimizing critical transmission path identification rules, adjusting frequency avoidance criteria, and adjusting path stiffness matching parameters and integrated frequency parameters of equipment, base, and hull. This forms a closed-loop iterative mechanism of "design, manufacturing, installation, operation and maintenance feedback optimization," enabling continuous evolution of the methodology.

[0077] The above-mentioned method for determining ship vibration suppression parameters first obtains a coupled finite element model of the equipment, base, and hull structure; then, it performs power flow analysis on the coupled finite element model to obtain the power flow analysis results; secondly, based on the power flow analysis results, it identifies the critical transmission paths; finally, based on the critical transmission paths, it determines the path stiffness matching parameters for path blocking and the integrated frequency parameters of the equipment, base, and hull. On the one hand, through the path stiffness matching parameters and integrated frequency parameters on the transmission path of vibration energy, the energy is reflected back to the equipment (vibration source) instead of heading towards the hull (receiving end); on the other hand, the integrated frequency parameters prevent resonance between the equipment and the base / hull, preventing the vibration energy from being amplified by the hull structure itself, thereby achieving the effect of suppressing vibration.

[0078] In one exemplary embodiment, such as Figure 2 As shown, the power flow analysis results include the energy transfer rate of each candidate path; power flow analysis is performed based on the coupled finite element model to obtain the power flow analysis results, including the following steps S202 to S206. Wherein:

[0079] Step S202: Obtain the discretized structural elements in the coupled finite element model and calculate the power flow density vector field of each structural element.

[0080] Among them, the power flow density vector field is a spatial distribution field that describes the direction and intensity of vibration energy flow in the structural unit, and it is the basis for subsequent analysis.

[0081] Optionally, the terminal acquires the discretized structural elements in the coupled finite element model, such as shell elements, beam elements, spring elements, bush elements, etc., or preset structural partitions, such as the base web region, panel region, elbow plate region, elastic base region, hull stiffener region, etc. Based on the stress tensor and velocity vector on the boundary of each structural element or structural partition, the terminal accurately locates the "flow direction" and "convergence / divergence region" of energy in space, thereby obtaining the power flow density vector field of each structural element.

[0082] Step S204: Perform divergence calculation on each power flux density vector field to determine key nodes, and determine candidate paths based on the key nodes.

[0083] Optionally, the terminal performs divergence calculations on each power flow density vector field to identify energy convergence nodes or diversion nodes as key nodes. Here, energy convergence nodes are those where "inflow > outflow," and diversion nodes are those where "outflow > inflow."

[0084] Furthermore, the terminal determines candidate paths based on key nodes. For example, if a node A is found to be an energy convergence node (input much greater than output) and is located at the connection between the equipment base and the hull, then node A can be considered a necessary point on the path. The candidate path might be, for example, equipment - energy convergence node - intermediate section of the path - energy distribution node - final receiving point (i.e., the hull).

[0085] It should be noted that there can be multiple candidate paths.

[0086] Step S206: Based on each power flux density vector field, calculate the energy flux along the candidate path to obtain the energy transfer rate of each candidate path.

[0087] Energy transfer rate, also known as energy transfer efficiency, is an inherent property of every path and is dimensionless or a percentage.

[0088] Optionally, the terminal performs line integration or summation of the energy flux along each candidate path (e.g., from the device to the top of the base, the web of the base, the bottom plate of the base, and the hull) based on the power flux density vector field of each candidate path, to obtain the energy transfer rate of each candidate path.

[0089] In this embodiment, the energy transfer rate of each candidate path is calculated by the power flow density vector field of each candidate path, which can serve as the basis for subsequent identification of critical transmission paths.

[0090] Following the previous embodiment, based on the power flow analysis results, critical transmission paths are identified, including: determining the energy transfer flux of each candidate path according to the energy transfer rate of each candidate path; and identifying candidate paths with energy transfer flux greater than the flux threshold as critical transmission paths.

[0091] Among them, energy transfer flux is the absolute value of energy transfer rate multiplied by total input energy (with dimensions, such as watts kW), representing the amount of vibration energy actually carried by the path.

[0092] Optionally, the terminal obtains the energy transfer rate of each candidate path, multiplies the energy transfer rate of each candidate path by the actual input energy spectrum of the device at the main disturbance frequency, and calculates the expected energy transfer flux of each path.

[0093] Furthermore, the terminal identifies candidate paths with energy transfer flux exceeding a flux threshold (e.g., 10% of total input energy) as critical transfer paths. Then, it traces back along these critical transfer paths to locate the component that bears the maximum energy flux density or generates the maximum impedance change, which is the critical vibration transmission component.

[0094] It should be noted that the reason for choosing energy transfer flux instead of energy transfer rate to determine the critical transfer path is that: a high-efficiency path with very low input energy has little actual harm; while a medium-efficiency path with a large input energy is the key to vibration suppression.

[0095] In this embodiment, the critical transmission path is determined by the energy transfer flux, which can accurately identify the critical transmission path and provide a basis for parameters for subsequent path blocking.

[0096] In one exemplary embodiment, such as Figure 3 As shown, based on the critical transmission path, the path stiffness matching parameters for path blocking are determined, including steps S302 to S308. Wherein:

[0097] Step S302: For each critical transmission path, the base in the critical transmission path is split into series impedance units.

[0098] Optionally, the terminal obtains the path segments inside the base and its supporting components from the identified key transmission paths, and divides the base into several series impedance units. Each impedance unit corresponds to a physical component or connection interface, such as the base top connecting plate (connecting to equipment), the base web plate, the base reinforcing elbow plate, the base bottom plate (connecting to the hull), and the bolt connection interface between the base and the hull.

[0099] Step S304: Determine the input and output terminals of the base based on the impedance unit.

[0100] Optionally, the terminal determines the input and output terminals of the base based on the physical components corresponding to the impedance units. For example, the impedance units corresponding to the top connecting plate of the base and part of the web plate are used as the input terminals of the base. The impedance units corresponding to the bottom plate of the base and the impedance units corresponding to the bolt connection interface between the base and the hull are used as the output terminals of the base.

[0101] Step S306: Based on the vibration transmission between the input and output terminals, establish a one-dimensional impedance transmission matrix model.

[0102] Optionally, when the vibration is transmitted from the input terminal A (connected to the equipment) to the output terminal B (connected to the hull), the terminal establishes a one-dimensional impedance transmission matrix model, as shown in formula (1). Under the action of simple harmonic excitation FeiωtFeiωt, the sum of the system's inertial force, damping force, and elastic force is equal to the external force.

[0103] Formula (1)

[0104] In the formula, m ij The elements of the mass matrix (i,j=1,2) represent the system's inertial distribution, in kg; c ijThese are damping matrix elements, representing energy dissipation characteristics, in N·s / m; stiffness matrix elements, k ij The elements represent the stiffness matrix, with units of N / m; X1 and X2 are the displacement responses of two degrees of freedom (input and output), with units of m; and It represents the velocity response (first derivative), and its unit is m / s; and Acceleration response (second derivative), unit is m / s²; F1 and F2 represent the magnitudes (complex numbers, including magnitude and phase) of the external forces acting on two degrees of freedom, unit is N; i represents the imaginary unit, t represents time. Represents the excitation angular frequency, with units of rad / s; The time-domain factor represents the harmonic excitation.

[0105] Assume that the above equation has a set of particular solutions in the time domain, and transform the time domain problem into a steady-state problem in the frequency domain, so as to solve the harmonic response as shown in formula (2).

[0106] Formula (2)

[0107] In the formula, X1 and X2 represent the complex amplitude of the displacement response (including amplitude and phase).

[0108] Substituting the particular solution into the equation, we obtain the following expression, and eliminate the non-zero values. Formula (3) is obtained, which is the frequency domain impedance matrix of the system, used to calculate the impedance value at different frequencies.

[0109] Formula (3)

[0110] In the formula, Represents the inertial term (mass × squared angular frequency); This represents the damping term (imaginary number).

[0111] The coefficients of the particular solution can be obtained from formula (3), and the steady-state response of the system is formula (4). The force and velocity response of the component at each frequency is actually solved, and then the input impedance is calculated.

[0112] Formula (4)

[0113] In the formula, Z ij It is called the impedance matrix and represents the inherent characteristics of the system.

[0114] Step S308: Adjust the parameters of the base in the one-dimensional impedance transfer matrix model until the objective function meets the preset conditions, and use the adjusted base parameters as path stiffness matching parameters.

[0115] The objective function is based on the input impedance Z of the base at the device excitation frequency (e.g., the pump's rotational frequency of 24.83 Hz and its main harmonic frequencies). in and the output impedance Z of the base out Determined. Objective function = Z in / Z out The preset condition is to maximize the ratio of input impedance to output impedance, i.e., max(Z). in / Z out ).

[0116] Optionally, the terminal adjusts the parameters of the base in the one-dimensional impedance transfer matrix model, including but not limited to the cross-sectional dimensions at the connection between the base and the equipment, the added mass, the material damping, the cross-sectional dimensions at the connection between the base and the hull, whether to introduce flexible connections or variable cross sections, the geometric parameters and material distribution of the intermediate impedance units, until the objective function satisfies the preset condition that the ratio of input impedance to output impedance is maximized, and the adjusted base parameters are used as path stiffness matching parameters.

[0117] In this embodiment, by maximizing the impedance mismatch of the path machinery (base), vibration energy can be reflected back to the equipment instead of heading towards the hull, thereby better suppressing vibration. Furthermore, in actual operation, adjusting the base parameters is more convenient than adjusting the hull and equipment parameters.

[0118] In one exemplary embodiment, such as Figure 4 As shown, the parameters of the base in the one-dimensional impedance transfer matrix model are adjusted until the objective function meets the preset conditions. The adjusted base parameters are then used as path stiffness matching parameters, including steps S402 to S406. Wherein:

[0119] Step S402: If the objective function does not meet the preset conditions, change the structural parameters of the base in the one-dimensional impedance transfer matrix model to obtain the updated structural parameters of the base.

[0120] Among them, updating the structural parameters is used to increase the input impedance of the base and / or decrease the output impedance of the base.

[0121] Optionally, if the objective function does not meet the preset conditions, the structural parameters of the base in the one-dimensional impedance transfer matrix model can be changed. Specifically, the structural parameters of the base can be adjusted with the goal of increasing the input impedance of the base and / or decreasing the output impedance of the base. For example, at the connection between the equipment and the base: increase the cross-section, increase the constraint mass, and use high-damping materials to improve local stiffness and energy dissipation capacity; at the connection between the base and the hull: reduce the cross-section, introduce flexible gaskets, design a variable cross-section transition zone, or use an elastic base to reduce local stiffness; in the middle section of the path: optimize the thickness and shape of the web and elbow plates, and introduce impedance abrupt change structures (such as material transition zones) to obtain updated structural parameters of the base.

[0122] Step S404: Determine the new objective function based on the updated structural parameters.

[0123] Step S406: If the new objective function meets the preset conditions, the updated structural parameters will be used as path stiffness matching parameters.

[0124] Optionally, the terminal determines a new objective function based on the updated structural parameters. Z is calculated for different combinations of design variables. in / Z out Until the preset target (e.g., maximization) is reached, the structural parameters will be updated as path stiffness matching parameters.

[0125] It should be noted that this process involves repeated calculations, and the objective function satisfies the preset condition, namely Z. in / Z out Maximization is the termination condition. Otherwise, new and updated structure parameters would be continuously determined.

[0126] In this embodiment, by maximizing impedance mismatch in a physical sense, the efficient transmission of vibration energy from the equipment to the hull is effectively blocked.

[0127] In one exemplary embodiment, such as Figure 5 As shown, based on the critical transmission path, the integrated frequency planning parameters for the equipment, base, and hull are determined, including steps S502 to S506. Wherein:

[0128] Step S502: Extract modal frequencies from the coupled finite element model based on the critical transmission path.

[0129] The modal frequencies include the equipment disturbance frequency, the first-order frequency of the base, and at least one first-order frequency of the base.

[0130] Optionally, the terminal extracts the device disturbance frequency (including the dominant frequency f) from the coupled finite element model based on the critical transmission path. excitation and the main harmonic frequencies), the first-order frequency f of the base. mount,1 The base has at least one first-order frequency f. mount,i (i is not equal to 1).

[0131] Step S504: Obtain the hull modal frequencies.

[0132] Among them, the hull modal frequency f hull,j This includes the local modal frequencies and the overall modal frequencies of the hull; the local modal frequencies are extracted from the coupled finite element model; the overall modal frequencies are determined by simplifying the ship as a continuous beam.

[0133] Optionally, the terminal simplifies the hull into a continuous beam, and its natural frequency is calculated as shown in formula (5).

[0134] Formula (5)

[0135] In the formula, EI represents the stiffness of the beam. Let be the mass of the beam, where A represents the density of the beam material; A represents the cross-sectional area of ​​the beam. Represents the nth natural angular frequency; The wavenumber constant (eigenvalue divided by length) represents the wavenumber constant related to the boundary conditions and modal order; L represents the length of the beam.

[0136] Step S506: Determine the first avoidance ratio based on the equipment disturbance frequency and the first-order frequency of the base.

[0137] Among them, the first avoidance ratio is greater than or equal to the first threshold, which is used to prevent the equipment from stimulating the base resonance.

[0138] Optionally, the terminal extracts the device's main frequency f from the device's disturbance frequency. excitation And calculate the device's main frequency f excitation With the first-order frequency f of the base mount,1 The ratio is used as the first avoidance ratio (requiring a value greater than or equal to 15%), and the first threshold ratio is, for example, 15%.

[0139] Step S508: Determine the second avoidance ratio based on at least one first-order frequency of the base and the modal frequency of the hull.

[0140] The second avoidance ratio is greater than or equal to the second threshold, used to prevent the base vibration from triggering overall hull resonance; the first threshold is less than the second threshold. The second threshold is, for example, 20%.

[0141] Optionally, the terminal computing base has at least one first-order frequency f. mount,i (i ≠ 1) and hull modal frequencies f hull,j The ratio of 20% to 20% is used as the second avoidance ratio, and the second avoidance ratio is greater than or equal to 20%.

[0142] The avoidance ratio is used to prevent the device's excitation frequency from approaching the natural frequency of any component or support structure along the path.

[0143] In this embodiment, by using the clearance ratio, resonance between the equipment and the base / hull can be avoided, preventing vibration energy from being amplified by the hull structure itself, thereby suppressing vibration.

[0144] It should be understood that although the steps in the flowcharts of the embodiments described above are shown sequentially according to the arrows, these steps are not necessarily executed in the order indicated by the arrows. Unless explicitly stated herein, there is no strict order restriction on the execution of these steps, and they can be executed in other orders. Moreover, at least some steps in the flowcharts of the embodiments described above may include multiple steps or multiple stages. These steps or stages are not necessarily completed at the same time, but can be executed at different times. The execution order of these steps or stages is not necessarily sequential, but can be performed alternately or in turn with other steps or at least some of the steps or stages of other steps.

[0145] Based on the same inventive concept, this application also provides a device for determining hull vibration suppression parameters to implement the above-mentioned method for determining hull vibration suppression parameters. The solution provided by this device is similar to the solution described in the above method. Therefore, the specific limitations of one or more embodiments of the device for determining hull vibration suppression parameters provided below can be found in the limitations of the method for determining hull vibration suppression parameters described above, and will not be repeated here.

[0146] In one exemplary embodiment, such as Figure 6 As shown, a device for determining ship vibration suppression parameters is provided, comprising: an acquisition module 601, an analysis module 602, an identification module 603, and a determination module 604, wherein:

[0147] The acquisition module 601 is used to acquire the coupled finite element model of the equipment, base, and hull structure.

[0148] Analysis module 602 is used to perform power flow analysis on the coupled finite element model and obtain the power flow analysis results.

[0149] The identification module 603 is used to identify critical transmission paths based on power flow analysis results.

[0150] The determination module 604 is used to determine the path stiffness matching parameters of the path blockage and the integrated frequency parameters of the equipment, base and hull based on the critical transmission path.

[0151] In an exemplary embodiment, the analysis module 602 is further configured to obtain the discretized structural elements in the coupled finite element model, calculate the power flux density vector field of each structural element, perform divergence calculation on each power flux density vector field, determine key nodes, determine candidate paths based on key nodes, and calculate the energy flux along the candidate paths based on each power flux density vector field to obtain the energy transfer rate of each candidate path.

[0152] In an exemplary embodiment, the identification module 603 is further configured to determine the energy transfer flux of each candidate path based on the energy transfer rate of each candidate path; and to identify candidate paths whose energy transfer flux is greater than the flux threshold as critical transfer paths.

[0153] In an exemplary embodiment, the determining module 604 is further configured to, for each critical transmission path, divide the base in the critical transmission path into series impedance units; determine the input and output terminals of the base based on the impedance units; establish a one-dimensional impedance transmission matrix model based on the vibration transmission between the input and output terminals; adjust the parameters of the base in the one-dimensional impedance transmission matrix model until the objective function meets the preset conditions, and use the adjusted base parameters as path stiffness matching parameters; wherein, the objective function is determined based on the input impedance and output impedance of the base at the device excitation frequency.

[0154] In an exemplary embodiment, the determining module 604 is further configured to change the structural parameters of the base in the one-dimensional impedance transfer matrix model to obtain updated structural parameters of the base when the objective function does not meet the preset conditions; wherein, the updated structural parameters are used to increase the input impedance of the base and / or decrease the output impedance of the base; a new objective function is determined based on the updated structural parameters; and if the new objective function meets the preset conditions, the updated structural parameters are used as path stiffness matching parameters.

[0155] In an exemplary embodiment, the determining module 604 is further configured to extract modal frequencies from the coupled finite element model based on the critical transmission path, wherein the modal frequencies include the equipment disturbance frequency, the first-order frequency of the base, and at least one first-order frequency of the base; obtain the hull modal frequencies; wherein the hull modal frequencies include hull local modal frequencies and hull global modal frequencies; the hull local modal frequencies are extracted from the coupled finite element model; the hull global modal frequencies are determined by simplifying the ship as a continuous beam; determine a first avoidance ratio based on the equipment disturbance frequency and the first-order frequency of the base; wherein the first avoidance ratio is greater than or equal to a first threshold, used to prevent equipment excitation from triggering base resonance; determine a second avoidance ratio based on at least one first-order frequency of the base and hull modal frequencies, wherein the second avoidance ratio is greater than or equal to a second threshold, used to prevent base vibration from triggering hull global resonance; the first threshold is less than the second threshold.

[0156] The modules in the aforementioned hull vibration suppression parameter determination device can be implemented entirely or partially through software, hardware, or a combination thereof. These modules can be embedded in the processor of a computer device in hardware form or independent of it, or stored in the memory of a computer device in software form, so that the processor can call and execute the corresponding operations of each module.

[0157] In one exemplary embodiment, a computer device is provided, which may be a server, and its internal structure diagram may be as follows: Figure 7 As shown, the computer device includes a processor, memory, input / output (I / O) interfaces, and a communication interface. The processor, memory, and I / O interfaces are connected via a system bus, and the communication interface is also connected to the system bus via the I / O interfaces. The processor provides computational and control capabilities. The memory includes non-volatile storage media and internal memory. The non-volatile storage media stores the operating system, computer programs, and a database. The internal memory provides the environment for the operating system and computer programs stored in the non-volatile storage media. The database stores ship vibration damping parameter data. The I / O interfaces are used for information exchange between the processor and external devices. The communication interface is used for communication with external terminals via a network connection. When the computer program is executed by the processor, it implements a method for determining ship vibration damping parameters.

[0158] Those skilled in the art will understand that Figure 7 The structure shown is merely a block diagram of a portion of the structure related to the present application and does not constitute a limitation on the computer device to which the present application is applied. Specific computer devices may include more or fewer components than those shown in the figure, or combine certain components, or have different component arrangements.

[0159] In one embodiment, a computer device is also provided, including a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to implement the steps in the above method embodiments.

[0160] In one embodiment, a computer-readable storage medium is provided having a computer program stored thereon that, when executed by a processor, implements the steps in the above method embodiments.

[0161] In one embodiment, a computer program product is provided, including a computer program that, when executed by a processor, implements the steps in the above method embodiments.

[0162] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium, and when executed, it can include the processes of the embodiments of the above methods. Any references to memory, databases, or other media used in the embodiments provided in this application can include at least one of non-volatile memory and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetic random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM can take many forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM). The databases involved in the embodiments provided in this application may include at least one type of relational database and non-relational database. Non-relational databases may include, but are not limited to, blockchain-based distributed databases. The processors involved in the embodiments provided in this application may be general-purpose processors, central processing units, graphics processing units, digital signal processors, programmable logic devices, quantum computing-based data processing logic devices, artificial intelligence (AI) processors, etc., and are not limited to these.

[0163] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this application.

[0164] The embodiments described above are merely illustrative of several implementation methods of this application, and while the descriptions are specific and detailed, they should not be construed as limiting the scope of this patent application. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of this application, and these all fall within the protection scope of this application. Therefore, the protection scope of this application should be determined by the appended claims.

Claims

1. A method for determining ship vibration suppression parameters, characterized in that, The method includes: Obtain a coupled finite element model of the equipment, base, and hull structure; Power flow analysis was performed on the coupled finite element model to obtain the power flow analysis results; Based on the power flow analysis results, key transmission paths are identified; Based on the key transmission path, the path stiffness matching parameters for path blocking and the integrated frequency parameters of the equipment, base, and hull are determined.

2. The method according to claim 1, characterized in that, The power flow analysis results include the energy transfer rate of each candidate path; The power flow analysis based on the coupled finite element model, to obtain the power flow analysis results, includes: Obtain the discretized structural elements in the coupled finite element model, and calculate the power flow density vector field of each structural element; Divergence calculations are performed on each of the power flux density vector fields to determine key nodes, and candidate paths are determined based on the key nodes. Based on the power flux density vector fields of each candidate path, the energy flux is calculated along the candidate paths to obtain the energy transfer rate of each candidate path.

3. The method according to claim 2, characterized in that, The identification of critical transmission paths based on power flow analysis results includes: The energy transfer flux of each candidate path is determined based on the energy transfer rate of each candidate path. The candidate paths whose energy transfer flux is greater than the flux threshold are identified as the critical transfer paths.

4. The method according to claim 1, characterized in that, The determination of path stiffness matching parameters for path blocking based on the critical transmission path includes: For each of the critical transmission paths, the base in the critical transmission path is divided into series impedance units; Based on the impedance unit, determine the input and output terminals of the base; A one-dimensional impedance transmission matrix model is established based on the vibration transmission between the input terminal and the output terminal. The parameters of the base in the one-dimensional impedance transfer matrix model are adjusted until the objective function meets the preset conditions. The adjusted base parameters are then used as path stiffness matching parameters. The objective function is determined based on the input impedance and output impedance of the base at the device excitation frequency.

5. The method according to claim 4, characterized in that, The process of adjusting the parameters of the base in the one-dimensional impedance transfer matrix model until the objective function satisfies a preset condition, and using the adjusted base parameters as path stiffness matching parameters, includes: If the objective function does not meet the preset conditions, the structural parameters of the base in the one-dimensional impedance transfer matrix model are changed to obtain the updated structural parameters of the base; wherein, the updated structural parameters are used to increase the input impedance of the base and / or decrease the output impedance of the base. Based on the updated structure parameters, determine the new objective function; If the new objective function satisfies the preset conditions, the updated structural parameters are used as the path stiffness matching parameters.

6. The method according to claim 1, characterized in that, Based on the aforementioned critical transmission path, the integrated frequency planning parameters for the equipment, base, and hull are determined, including: Modal frequencies are extracted from the coupled finite element model based on the key transmission path, wherein the modal frequencies include the device disturbance frequency, the first-order frequency of the base, and at least one first-order frequency of the base; Obtain the ship's hull modal frequencies; wherein, the ship's hull modal frequencies include the ship's local modal frequencies and the ship's overall modal frequencies; the ship's local modal frequencies are extracted from the coupled finite element model; the ship's overall modal frequencies are determined by simplifying the ship as a continuous beam; A first avoidance ratio is determined based on the device disturbance frequency and the first-order frequency of the base; wherein the first avoidance ratio is greater than or equal to a first threshold, which is used to prevent the device excitation from triggering resonance in the base; A second avoidance ratio is determined based on at least one first-order frequency of the base and the modal frequency of the hull, wherein the second avoidance ratio is greater than or equal to a second threshold, used to prevent the base vibration from exciting overall hull resonance; the first threshold is less than the second threshold.

7. A device for determining ship vibration suppression parameters, characterized in that, The device includes: The acquisition module is used to acquire coupled finite element models of the equipment, base, and hull structure. The analysis module is used to perform power flow analysis on the coupled finite element model and obtain the power flow analysis results; The identification module is used to identify critical transmission paths based on power flow analysis results; The determination module is used to determine the path stiffness matching parameters of the path blockage and the integrated frequency parameters of the equipment, base, and hull based on the key transmission path.

8. A computer device comprising a memory and a processor, wherein the memory stores a computer program, characterized in that, When the processor executes the computer program, it implements the steps of the method according to any one of claims 1 to 6.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by a processor, it implements the steps of the method according to any one of claims 1 to 6.

10. A computer program product, comprising a computer program, characterized in that, When the computer program is executed by a processor, it implements the steps of the method according to any one of claims 1 to 6.

Images