CN110309521A - A kind of hard sail wind vibration response calculation method based on fluid structurecoupling simulation - Google Patents
A kind of hard sail wind vibration response calculation method based on fluid structurecoupling simulation Download PDFInfo
- Publication number
- CN110309521A CN110309521A CN201811607705.4A CN201811607705A CN110309521A CN 110309521 A CN110309521 A CN 110309521A CN 201811607705 A CN201811607705 A CN 201811607705A CN 110309521 A CN110309521 A CN 110309521A
- Authority
- CN
- China
- Prior art keywords
- sail
- model
- wind
- fluid
- domain
- 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.)
- Withdrawn
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Automation & Control Theory (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
The invention discloses a kind of hard sail wind vibration response calculation methods based on fluid structurecoupling simulation, comprising the following steps: 1, establish using gambit the fluid domain grid model of ship sail；2, suitable turbulence model is selected in CFX, acts on fluid domain grid model and surface pressing file suffered by sail is obtained by numerical simulation calculation；3, the solid domain structural model of ship sail is established；4, the surface pressing numeric file that CFX is obtained is imported by interpolation method and is acted on sail solid domain structural model；5, it completes to couple and carry out numerical simulation sail wind shake situation, calculates the stress value that sail is subject to；6, judge whether sail structural strength meets code requirement under each operating condition according to calculated result, and then instruct the design and construction of ship sail.The present invention can relatively accurately calculate the stress value that sail structure is subject to, and then compared with the allowable stress of structural material, and judgement meets specification for structure, therefore the design and construction cost of ship sail can be effectively reduced.
Description
Technical field
The present invention relates to a kind of Computer Numerical Simulation and the fields CAE, are based on fluid structurecoupling more specifically to one kind
The hard sail wind vibration response calculation method of simulation.
Background technique
Ship will cause a large amount of water, exhaust emission during construction and operation.And sail-assistant ship utilizes ocean
Upper natural wind abundant provides power for sail, meets national energy conservation and emission reduction policy.Modern sail is hard sail, and rigidity rises
Lifan, since the streamline shape of sail two sides is different, when the pressure difference that natural airflow is generated through sail provides lift for sail, work
For the power for helping ship to advance.While wind field provides power for sail, the dynamic load generated will also be produced structure
It is raw to influence.Ship's navigation can generate rolling when meeting with wave, and rolling will drive towering sail structure and move together, and work as
When rolling period is close with the intrinsic frequency of sail, the movement of sail can be made to aggravate, so as to cause the increase of structural stress value.Cause
This needs in the design process to analyze this unsafe condition.
Summary of the invention
The dangerous situation that the occurs when present invention meets with wave for ship's navigation provides a kind of based on fluid structurecoupling simulation
Hard sail wind vibration response calculation method, specific steps of the invention are as follows:
S1, hull coordinate system and sail coordinate system are established in gambit, ship is established according to the actual parameter of sail
The fluid domain grid model of hard sail, and the fluid domain grid model is imported in CFX；
The initial value wind speed and boundary condition of ship hard sail are established in S2, CFX, select turbulence model and are acted on described
On fluid domain grid model, surface pressing suffered by ship hard sail is calculated using finite volume method, and obtain surface
Pressure value file；The turbulence model is SST-SAS turbulence model, the equation group of the turbulence model are as follows:
In formula, ρ is density, and t is the time, and k is Turbulent Kinetic, PkItem, μ are generated for kinetic energytIt is turbulent flow stickiness, ω is specific consumption
The rate of dissipating, μiIt is time averaged velocity, F1For mixed function, σkAnd σεIt is k, turbulent flow Pood's number of ε, β*It is constant, S with βijAveragely to answer
Variability；
The SST-SAS turbulence model equation is in view of turbulent viscosity dissipation in Disengagement zone, in the ratio dissipative shock wave of standard SST
Q is added in ω equationSAS, wherein SST turbulence model is the deformation of standard k-w model；It is using mixed function that standard k-e is rapid
Flow model combines with k-w turbulence model, and k-w is used to solve in boundary layer, and k-e outside boundary layer for solving, and QSAS
By stationary zones, i.e. QSAS=0, with unstable state region, i.e. QSAS> 0, it distinguishes；QSASL invkFor von Karman ruler
Degree, L are modelling rapids stress scale, LvkIn unstable state area according to currently known flow field analysis whirlpool dynamic, turbulent flow mould is adjusted in real time
Length dimension in type；Model constants are as follows:
C=2.0, σΦ=2/3,κ=0.41, Cμ=0.09,
S3, sail sail face is established in ANSYS software, mast, the solid domain structural model on pedestal and neighbouring deck is simultaneously
Edge-restraint condition is added to the solid domain structural model, is constrained to the freedom degree constraint of solid domain structural model fringe node
And torsional restraint；
S4, the surface pressing numeric file is imported by interpolation method and acts on the solid domain of ANSYS ship sail
On structural model, at the fluid structurecoupling interface of fluid domain grid model and solid domain structural model, fluid should be met and consolidated
Body stress τ, displacement d, heat flow q, the variables such as temperature T are equal or conservation, that is, meet following equation:
In formula, subscript s and f respectively represent fluid domain and solid domain；
S5, numerical simulation calculation hard sail wind shake situation is carried out using Finite Element, calculates what hard sail was subject to
Stress value；
S6, judge that hard sail structural strength is under each operating condition according to the stress value result that the hard sail of calculating is subject to
It is no to meet code requirement, and then instruct the design and construction of ship hard sail.
Under preferred embodiment, the wind type that the fluid domain grid model established in the S1 is selected is fluctuating wind, wind speed when
Journey formula is as follows:
V=V0+Asinωt
Wherein V0For mean wind speed, ω is wind speed change frequency, A 0.4V0；
Boundary condition are as follows: the fluid domain grid model of foundation is a hexahedron cuboid, and right side is speed entrance boundary
Condition sets up corresponding wind inflow velocity, and left side sets up pressure export boundary condition, for defining the quiet of flowing exit
Pressure, hexahedral front and back and upside are set as symmetrical boundary condition, overflow wind will not from these three faces, downside is set as Gu Bi
Without slip boundary condition, for simulating the deck for placing sail.
Under preferred embodiment, vonKarman scale L is improvedvkFor
Wherein ΩcvFor control volume unit volume, constant CsIt can be demarcated by the simulation to every uniform same sex turbulent flow,
Finally take Cs=0.26, LvkIt is the minimum resolution for characterizing flow field whirlpool.
The invention has the benefit that calculation method of the invention passes through mould based on Fluid structure interaction model
Intend the sail wind shake under each operating condition to carry out the stress value that sail is subject to and calculate, since analog case and the practical wind of ship sail connect
Closely, the stress value that sail structure is subject to can be relatively accurately calculated, and then compared with the allowable stress of structural material, judges to accord with
Specification for structure is closed, therefore the design and construction cost of ship sail can be effectively reduced.
Detailed description of the invention
Fig. 1 is calculation flow chart of the invention；
Fig. 2 is hull coordinate system, sail coordinate system in one embodiment of the present of invention；
Fig. 3 is that the sail in one embodiment of the present of invention rotates angle；
Fig. 4 is the fluid domain grid model that the gambit in one embodiment of the present of invention is established；
Sail surface pressure distribution when Fig. 5 is 0 ° of angle of attack in one embodiment of the present of invention；
Fig. 6 is the finite element model of sail main body in one embodiment of the present of invention；
Fig. 7 is sail pedestal finite element model in one embodiment of the present of invention；
Fig. 8 is that one embodiment of the present of invention middle deck constrains schematic diagram；
Fig. 9 is the observation point top view that the response condition of sail is observed in one embodiment of the present of invention；
Figure 10 is the observation point front view that the response condition of sail is observed in one embodiment of the present of invention；
Figure 11 is the curve that 0 ° of angle of attack sail face maximum moving stress changes with frequency domain.
Specific embodiment
With reference to embodiment to i.e. a kind of hard sail wind vibration response based on fluid structurecoupling simulation of the present invention
Calculation method is further described.
In one embodiment, the invention discloses a kind of hard sail wind vibration response calculating based on fluid structurecoupling simulation
Method includes the following steps:
S1, the fluid domain numerical model that ship sail is established using gambit；
Suitable turbulence model is selected in S2, CFX, acts on fluid domain grid model and is obtained by numerical simulation calculation
Surface pressing file suffered by sail out；
S3, the solid domain structural model for establishing ship sail；
S4, the surface pressing numeric file that CFX is obtained is imported by interpolation method and is acted on sail structural model, it is complete
At coupling and numerical simulation calculation response；
S5, coupling and numerical simulation sail wind shake situation are completed, calculates the stress value that sail is subject to；
S6, judge that sail structural strength meets code requirement under each operating condition according to calculated result, and then instruct ship wind
The design and construction of sail.
Calculation method disclosed in the present embodiment, it is first determined hull coordinate system and sail coordinate system establish fluid
The geometrical model of domain and solid domain and given initial value and boundary condition, initial value is wind speed, boundary condition: the fluid domain grid of foundation
Model is a hexahedron cuboid, and right side is speed entrance boundary condition, sets up corresponding wind inflow velocity, left side is set up
Pressure export boundary condition, for defining the static pressure in flowing exit, hexahedral front and back and upside are set as symmetrical border item
Part overflow wind will not from these three faces, and downside is set as Gu Bi without slip boundary condition, for simulating the first for placing sail
Plate.
Numerical model is established later, and the stress value meter that sail is subject to is realized by simulating the sail wind shake under each operating condition
It calculates, which has comprehensively considered ship sail in practice based on the complete theoretical model and numerical model
Factor, it is practical close by wind environment with ship sail, can relatively accurately calculate the stress value that sail structure is subject to, so with knot
The allowable stress of structure material compares, and judgement meets specification for structure, therefore the design and construction of ship sail can be effectively reduced
Cost.
Described in the present embodiment using gambit establish fluid domain model selection wind parameter type be fluctuating wind；Wind
Load is towering, high-level structure primary load, even plays conclusive effect in some cases.According to the reality of a large amount of wind
Survey data can be seen that in the time-history curves of wind, including two kinds of ingredients: long period part, value is often at 10 minutes or more；It is short
Cycle portions are typically only several seconds (frequency of fluctuating wind is 0.1-1Hz, i.e. the period is 1-10s).According to it is above two at
Point, practical upper often keep watch is divided into average wind (i.e. constant wind) and fluctuating wind to be analyzed.Fluctuating wind is irregular due to wind
Caused by property, its intensity changes at random at any time, and since its period is shorter, thus its interaction property is power,
It will cause the vibration of structure.Take the time-histories formula of fluctuating wind wind speed as follows:
V=V0+Asinωt
Wherein V0For mean wind speed, ω is wind speed change frequency, A 0.4V0。
SST-SAS turbulence model equation is used using the ship sail fluid domain that CFX is established described in the present embodiment.
SST-SAS turbulence model is based on SST-RANS, it is contemplated that turbulent viscosity dissipates in Disengagement zone, in standard SST
Than adding Q in dissipative shock wave ω equationSAS, wherein SST turbulence model is the deformation of standard k-w model；It will using mixed function
Standard k-e turbulence model combines with k-w turbulence model, and k-w is used to solve in boundary layer, and k-e is for outside boundary layer
It solves, and QSASBy stationary zones (QSAS=0) with unstable state region (QSAS> 0) it distinguishes；Specific equation are as follows:
In formula, ρ is density, and t is the time, and k is Turbulent Kinetic, PkItem, μ are generated for kinetic energytIt is turbulent flow stickiness, ω is specific consumption
The rate of dissipating, μiIt is time averaged velocity, F1For mixed function, σkAnd σεIt is k, turbulent flow Pood's number of ε, β*It is constant, S with βijAveragely to answer
Variability.
QSASL invkFor von Karman scale, L is modelling rapids stress scale, LvkIt can be in unstable state area according to current
Known flow field analysis whirlpool dynamic, adjusts length dimension in turbulence model in real time, quasi- on the basis of adaptive whole flow field grid
Local flow scale really is embodied, when flowing into unstable state region, causes to increase than dissipative shock wave ω, tubulence energy k is reduced, and is added
Fast model separates the generation of stream in unstable state region；Wherein model constants are C=2.0, σΦ=2/3,κ=
0.41,Cμ=0.09,
QSASBy stationary zones (QSAS=0) with unstable state region (QSAS> 0) it distinguishes, successfully remains RANS
Model near wall boundary layer (stationary zones) excellent performance.In unstable state region, LvkUnstable state region can be differentiated
It interior whirlpool situation and realizes dynamic change, causes to increase than dissipative shock wave ω, tubulence energy k is reduced, so that turbulent viscosity is dissipated in
It is reduced in distinguishable pulsation, develops the equation of motion from stable state to unstable state, and turbulent viscosity is promoted to dissipate, exacerbate SAS
Model separates the generation of stream in unstable state region.
But the turbulent viscosity dissipation deficiency of higher wave number may result in SAS model performance variation, only when grid and when
Spacer step long enough hours, its excellent performance of SAS model competence exertion.The emphasis that turbulent dissipation appropriate becomes SAS model is provided,
Now commonly used following methods solve:
Improve von karman scale Lvk:
Wherein ΩcvFor control volume unit volume, constant CsIt can be demarcated by the simulation to every uniform same sex turbulent flow,
Finally take Cs=0.26.LvkIt is the minimum resolution for characterizing flow field whirlpool, is transmitted when flow field whirlpool is descending, whirlpool ruler
When degree is less than the grid limit of resolution in the region RANS, LvkReduction leads to QSASDissipative term increases to which whirlpool separates in advance, serious shadow
RANS model performance is rung.And L in above formulavkIt can be with basin grid dynamic change, when coarse grid occurs, LvkIncrease to reduce
QSASInfluence, ensure that the stability in the region RANS.Turbulent flow mould according to existing experimental study, as vortex-induced vibration research
Type SAS is highly effective, so will be studied using the SAS turbulence model.
After establishing numerical model, then pass through the sail stress of (0 °, 15 °, 30 °, 45 °) under each angle of attack of numerical simulation calculation
Situation obtains surface stress file.Numerical simulation calculation uses finite volume method.Finite volume method, this method is exactly will stream
Domain space is divided into limited volume, and each volume is indicated with one group of physical quantity come the continuity equation in discrete flow field, according to control
Equation integrates each limited bulk to obtain the algebraic equation to each volume, is then iterated solution.
Establish the solid domain structural model of ship sail.The surface pressing numeric file that CFX is obtained is led by interpolation method
Enter and acts on sail structural model.It is because can be incited somebody to action when carrying out multiple physical field wind-structure interaction using interpolation method
The physical quantity of different physical field is contacted, this is related to the accuracy of subsequent solution procedure.It is flowed in sail vortex-induced vibration problem
Physical quantity between body field and solid domain is mutually transmitted by the interface of the two, but due to structured grid node
It is not correspondingly that load can not be transmitted directly, so must be according to the location information of node with flow field grid node
Interpolation processing is carried out to respective physical amount and completes data exchange.
Interpolation is carried out using the node pressure in interpolation method stream field result, applies the calculated result in flow field to realize
Onto coupling interface, completes one-way flow and consolidate coupling operational.For the accuracy for guaranteeing data transmitting, at structure and fluid interface face
It must satisfy the identical condition of boundary condition, while under the same coordinate system, locations of structures and structure in fluid calculation model
The position consistency of finite element model.ANSYS Workbench built-in interpolating function, principle are to create on the coupling surface of structure
Node is made, interpolation is carried out according to CFX calculated result, then replicates node along the practical length of sail, when the node of duplication
When consistent with actual sail, each node pressure that interpolation goes out is applied on duplication posterior nodal point.Therefore interpolation method can be accurate
Load node all directions on pressure changing keep result relatively accurate.
It completes coupling and numerical simulation calculation uses Finite Element Method, applied program is ANSYS
Workbench.Its basic thought be surface pressing suffered by sail is calculated using finite volume method by CFX software, and
Surface pressing in the destination file of CFX is imported in Stastic Structure and made by ANSYS Workbench platform
For sail structure, the response of coupling joint account sail is completed.The visible Fig. 1 of main-process stream.
Wherein solid portion conservation equation can be exported by Newton's second law, it may be assumed that
In formula, ρsFor density of solid；σsFor cauchy stress tensor；fsFor volume vector；For solid domain local acceleration arrow
Amount.
At fluid structurecoupling interface, the variables phase such as fluid and solid stress τ, displacement d, heat flow q, temperature T should be met
Deng or conservation, that is, meet following equation:
In formula, subscript s and f respectively represent fluid domain and solid domain.
Coupling and numerical simulation sail wind shake situation are completed, the stress value that sail is subject to is calculated.
Numerical simulation calculation use Finite Element, this method by the discrete subdivision of computational domain be it is limited do not overlap and
Unit interconnected selects basic function in each unit, with the linear combination of unit basic function come true in approximation unit
Solution, on entire computational domain overall basic function can be seen as by each unit set of basis function at, then the solution in entire computational domain
It can be regarded as being deconstructed by the approximation on all units.
According to the stress response value for determining hard sail wind shake is calculated, it is compared with the structure of ship sail with material,
And then determine whether the design of the sail meets the requirement of structural strength.Play the work for instructing design and the construction of ship sail
With.Greatly reduce the cost of ship sail design.
In one embodiment, the step of a kind of hard sail wind vibration response calculates is disclosed:
(1) the fluid domain grid model of ship hard sail is established using gambit.Hull coordinate system and wind are determined first
Sail coordinate system, as shown in Figure 2.
Hull coordinate system is denoted as Coor.0 as right-handed coordinate system.Wherein origin: Frame0, central fore-and-aft vertical plane and baseline
Intersection；X-axis: longitudinal axis, positive direction are to be directed toward bow from stern；Y-axis: lateral shaft, positive direction are to be directed toward a left side from center line of boat
The side of a ship；Z axis: vertical axis, positive direction are that vertical XY faces upward (right-handed coordinate system).
Two local fixed coordinate system Coor.P and Coor.S, right-handed coordinate system, for showing that sail is located at ship are set again
Spatial position on body.Origin: sail pedestal upper flange bottom center point is located at the spatial position on hull, and coordinate is
(Fr214, ± 21250,31225A/B)@Coor.0.L axis is positive: with Coor.0X axis；T axis is positive: with Coor.0Y axis；V axis is just
To: with Coor.0Z axis.
Sail coordinate system Coor.P ' and Coor.S ', right-handed coordinate system, for calculating sail rotation angle, showing are set again
The angle of attack of wind.Wherein origin: sail pedestal upper flange bottom center point is located at the spatial position on hull, coordinate be (Fr214,
±21250,31225A/B)@Coor.0。
X-axis is positive: being parallel to sail face direction；Y-axis is positive: perpendicular to sail face direction；Z axis is positive: with Coor.0z axis.
Provide that sail rotation angle is x-axis and L axle clamp angle, dextrorotation is positive；As shown in Figure 3.
It is assumed that sail is rigid model, regards sail as cross section identical structural body and study, the height of sail is
39.68m, width 14.80m, with a thickness of 4.172m, sail body height is 33.92m, mast height 37.415m, base height
Model such as Fig. 4 for 2.265m, after establishing fluid domain geometrical model and grid division with gambit, it is seen that size of mesh opening is by wing sail
Shape express than more complete, the transition in boundary layer is also relatively uniform, and mesh quality is totally good；Model scale ratio is 1:
16。
(2) boundary condition is established in CFX, selects suitable turbulence model, and sail institute is obtained by numerical simulation calculation
The surface pressing received.
Wind speed formula is V=V0+ Asin ω t, wherein V0It is taken as 25m/s, ω value is 0.1-1.0Hz, A 0.4V0.The angle of attack
Respectively 0 °, 15 °, 30 °, 45 ° four kinds substitute into be calculated respectively.The size of its computational domain chooses entrance away from sail cloth leading edge
1.5 times of sails are wide, and outlet is 1 times of sail height away from sail cloth rear.Boundary condition is arranged as follows: the fluid domain grid model of foundation is one
Hexahedron cuboid, right side are speed entrance boundary condition, set up corresponding wind inflow velocity, left side sets up pressure export side
Boundary's condition, for defining the static pressure in flowing exit, hexahedral front and back and upside are set as symmetrical boundary condition, make wind will not
It is overflowed from these three faces, downside is set as Gu Bi without slip boundary condition, for simulating the deck for placing sail.Turbulence model choosing
With SST-SAS turbulence model.The sail force-bearing situation under each angle of attack is calculated using finite volume method, wherein sail table when 0 ° of angle of attack
The distribution situation of surface pressure；Such as Fig. 5.
(3) the solid domain structural model of ship sail is established.It is influenced to reduce boundary to the calculating of sail inherent characteristic,
To near sail pedestal and deck etc. also modeled.The sail face of sail, mast, pedestal plate face use shell unit
(shell63) it simulates, the reinforcing rib of sail, pedestal etc. is simulated with beam element (beam188).Fig. 6 and Fig. 7 is respectively that sail is limited
Meta-model and pedestal finite element model.
Boundary constraint is applied to structural model, the visible Fig. 8 of the boundary constraint on deck is constrained to 1 partial bulkhead of arrow in figure
The freedom degree in the x, y, z direction of fringe node and the torsion in the direction y；The x, y, z direction of 2 partial bulkhead fringe node of arrow from
By the torsion in degree and the direction x.
(4) the surface pressing numeric file that CFX is obtained is imported by interpolation method and is acted on sail structural model.Wind
The approximately periodic variation of pressure caused when sail surface is flowed through, sail surface normal pressure maximum is chosen and negative pressure is maximum
When corresponding time step sail pressure loading, subtract each other to obtain the distribution of sail surface fluctuation pressure range, will after load transfer plus
To finite element model.
(5) sail wind shake situation under coupling and each operating condition of numerical simulation is completed, the stress value that sail is subject to is calculated and sets up 12
The response condition of a observation point observation sail, 12 observation points are respectively A1, A2, A3, A4, B1, B2, B3, B4, C1, C2,
C3, C4, observation point schematic diagram are shown in Fig. 9 and Figure 10.It can be calculated under each angle of attack of each frequency, each observation point dynamic respond value, and
Dynamic stress distribution, and then draw out the curve that each angle of attack sail face maximum moving stress changes with frequency domain.Figure 11 be 0 ° of angle of attack sail face most
The curve that big dynamic stress changes with frequency domain.
(6) judge that sail structural strength meets code requirement under each operating condition according to calculated result, and then instruct ship wind
The design and construction of sail.
The covering and skeleton of known sail select stainless steel material (316L, allowable stress=195/1.33=146.61MPa
With allowable stress=205/1.33=154.13MPa), mast selects high intensity for ship steel (DH36, allowable stress=355/1.77
=200.56MPa).And for being in the fluctuating wind of sinusoidal variations (frequency 0.1-1.0Hz), and consider upper 0.2 times of wind shake
Coefficient, comparison between calculation results show that the stress value in sail face and sail skeleton under each operating condition is both less than corresponding allowable stress,
Deformation also meets related specification requirement, and the sail wind shake of the condition not will lead to sail destruction.This structure and material can be used for
The design of sail reduces the cost of ship sail design.
The foregoing is only a preferred embodiment of the present invention, but scope of protection of the present invention is not limited thereto,
Anyone skilled in the art within the technical scope of the present disclosure, according to the technique and scheme of the present invention and its
Inventive concept is subject to equivalent substitution or change, should be covered by the protection scope of the present invention.
Claims (3)
1. a kind of hard sail wind vibration response calculation method based on fluid structurecoupling simulation, which is characterized in that specific steps are as follows:
S1, hull coordinate system and sail coordinate system are established in gambit, ship hard is established according to the actual parameter of sail
The fluid domain grid model of sail, and the fluid domain grid model is imported in CFX；
The initial value wind speed and boundary condition of ship hard sail are established in S2, CFX, select turbulence model and act on the fluid
On the grid model of domain, surface pressing suffered by ship hard sail is calculated using finite volume method, and obtain surface pressing
Numeric file；The turbulence model is SST-SAS turbulence model, the equation group of the turbulence model are as follows:
In formula, ρ is density, and t is the time, and k is Turbulent Kinetic, PkItem, μ are generated for kinetic energytTurbulent flow stickiness, ω be than dissipative shock wave,
μiIt is time averaged velocity, F1For mixed function, σkAnd σεIt is k, turbulent flow Pood's number of ε, β*It is constant, S with βijFor mean strain rate；
The SST-SAS turbulence model equation is in view of turbulent viscosity dissipation in Disengagement zone, in the ratio dissipative shock wave side ω of standard SST
Q is added in journeySAS?；SST turbulence model is the deformation of standard k-w model；Using mixed function by standard k-e turbulence model with
K-w turbulence model combines, and k-w is used to solve in boundary layer, and k-e outside boundary layer for solving, and QSASBy stable state
Region, i.e. QSAS=0, with unstable state region, i.e. QSAS> 0, it distinguishes；QSASL invkFor von Karman scale, L is modelling
Rapid stress scale, LvkIn unstable state area according to currently known flow field analysis whirlpool dynamic, length ruler in turbulence model is adjusted in real time
Degree；Model constants are as follows:
C=2.0, σΦ=2/3,κ=0.41, Cμ=0.09,
S3, sail sail face is established in ANSYS software, mast, the solid domain structural model on pedestal and neighbouring deck and to institute
Solid domain structural model addition edge-restraint condition is stated, the constraint condition is the freedom degree of solid domain structural model fringe node
Constraint and torsional restraint；
S4, the surface pressing numeric file is imported by interpolation method and acts on the solid domain structure of ANSYS ship sail
On model, at the fluid structurecoupling interface of fluid domain grid model and solid domain structural model, fluid should be met and answered with solid
Power τ, displacement d, heat flow q, the variables such as temperature T are equal or conservation, that is, meet following equation:
In formula, subscript s and f respectively represent fluid domain and solid domain；
S5, numerical simulation calculation hard sail wind shake situation is carried out using Finite Element, calculates the stress that hard sail is subject to
Value；
S6, judge whether hard sail structural strength accords under each operating condition according to the stress value result that the hard sail of calculating is subject to
Code requirement is closed, and then instructs the design and construction of ship hard sail.
2. the hard sail wind vibration response calculation method according to claim 1 based on fluid structurecoupling simulation, feature exist
In: the wind type that the fluid domain grid model established in the S1 is selected is fluctuating wind, and the time-histories formula of wind speed is as follows:
V=V0+Asinωt
Wherein V0For mean wind speed, ω is wind speed change frequency, A 0.4V0；
Boundary condition are as follows: the fluid domain grid model of foundation is a hexahedron cuboid, and right side is speed entrance boundary condition,
Corresponding wind inflow velocity is set up, left side sets up pressure export boundary condition, for defining the static pressure in flowing exit, six faces
The front and back and upside of body are set as symmetrical boundary condition, overflow wind will not from these three faces, downside is set as Gu Bi without sliding
Boundary condition, for simulating the deck for placing sail.
3. the hard sail wind vibration response calculation method according to claim 1 based on fluid structurecoupling simulation, feature exist
In: improve von Karman scale LvkFor
Wherein ΩcvFor control volume unit volume, constant CsIt is demarcated by the simulation to every uniform same sex turbulent flow, is finally taken
Cs=0.26, LvkIt is the minimum resolution for characterizing flow field whirlpool.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811607705.4A CN110309521A (en) | 2018-12-27 | 2018-12-27 | A kind of hard sail wind vibration response calculation method based on fluid structurecoupling simulation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811607705.4A CN110309521A (en) | 2018-12-27 | 2018-12-27 | A kind of hard sail wind vibration response calculation method based on fluid structurecoupling simulation |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110309521A true CN110309521A (en) | 2019-10-08 |
Family
ID=68074182
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811607705.4A Withdrawn CN110309521A (en) | 2018-12-27 | 2018-12-27 | A kind of hard sail wind vibration response calculation method based on fluid structurecoupling simulation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110309521A (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110955991A (en) * | 2019-11-18 | 2020-04-03 | 华北水利水电大学 | Fluid-solid coupling calculation method for interface bidirectional data exchange |
| CN113378489A (en) * | 2021-05-27 | 2021-09-10 | 杭州电子科技大学 | OpenFOAM platform-based low-magnetic Reynolds number magnetic fluid turbulence analysis system and analysis method |
| CN115034162A (en) * | 2022-07-12 | 2022-09-09 | 北京航空航天大学 | Grid self-adaptive turbulence simulation method based on turbulence energy spectrum coupling k-epsilon series model |
| CN117669277A (en) * | 2024-02-02 | 2024-03-08 | 中国科学院大学 | Heat exchange flow channel optimization method based on fluid-solid coupling effect and heat exchange flow channel |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104508422A (en) * | 2012-05-30 | 2015-04-08 | 赛创尼克株式会社 | System and method for providing information on fuel savings, safe operation, and maintenance by real-time predictive monitoring and predictive controlling of aerodynamic and hydrodynamic environmental internal/external forces, hull stresses, motion with six degrees of freedom, and the location of marine structure |
| CN106383930A (en) * | 2016-08-31 | 2017-02-08 | 中国石油大学(华东) | Multiple fluid-solid coupling calculation method for tail bearing-rotor system |
| CN107605666A (en) * | 2017-09-11 | 2018-01-19 | 北京金风科创风电设备有限公司 | Envelope structure with vortex-induced vibration inhibiting function and method for inhibiting vortex-induced vibration |
-
2018
- 2018-12-27 CN CN201811607705.4A patent/CN110309521A/en not_active Withdrawn
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104508422A (en) * | 2012-05-30 | 2015-04-08 | 赛创尼克株式会社 | System and method for providing information on fuel savings, safe operation, and maintenance by real-time predictive monitoring and predictive controlling of aerodynamic and hydrodynamic environmental internal/external forces, hull stresses, motion with six degrees of freedom, and the location of marine structure |
| CN106383930A (en) * | 2016-08-31 | 2017-02-08 | 中国石油大学(华东) | Multiple fluid-solid coupling calculation method for tail bearing-rotor system |
| CN107605666A (en) * | 2017-09-11 | 2018-01-19 | 北京金风科创风电设备有限公司 | Envelope structure with vortex-induced vibration inhibiting function and method for inhibiting vortex-induced vibration |
Non-Patent Citations (2)
| Title |
|---|
| 李思洋: "《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》", 15 April 2018 * |
| 陈东阳等: "旋转弹丸气动特性的尺度自适应模拟", 《哈尔滨工程大学学报》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110955991A (en) * | 2019-11-18 | 2020-04-03 | 华北水利水电大学 | Fluid-solid coupling calculation method for interface bidirectional data exchange |
| CN110955991B (en) * | 2019-11-18 | 2023-03-21 | 华北水利水电大学 | Fluid-solid coupling calculation method for interface bidirectional data exchange |
| CN113378489A (en) * | 2021-05-27 | 2021-09-10 | 杭州电子科技大学 | OpenFOAM platform-based low-magnetic Reynolds number magnetic fluid turbulence analysis system and analysis method |
| CN115034162A (en) * | 2022-07-12 | 2022-09-09 | 北京航空航天大学 | Grid self-adaptive turbulence simulation method based on turbulence energy spectrum coupling k-epsilon series model |
| CN115034162B (en) * | 2022-07-12 | 2024-06-07 | 北京航空航天大学 | Grid self-adaptive turbulence simulation method based on turbulence energy spectrum coupling k-epsilon series model |
| CN117669277A (en) * | 2024-02-02 | 2024-03-08 | 中国科学院大学 | Heat exchange flow channel optimization method based on fluid-solid coupling effect and heat exchange flow channel |
| CN117669277B (en) * | 2024-02-02 | 2024-04-26 | 中国科学院大学 | Heat exchange flow channel optimization method based on fluid-solid coupling effect and heat exchange flow channel |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110309521A (en) | A kind of hard sail wind vibration response calculation method based on fluid structurecoupling simulation | |
| Vukčević et al. | Implementation of the ghost fluid method for free surface flows in polyhedral finite volume framework | |
| Bhalla et al. | Simulating water-entry/exit problems using Eulerian–Lagrangian and fully-Eulerian fictitious domain methods within the open-source IBAMR library | |
| Feng et al. | Numerical investigation of viscous effects on the gap resonance between side-by-side barges | |
| Kristiansen et al. | Gap resonance analyzed by a new domain-decomposition method combining potential and viscous flow DRAFT | |
| Yang et al. | A simple and efficient direct forcing immersed boundary framework for fluid–structure interactions | |
| KR101267627B1 (en) | SPH Fluid simulation method and system for Multi-Level Vorticity, recording medium for the same | |
| Westphalen et al. | Investigation of wave-structure interaction using state of the art CFD techniques | |
| Xie | Numerical study of breaking waves by a two‐phase flow model | |
| Yue et al. | Large eddy simulation of turbulent open-channel flow with free surface simulated by level set method | |
| Kumar et al. | Implicit boundary method for finite element analysis using non‐conforming mesh or grid | |
| Apsley | CFD calculation of turbulent flow with arbitrary wall roughness | |
| Nasif et al. | Supercritical flow characteristics in smooth open channels with different aspect ratios | |
| Lee et al. | Comparative study on the breaking waves by a piston-type wavemaker in experiments and SPH simulations | |
| Kim et al. | Nonlinear time-domain NWT simulations for two types of a backward bent duct buoy (BBDB) compared with 2D wave-tank experiments | |
| Martin et al. | Accurate modeling of the interaction of constrained floating structures and complex free surfaces using a new quasistatic mooring model | |
| Xu et al. | An immersed boundary method with y+‐adaptive wall function for smooth wall shear | |
| Yue et al. | Coherent structures in open-channel flows over a fixed dune | |
| Rapaka et al. | An immersed boundary method for direct and large eddy simulation of stratified flows in complex geometry | |
| CN109740233A (en) | Hard sail vortex-induced vibration based on fluid structurecoupling simulation responds preparation method | |
| Ghoreishi et al. | Vorticity-based polynomial adaptation for moving and deforming domains | |
| Li et al. | Numerical investigation of swimmer’s gliding stage with 6-DOF movement | |
| Kimura et al. | Computations of suspended sediment transport in a shallow side-cavity using depth-averaged 2D models with effects of secondary currents | |
| Liu | A CFD study of fluid-structure interaction problems for floating offshore wind turbines | |
| Garzon et al. | A coupled level set-boundary integral method for moving boundary simulations |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WW01 | Invention patent application withdrawn after publication |
Application publication date: 20191008 |
|
| WW01 | Invention patent application withdrawn after publication |