Modal analysis method and device of Coriolis mass flowmeter and storage medium

Abstract

The invention discloses a modal analysis method and device of a Coriolis mass flowmeter and a storage medium, relates to the technical field of flow measurement, and aims to solve the problem that the accuracy of a modal analysis result of the Coriolis mass flowmeter is low. The modal analysis method of the Coriolis mass flowmeter comprises the steps that the real-time relation between fluid pressure and structural rigidity is determined through a bidirectional fluid-solid coupling model; dynamically correcting the flow data of the Coriolis mass flow meter based on the real-time relationship; and determining a modal analysis result of the Coriolis mass flowmeter based on the dynamically corrected flow data and a nonlinear relationship between the straight pipe length and the phase difference.

Description

Technical Field

[0001] The present application relates to the field of flow measurement technology, and in particular to a modal analysis method, device, and storage medium for a Coriolis mass flowmeter. Background Art

[0002] As an instrument that directly measures fluid mass flow, Coriolis mass flowmeter is widely used in chemical, petroleum, pharmaceutical and other fields.

[0003] Currently, the modal analysis of Coriolis mass flowmeters is usually performed based on an empty pipe structure analysis method (such as an empty pipe second-order frequency of 101.73 Hz).

[0004] However, with the influence of the added mass effect of the fluid, the prediction error of the actual natural frequency of the flow-carrying tube may increase, thereby reducing the accuracy of the modal analysis results of the Coriolis mass flowmeter. Summary of the Invention

[0005] The purpose of this application is to provide a modal analysis method, device and storage medium for a Coriolis mass flowmeter, aiming to solve the problem of low accuracy of modal analysis results of Coriolis mass flowmeters.

[0006] To achieve the above objectives, this application adopts the following technical solutions:

[0007] In the first aspect, the present application provides a modal analysis method for a Coriolis mass flowmeter, comprising: determining the real-time relationship between fluid pressure and structural stiffness through a bidirectional fluid-solid coupling model; dynamically correcting the flow data of the Coriolis mass flowmeter based on the real-time relationship; and determining the modal analysis results of the Coriolis mass flowmeter based on the dynamically corrected flow data and the nonlinear relationship between the straight pipe length and the phase difference.

[0008] The modal analysis method for a Coriolis mass flowmeter provided in an embodiment of the present application can determine the fluid-solid coupling condition by correlating the fluid pressure with the structural stiffness in real time, and compensate and correct the flow data of the Coriolis mass flowmeter based on the dynamic fluid-solid coupling condition. Then, based on the corrected flow data, the modal analysis results of the Coriolis mass flowmeter are determined in combination with the nonlinear relationship between the straight pipe length and the phase difference. In this way, by referring to the fluid-solid coupling condition and quantifying the nonlinear effects of the structural parameters, the accuracy of the modal analysis results of the Coriolis mass flowmeter can be provided.

[0009] In some embodiments, the parameter settings of the bidirectional fluid-structure coupling model satisfy at least one of the following:

[0010] The pipe diameter and wall thickness of the measuring pipe are 15 mm and 1.5 mm, respectively.

[0011] The structural parameters of the vibration damping plate are set according to the actual size;

[0012] The material of the measuring tube is 316 stainless steel;

[0013] The elastic modulus of the measuring tube is Gigapascal;

[0014] Poisson's ratio is 0.275;

[0015] Density is 7980 kg / m3;

[0016] The fluid medium is methane at 300 Kelvin;

[0017] The fluid domain mesh adopts a hexahedral structured network;

[0018] An expansion layer is provided at the fluid-solid interface;

[0019] The solid domain mesh is a hexahedral mesh.

[0020] In some embodiments, the multiphysics coupling settings of the bidirectional fluid-structure interaction model satisfy at least one of the following:

[0021] A pressure-based solver and Realizable k-ε turbulence model are used for the fluid domain;

[0022] Set up the gravity field for the fluid domain;

[0023] Activate the energy equation for the fluid domain;

[0024] A transient dynamics analysis method is used for the solid domain;

[0025] For the solid domain, boundary conditions are set, including fixed supports at the inlet and outlet, and the inner wall is set as the fluid-solid interface;

[0026] For the solid domain, a loading condition is set up that includes a periodic pressure applied at the shaker locations.

[0027] In some embodiments, the above-mentioned method of "dynamically correcting the flow data of the Coriolis mass flowmeter based on the real-time relationship" includes: dynamically correcting the flow data of the Coriolis mass flowmeter based on the real-time relationship using a temperature expansion compensation algorithm.

[0028] In some embodiments, the above-mentioned method of "determining the modal analysis results of the Coriolis mass flowmeter based on the dynamically corrected flow data and the nonlinear relationship between the straight pipe length and the phase difference" includes: obtaining each order of natural frequencies based on the dynamically corrected flow data and the nonlinear relationship between the straight pipe length and the phase difference; analyzing the modal vibration characteristics based on each order of natural frequencies; selecting the second-order modal frequency as the excitation frequency for matching verification; and determining the modal analysis results based on the matching verification results of the second-order mode.

[0029] In some embodiments, the method further includes: obtaining an initial flow field distribution based on dynamically corrected flow data and the nonlinear relationship between the straight pipe length and the phase difference; determining field quantity information based on the initial flow field distribution; calculating the stress field under the action of the fluid based on the field quantity information; obtaining the initial deformation of the structure based on the stress field; and determining the steady-state analysis results based on the initial deformation of the structure.

[0030] In some embodiments, the method further includes: performing transient coupling calculations based on dynamically corrected flow data and the nonlinear relationship between straight pipe length and phase difference; optimizing transient coupling calculation results using cubic spline interpolation; extracting time domain characteristics of the optimized transient coupling calculation results; determining phase difference analysis results and frequency domain energy analysis results based on the time domain characteristics; and determining vibration data analysis results of the transient coupling calculation based on the phase difference analysis results and the frequency domain energy analysis results.

[0031] In some embodiments, the time step size in the transient coupling calculation is set to between one tenth and one twentieth of the modal frequency.

[0032] In a second aspect, the present application provides a modal analysis device for a Coriolis mass flowmeter, the device comprising: a processing module.

[0033] The processing module is used to determine the real-time relationship between fluid pressure and structural stiffness through a bidirectional fluid-solid coupling model; the processing module is also used to dynamically correct the flow data of the Coriolis mass flowmeter based on the real-time relationship; the processing module is also used to determine the modal analysis results of the Coriolis mass flowmeter based on the dynamically corrected flow data and the nonlinear relationship between the straight pipe length and the phase difference.

[0034] In some embodiments, the parameter settings of the bidirectional fluid-structure coupling model satisfy at least one of the following:

[0035] The pipe diameter and wall thickness of the measuring pipe are 15 mm and 1.5 mm, respectively.

[0036] The structural parameters of the vibration damping plate are set according to the actual size;

[0037] The material of the measuring tube is 316 stainless steel;

[0038] The elastic modulus of the measuring tube is Gigapascal;

[0039] Poisson's ratio is 0.275;

[0040] Density is 7980 kg / m3;

[0041] The fluid medium is methane at 300 Kelvin;

[0042] The fluid domain mesh adopts a hexahedral structured network;

[0043] An expansion layer is provided at the fluid-solid interface;

[0044] The solid domain mesh is a hexahedral mesh.

[0045] In some embodiments, the multiphysics coupling settings of the bidirectional fluid-structure interaction model satisfy at least one of the following:

[0046] A pressure-based solver and Realizable k-ε turbulence model are used for the fluid domain;

[0047] Set up the gravity field for the fluid domain;

[0048] Activate the energy equation for the fluid domain;

[0049] A transient dynamics analysis method is used for the solid domain;

[0050] For the solid domain, boundary conditions are set, including fixed supports at the inlet and outlet, and the inner wall is set as the fluid-solid interface;

[0051] For the solid domain, a loading condition is set up that includes a periodic pressure applied at the shaker locations.

[0052] In some embodiments, the processing module is specifically configured to dynamically correct the flow data of the Coriolis mass flowmeter using a temperature expansion compensation algorithm based on a real-time relationship.

[0053] In some embodiments, the device further includes: an acquisition module for acquiring each order of natural frequencies based on the dynamically corrected flow rate data and the nonlinear relationship between the straight pipe length and the phase difference; a processing module for analyzing modal vibration characteristics based on each order of natural frequencies; the processing module for selecting a second-order modal frequency as an excitation frequency for matching verification; and the processing module for determining a modal analysis result based on the matching verification result of the second-order mode.

[0054] In some embodiments, the acquisition module is also used to obtain the initial flow field distribution based on the dynamically corrected flow data and the nonlinear relationship between the straight pipe length and the phase difference; the processing module is also used to determine the field quantity information based on the initial flow field distribution; the processing module is also used to calculate the stress field under the action of the fluid based on the field quantity information; the acquisition module is also used to obtain the initial deformation of the structure based on the stress field; the processing module is also used to determine the steady-state analysis results based on the initial deformation of the structure.

[0055] In some embodiments, the processing module is also used to perform transient coupling calculations based on the dynamically corrected flow data and the nonlinear relationship between the straight pipe length and the phase difference; the processing module is also used to optimize the transient coupling calculation results using cubic spline interpolation; the processing module is also used to extract the time domain characteristics of the optimized transient coupling calculation results; the processing module is also used to determine the phase difference analysis results and the frequency domain energy analysis results based on the time domain characteristics; the processing module is also used to determine the vibration data analysis results of the transient coupling calculation based on the phase difference analysis results and the frequency domain energy analysis results.

[0056] In some embodiments, the time step size in the transient coupling calculation is set to between one tenth and one twentieth of the modal frequency.

[0057] In a third aspect, the present application provides a modal analysis device for a Coriolis mass flowmeter, the device comprising: a processor and a memory, the processor and the memory being coupled, the memory being used to store one or more programs, the one or more programs comprising computer-executable instructions, and when the modal analysis device for the Coriolis mass flowmeter is running, the processor executes the computer-executable instructions stored in the memory to implement any of the modal analysis methods for the Coriolis mass flowmeter described in the first aspect above.

[0058] In a fourth aspect, the present application provides a computer-readable storage medium having instructions stored therein. When the instructions are executed on a computer, the computer executes the modal analysis method for a Coriolis mass flowmeter described in any one of the first aspects above.

[0059] In a fifth aspect, the present application provides a computer program product applied to a modal analysis device of a Coriolis mass flowmeter, the computer program product including computer instructions. When the computer instructions are executed on the modal analysis device of the Coriolis mass flowmeter, the modal analysis device of the Coriolis mass flowmeter implements any modal analysis method of the Coriolis mass flowmeter described in the first aspect above.

[0060] In the above scheme, the technical problems that can be solved and the technical effects achieved by the modal analysis device, equipment, computer storage medium or computer program product of the Coriolis mass flowmeter can be referred to the technical problems and technical effects solved by the first aspect above, and will not be repeated here. BRIEF DESCRIPTION OF THE DRAWINGS

[0061] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the following briefly introduces the drawings required for use in the description of the embodiments. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without any creative work.

[0062] Figure 1 A schematic structural diagram of a modal analysis system for a Coriolis mass flowmeter provided in an embodiment of the present application;

[0063] Figure 2 A schematic flow chart of a modal analysis method for a Coriolis mass flowmeter provided in an embodiment of the present application;

[0064] Figure 3 A schematic structural diagram of a modal analysis device for a Coriolis mass flowmeter provided in an embodiment of the present application;

[0065] Figure 4 A schematic structural diagram of a modal analysis device for a Coriolis mass flowmeter provided in an embodiment of the present application;

[0066] Figure 5 A conceptual partial view of a computer program product provided in an embodiment of the present application. DETAILED DESCRIPTION

[0067] The following will be combined with the drawings in the embodiments of this application to clearly and completely describe the technical solutions in the embodiments of this application. Obviously, the embodiments described are only part of the embodiments of this application, not all of the embodiments. Based on the embodiments in this application, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of this application.

[0068] In the description of this application, it should be understood that the terms "upper," "lower," "left," "right," "front," "back," "inner," "outer," and the like, indicating directions or positional relationships, are based on the directions or relative positional relationships shown in the accompanying drawings and are intended solely to facilitate the description of this application and simplify the description. They do not indicate or imply that the devices or components referred to must have a specific direction, be constructed, or operate in a specific direction. Therefore, they should not be construed as limitations on this application. Unless otherwise specified, the above-mentioned directionality descriptions may be flexibly set in actual application, provided that the relative positional relationships shown in the accompanying drawings are met.

[0069] The terms "first" and "second" are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of the technical features being referred to. Thus, a feature specified as "first" or "second" may explicitly or implicitly include one or more of such features. Throughout this application, unless otherwise specified, "plurality" means two or more.

[0070] In the description of this application, it should be noted that, unless otherwise expressly specified or limited, the terms "mounted," "connected," "connected," and "connected" should be understood broadly. For example, they may refer to fixed connections, detachable connections, or integral connections. They may be directly connected, indirectly connected through an intermediary, or internally connected between two components. Those skilled in the art will understand the specific meanings of these terms in this application based on the specific circumstances.

[0071] In the embodiments of the present application, the terms "comprises," "comprising," or any other variations thereof are intended to encompass non-exclusive inclusion, such that a process, article, or device comprising a series of elements includes not only those elements but also other elements not explicitly listed, or elements inherent to such process, article, or device. In the absence of further limitations, an element defined by the phrase "comprising a ..." does not preclude the presence of other identical elements in the process, article, or device comprising the element.

[0072] In the embodiments of this application, words such as "exemplary" or "for example" are used to indicate examples, illustrations, or descriptions. Any embodiment or design described as "exemplary" or "for example" in the embodiments of this application should not be interpreted as being preferred or advantageous over other embodiments or designs. Rather, the use of words such as "exemplary" or "for example" is intended to present the relevant concepts in a concrete manner.

[0073] In the description of this specification, specific features, structures, materials or characteristics may be combined in an appropriate manner in any one or more embodiments or examples.

[0074] A Coriolis mass flowmeter (CMF) operates based on the Coriolis effect. The core principle is that when a fluid flows through a vibrating pipe, the Coriolis force creates a phase or amplitude difference proportional to the fluid's mass flow rate at both ends of the pipe. By detecting this difference and combining it with parameters such as the pipe's natural frequency and mode shape, the fluid's mass flow rate can be calculated.

[0075] Modal analysis: Through theoretical calculation or experimental analysis, the natural frequency and corresponding vibration mode of the measuring pipeline are obtained, laying the foundation for fluid-structure coupling vibration analysis.

[0076] Fluid-Solid Interaction Analysis: Considering the interaction between fluid and solid structure, computational fluid dynamics (CFD) and finite element analysis (FEA) methods are used to perform fluid-solid interaction analysis to simulate and measure the vibration characteristics of the pipeline under the action of fluid flow.

[0077] Signal detection and processing: Install a vibration pickup on the measuring pipe to detect the vibration signal caused by fluid flow, and extract phase difference or time difference information through signal processing algorithms.

[0078] Mass flow calculation: Calculate the mass flow of the fluid based on the detected phase difference or time difference and the flow coefficient obtained from the calibration experiment.

[0079] In the existing technology, the accuracy and adaptability of Coriolis mass flowmeters are usually improved based on the following methods:

[0080] 1. Flow field modeling: Establish a steady-state flow field model in Fluent (a computational fluid dynamics software), set the fluid properties (such as natural gas density ρ = 1.13 kg / m 3 ) and boundary conditions.

[0081] 2. Prestress loading: Map the flow field pressure distribution to the solid domain and generate a prestress file.

[0082] 3. Modal calculation: Apply fixed constraints in ANSYS Mechanical (a finite element analysis software), calculate the first six modes, and extract the second-order frequency.

[0083] 4. Phase difference estimation: Obtain the empirical relationship table between structural parameters and ΔT through orthogonal experiments.

[0084] However, limitations of existing technologies include:

[0085] Static coupling: Dynamic stiffness changes caused by pressure and temperature changes are not considered (the phase difference shift is as much as 3.2 microseconds (μs) under a pressure of 5 megapascals (MPa)).

[0086] Deviation of the empirical model: The linear assumption error between the structural parameters and ΔT is large (L7 every ±10 millimeters (mm) causes ΔT to change by ±1.2μs, but it is actually nonlinear).

[0087] Multi-field decoupling: Fluid density and viscosity changes are not fed back to modal analysis in real time.

[0088] That is to say, the existing technology has the following problems:

[0089] 1. Modal prediction deviation problem:

[0090] Traditional methods, based on empty tube structural analysis (e.g., the second-order frequency of an empty tube is 101.73 Hz), fail to account for the added mass effect of the fluid. This results in a prediction error of over 6% in the actual natural frequency of the flow-carrying tube (approximately 95 Hz). This causes the excitation frequency to deviate from the resonant point, leading to a phase difference measurement error exceeding 5%, directly impacting the accuracy of mass flow calculations.

[0091] 2. Poor adaptability to dynamic working conditions:

[0092] Pressure and temperature changes significantly affect the stiffness of the measuring tube (for example, a phase shift of 3.2 μs occurs when the pressure increases from 2 MPa to 5 MPa). However, existing static models cannot correct these parameters in real time, resulting in a mass flow calculation error of up to 1.8%. The lack of a real-time compensation mechanism under dynamic conditions is a key bottleneck restricting measurement accuracy.

[0093] 3. Problems of unmodeled nonlinear effects of structural parameters:

[0094] The nonlinear effects of parameters such as straight pipe length (e.g., every ±10mm of L7 results in a ±1.2μs change in phase difference) and pipe wall thickness are not precisely quantified, resulting in poor consistency in phase difference predictions across multiple flowmeter specifications (error >1.5%). Traditional models, due to simplified assumptions, cannot account for the impact of complex structural parameters.

[0095] Among them, when designing and analyzing traditional Coriolis mass flowmeter models, in order to simplify calculations or improve efficiency, some complex physical phenomena or structural parameters are usually simplified or ignored, resulting in the model being unable to fully reflect the measurement accuracy issues under actual working conditions.

[0096] In summary, the causal reasoning of the shortcomings of existing technologies includes:

[0097] 1. The root cause of modal prediction deviation:

[0098] Cause and effect chain: empty tube modal analysis → ignoring the added mass of the fluid → the calculated value of the natural frequency is too high → the excitation frequency deviates from the resonance point → the phase difference measurement error increases.

[0099] Data support: The empty tube frequency is 101.73Hz relative to the actual frequency of the current-carrying tube 95Hz (error 6%), and the phase difference error is greater than 5%.

[0100] 2. Causes of poor adaptability to dynamic working conditions:

[0101] Causal chain: static pressure assumption → uncorrelated pressure-stiffness dynamics → increased pressure causing increased pipe wall stress → nonlinear phase shift → accumulated flow calculation errors.

[0102] Data support: When the pressure increases from 2MPa to 5MPa, the phase difference shifts by 3.2μs, with an error of 1.8%.

[0103] 3. Consequences of not modeling the nonlinear effects of structural parameters:

[0104] Cause and effect chain: Empirical formula linear assumption → Ignoring the nonlinear relationship between straight pipe length and phase difference → Parameter changes lead to prediction deviation → Poor consistency of multi-specification flow meters.

[0105] Data support: Every ±10mm of L7 causes a ΔT change of ±1.2μs, the model error is ±0.7μs, and the prediction error of multiple specifications is >1.5%.

[0106] To address these issues, this application proposes a new modal analysis method for Coriolis mass flowmeters. This method determines fluid-structure coupling by correlating fluid pressure with structural stiffness in real time, and then compensates and corrects the Coriolis mass flowmeter's flow data based on the dynamic fluid-structure coupling. Next, based on the corrected flow data, the modal analysis results of the Coriolis mass flowmeter are determined by combining the nonlinear relationship between straight pipe length and phase difference. This method, by referencing fluid-structure coupling and quantifying the nonlinear effects of structural parameters, can provide accurate results for the modal analysis of the Coriolis mass flowmeter.

[0107] In this regard, the present application provides a modal analysis system for a Coriolis mass flowmeter, which is applied to a modal analysis device for a Coriolis mass flowmeter, such as Figure 1 As shown, it includes: model establishment and parameter setting module, calculation method setting module, fluid-solid coupling analysis module, and data processing module.

[0108] 1. The model building and parameter setting module is used to realize geometric model construction and mesh division.

[0109] The geometric model construction includes the setting of modeling parameters and material parameters, and the meshing includes the division of fluid domain meshes and solid domain meshes.

[0110] 2. The calculation method setting module is used to implement fluid domain calculation settings and solid domain calculation settings.

[0111] The fluid domain calculation settings include the solver parameter settings and the turbulence model parameter settings, and the solid domain calculation settings include the boundary condition settings and the analysis type settings.

[0112] 3. The fluid-solid coupling analysis module is used to realize steady-state analysis, modal analysis and transient coupling calculation.

[0113] Among them, steady-state analysis includes the flow field calculation process and the prestress analysis process, modal analysis includes the modal calculation process combined with fluid action, and transient coupling calculation includes the setting of time step and data transmission.

[0114] 4. The data processing module is used to realize vibration data processing.

[0115] Among them, vibration data processing includes cubic spline interpolation process, time difference calculation process, and phase difference calculation process.

[0116] After introducing the application scenario and implementation environment of the embodiment of the present application, the modal analysis method of the Coriolis mass flowmeter provided by the embodiment of the present application is introduced in detail in combination with the above-mentioned implementation environment.

[0117] The methods in the following embodiments can all be implemented in the above application scenarios and implementation environments. The embodiments of the present application are described in detail below with reference to the accompanying drawings.

[0118] Figure 2 The following is a flow chart of a modal analysis method for a Coriolis mass flowmeter provided in an embodiment of the present application. Figure 2 As shown, the method is applied to a modal analysis device of a Coriolis mass flowmeter, including:

[0119] S201. Determine the real-time relationship between fluid pressure and structural stiffness through a bidirectional fluid-structure coupling model.

[0120] In some embodiments, the parameter settings of the bidirectional fluid-structure coupling model may satisfy at least one of the following:

[0121] The pipe diameter and wall thickness of the measuring pipe are 15 mm and 1.5 mm, respectively.

[0122] The structural parameters of the vibration damping plate are set according to the actual size;

[0123] The material of the measuring tube is 316 stainless steel;

[0124] The elastic modulus of the measuring tube is Gigapascal;

[0125] Poisson's ratio is 0.275;

[0126] Density is 7980 kg / m3;

[0127] The fluid medium is methane at 300 Kelvin;

[0128] The fluid domain mesh adopts a hexahedral structured network;

[0129] An expansion layer is provided at the fluid-solid interface;

[0130] The solid domain mesh is a hexahedral mesh.

[0131] As a possible implementation method, in the process of constructing a geometric model of a modal analysis device for a Coriolis mass flowmeter, a three-dimensional model of a double-U-type Coriolis flowmeter can be established using SolidWorks (a three-dimensional solid model creation tool).

[0132] Among them, the structural parameters of the measuring tube are: tube diameter 15mm, wall thickness 1.5mm;

[0133] Straight pipe section L7: 410-430mm adjustable;

[0134] Straight pipe section L8: 191-211mm adjustable;

[0135] The structural parameters of the vibration damping plate are set according to the actual size.

[0136] In addition, the modal analysis equipment of the Coriolis mass flowmeter can also set the material parameters:

[0137] Measuring tube material: 316 stainless steel;

[0138] Elastic modulus: 193 GPa;

[0139] Poisson's ratio: 0.275;

[0140] Density: 7980 kg / m3;

[0141] Fluid medium: 300 Kelvin (K) methane.

[0142] In addition, the modal analysis equipment of the Coriolis mass flowmeter can also divide the grid:

[0143] 1. Fluid domain mesh:

[0144] 1.1, using hexahedral structured grid;

[0145] 1.2. Set up an expansion layer at the fluid-solid interface;

[0146] Among them, adding an expansion layer is mainly used in fluid domain meshing, with the aim of improving the accuracy of fluid flow simulation, especially the flow details at the fluid-solid interface.

[0147] 1.3. The inlet and outlet surfaces are divided using the sweeping method.

[0148] 2. Solid domain mesh:

[0149] 2.1 Using hexahedral grid;

[0150] 2.2. Carry out expansion treatment on the inner wall of the pipeline;

[0151] Among them, the expansion treatment of the inner wall of the pipe is mainly used for solid domain meshing, with the aim of improving the calculation accuracy of the stress distribution in the solid domain.

[0152] 2.3. Ensure that the mesh quality meets the calculation requirements.

[0153] Exemplarily, the specific steps include:

[0154] Step 1. Import the 3D model: Use Meshing software to import the 3D model of the Coriolis flowmeter to ensure that the geometric structure is complete and correct.

[0155] This provides the basic conditions for subsequent meshing and other operations.

[0156] Step 2: Set geometric parameters: Determine the geometric dimensions of the fluid domain, including the straight section length of the measuring tube (L7, L8), bend radius, tube diameter, wall thickness, etc.

[0157] In this way, setting geometric parameters based on the imported model can ensure that the size and structure of the model meet actual requirements.

[0158] Step 3. Select the mesh type: Select hexahedral structured mesh and set an expansion layer at the fluid-solid interface.

[0159] In this way, choosing the appropriate mesh type (such as hexahedral structured mesh) and setting the expansion layer determine the quality and accuracy of the mesh division.

[0160] Step 4. Grid division: Grid the fluid domain to ensure grid quality.

[0161] In this way, the fluid domain is meshed according to the selected mesh type and settings. This step is the core operation of meshing.

[0162] Step 5. Set the expansion layer: Set the expansion layer at the fluid-solid interface to capture the flow details of the fluid at the interface.

[0163] In this way, an expansion layer is set at the fluid-solid interface to capture the flow details of the fluid at the interface. This step needs to be explicitly set before meshing.

[0164] Step 6. Divide the inlet and outlet surfaces: Use the sweeping method to divide the fluid inlet and outlet surfaces to ensure the continuity and quality of the grid.

[0165] In this way, the fluid inlet and outlet surfaces are meshed using a sweeping method to ensure the continuity and quality of the mesh. This step depends on the preliminary results of meshing.

[0166] Step 7. Check the mesh quality: Check the element quality of the mesh to ensure that there are no distorted or undersized elements.

[0167] In this way, the element quality of the mesh is checked to ensure that there are no distorted or undersized elements. This step verifies the correctness of all the previous steps.

[0168] Step 8. Export the mesh file: After completing the mesh division, export the mesh file for subsequent fluid dynamics simulation calculations.

[0169] In this way, after the meshing is completed, the mesh file is exported for subsequent fluid dynamics simulation calculations. This step is the final output of the entire process.

[0170] In some embodiments, the multiphysics coupling settings of the bidirectional fluid-structure interaction model may satisfy at least one of the following:

[0171] A pressure-based solver and Realizable k-ε turbulence model are used for the fluid domain;

[0172] Set up the gravity field for the fluid domain;

[0173] Activate the energy equation for the fluid domain;

[0174] A transient dynamics analysis method is used for the solid domain;

[0175] For the solid domain, boundary conditions are set, including fixed supports at the inlet and outlet, and the inner wall is set as the fluid-solid interface;

[0176] For the solid domain, a loading condition is set up that includes a periodic pressure applied at the shaker locations.

[0177] As a possible implementation method, the modal analysis equipment of the Coriolis mass flowmeter can set the calculation method during the geometric model construction process, including:

[0178] 1. Fluid domain calculation settings, including:

[0179] 1.1 Solver settings:

[0180] 1.1.1. Select the pressure-based solver;

[0181] 1.1.2. Turn on the gravity option;

[0182] 1.1.3. Operation density is set to zero;

[0183] 1.1.4. Enable the energy equation.

[0184] 1.2 Turbulence model settings:

[0185] 1.2.1. Select the k-εRealizable model;

[0186] 1.2.2, using scalable wall functions;

[0187] 1.2.3 Turn on the solution stability option.

[0188] 2. Solid domain calculation settings, including:

[0189] 2.1 Boundary conditions:

[0190] 2.1.1、Add fixed support to the entrance and exit;

[0191] Among them, when setting boundary conditions, applying fixed supports to the entrances and exits in the solid domain model usually means restricting the displacement of these areas to simulate the situation in actual engineering applications where certain parts are rigidly fixed.

[0192] 2.1.2. The inner wall is set as the fluid-solid interface;

[0193] 2.1.3. Apply periodic pressure at the vibrator position.

[0194] 2.2 Analysis Type:

[0195] 2.2.1, Steady-state structural analysis;

[0196] 2.2.2 Modal analysis;

[0197] 2.2.3. Transient dynamic analysis.

[0198] Among them, transient dynamics analysis mainly studies the dynamic response of the structure under time-varying loads, focusing on the changes in the acceleration, velocity and displacement of the structure over time.

[0199] It should be noted that the multiphysics coupling settings described above are a collection of configuration items, describing the parameters and conditions that need to be configured when performing calculations in the fluid and solid domains. There is no strict order of precedence and they can be viewed as a parallel configuration process.

[0200] In other words, the interaction between fluid and solid structure: when the fluid flows, it exerts forces on the surrounding solid structures, causing them to deform or move; in turn, the deformation or movement of the solid structures changes the flow pattern of the fluid. The present embodiment addresses the lack of fluid-solid coupling by introducing a bidirectional fluid-solid coupling model that correlates fluid pressure with structural stiffness in real time, reducing the modal frequency prediction error from 6% to ≤3% and improving the excitation frequency matching by 60%.

[0201] S202. Dynamically correct the flow data of the Coriolis mass flowmeter based on the real-time relationship.

[0202] As a possible implementation method, the modal analysis equipment of the Coriolis mass flowmeter can dynamically correct the flow data of the Coriolis mass flowmeter based on the real-time relationship and the temperature expansion compensation algorithm.

[0203] In other words, by establishing a pressure-stiffness correlation model (σ=β·P) and a temperature expansion compensation algorithm, dynamic parameter correction is achieved, so that the phase difference offset compensation rate caused by pressure changes is greater than 90%, and the working condition adaptability error is compressed from 1.8% to 0.3%.

[0204] S203 : Determine a modal analysis result of the Coriolis mass flowmeter based on the dynamically corrected flow data and the nonlinear relationship between the straight pipe length and the phase difference.

[0205] The nonlinear relationship between the straight tube length and the phase difference is determined by establishing a weight coefficient library (k1 = 0.12 μs / mm, k2 = 0.09 μs / mm) of L7, L8, and t through parametric scanning.

[0206] It should be noted that L7 and L8 represent the lengths of the two straight sections of the measuring tube, respectively. In Coriolis mass flowmeters, the measuring tube is typically designed in a U-shape or other configuration. The length of the straight section significantly influences the vibration characteristics of the measuring tube. Variations in L7 and L8 affect the natural frequency and mode shape of the measuring tube, thus affecting the accuracy of phase and time difference measurements.

[0207] t represents the wall thickness of the measuring tube. Wall thickness affects the stiffness and natural frequency of the measuring tube. Thicker walls increase stiffness, thereby increasing its natural frequency; thinner walls decrease stiffness and lower natural frequency. Variations in wall thickness also affect the accuracy of phase and time difference measurements.

[0208] In other words, by combining the nonlinear effects of quantified structural parameters, the phase difference prediction accuracy can reach ±0.5μs, which is three times higher than the traditional model, and the consistency error of multiple specifications is less than 0.5%, thereby improving the accuracy of the modal analysis results of the Coriolis mass flowmeter.

[0209] In some embodiments, during the process (i.e., S203 above) in which the modal analysis device of the Coriolis mass flowmeter determines the modal analysis results of the Coriolis mass flowmeter based on the dynamically corrected flow data and the nonlinear relationship between the straight pipe length and the phase difference, the modal analysis device of the Coriolis mass flowmeter can obtain each order of natural frequencies based on the dynamically corrected flow data and the nonlinear relationship between the straight pipe length and the phase difference, and analyze the modal vibration characteristics based on each order of natural frequencies. The modal analysis device of the Coriolis mass flowmeter can then select the second-order modal frequency as the excitation frequency for matching verification, and determine the modal analysis results based on the matching verification results of the second-order mode.

[0210] For example, in the process of determining the modal analysis results through fluid-structure coupling analysis, the modal calculation of fluid action can be considered, including:

[0211] 1. Obtaining the natural frequencies of each order: Performing modal decomposition under fluid load is the basis of modal analysis and must be completed before subsequent analysis can be performed.

[0212] 2. Analyze the characteristics of the modal vibration mode: Based on the natural frequency, further analyze the characteristics of the modal vibration mode to provide a basis for selecting the excitation frequency.

[0213] 3. Select the second-order modal frequency as the excitation frequency: Based on the analysis results of the modal vibration shape, select the second-order modal frequency that best matches the actual vibration state as the excitation frequency.

[0214] As can be understood, by correlating fluid pressure with structural stiffness in real time, the fluid-structure interaction is determined, and the flow data of the Coriolis mass flowmeter is compensated based on the dynamic fluid-structure interaction. Next, based on this corrected flow data, the modal analysis results of the Coriolis mass flowmeter are determined by combining the nonlinear relationship between straight pipe length and phase difference. In this way, by considering the fluid-structure interaction and quantifying the nonlinear effects of structural parameters, the accuracy of the modal analysis results of the Coriolis mass flowmeter can be improved.

[0215] In some embodiments, the modal analysis device of the Coriolis mass flowmeter can also obtain an initial flow field distribution based on the dynamically corrected flow data and the nonlinear relationship between straight pipe length and phase difference, and determine field quantity information based on the initial flow field distribution. The modal analysis device of the Coriolis mass flowmeter can then calculate the stress field under the action of the fluid based on the field quantity information, obtain the initial structural deformation based on the stress field, and further determine the steady-state analysis results based on the initial structural deformation.

[0216] For example, the process of determining steady-state analysis results through fluid-structure coupling analysis may include:

[0217] 1. Flow field calculation:

[0218] 1.1. Obtain the initial flow field distribution;

[0219] 1.2. Output pressure, velocity and other field quantities;

[0220] 1.3. As the initial condition for transient analysis.

[0221] 2. Prestress analysis:

[0222] 2.1. Calculate the stress field under the action of fluid;

[0223] Among them, this is the first step of prestress analysis, which obtains the stress distribution by calculating the force exerted by the fluid on the solid structure.

[0224] 2.2. Obtain the initial deformation of the structure.

[0225] Based on the calculated stress field, the initial deformation of the solid structure under the action of the fluid is further analyzed. The initial deformation is the structural response caused by the fluid action and provides a basis for subsequent modal analysis and transient coupling calculations.

[0226] In some embodiments, the modal analysis device of the Coriolis mass flowmeter can also perform transient coupling calculations based on the dynamically corrected flow data and the nonlinear relationship between straight pipe length and phase difference, and optimize the transient coupling calculation results using cubic spline interpolation. The modal analysis device can then extract the time domain characteristics of the optimized transient coupling calculation results, determine phase difference analysis results and frequency domain energy analysis results based on the time domain characteristics, and further determine the vibration data analysis results of the transient coupling calculation based on the phase difference analysis results and frequency domain energy analysis results.

[0227] The time step in the transient coupling calculation is set between one tenth and one twentieth of the modal frequency.

[0228] Exemplarily, the process of determining the vibration data analysis results of the transient coupling calculation through the fluid-structure coupling analysis may include:

[0229] 1. Transient coupling calculation:

[0230] 1.1, Time step setting:

[0231] 1.1.1. Take 1 / 10 to 1 / 20 of the modal frequency;

[0232] 1.1.2. Ensure calculation stability.

[0233] It should be noted that transient coupling calculation is a dynamic process that determines and studies the interaction between fluid and solid over time. It focuses not only on the dynamic response of the structure, but also on the dynamic characteristics of the fluid flow and the interaction between the two.

[0234] 1.2 Data transmission settings:

[0235] 1.2.1、Transmission of force;

[0236] 1.2.2, Transmission of displacement;

[0237] 1.2.3. Ensure convergence.

[0238] 2. Data processing method:

[0239] 2.1. Vibration data processing:

[0240] 2.1.1, using cubic spline interpolation:

[0241] 2.1.1.1. Process the vibration curve of the monitoring point;

[0242] The method is realized through cubic spline interpolation, and the specific steps include data acquisition, preprocessing, cubic spline interpolation and curve generation.

[0243] 2.1.1.2. Improve the accuracy of time difference calculation.

[0244] It should be noted that the use of cubic spline interpolation is the first step in data processing. The vibration curve of the monitoring point is processed by the interpolation method to make the data smoother and more accurate, providing a more accurate basis for the subsequent time difference calculation.

[0245] 2.1.2, Time difference calculation:

[0246] 2.1.2.1. Take two adjacent points with zero displacement;

[0247] Among them, in the process of calculating the time difference, it is necessary to first find two adjacent points with zero displacement, which is the basis for calculating the time difference.

[0248] 2.1.2.2. Calculate the time difference Δt.

[0249] Here, based on the two points found, the time difference Δt between them is calculated. This step depends on the results of the previous step.

[0250] It should be noted that the time difference calculation can be based on the interpolated vibration curve, finding two adjacent points with zero displacement and calculating the time difference Δt. This step directly relies on the interpolated data.

[0251] 2.1.3 Phase difference calculation:

[0252] 2.1.3.1, ΔΘ=ω·Δt=2πf·Δt;

[0253] 2.1.3.2. Analyze the changing pattern of phase difference.

[0254] The phase difference variation law is the phase difference variation law of the vibration signal of the measuring tube of the Coriolis mass flowmeter.

[0255] It should be noted that the phase difference calculation can use the calculated time difference Δt, combined with the angular velocity ω or frequency f, to calculate the phase difference ΔΘ. This step depends on the result of the time difference calculation.

[0256] In conjunction with the above embodiment, during the steady-state analysis process, a steady-state analysis of the fluid is first performed to obtain the initial flow field distribution, which provides initial conditions for the subsequent transient analysis.

[0257] During prestress analysis, the fluid pressure distribution obtained from steady-state analysis is mapped to the solid domain for prestress analysis. This step converts fluid action into prestress, providing data that is more consistent with actual working conditions for modal analysis.

[0258] During the modal analysis process, modal analysis is performed on the basis of prestressing to obtain the natural frequencies and vibration modes of each order. This step determines the vibration characteristics of the structure and provides the excitation frequency for subsequent transient coupling calculations.

[0259] During the transient coupling calculation process, the time step and data transmission are set to perform the transient coupling calculation. This step requires ensuring that the time step and total calculation time are set appropriately to ensure the stability and convergence of the calculation.

[0260] During data processing, the time and phase differences are extracted from the vibration curves at the monitoring points through interpolation and phase difference calculation. This step processes the results of the transient coupling calculation and is used to analyze the measurement performance of the flowmeter.

[0261] In summary, the modal analysis method for the Coriolis mass flowmeter provided in the embodiments of the present application may include:

[0262] 1. Fluid-solid coupling analysis method:

[0263] 1.1. Bidirectional fluid-structure coupling using the separation solution method;

[0264] 1.2. Dynamic mesh technology and fluid-solid interface treatment;

[0265] 1.3. Bidirectional data transmission mechanism of force and displacement;

[0266] 1.4. Convergence control strategy.

[0267] 2. Prestressed modal analysis technology:

[0268] 2.1. Obtain the initial field by steady-state flow field analysis;

[0269] 2.2, fluid action is converted into prestress;

[0270] 2.3. Modal calculation considering fluid action;

[0271] 2.4. Natural frequency and vibration mode prediction.

[0272] 3. Vibration data processing technology:

[0273] 3.1, cubic spline interpolation algorithm;

[0274] 3.2, High-precision time difference calculation;

[0275] 3.3. Phase difference analysis method.

[0276] In other words, in the modal analysis method of multi-physics field coupling, it is possible to build a geometric model based on shared topology, mesh the expansion layer of the fluid-solid interface, complete the analysis chain of steady-state flow field → prestressed analysis → modal analysis → transient coupling, set the association between time step and modal frequency, and establish a bidirectional fluid-solid coupling data transmission mechanism.

[0277] The vibration data processing method may be based on vibration curve processing of cubic spline interpolation, time difference calculation method, and phase difference analysis method.

[0278] In the calculation parameter setting method, pressure-based solver configuration, k-εRealizable turbulence model parameter setting, dynamic mesh parameter setting, and convergence control parameter setting can be realized.

[0279] In the structural parameter optimization method, the influence analysis of straight pipe length, pipe wall thickness and vibration pickup position optimization can be realized.

[0280] In the operating parameter influence law analysis method, flow-phase difference relationship prediction, pressure influence analysis, and temperature effect compensation can be achieved.

[0281] Optionally, in terms of geometric model construction, the measuring tube structure can be a single U-shaped, μ-shaped, or straight tube structure. The vibration damping structure can also adopt a variety of forms, such as single plate, multi-plate, frame, or composite. Meshing can be done using tetrahedral, hybrid, or unstructured meshes. Interface processing can be done using either a consistent mesh or a non-conformant mesh with interface interpolation, or using dynamic mesh reconstruction technology.

[0282] Form provides the design framework for sizing, which in turn influences its implementation and performance. The diverse range of vibration-damping structures emphasizes the flexibility of form, while the actual sizing requires detailed planning and optimization based on specific design needs and performance requirements.

[0283] Interface treatment is equivalent to the treatment method of the fluid-solid interface in fluid-solid coupling analysis.

[0284] Optionally, for fluid calculations, turbulence models can include the standard k-ε model, RNG k-ε model, k-ω model, SST model, or large eddy simulation (LES). Solution algorithms can include SIMPLE, SIMPLEC, PISO, or a segregated solver. For solid structure analysis, explicit dynamics, implicit dynamics, harmonic response, or random vibration analysis can be used. Modal extraction can utilize the Block Lanczos method, PCG Lanczos method, subspace method, or QR damping method.

[0285] Among them, the solution algorithm is a specific numerical method used to solve the fluid dynamics equations, which can be pressure-based or density-based.

[0286] A segregated solver is a solution strategy that solves the fluid and solid equations separately and transfers data through the fluid-solid interface. It is a more general solution method that can be combined with a pressure-based solver or other types of solvers.

[0287] Solid structure analysis is equivalent to the fluid-solid coupling analysis method mentioned above.

[0288] Alternatively, for fluid-structure interaction analysis, direct solutions can be used to simultaneously solve the coupled equations, or variants of segregated solutions such as strong coupling iteration, weak coupling iteration, explicit coupling, and semi-implicit coupling can be used. Data transfer can be performed using methods such as radial basis function interpolation, polynomial interpolation, inverse distance weighted interpolation, or least squares interpolation. Interface mapping can be performed using node-to-node mapping, node-to-element mapping, Gaussian point mapping, or conservative interpolation.

[0289] Among them, node-to-node mapping refers to mapping the nodes of the fluid domain directly to the nodes of the solid domain at the fluid-solid interface to ensure accurate transmission of data between nodes.

[0290] Optionally, vibration data processing can employ curve fitting methods such as polynomial fitting, B-spline fitting, Bessel curve fitting, Fourier series fitting, or wavelet analysis; signal processing can utilize techniques such as Hilbert transform, wavelet transform, empirical mode decomposition, adaptive filtering, or neural network processing. Parameter optimization can utilize methods such as genetic algorithms, particle swarm optimization, response surface methodology, neural network optimization, or grey system theory. Sensitivity analysis can utilize methods such as local sensitivity analysis, global sensitivity analysis, analysis of variance, or orthogonal experimental design.

[0291] Among them, signal processing is used to analyze and process the vibration signal of the Coriolis flowmeter to extract key information such as phase difference and time difference.

[0292] Parameter optimization is used to optimize the design parameters of the Coriolis flowmeter, such as structural dimensions and material properties, to improve the performance of the flowmeter.

[0293] Sensitivity analysis is used to evaluate the sensitivity of the Coriolis flowmeter performance to changes in different parameters to determine which parameters have the greatest impact on measurement accuracy.

[0294] The above mainly introduces the solution provided by the embodiment of the present application from the perspective of computer equipment. It is understandable that, in order to realize the above functions, the computer equipment includes hardware structures and / or software modules corresponding to the execution of each function. Those skilled in the art should easily realize that, in combination with the modal analysis method steps of the Coriolis mass flowmeter described in each example of the embodiment disclosed in this application, the present application can be implemented in the form of hardware or a combination of hardware and computer software. Whether a function is executed in the form of hardware or computer software driving hardware depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of this application.

[0295] The present application also provides a modal analysis device for a Coriolis mass flowmeter. The modal analysis device for a Coriolis mass flowmeter can be a computer device, a CPU in the computer device, a processing module in the computer device for analyzing the modal state of the Coriolis mass flowmeter, or a client in the computer device for analyzing the modal state of the Coriolis mass flowmeter.

[0296] In the embodiment of the present application, the modal analysis device of the Coriolis mass flowmeter can be divided into functional modules or functional units according to the above-mentioned method example. For example, each functional module or functional unit can be divided according to each function, or two or more functions can be integrated into one processing module. The above-mentioned integrated module can be implemented in the form of hardware or in the form of software functional modules or functional units. Among them, the division of modules or units in the embodiment of the present application is schematic and is only a logical functional division. There may be other division methods in actual implementation.

[0297] like Figure 3 FIG. 3 is a schematic diagram of a modal analysis device for a Coriolis mass flowmeter according to an embodiment of the present invention. The modal analysis device 300 for a Coriolis mass flowmeter is applied to a modal analysis device for a Coriolis mass flowmeter. The modal analysis device 300 for a Coriolis mass flowmeter may include: a processing module 301.

[0298] Processing module 301 is used to determine the real-time relationship between fluid pressure and structural stiffness through a bidirectional fluid-solid coupling model; processing module 301 is also used to dynamically correct the flow data of the Coriolis mass flowmeter based on the real-time relationship; processing module 301 is also used to determine the modal analysis results of the Coriolis mass flowmeter based on the dynamically corrected flow data and the nonlinear relationship between the straight pipe length and the phase difference.

[0299] Optionally, the parameter setting of the bidirectional fluid-structure coupling model satisfies at least one of the following:

[0300] The pipe diameter and wall thickness of the measuring pipe are 15 mm and 1.5 mm, respectively.

[0301] The structural parameters of the vibration damping plate are set according to the actual size;

[0302] The material of the measuring tube is 316 stainless steel;

[0303] The elastic modulus of the measuring tube is Gigapascal;

[0304] Poisson's ratio is 0.275;

[0305] Density is 7980 kg / m3;

[0306] The fluid medium is methane at 300 Kelvin;

[0307] The fluid domain mesh adopts a hexahedral structured network;

[0308] An expansion layer is provided at the fluid-solid interface;

[0309] The solid domain mesh is a hexahedral mesh.

[0310] Optionally, the multi-physics coupling setting of the bidirectional fluid-structure coupling model satisfies at least one of the following:

[0311] A pressure-based solver and Realizable k-ε turbulence model are used for the fluid domain;

[0312] Set up the gravity field for the fluid domain;

[0313] Activate the energy equation for the fluid domain;

[0314] A transient dynamics analysis method is used for the solid domain;

[0315] For the solid domain, boundary conditions are set, including fixed supports at the inlet and outlet, and the inner wall is set as the fluid-solid interface;

[0316] For the solid domain, a loading condition is set up that includes a periodic pressure applied at the shaker locations.

[0317] Optionally, the processing module 301 is specifically configured to dynamically correct the flow data of the Coriolis mass flowmeter based on the real-time relationship using a temperature expansion compensation algorithm.

[0318] Optionally, the device further includes: an acquisition module 302. The acquisition module 302 is configured to acquire each order of natural frequency based on the dynamically corrected flow data and the nonlinear relationship between the straight pipe length and the phase difference; a processing module 301 is specifically configured to analyze modal vibration characteristics based on each order of natural frequency; the processing module 301 is further configured to select a second-order modal frequency as an excitation frequency for matching verification; and the processing module 301 is further configured to determine a modal analysis result based on the matching verification result of the second-order mode.

[0319] Optionally, the acquisition module 302 is also used to obtain the initial flow field distribution based on the dynamically corrected flow data and the nonlinear relationship between the straight pipe length and the phase difference; the processing module 301 is also used to determine the field quantity information based on the initial flow field distribution; the processing module 301 is also used to calculate the stress field under the action of the fluid based on the field quantity information; the acquisition module 302 is also used to obtain the initial deformation of the structure based on the stress field; the processing module 301 is also used to determine the steady-state analysis results based on the initial deformation of the structure.

[0320] Optionally, the processing module 301 is also used to perform transient coupling calculations based on the dynamically corrected flow data and the nonlinear relationship between the straight pipe length and the phase difference; the processing module 301 is also used to optimize the transient coupling calculation results using cubic spline interpolation; the processing module 301 is also used to extract the time domain characteristics of the optimized transient coupling calculation results; the processing module 301 is also used to determine the phase difference analysis results and the frequency domain energy analysis results based on the time domain characteristics; the processing module 301 is also used to determine the vibration data analysis results of the transient coupling calculation based on the phase difference analysis results and the frequency domain energy analysis results.

[0321] Optionally, the time step in the transient coupled calculation is set to between one tenth and one twentieth of the modal frequency.

[0322] Figure 4 FIG4 is a schematic diagram illustrating a modal analysis device for a Coriolis mass flowmeter according to an exemplary embodiment. The modal analysis device for a Coriolis mass flowmeter may include a processor 402 configured to execute application code to implement the modal analysis method for a Coriolis mass flowmeter disclosed herein.

[0323] The processor 402 may be a CPU, a microprocessor, an application-specific integrated circuit (ASIC), or one or more integrated circuits for controlling the execution of the program of the present application.

[0324] like Figure 4As shown, the modal analysis device of the Coriolis mass flowmeter may further include a memory 403. The memory 403 is used to store application code for executing the solution of the present application, and is controlled by the processor 402 to execute the code.

[0325] The memory 403 may be a read-only memory (ROM) or other type of static storage device that can store static information and instructions, a random access memory (RAM) or other type of dynamic storage device that can store information and instructions, or an electrically erasable programmable read-only memory (EEPROM), a compact disc read-only memory (CD-ROM) or other optical disc storage, an optical disc storage (including a compact disc, laser disc, optical disc, digital versatile disc, Blu-ray disc, etc.), a magnetic disk storage medium or other magnetic storage device, or any other medium that can be used to carry or store desired program code in the form of instructions or data structures and can be accessed by a computer, but is not limited thereto. The memory 403 may exist independently and be connected to the processor 402 via the bus 404. The memory 403 may also be integrated with the processor 402.

[0326] like Figure 4 As shown, the modal analysis device for a Coriolis mass flowmeter may further include a communication interface 401, wherein the communication interface 401, the processor 402, and the memory 403 may be coupled to each other, for example, via a bus 404. The communication interface 401 is used to exchange information with other devices, for example, to support information exchange between the modal analysis device for the Coriolis mass flowmeter and other devices.

[0327] It should be pointed out that Figure 4 The device structure shown in the figure does not constitute a limitation on the modal analysis device of the Coriolis mass flowmeter, except Figure 4 In addition to the components shown, the modal analysis apparatus for the Coriolis mass flowmeter may include more or fewer components than shown, or may combine certain components, or have a different arrangement of components.

[0328] In actual implementation, the functions implemented by the processing module 301 can be Figure 4 The processor 402 shown calls the program code in the memory 403 to implement it.

[0329] The present application also provides a computer-readable storage medium having instructions stored thereon. When the instructions in the computer-readable storage medium are executed by a processor of a computer device, the computer device is enabled to perform the modal analysis of the Coriolis mass flowmeter provided in the above-described embodiment. For example, the computer-readable storage medium may be a memory 403 including instructions, and the instructions may be executed by the processor 402 of the computer device to perform the above-described method. Alternatively, the computer-readable storage medium may be a non-transitory computer-readable storage medium, such as a ROM, RAM, CD-ROM, magnetic tape, floppy disk, or optical data storage device.

[0330] Figure 5 A conceptual partial view of a computer program product provided by an embodiment of the present application is exemplarily shown, where the computer program product includes a computer program for executing a computer process on a computing device.

[0331] In one embodiment, the computer program product is provided using a signal bearing medium 500. The signal bearing medium 500 may include one or more program instructions that, when executed by one or more processors, may provide the above-described Figure 2 Thus, for example, reference to Figure 2 In the embodiment shown in , one or more features of S201 to S203 may be undertaken by one or more instructions associated with the signal bearing medium 500. In addition, Figure 5 The program instructions in also describe example instructions.

[0332] In some examples, the signal-bearing medium 500 may include a computer-readable medium 501, such as, but not limited to, a hard drive, a compact disk (CD), a digital video disk (DVD), a digital tape, a memory, a read-only memory (ROM), or a random access memory (RAM), and the like.

[0333] In some implementations, signal bearing medium 500 may include computer recordable medium 502 such as, but not limited to, memory, read / write (R / W) CD, R / W, DVD, and the like.

[0334] In some embodiments, signal bearing medium 500 may include communication medium 503 such as, but not limited to, digital and / or analog communication media (eg, fiber optic cables, waveguides, wired communication links, wireless communication links, etc.).

[0335] The signal bearing medium 500 may be conveyed by a wireless form of communication medium 503. The one or more program instructions may be, for example, computer executable instructions or logic implemented instructions.

[0336] In some examples, such as for Figure 5 The described modal analysis apparatus of the Coriolis mass flowmeter may be configured to provide various operations, functions, or actions in response to one or more program instructions via the computer-readable medium 501 , the computer-recordable medium 502 , and / or the communication medium 503 .

[0337] Through the description of the above implementation methods, technical personnel in the relevant field can clearly understand that for the convenience and simplicity of description, only the division of the above-mentioned functional modules is used as an example. In actual applications, the above-mentioned functions can be distributed and completed by different functional modules as needed, that is, the internal structure of the device can be divided into different functional modules to complete the full classification or partial functions described above.

[0338] In the several embodiments provided in this application, it should be understood that the disclosed devices and methods can be implemented in other ways. For example, the device embodiments described above are merely schematic. For example, the division of modules or units is only a logical function division. In actual implementation, there may be other division methods, such as multiple units or components can be combined or integrated into another device, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some interfaces, indirect coupling or communication connection of devices or units, which can be electrical, mechanical or other forms.

[0339] The units described as separate components may or may not be physically separate, and the components shown as units may be one physical unit or multiple physical units, that is, they may be in one location or distributed in multiple different locations. Depending on actual needs, some or all of the units may be selected to achieve the purpose of this embodiment.

[0340] In addition, the functional units in the various embodiments of the present application may be integrated into a single processing unit, or each unit may exist physically separately, or two or more units may be integrated into a single unit. The aforementioned integrated units may be implemented in the form of hardware or software functional units.

[0341] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a readable storage medium. Based on this understanding, the technical solution of the embodiment of the present application is essentially or the part that contributes to the prior art or the full classification part or part of the technical solution can be embodied in the form of a software product, which is stored in a storage medium and includes several instructions to enable a device (which can be a single-chip microcomputer, chip, etc.) or a processor (processor) to execute the full classification part or part of the steps of the various embodiments of the present application. The aforementioned storage medium includes various media that can store program codes, such as a USB flash drive, a mobile hard disk, a ROM, a RAM, a magnetic disk or an optical disk.

[0342] The above are only specific embodiments of the present application, but the scope of protection of this application is not limited thereto. Any changes or substitutions that can be easily conceived by a person skilled in the art within the technical scope disclosed in this application should be included in the scope of protection of this application. Therefore, the scope of protection of this application should be based on the scope of protection of the claims.

Claims

1. A modal analysis method for a Coriolis mass flowmeter, characterized in that: The method comprises: Determine the real-time relationship between fluid pressure and structural stiffness through a two-way fluid-structure coupling model; Based on the real-time relationship, dynamically correcting the flow data of the Coriolis mass flowmeter; Based on the dynamically corrected flow data and the nonlinear relationship between the straight pipe length and the phase difference, a modal analysis result of the Coriolis mass flowmeter is determined.

2. The method according to claim 1, characterized in that The parameter setting of the bidirectional fluid-structure coupling model satisfies at least one of the following: The pipe diameter and wall thickness of the measuring pipe are 15 mm and 1.5 mm, respectively. The structural parameters of the vibration damping plate are set according to the actual size; The material of the measuring tube is 316 stainless steel; The elastic modulus of the measuring tube is Gigapascal; Poisson's ratio is 0.275; Density is 7980 kg / m3; The fluid medium is methane at 300 Kelvin; The fluid domain mesh adopts a hexahedral structured network; An expansion layer is provided at the fluid-solid interface; The solid domain mesh is a hexahedral mesh.

3. The method according to claim 1, characterized in that The multi-physics coupling setting of the bidirectional fluid-structure coupling model satisfies at least one of the following: A pressure-based solver and Realizable k-ε turbulence model are used for the fluid domain; Set up the gravity field for the fluid domain; Activate the energy equation for the fluid domain; A transient dynamics analysis method is used for the solid domain; For the solid domain, boundary conditions are set, including fixed supports at the inlet and outlet, and the inner wall is set as the fluid-solid interface; For the solid domain, a loading condition is set up that includes a periodic pressure applied at the shaker locations.

4. The method according to claim 1, wherein The method of dynamically correcting the flow data of the Coriolis mass flowmeter based on the real-time relationship includes: Based on the real-time relationship, a temperature expansion compensation algorithm is used to dynamically correct the flow data of the Coriolis mass flowmeter.

5. The method according to claim 1, wherein Determining the modal analysis result of the Coriolis mass flowmeter based on the dynamically corrected flow data and the nonlinear relationship between the straight pipe length and the phase difference includes: Obtaining each order of natural frequency based on the dynamically corrected flow data and the nonlinear relationship between the straight pipe length and the phase difference; Analyze the modal vibration characteristics based on the natural frequencies of each order; The second-order modal frequency is selected as the excitation frequency for matching verification; The modal analysis result is determined according to the matching verification result of the second-order mode.

6. The method according to claim 1, characterized in that The method further comprises: Obtaining an initial flow field distribution based on the dynamically corrected flow data and a nonlinear relationship between the straight pipe length and the phase difference; determining field quantity information based on the initial flow field distribution; Calculating the stress field under the action of the fluid based on the field quantity information; Based on the stress field, obtaining an initial deformation of the structure; Based on the initial deformation of the structure, steady-state analysis results are determined.

7. The method according to claim 1, characterized in that The method further comprises: Performing transient coupling calculation based on the dynamically corrected flow data and the nonlinear relationship between the straight pipe length and the phase difference; Cubic spline interpolation is used to optimize transient coupling calculation results; Extracting the time domain characteristics of the optimized transient coupling calculation results; Based on the time domain characteristics, determining a phase difference analysis result and a frequency domain energy analysis result; The vibration data analysis result of the transient coupling calculation is determined based on the phase difference analysis result and the frequency domain energy analysis result.

8. The method according to claim 7, characterized in that The time step in the transient coupling calculation is set to between one tenth and one twentieth of the modal frequency.

9. A modal analysis device for a Coriolis mass flowmeter, characterized in that: include: processor and memory; The processor is coupled to the memory; The memory is used to store one or more programs, which include computer-executable instructions. When the modal analysis device of the Coriolis mass flowmeter is running, the processor executes the computer-executable instructions stored in the memory to enable the modal analysis device of the Coriolis mass flowmeter to perform the modal analysis method of the Coriolis mass flowmeter as described in any one of claims 1-8.

10. A computer-readable storage medium storing instructions, characterized in that: When a computer executes the instructions, the computer performs the modal analysis method for the Coriolis mass flowmeter according to any one of claims 1 to 8.

Images