CN106383930B - A kind of multiple fluid and structural simulation method of tail bearing-rotor-support-foundation system - Google Patents
A kind of multiple fluid and structural simulation method of tail bearing-rotor-support-foundation system Download PDFInfo
- Publication number
- CN106383930B CN106383930B CN201610785265.6A CN201610785265A CN106383930B CN 106383930 B CN106383930 B CN 106383930B CN 201610785265 A CN201610785265 A CN 201610785265A CN 106383930 B CN106383930 B CN 106383930B
- Authority
- CN
- China
- Prior art keywords
- fluid
- domain
- discretization
- rotor
- propeller
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
A kind of multiple fluid and structural simulation method of tail bearing-rotor-support-foundation system, belong to ocean platform fluid structurecoupling technical field, including multiple tail bearings, a propeller shaft, total time setting, fluid domain cell discretization, solid domain cell discretization, the calculating of discretization fluid domain, the calculating of discretization solid domain, the data coupling for being arranged including unit time step-length and calculating calculate, calculate end step；The multiple fluid and structural simulation method of the tail bearing-rotor-support-foundation system realizes the accurate simulation of tail bearing Yu the rotor-support-foundation system course of work, solves existing Dynamic mesh, the area of space occupied applied to multiple tail bearings, it is excessive to carry out mesh distortion present in numerical simulation, so as to cause the low problem of fluid domain computational accuracy；The multiple coupling calculating for realizing fluid-solid overcomes the drawbacks of symmetrical Rotor-Bearing System coupling calculates, being capable of the multiple coupling effect of accurate simulation propeller shaft and multiple tail bearings between vibration, lubrication.
Description
Technical field
A kind of multiple fluid and structural simulation method of tail bearing-rotor-support-foundation system, belongs to ocean platform fluid structurecoupling technical field.
Background technique
As petroleum investigation, investigation and prospecting center of gravity gradually turn to deep-sea from land, shallow sea, the application of ocean platform is constantly expanded
Greatly.Ocean platform is divided into fixed, movable and semi-fixed type by structure feature and working condition.Propeller is that ocean is floating type
One of critical component of platform, ocean floating platform support its platform weight by the buoyancy of itself, and by propeller
Dynamic positioning platform, propeller shaft are the tailing axle of the argosy system with complicated shafting, support fortune by the tail bearing of ship tail portion
Row runs the effects of reaction force generated resists wind, wave, stream power, reaches balance and position purpose, it is ensured that ocean floating platform
It can long period safety and stability operation.Therefore, the tail bearing-rotor-support-foundation system fluid and structural simulation method of research asymmetric drift off field
It is of great significance to ocean platform long period trouble free service；Currently, the tail bearing that China's propeller shafting uses is largely
By introduction of foreign technology or external product is copied to be produced, autonomous Design power is not strong, multi- scenarios method operation process simulation
Weak foundation and there are various problems, can not develop the boat with the deep sea drilling floating platform of autonomous property right for China's depth
It is capable to provide technical support with positioning power system.
In realizing process of the present invention, inventor is had found in the prior art the prior art has at least the following problems: 1, tail bearing and spiral shell
Revolve the Working Process inaccuracy of paddle shaft system；2, tail bearing-rotor-support-foundation system carries out the fluid during fluid structurecoupling
Domain computational accuracy is low；3, tail bearing-rotor-support-foundation system carries out accurately disclosing propeller shaft and multiple tails during fluid structurecoupling
Coupling mechanism of the bearing between vibration, lubrication.
Found through analysis, occur the main reason for above problem be due to:
1, the CAE software of existing fluid structurecoupling technical application large size, to the discretization data of fluid domain and solid domain into
Row unified calculation processing leads to tail bearing and the Working Process of propeller shaft system inaccuracy；
2, existing fluid domain calculate be applied to tail bearing numerical simulation using Dynamic mesh there are mesh distortions it is excessive
Problem；
3, in propeller shafting, propeller shaft is supported jointly by multiple tail bearings, is had between each tail bearing very strong
Coupled relation, the variation of working condition influences whether the working condition of adjacent tail bearing between rotor-support-foundation system and a tail bearing,
The lubrication working condition interaction of each tail bearing influences each other；Therefore, the coupling calculating process of tail bearing and propeller shaft will
It different from the calculating of symmetrical Rotor-Bearing System, needs to consider the asymmetry between each tail bearing, realizes fluid-fluid-
The multiple coupling of solid calculates.
Summary of the invention
The technical problem to be solved by the present invention is overcome the deficiencies of the prior art and provide a kind of accurate simulation tail bearing with
The rotor-support-foundation system course of work improves fluid domain computational accuracy, accurate simulation propeller shaft and multiple tail bearing coupling mechanism
The multiple fluid and structural simulation method of tail bearing-rotor-support-foundation system.
The technical solution adopted by the present invention to solve the technical problems is: the multiple fluid structurecoupling of tail bearing-rotor-support-foundation system
Calculation method, including multiple tail bearings, a propeller shaft, further comprising the steps of:
The setting of unit time step-length and the total time setting calculated；
Fluid domain cell discretization, the cell discretization for the area of space that multiple tail bearings occupy；
Solid domain cell discretization, the cell discretization of propeller shaft rotor dynamics equation；
The calculating of discretization fluid domain updates discretization fluid domain grid node using structuring dynamic mesh update method,
And fluid matasomatism is obtained in the borderline Nonlinear Oil-Film Forces of Sliding of propeller shaft；
The calculating of discretization solid domain solves propeller shaft rotor dynamics equation using time-domain integration method, obtains spiral
Displacement and center at paddle shaft neck；
Data coupling calculates, in unit time step-length, using multiple Nonlinear Oil-Film Forces of Sliding as discretization solid domain
Boundary condition is input to discretization solid domain and is calculated；In same time step, by displacement and centre bit at propeller axle journal
It sets, as the grid updating coordinate of discretization fluid domain, is input to multiple discretization fluid domains and is calculated；
Calculating terminates, and the total time of calculating arrives.
Preferably, the discretization fluid domain calculating the following steps are included:
S201 establishes multiple Nonlinear Oil-Film Forces of Sliding equations；
S202 calls rigid boundary condition's database file, determines that current time walks the grid section at lower propeller axle journal
Point coordinate；
S203 updates discretization fluid domain using structuring dynamic mesh update method in computation fluid dynamics software
Grid node, obtain multiple Nonlinear Oil-Film Forces of Sliding under newest grid node；
Multiple Nonlinear Film force boundary condition database files are written in multiple Nonlinear Oil-Film Forces of Sliding by S204；
S205, rigid boundary condition's database file exist, return step S202；Otherwise, the calculating of discretization fluid domain
Terminate.
Preferably, the discretization solid domain calculating the following steps are included:
S301 establishes propeller shaft rotor dynamics equation；
S302 calls multiple Nonlinear Film force boundary condition database files, determines integral domain；
S303 couples multiple Nonlinear Oil-Film Forces of Sliding with propeller shaft rotor dynamics equation；
S304 solves propeller shaft rotor dynamics equation using time-domain integration method, obtains at newest propeller axle journal
Displacement and center；
Rigid boundary condition's database file is written in displacement at the propeller axle journal and center by S305；
S306, multiple Nonlinear Film force boundary condition database files exist, return step S302；Otherwise, discretization
The calculating of solid domain terminates.
Preferably, the structuring dynamic mesh update method the following steps are included:
S401 calls rigid boundary condition's database file, obtains displacement and center at propeller axle journal；
S402, mesh point coordinate identification；
S403 determines grid node region according to bearing shell line of demarcation equation；
S404, grid node moving distance assignment；
S405 determines newest mesh point coordinate.
Preferably, the method for mesh point coordinate moving distance assignment described in step S404 are as follows: with propeller axle journal weight
The grid node and propeller axle journal synchronizing moving assignment of conjunction, the grid node being overlapped with bearing shell keep assignment constant, propeller
Grid node between axle journal and bearing shell moves assignment as follows:
Wherein, (xi, yi) it is the coordinate that propeller axle journal moves front nodal point i, (x'i, y'i) it is after propeller axle journal is mobile
The coordinate of node i, NiIndicate the radial grid number of plies where node i, reticulate layer number node from bearing shell is saved at axle journal
O'clock incremented by successively from 0 to N, N indicates total radial grid number of plies, and N≤10, q indicate the growth ratio of grid height, 0.96≤q≤
1.06, and q ≠ 1.
Preferably, the expression formula of the propeller shaft rotor dynamics equation are as follows:
Wherein, M is the mass matrix of propeller shaft rotor-support-foundation system, and s is the displacement of propeller shaft rotor-support-foundation system, and C is spiral shell
The damping matrix of paddle shaft rotor-support-foundation system is revolved, K is the stiffness matrix of propeller shaft rotor-support-foundation system, and Q is that the quality of propeller shaft is uneven
It measures, F2For propeller shaft gravity, F3For Nonlinear Oil-Film Forces of Sliding；Kr0ejωtFor the exciting force as caused by initial bending, r0It indicates just
Begin to be bent, ω is propeller shaft angular velocity of rotation；BiFor the location matrix of active force, i≤4.
Preferably, bearing shell line of demarcation equation described in step S403 are as follows:
When tail bearing is elliptic bearing, the general expression of line of demarcation equation are as follows:
Wherein: R is the radius of tail bearing, and r is the radius of propeller axle journal, (a0, b0) it is that propeller axle journal center is sat
Mark；
When tail bearing is that four oily rachises are held, general expression are as follows:
Wherein: R is the radius of tail bearing, and r is the radius of propeller axle journal, (a0, b0) it is that propeller axle journal center is sat
Mark,
Preferably, the time-domain integration method is any one of Newmark method, Runge-Kutta method.
Preferably, the unit time step-length is according to formula: unit time step-length=L × minimum grid length/fluid is special
Sign speed is set, 0.1≤L≤1.
Preferably, the total time of the calculating is according to formula: total time=calculating total step number × unit time step is long to carry out
Setting.
The technical solution adopted by the present invention to solve the technical problems is described as follows:
By fluid domain cell discretization and solid domain cell discretization step, fluid domain and solid domain are carried out respectively single
First discretization overcomes the drawbacks of traditional mainframe computer software unifies sliding-model control, simplifies subsequent calculating step, after saving
Continuous computing resource.
It is coupled by the calculating of discretization fluid domain, the calculating of discretization solid domain and data and calculates step, in unit
In spacer step, the calculating process that Nonlinear Film force boundary condition is discretized solid domain is called, and determines its time-domain integration region, and
Using in Newmark method or Runge-Kutta method any time-domain integration method calculate output propeller axle journal displacement and
Center, propeller axle journal displacement and center are written into rigid boundary condition's database file；In same time step
Interior, the calculating process which is discretized fluid domain is called, and more using structuring dynamic mesh
New method updates grid node to obtain multiple Nonlinear Oil-Film Forces of Sliding under new grid node, and multiple Nonlinear Oil-Film Forces of Sliding is same
When be written Nonlinear Film force boundary condition database file, wait the calculating process of discretization solid domain in next time step
It calls, so on circulate, until what is calculated arrives total time, calculating terminates；Realize discretization fluid domain in same time step
With the multiple fluid and structural simulation of discretization solid domain, the accurate simulation course of work of tail bearing and propeller shaft system.
Structuring dynamic mesh update method accurately determines grid node region using bearing shell line of demarcation equation, passes through
The assignment method for following propeller axle journal mobile using mesh point coordinate, i.e., the grid node and spiral shell being overlapped with propeller axle journal
Paddle shaft neck synchronizing moving assignment is revolved, the grid node being overlapped with bearing shell keeps assignment constant, between propeller axle journal and bearing shell
Grid node moves assignment as follows:
Overcome caused by mesh distortion is excessive in traditional dynamic mesh update method that fluid domain computational accuracy is too low to ask
Topic.
The propeller shaft rotor dynamics equation that discretization solid domain is established has fully considered that the quality of propeller shaft is uneven
It measures, propeller shaft gravity, the various mechanism of action of Nonlinear Oil-Film Forces of Sliding and the exciting force as caused by initial bending complete more
The coupling of a Nonlinear Oil-Film Forces of Sliding and propeller shaft rotor dynamics equation overcomes symmetrical Rotor-Bearing System coupling
The drawbacks of calculating, the coupling mechanism of more acurrate simulation and announcement propeller shaft and tail bearing between vibration, lubrication.
Compared with prior art, the present invention has the beneficial effects that
1, the multiple fluid and structural simulation method of the tail bearing-rotor-support-foundation system carries out unit for fluid domain and solid domain respectively
Discretization, and be respectively completed by the calculating step of discretization fluid domain and the calculating step of discretization solid domain for fluid domain
With the calculating of solid domain；The feature modeling for calculating step coupling discretization fluid domain and solid domain is coupled by data as a result, complete
At the multiple fluid and structural simulation of tail bearing-rotor-support-foundation system, the accurate of tail bearing and the propeller shaft system course of work is realized
Simulation；
2, using the structuring dynamic mesh update method for following propeller axle journal moving method based on mesh point coordinate, solution
Existing Dynamic mesh of having determined is excessive applied to mesh distortion existing for tail bearing numerical simulation, calculates essence so as to cause fluid domain
Spend low problem；
3, the multiple fluid and structural simulation method of the tail bearing-rotor-support-foundation system fully considers the asymmetry between multiple tail bearings
Property, it realizes that fluid-fluid-solid multiple coupling calculates, overcomes the drawbacks of symmetrical Rotor-Bearing System coupling calculates, more
The coupling mechanism of accurate simulation and announcement propeller shaft and tail bearing between vibration, lubrication, develops for China's depth and has
The navigation and positioning power system for having the deep sea drilling floating platform of autonomous property right provide important technical support.
Detailed description of the invention
The multiple fluid and structural simulation method and step flow diagram of Fig. 1 tail bearing-rotor-support-foundation system.
The calculating steps flow chart block diagram of Fig. 2 discretization fluid domain.
The calculating steps flow chart block diagram of Fig. 3 discretization solid domain.
Fig. 4 structuring dynamic mesh update method steps flow chart block diagram.
Specific embodiment
1~4 pair of multiple fluid and structural simulation method of tail bearing-rotor-support-foundation system of the present invention is done furtherly with reference to the accompanying drawing
It is bright.
Fig. 1 be the multiple fluid and structural simulation method and step flow diagram of tail bearing-rotor-support-foundation system, including multiple tail bearings,
One propeller shaft, further comprising the steps of:
Step S101, unit time step-length installation warrants formula: unit time step-length=L × minimum grid length/fluid
Characteristic velocity is set, wherein 0.1≤L≤1；The total time installation warrants formula of calculating: total time=calculating total step number ×
Unit time step length is set, and is calculated step number and is preset；
Step S102, fluid domain cell discretization；That is, the cell discretization for the area of space that multiple tail bearings occupy；It adopts
It is discrete that fluid domain is carried out with Gambit software, and grid cell size is set, using Cooper gridding method, by tail bearing oil film
Gap is divided into structured grid；
Step S103, solid domain cell discretization；That is, the cell discretization of propeller shaft rotor dynamics equation；Using
The cubic polynomial to overlap each other is as Rayleigh function, and using propeller shaft as continuous beam processing, whole axis is not only displaced and speed
It spends continuously, and moment of flexure and shearing are also kept continuously, the corner of shaft is then indicated by the derivative of shaft displacement, is being examined in this way
Under the premise of considering the influence of system gyroscope torque, number of degrees of freedom, reduces half, propeller than conventional Finite Element Discretization Technique
The equation of motion of axis be it is linear, other non-linear factors are described by acting on the generalized force on propeller shaft, and are used
Nonlinear fault power is coupled in system differential equation by Modal Synthesis Technique and position function；
Step S104, the calculating of discretization fluid domain；Discretization fluid domain is updated using structuring dynamic mesh update method
Grid node, and fluid matasomatism is obtained in the borderline Nonlinear Oil-Film Forces of Sliding of propeller shaft；
Step S105, the calculating of discretization solid domain；Propeller shaft rotor dynamics equation is solved using time-domain integration method,
Obtain displacement and center at propeller axle journal；Newmark method or Runge-Kutta method can be used in time-domain integration method；
Step S106, data coupling calculate；In unit time step-length, using multiple Nonlinear Oil-Film Forces of Sliding as discretization
The boundary condition of solid domain is input to discretization solid domain and is calculated；In same time step, it will be displaced at propeller axle journal
And center is input to multiple discretization fluid domains and is calculated as the grid updating coordinate of discretization fluid domain；
The total time of step S107, calculating arrive, and calculating terminates, and otherwise continue discretization solid domain or discretization fluid domain
It calculates.
Fig. 2 is the calculating steps flow chart block diagram of discretization fluid domain, and the calculating of the discretization fluid domain includes following step
It is rapid:
Step S201 establishes multiple Nonlinear Oil-Film Forces of Sliding equations, equation expression formula are as follows:
Wherein: R is tail bearing radius, and L is tail bearing length, and θ is tail bearing angle of circumference, and p is oil film pressure；Pass through calculating
Hydrodynamics software, setting unit time step, after calculating total step number, boundary condition, the parameters such as propeller shaft revolving speed, calculating can
Obtain multiple Nonlinear Oil-Film Forces of Sliding；
Step S202 calls rigid boundary condition's database file, determines that current time walks grid at lower propeller axle journal
Node coordinate, rigid boundary condition's database file extend entitled .dat；
Step S203 updates the grid node of discretization fluid domain using structuring dynamic mesh update method, obtains newest
Multiple Nonlinear Oil-Film Forces of Sliding under grid node；
Multiple Nonlinear Film force boundary condition database files are written in multiple Nonlinear Oil-Film Forces of Sliding by step S204；It is more
A Nonlinear Film force boundary condition database file extends entitled .dat；
Step S205, rigid boundary condition's database file exist, return step S202；Otherwise, discretization fluid domain
Calculating terminates.
Nonlinear Film force boundary condition database file and rigid boundary condition's database file are being calculated accordingly
It is automatically deleted after the invocation of procedure.
Fig. 3 is the calculating steps flow chart block diagram of discretization solid domain, and the calculating of the discretization solid domain includes following step
It is rapid:
Step S301 establishes propeller shaft rotor dynamics equation:
Wherein, M is the mass matrix of propeller shaft rotor-support-foundation system, and s is the displacement of propeller shaft rotor-support-foundation system, and C is spiral shell
The damping matrix of paddle shaft rotor-support-foundation system is revolved, K is the stiffness matrix of propeller shaft rotor-support-foundation system, and Q is that the quality of propeller shaft is uneven
It measures, F2For propeller shaft gravity, F3For Nonlinear Oil-Film Forces of Sliding；Kr0ejωtFor the exciting force as caused by initial bending, r0It indicates just
Begin to be bent, ω is propeller shaft angular velocity of rotation；BiFor the location matrix of active force, i≤4；
Step S302 calls multiple Nonlinear Film force boundary condition database files, determines time-domain integration region；
Step S303 couples multiple Nonlinear Oil-Film Forces of Sliding with propeller shaft rotor dynamics equation, coupled wave equation are as follows:
Wherein: F is generalized force, herein refers to Nonlinear Oil-Film Forces of Sliding, this equation can be used for other non-linear exciting-vibration forces
Coupling,For modal transfer matrix, Respectively represent propeller
1 rank of axis rotor dynamics equation, 2 ranks ..., n rank characteristic vector, JTAnd JRThe respectively translation and rotation of Nonlinear Oil-Film Forces of Sliding
Location matrix,W is the axial position function of propeller shaft, fi eAnd li eRespectively act on
The concentrated force and concentrated moment of i-th of Nonlinear Oil-Film Forces of Sliding on propeller shaft, active position areThe direction z is parallel
In propeller shaft axially direction.
Step S304 solves propeller shaft rotor dynamics equation using time-domain integration method, obtains newest propeller shaft
Displacement and center at neck；
Rigid boundary condition's database file is written in displacement at the propeller axle journal and center by step S305；
Step S306, multiple Nonlinear Film force boundary condition database files exist, return step S302；Otherwise, from
The calculating of dispersion solid domain terminates.
Nonlinear Film force boundary condition database file and rigid boundary condition's database file are being calculated accordingly
It is automatically deleted after the invocation of procedure.
Fig. 4 be structuring dynamic mesh update method steps flow chart block diagram, the structuring dynamic mesh update method include with
Lower step:
Step S401 calls rigid boundary condition's database file, obtains displacement and center at propeller axle journal；
Step S402 identifies mesh point coordinate according to displacement at propeller axle journal and center；
Step S403 determines grid node region according to bearing shell line of demarcation equation；The bearing shell line of demarcation of different structure bearing
Equation is also different, when tail bearing is elliptic bearing, the general expression of line of demarcation equation are as follows:
Wherein: R is the radius of tail bearing, and r is the radius of propeller axle journal, (a0, b0) it is journal centre's position coordinates；
When tail bearing is that four oily rachises are held, general expression are as follows:
Wherein: R is the radius of tail bearing, and r is the radius of propeller axle journal, (a0, b0) it is journal centre's position coordinates,
Step S404, grid node moving distance assignment；The grid node and propeller axle journal being overlapped with propeller axle journal
Synchronizing moving assignment, the grid node being overlapped with bearing shell keep assignment constant, the grid node between propeller axle journal and bearing shell
Assignment is moved as follows:
Wherein, (xi, yi) it is the coordinate that axle journal moves front nodal point i, (x'i, y'i) it is the coordinate that axle journal moves posterior nodal point i, Ni
Indicate the radial grid number of plies where node i, reticulate layer number node node at axle journal from bearing shell is successively passed from 0 to N
Increasing, N indicates the total radial grid number of plies, and N is positive integer, and N≤10, q indicate the growth ratio of grid height, 0.96≤q≤1.06, and
q≠1；
Step S405 is determined and is updated mesh point coordinate, i.e., newest mesh point coordinate determines.
The initial fluid domain grid model of Examples 1 to 2 is divided using structured grid, and radial grid gap size is consistent,
The grid node center of circle having the same in identical tile fragment region, propeller shaft system have multiple thin circles as continuous beam processing
Disk rigid body composition, propeller shaft entirety are support by two tail bearings, and propeller shaft quality is divided in the width range of itself
The corner vibration of cloth quality treatment, propeller shaft is indicated by the displacement derivatives of propeller shaft.
1: two four multiple fluid and structural simulation process of oily leaf tail Rotor-Bearing System of embodiment, it is best real for the present invention
Apply example.
1.1, setting unit time step Δ t1, the total time T of calculating1；
1.2, numerical value calculating is carried out using computation fluid dynamics software, while to two four oily leaf tail bearings, obtains two
Nonlinear Oil-Film Forces of Sliding (the Fx of a four oily leaf tail bearings1, Fy1)、(Fx2, Fy2)；
1.3, by Δ t1, current calculate time Td1、(Fx1, Fy1)、(Fx2, Fy2) as the non-of discretization solid domain calculating
Two four oily leaf tail bearing boundary condition database files are written in linear oil-film force boundary condition；
1.4, propeller shaft rotor dynamics equation is established in the calculating of discretization solid domain:
Wherein, M is the mass matrix of propeller shaft rotor-support-foundation system, and s is the displacement of propeller shaft rotor-support-foundation system, and C is spiral shell
The damping matrix of paddle shaft rotor-support-foundation system is revolved, K is the stiffness matrix of propeller shaft rotor-support-foundation system, and Q is that the quality of propeller shaft is uneven
It measures, F2For propeller shaft gravity, F3For Nonlinear Oil-Film Forces of Sliding, Kr0ejωtFor propeller shaft in exciting, touch and rub, misalign situation
Under, the exciting force as caused by initial bending, r0Indicate that initial bending, ω are propeller shaft angular velocity of rotation；BiFor above-mentioned active force
Location matrix, i≤4, i be rounded numerical value；
1.5, discretization solid domain, which calculates, calls above-mentioned two four oily leaf tail bearing boundary condition database files, determines
Newmark method integral domain；
1.6, pass through formula:
By (Fx1, Fy1)、(Fx2, Fy2) coupled with propeller shaft rotor dynamics equation；Wherein, XTFor modal transformation square
Battle array,
Respectively represent propeller shaft rotor dynamics equation 1 rank, 2
Rank ..., n rank characteristic vector, JTAnd JRThe respectively translation of Nonlinear Oil-Film Forces of Sliding and turned position matrix,
W is the axial position function of propeller shaft, fi eAnd li eI-th respectively acted on propeller shaft is non-linear
The concentrated force and concentrated moment of oil-film force；
1.7, above-mentioned propeller shaft rotor dynamics equation is solved using Newmark integration method；
1.8, by the propeller axle journal being calculated displacement (the Δ x at two four oily leaf tail bearings respectively1、Δy1)、
(Δx2、Δy2) and center (x1、y1)、(x2、y2) deposit rigid boundary condition's database file；
1.9, the total time T of calculating1It does not arrive and Nonlinear Film force boundary condition database file exists, then return step
1.5；Otherwise, if the total time T calculated1It does not arrive, enters step 1.10, the total time T of calculating1Terminate to then calculating；
1.10, using structuring dynamic mesh update method, rigid boundary condition's database file is called, current spiral is read
Displacement (Δ x of the paddle shaft neck at two four oily leaf tail bearings1、Δy1)、(Δx2、Δy2) and center (x1、y1)、(x2、
y2)；Identify that current time walks Δ t by calculating fluid software CFD1Interior mesh point coordinate is (m1, n1), two four oily leaves
Region occupied by the oily leaf tail bearing of single the four of tail bearing is divided into four parts by four tile fragment lines of demarcation, and line of demarcation with
The movement of propeller axle journal and change, the equation in four lines of demarcation is respectively as follows:
The bearing shell line of demarcation equation of another four oily leaf tail bearing is
(m is judged according to above-mentioned bearing shell line of demarcation equation1, n1) where flow field domain, it may be assumed that
Work as n1≤b1And n1< b4When, mesh point coordinate (m1, n1) it is in a1b1And a4b4The fluid domain that line of demarcation determines
It is interior；
Work as n1≤b2And n1> b1When, mesh point coordinate (m1, n1) it is in a1b1And a2b2The fluid domain that line of demarcation determines
It is interior；
Work as n1> b2And n1≥b3When, mesh point coordinate (m1, n1) it is in a2b2And a3b3The fluid domain that line of demarcation determines
It is interior；
Work as n1< b3And n1≥b4When, mesh point coordinate (m1, n1) it is in a3b3And a4b4The fluid domain that line of demarcation determines
It is interior；
Work as n1≤b5And n1< b8When, mesh point coordinate (m1, n1) it is in a5b5And a8b8The fluid domain that line of demarcation determines
It is interior；
Work as n1≤b6And n1> b5When, mesh point coordinate (m1, n1) it is in a5b5And a6b6The fluid domain that line of demarcation determines
It is interior；
Work as n1> b6And n1≥b7When, mesh point coordinate (m1, n1) it is in a6b6And a7b7The fluid domain that line of demarcation determines
It is interior；
Work as n1< b7And n1≥b8When, mesh point coordinate (m1, n1) it is in a7b7And a8b8The fluid domain that line of demarcation determines
It is interior；
1.11 for the grid node moving distance in each flow field domain assignment, the grid node being overlapped with propeller axle journal
With propeller axle journal synchronizing moving assignment, the grid node that is overlapped with bearing shell keeps assignment constant, propeller axle journal and bearing shell it
Between grid node move assignment as follows:
I.e. for time step Δ t1Interior mesh point coordinate (m1, n1), when it is overlapped with propeller axle journal, net at this time
Lattice node coordinate is (m1+ Δ x, n1+Δy)；When it is overlapped with bearing shell, mesh point coordinate is (m at this time1, n1)；When it is in spiral shell
Revolve (N when between paddle shaft neck and bearing shelli=2, N=7, q=0.96), mesh point coordinate is at this time
Wherein, as mesh point coordinate (m1, n1) determine in a1b1~a4b4When flow field domain where bearing shell line of demarcation, Δ x
=Δ x1, Δ y=Δ y1；Mesh point coordinate (m1, n1) determine in a5b5~a8b8Δ x when flow field domain where bearing shell line of demarcation
=Δ x2, Δ y=Δ y2；
1.12, grid node, which updates, to be completed, and obtains newest grid node (m'1, n'1) under two four oily leaf tail bearings
Nonlinear Oil-Film Forces of Sliding (Fx3, Fy3)、(Fx4, Fy4)；
1.13, the total time T of calculating1It does not arrive and rigid boundary condition's database file exists, return step 1.4；Otherwise,
If the total time T calculated1It does not arrive, enters step 1.3, the total time T of calculating1Terminate to then calculating.
The multiple fluid and structural simulation process of the oval tail bearing-rotor-support-foundation system of embodiment 2: three.
2.1, setting unit time step Δ t2, the total time T of calculating2；
2.2, numerical value calculating is carried out using computation fluid dynamics software, while to three oval tail bearings, obtains three
Nonlinear Oil-Film Forces of Sliding (the Fx of oval tail bearing5, Fy5)、(Fx6, Fy6)、(Fx7, Fy7)；
2.3, by Δ t2, current calculate time Td2、(Fx5, Fy5)、(Fx6, Fy6)、(Fx7, Fy7) it is used as discretization solid domain
Three oval tail bearing boundary condition database files are written in the Nonlinear Film force boundary condition of calculating；
2.4, propeller shaft rotor dynamics equation is established in the calculating of discretization solid domain:
2.5, discretization solid domain, which calculates, calls above three ellipse tail bearing boundary condition database file, determines
Runge-Kutta method integral domain；
2.6, Nonlinear Oil-Film Forces of Sliding couples formula for (Fx with propeller shaft rotor dynamics equation through the invention5, Fy5)、
(Fx6, Fy6)、(Fx7, Fy7) coupled with propeller shaft rotor dynamics equation；
2.7, propeller shaft rotor dynamics equation is solved using Runge-Kutta integration method；
2.8, by the propeller axle journal being calculated displacement (the Δ x at three oval tail bearings respectively3、Δy3)、(Δ
x4、Δy4)、(Δx5、Δy5) and center (x3、y3)、(x4、y4)、(x5、y5) deposit rigid boundary condition's database file；
2.9, the total time T of calculating2It does not arrive and Nonlinear Film force boundary condition database file exists, then return step
2.5；Otherwise, if the total time T calculated2It does not arrive, enters step 2.10, the total time T of calculating2Terminate to then calculating；
2.10, using structuring dynamic mesh update method, rigid boundary condition's database file is called, is read current
(Δx3、Δy3)、(Δx4、Δy4)、(Δx5、Δy5) and (x3、y3)、(x4、y4)、(x5、y5), by calculating fluid software CFD
Identify that current time walks Δ t2Interior mesh point coordinate is (m2, n2), shared by the single oval tail bearing of three oval tail bearings
According to region be divided into two parts by two tile fragment lines of demarcation, and line of demarcation changes, two with the movement of propeller axle journal
The equation in line of demarcation is respectively as follows:
The bearing shell line of demarcation equation of other two oval tail bearings are as follows:
(m is judged according to above-mentioned bearing shell line of demarcation equation2, n2) where flow field domain, it may be assumed that
Work as m2>=0 and n2≤b1When or m2< 0 and n2≤b2When, mesh point coordinate (m2, n2) it is in a1b1Line of demarcation is along suitable
Clockwise is to a2b2In fluid domain where line of demarcation；
Work as m2>=0 and n2> b1When or m2< 0 and n2> b2When, mesh point coordinate (m2, n2) it is in a1b1Line of demarcation is along inverse
Clockwise is to a2b2In fluid domain where line of demarcation；
Work as m2>=0 and n2≤b3When or m2< 0 and n2≤b4When, mesh point coordinate (m2, n2) it is in a3b3Line of demarcation is along suitable
Clockwise is to a4b4In fluid domain where line of demarcation；
Work as m2>=0 and n2> b3When or m2< 0 and n2> b4When, mesh point coordinate (m2, n2) it is in a3b3Line of demarcation is along inverse
Clockwise is to a4b4In fluid domain where line of demarcation；
Work as m2>=0 and n2≤b5When or m2< 0 and n2≤b6When, mesh point coordinate (m2, n2) it is in a5b5Line of demarcation is along suitable
Clockwise is to a6b6In fluid domain where line of demarcation；
Work as m2>=0 and n2> b5When or m2< 0 and n2> b6When, mesh point coordinate (m2, n2) it is in a5b5Line of demarcation is along inverse
Clockwise is to a6b6In fluid domain where line of demarcation；
2.11, grid node moving distance assignment moves assignment method, and root also in compliance with grid node of the present invention
According to propeller axle journal displacement (the Δ x at three oval tail bearings respectively3、Δy3)、(Δx4、Δy4)、(Δx5、Δy5) assign
It is worth and accurately determines current mesh node region；I.e. for time step Δ t2Interior mesh point coordinate (m2, n2), when it
When being overlapped with propeller axle journal, mesh point coordinate is (m at this time2+ Δ x, n2+Δy)；When it is overlapped with bearing shell, grid section at this time
Point coordinate is (m2, n2)；(the N when it is between propeller axle journal and bearing shelli=9, N=10, q=1.01), grid node at this time
Coordinate is
Wherein, as mesh point coordinate (m2, n2) determine in a1b1、a2b2When flow field domain where bearing shell line of demarcation, Δ x=
Δx3, Δ y=Δ y3；Mesh point coordinate (m2, n2) determine in a3b3、a4b4When flow field domain where bearing shell line of demarcation, Δ x=
Δx4, Δ y=Δ y4；Mesh point coordinate (m2, n2) determine in a5b5、a6b6When flow field domain where bearing shell line of demarcation, Δ x=
Δx5, Δ y=Δ y5；
2.12, grid node, which updates, to be completed, and obtains the Nonlinear Film of three oval tail bearings under newest grid node
Power (Fx8, Fy8)、(Fx9, Fy9)、(Fx10, Fy10)；
2.13, the total time T of calculating2It does not arrive and rigid boundary condition's database file exists, return step 2.4；Otherwise,
If the total time T calculated2It does not arrive, enters step 2.3, the total time T of calculating2Terminate to then calculating.
The above described is only a preferred embodiment of the present invention, being not that the invention has other forms of limitations, appoint
What those skilled in the art changed or be modified as possibly also with the technology contents of the disclosure above equivalent variations etc.
Imitate embodiment.But without departing from the technical solutions of the present invention, according to the technical essence of the invention to above embodiments institute
Any simple modification, equivalent variations and the remodeling made, still fall within the protection scope of technical solution of the present invention.
Claims (8)
1. a kind of multiple fluid and structural simulation method of tail bearing-rotor-support-foundation system, including multiple tail bearings, a propeller shaft,
It is characterized in that: further comprising the steps of:
The setting of unit time step-length and the total time setting calculated；
Fluid domain cell discretization, the cell discretization for the area of space that multiple tail bearings occupy；
Solid domain cell discretization, the cell discretization of propeller shaft rotor dynamics equation；
The calculating of discretization fluid domain updates discretization fluid domain grid node using structuring dynamic mesh update method, and obtains
Fluid matasomatism is obtained in the borderline Nonlinear Oil-Film Forces of Sliding of propeller shaft；
The calculating of discretization solid domain solves propeller shaft rotor dynamics equation using time-domain integration method, obtains propeller shaft
Displacement and center at neck；
Data coupling calculates, in unit time step-length, using multiple Nonlinear Oil-Film Forces of Sliding as the boundary of discretization solid domain
Condition is input to discretization solid domain and is calculated；In same time step, will at propeller axle journal displacement and center,
As the grid updating coordinate of discretization fluid domain, it is input to multiple discretization fluid domains and is calculated；
Calculating terminates, and the total time of calculating arrives；
The calculating of the discretization fluid domain the following steps are included:
S201 establishes multiple Nonlinear Oil-Film Forces of Sliding equations；
S202 calls rigid boundary condition's database file, and the grid node for determining that current time walks at lower propeller axle journal is sat
Mark；
S203 updates the net of discretization fluid domain using structuring dynamic mesh update method in computation fluid dynamics software
Lattice node obtains multiple Nonlinear Oil-Film Forces of Sliding under newest grid node；
Multiple Nonlinear Film force boundary condition database files are written in multiple Nonlinear Oil-Film Forces of Sliding by S204；
S205, rigid boundary condition's database file exist, return step S202；Otherwise, the calculating of discretization fluid domain terminates；
The calculating of the discretization solid domain the following steps are included:
S301 establishes propeller shaft rotor dynamics equation；
S302 calls multiple Nonlinear Film force boundary condition database files, determines integral domain；
S303 couples multiple Nonlinear Oil-Film Forces of Sliding with propeller shaft rotor dynamics equation；
S304 solves propeller shaft rotor dynamics equation using time-domain integration method, obtains and is displaced at newest propeller axle journal
The center and；
Rigid boundary condition's database file is written in displacement at the propeller axle journal and center by S305；
S306, multiple Nonlinear Film force boundary condition database files exist, return step S302；Otherwise, discretization solid
The calculating in domain terminates.
2. the multiple fluid and structural simulation method of tail bearing-rotor-support-foundation system according to claim 1, it is characterised in that: described
Structuring dynamic mesh update method the following steps are included:
S401 calls rigid boundary condition's database file, obtains displacement and center at propeller axle journal；
S402, mesh point coordinate identification；
S403 determines grid node region according to bearing shell line of demarcation equation；
S404, grid node moving distance assignment；
S405 determines newest mesh point coordinate.
3. the multiple fluid and structural simulation method of tail bearing-rotor-support-foundation system according to claim 2, it is characterised in that: step
The method of the assignment of mesh point coordinate moving distance described in S404 are as follows: the grid node and propeller being overlapped with propeller axle journal
Axle journal synchronizing moving assignment, the grid node being overlapped with bearing shell keep assignment constant, the grid between propeller axle journal and bearing shell
Node moves assignment as follows:
Wherein, (xi, yi) it is the coordinate that propeller axle journal moves front nodal point i, (x'i, y'i) it is the mobile posterior nodal point i of propeller axle journal
Coordinate, NiIndicate the radial grid number of plies where node i, reticulate layer number from bearing shell node at axle journal node from 0
Incremented by successively to N, N indicates total radial grid number of plies, and N≤10, q indicate the growth ratio of grid height, 0.96≤q≤1.06, and q
≠1。
4. the multiple fluid and structural simulation method of tail bearing-rotor-support-foundation system according to claim 1, it is characterised in that: described
The expression formula of propeller shaft rotor dynamics equation are as follows:
Wherein, M is the mass matrix of propeller shaft rotor-support-foundation system, and s is the displacement of propeller shaft rotor-support-foundation system, and C is propeller
The damping matrix of bearing-rotor system, K are the stiffness matrix of propeller shaft rotor-support-foundation system, and Q is the mass unbalance amount of propeller shaft,
F2For propeller shaft gravity, F3For Nonlinear Oil-Film Forces of Sliding；Kr0ejωtFor the exciting force as caused by initial bending, r0Indicate initial curved
Song, ω are propeller shaft angular velocity of rotation；BiFor the location matrix of active force, i≤4.
5. the multiple fluid and structural simulation method of tail bearing-rotor-support-foundation system according to claim 2, it is characterised in that: step
The equation of bearing shell line of demarcation described in S403 are as follows:
When tail bearing is elliptic bearing, the general expression of line of demarcation equation are as follows:
Wherein: R is the radius of tail bearing, and r is the radius of propeller axle journal, (a0, b0) it is propeller axle journal center position coordinates；
When tail bearing is that four oily rachises are held, general expression are as follows:
Wherein: R is the radius of tail bearing, and r is the radius of propeller axle journal, (a0, b0) it is propeller axle journal center position coordinates,
6. the multiple fluid and structural simulation method of tail bearing-rotor-support-foundation system according to claim 1, it is characterised in that: described
Time-domain integration method is any one of Newmark method, Runge-Kutta method.
7. the multiple fluid and structural simulation method of tail bearing-rotor-support-foundation system according to claim 1, it is characterised in that: described
Unit time step-length is according to formula: unit time step-length=L × minimum grid length/characteristic of fluid speed is set, 0.1≤
L≤1。
8. the multiple fluid and structural simulation method of tail bearing-rotor-support-foundation system according to claim 1, it is characterised in that: described
The total time of calculating is according to formula: total time=calculating total step number × unit time step length is set.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610785265.6A CN106383930B (en) | 2016-08-31 | 2016-08-31 | A kind of multiple fluid and structural simulation method of tail bearing-rotor-support-foundation system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610785265.6A CN106383930B (en) | 2016-08-31 | 2016-08-31 | A kind of multiple fluid and structural simulation method of tail bearing-rotor-support-foundation system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106383930A CN106383930A (en) | 2017-02-08 |
| CN106383930B true CN106383930B (en) | 2019-06-28 |
Family
ID=57939408
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610785265.6A Expired - Fee Related CN106383930B (en) | 2016-08-31 | 2016-08-31 | A kind of multiple fluid and structural simulation method of tail bearing-rotor-support-foundation system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106383930B (en) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107895077B (en) * | 2017-11-10 | 2020-03-17 | 西安交通大学 | Gas turbine pull rod rotor assembly parameter optimization method considering multiple disk and drum manufacturing factors |
| CN108197345B (en) * | 2017-12-06 | 2021-03-30 | 华北水利水电大学 | Long-line water delivery structure power time course analysis modeling method based on coupling theory |
| CN108205516B (en) * | 2017-12-26 | 2020-10-16 | 中国海洋大学 | Ocean structure fluid-solid coupling dynamic response method based on condition transfer |
| CN109740233A (en) * | 2018-12-27 | 2019-05-10 | 大连船舶重工集团有限公司 | Hard sail vortex-induced vibration based on fluid structurecoupling simulation responds preparation method |
| CN110309521A (en) * | 2018-12-27 | 2019-10-08 | 大连船舶重工集团有限公司 | A kind of hard sail wind vibration response calculation method based on fluid structurecoupling simulation |
| CN109977475B (en) * | 2019-02-28 | 2021-03-30 | 浙江大学 | Dynamic grid updating method for three-dimensional complex pipeline fluid-solid coupling calculation |
| CN110309529B (en) * | 2019-04-17 | 2022-12-09 | 中国舰船研究设计中心 | Rapid algorithm for calculating ship propeller flow field grid motion |
| CN111222291B (en) * | 2019-12-30 | 2023-03-10 | 太原重工股份有限公司 | Oil film fluid domain discretization method |
| CN112949112B (en) * | 2021-01-29 | 2022-05-20 | 中国石油大学(华东) | Rotor-sliding bearing system lubrication basin dynamic grid parallel computing method |
| CN116341419B (en) * | 2023-05-17 | 2023-08-01 | 中国科学院、水利部成都山地灾害与环境研究所 | Numerical value determining method and system for fluid-solid coupling and electronic equipment |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102902859A (en) * | 2012-10-09 | 2013-01-30 | 西安交通大学 | Hybrid sliding bearing fluid-solid coupling design method based on phase change model |
| CN102913464A (en) * | 2012-09-05 | 2013-02-06 | 江苏大学 | Method for predicting transient fluid-solid coupling characteristic of centrifugal pump rotor |
| CN104091022A (en) * | 2014-07-11 | 2014-10-08 | 北京工业大学 | Fluid-solid coupling simulation analysis method for static pressure main shaft system under working condition of velocity slip |
| CN105473898A (en) * | 2013-08-23 | 2016-04-06 | 美国轮轴制造公司 | Power transmitting component with torque transfer device configured with drag reduction system |
| CN105699074A (en) * | 2016-04-01 | 2016-06-22 | 山东大学 | Fluid-structure interaction dynamic property testing device for blade-rotor system |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10423730B2 (en) * | 2014-10-02 | 2019-09-24 | Siemens Industry Software Nv | Contact modeling between objects |
-
2016
- 2016-08-31 CN CN201610785265.6A patent/CN106383930B/en not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102913464A (en) * | 2012-09-05 | 2013-02-06 | 江苏大学 | Method for predicting transient fluid-solid coupling characteristic of centrifugal pump rotor |
| CN102902859A (en) * | 2012-10-09 | 2013-01-30 | 西安交通大学 | Hybrid sliding bearing fluid-solid coupling design method based on phase change model |
| CN105473898A (en) * | 2013-08-23 | 2016-04-06 | 美国轮轴制造公司 | Power transmitting component with torque transfer device configured with drag reduction system |
| CN104091022A (en) * | 2014-07-11 | 2014-10-08 | 北京工业大学 | Fluid-solid coupling simulation analysis method for static pressure main shaft system under working condition of velocity slip |
| CN105699074A (en) * | 2016-04-01 | 2016-06-22 | 山东大学 | Fluid-structure interaction dynamic property testing device for blade-rotor system |
Non-Patent Citations (1)
| Title |
|---|
| 非线性转子－轴承耦合系统润滑及稳定性分析;李强 等;《浙江大学学报（工学版）》;20121031;第46卷(第10期);1729-1736 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106383930A (en) | 2017-02-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106383930B (en) | A kind of multiple fluid and structural simulation method of tail bearing-rotor-support-foundation system | |
| Gallino et al. | Physics of Bubble‐Propelled Microrockets | |
| KR20100041510A (en) | An effective method to solve rigid body interactions in particle based fluid simulations | |
| Thai et al. | Finite-element modeling for static bending analysis of rotating two-layer FGM beams with shear connectors resting on imperfect elastic foundations | |
| Cai et al. | Influence of the mass conservation cavitation boundary on the tribo-dynamic responses of the micro-groove water-lubricated bearing | |
| Zan et al. | Numerical simulations of dynamic pipeline-vessel response on a deepwater S-laying vessel | |
| CN101294811A (en) | Strapdown inertial navigation system adopting strange perturbation method for taper cone error and rowing error compensation | |
| CN109614707B (en) | Dynamic modeling method based on stepped shaft-flexible disk coupling system | |
| Sun et al. | Dynamics modeling and hydrodynamic coefficients identification of the wave glider | |
| CN116702648B (en) | Vortex-induced vibration calculation method and device for three-dimensional drill string | |
| Wang et al. | Dynamic behavior of the deepwater flexible pipeline during pipe laying process | |
| Xu et al. | The stability analysis of tension-leg platforms under marine environmental loads via altering the connection angle of tension legs | |
| Wang et al. | Transient flow around an impulsively started cylinder using a dynamic mesh method | |
| WO2020174536A1 (en) | Numerical value calculation method and numerical value calculation program | |
| Qiu et al. | Computation of wave-body interactions using the panel-free method and exact geometry | |
| Felici et al. | Analysis of vortex induced vibration of a thermowell by high fidelity FSI numerical analysis based on RBF structural modes embedding | |
| Liu et al. | Mesh-free GPU simulation of complex flows in three dimensional distorting domain | |
| Yang et al. | On some azimuthal equatorial flows | |
| JP6167554B2 (en) | Channel shape optimization method and channel shape optimization apparatus | |
| Vu Nguyen et al. | A Moving Element Method for Hydroelastic Response of a Floating Thin Plate Due to a Moving Load | |
| Chen et al. | A Vortex-Lattice Modeling Approach for Free-Surface Effects on Submerged Bodies and Propellers | |
| CN110879553B (en) | Control method and system of micro gyroscope available based on output state | |
| Yang | High fidelity fluid-structure turbulence modeling using an immersed-body method | |
| Ma et al. | Simulation of droplet impacting on elastic solid with the SPH method | |
| Lighthill | What Is Mechanics? |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190628 Termination date: 20210831 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |