10/22/24, 7:18 PM PDF translate: -----A numerical method for predicting the hydroelastic response ofmarine propeller.

pdf | Sider

Recherche océanique appliquée 74 (2018) 188-204

Listes de matières disponibles sur ScienceDirect

Recherche océanique appliquée

page d'accueil du journal : www.elsevier.com/locate/apor

Une méthode numérique pour prédire la réponse hydroélastique des hélices

marines
Jiasheng Li a,b , Zhiqiang Rao c , Jinpeng Su a,b , Yegao Qu a , Hongxing Hua a,b,∗
a Laboratoire clé d'État des systèmes mécaniques et des vibrations, Université Jiao Tong de Shanghai, Shanghai, 200240, RP
Chine b Centre d'innovation collaborative pour l'exploration avancée des navires et des eaux profondes (CISSE), Shanghai,
200240, RP Chine c Centre de recherche scientifique sur les navires de Chine ( CSSRC), Wuxi, 214082, Jiangsu, RP Chine

article infos abstrait

Historique de l'article : Reçu le 25 juin Cet article se concentre sur le développement d'un modèle numérique pour prédire les réponses hydroélastiques des
2017 Reçu sous sa forme révisée le 25 hélices marines oscillant dans le sillage d'un sous-marin. Les matrices de masse ajoutée et d'amortissement dues à
janvier 2018 Accepté le 14 février 2018 une interaction fluide-structure fortement couplée ont été prises en compte. Des méthodes de panneaux
tridimensionnels dans les domaines temporel et fréquentiel combinées à la méthode des éléments finis ont été
utilisées pour étudier les performances hydrodynamiques de l'hélice. Les méthodes des panneaux ont été utilisées
Mots-clés : Interaction pour évaluer les forces hydrodynamiques générées par le sillage, et la méthode des éléments finis a été appliquée
fluide-structure pour déterminer la réponse hydroélastique de l'hélice due à la fluctuation de pression. Pour prédire les réponses
Hydroélastique structurelles, une méthode de superposition de modes combinée à la méthode de Wilson a été utilisée pour
Méthode du panneau d'hélice surmonter la faible efficacité numérique causée par la masse ajoutée asymétrique et les matrices d'amortissement.
Méthode des éléments finis Le modèle numérique proposé a été validé par comparaison avec d'autres solutions numériques. Les performances
Méthode de superposition de des hélices ont été examinées en considérant les interactions fluide-structure unidirectionnelles et bidirectionnelles.
modes
© 2018 Elsevier Ltd. Tous droits réservés.

1. Présentation ont été utilisés dans les industries maritimes pratiques, et il est d’une
importance cruciale pour les concepteurs d’avoir une connaissance
Lors de la rotation d'une hélice dans le sillage spatialement non uniforme approfondie du comportement hydrodynamique des hélices élastiques.
d'un sous-marin, des forces fluctuantes instables peuvent induire des Pour analyser le problème d'interaction fluide-structure des hélices sous-
vibrations de l'hélice elle-même ou être transmises à la structure de la marines, deux méthodes peuvent être utilisées, à savoir la méthode
coque du sous-marin à travers l'arbre et les roulements. Ces forces unidirectionnelle et la méthode bidirectionnelle. Dans la première
comprennent les forces et moments de poussée, de couple, verticaux et méthode, les forces hydrodynamiques et les réponses structurelles de
transversaux. Étant donné que ces forces et moments peuvent entraîner de l'hélice sont calculées l'une après l'autre. Les équations régissant le fluide
fortes vibrations de la coque et des bruits sous-marins, une méthode fiable entourant l'hélice sont d'abord résolues, puis les forces hydrodynamiques
pour prédire la réponse hydroélastique des hélices marines est d'une grande résultantes agissant sur l'hélice sont calculées, qui sont utilisées pour
importance. calculer les réponses structurelles de l'hélice. Dans les interactions fluide-
Au cours des dernières décennies, un certain nombre de méthodes structure bidirectionnelles, le fluide et la structure sont modélisés comme
analytiques et numériques ont été développées pour prédire les forces
un système couplé. Le mouvement structurel de l'hélice provoque une
hydrodynamiques des hélices sous-marines, notamment la méthode du
modification des forces hydrodynamiques agissant sur la surface
moment ou des éléments de pale, la méthode de la ligne de levage, la
méthode de la surface de levage et la méthode des panneaux. Ces méthodes
mouillée de l'hélice, ce qui entraîne à son tour une modification du
supposent que les pales de l’hélice sont rigides et négligent l’effet de la
mouvement de la structure. Deux grandes classes de formulations
déformation des pales sur le champ d’écoulement environnant. En général, d'interaction fluide-structure (FSI) bidirectionnelles existent dans la
ces méthodes sont efficaces sur le plan informatique, mais dans certains cas, littérature [1,2], à savoir les méthodes faiblement et fortement couplées,
la précision de ces méthodes pour la prédiction des forces hydrodynamiques qui peuvent être utilisées pour l'analyse des problèmes FSI des hélices
des hélices n'est pas satisfaisante. Ces dernières années, les hélices marines sous-marines. Dans les approches faiblement couplées (également
flexibles ont appelées approches itératives), les équations fluides et structurelles sont
résolues de manière itérative de manière découplée, ce qui permet
l’utilisation potentielle de solveurs fluides et structurels existants et bien
validés. Cependant, ces méthodes peuvent souffrir de certains manques
∗ Auteur correspondant à : State Key Laboratory of Mechanical System and de convergence, qui ont été constatés dans de nombreuses situations, et
Vibration, Shanghai Jiao Tong University, Shanghai 200240, RP Chine. une littérature considérable a été consacrée à ce sujet.
Adresse e-mail : [email protected] (H. Hua).

https://doi.org/10.1016/j.apor.2018.02.012 0141-1187/©
2018 Elsevier Ltd. Tous droits réservés.

https://sider.ai/pdf-translate/E0G5UEKEG8A?fid=0 1/17
10/22/24, 7:18 PM PDF translate: -----A numerical method for predicting the hydroelastic response ofmarine propeller.pdf | Sider

J. Li et coll. / Recherche océanique appliquée 74 (2018) 188-204 189

Cette théorie a été employée par Conolly [4]. Cette méthode a été
Nomenclature appliquée avec succès pour les cas de pales larges, mais l'hypothèse de
formes symétriques et de déflexion normale invariante d'une section limite
Coefficient de Poisson
l'application de l'approche. Les méthodes par éléments finis (FEM) ont
E Module de Young également été développées pour l'analyse de la résistance des pales, car
s Densité structurelle aucune hypothèse supplémentaire sur la géométrie des pales n'est
M Matrice de masse globale nécessaire dans de telles méthodes. Des éléments plans triangulaires ont
C Matrice d'amortissement globale été utilisés par Genalis [5] pour la modélisation structurelle de la pale.
K Matrice de rigidité globale Atkinson [6] a appliqué des éléments de coque épaisse super-
C̋ Matrice d'amortissement provoquée par les forces de Coriolis paramétriques dans l'analyse de la pale. Ma [7] a utilisé des éléments de
brique isopérimétriques quadratiques 3D pour analyser une hélice très
K̋ Matrice de rigidité causée par les forces centrifuges
asymétrique. Des éléments triangulaires à coque mince ont été utilisés par
̈u Vecteur d'accélération nodale
Sontvedt [8] pour analyser des hélices modérément et fortement inclinées.
u̇ Vecteur de vitesse nodale
Il est à noter que dans toutes ces méthodes, les charges hydrodynamiques
toi Vecteur de déplacement nodal
agissant sur la surface de la pale ont été obtenues par la méthode quasi-
Fw Vecteur de force hydrodynamique
stationnaire, la méthode des lignes de levage ou la méthode des surfaces
F̋ Vecteur de force centrifuge portantes.
V total Vitesse totale du flux
Pour les interactions fluide-structure bidirectionnelles, Kuo et Vorus [9]
V Sillage d'afflux non uniforme
ont introduit une approche structurelle et hydrodynamique entièrement
Potentiel de vitesse de perturbation correspondant au champ combinée pour l'analyse des contraintes des pales. Un ensemble de
d'écoulement induit par l'hélice conditions aux limites ont été établies sur l'interaction fluide-structure, et
n Vecteur normal de l'unité extérieure sur la surface de la lame
une méthode de panneaux sur la surface de cambrure a été appliquée pour
∂ ∂n
La dérivée de par rapport à la direction normale déterminer les charges hydrodynamiques de la pale, la masse ajoutée et
les matrices d'amortissement dues au mouvement élastique de la pale. Le
je Vecteur de déplacement des points nodaux sur la surface de la problème couplé a été résolu dans le domaine fréquentiel en utilisant le
pale FEM basé sur des éléments de briques isoparamétriques linéaires 3D. Il
V0 Partie uniforme du sillage s'avère que l'inclinaison de la pale a un effet sur les matrices
w (Rwake , t) Le saut potentiel à travers les feuilles de sillage (Rre , d'amortissement et de masse ajoutées et sur les modes humides. Atkinson
t − t′) Le saut potentiel à travers la surface de la pale au et Glover [10] ont étudié les comportements dynamiques des hélices
bord de fuite élastiques en utilisant une méthode de surface de portance instable, qui a
t' Temps nécessaire au fluide pour parcourir la surface du été couplée à un modèle par éléments finis de l'hélice basé sur des
éléments de coque épais. Les résultats ont montré que la taille et la forme
sillage depuis le bord de fuite de la pale Rre jusqu'au
des cavités produites peuvent être grandement affectées par les
point de sillage Rwake
mouvements hydroélastiques de l'hélice. Cependant, leur analyse s’est
Densité du fluide
limitée aux problèmes d’écoulement constant et, par conséquent, les
r Le potentiel de perturbation dû à l'hélice rigide avançant dans forces dynamiques des pales ne peuvent pas être capturées. Lin et Lin
le sillage d'afflux non uniforme
[11] ont analysé l'analyse hydroélastique des hélices marines basée sur un
v Le potentiel de perturbation dû à l’hélice élastique vibrant FEM non linéaire 3D couplé et une méthode de surface de levage stable
dans la partie uniforme du sillage sans cavitation. Les coefficients de poussée, d'efficacité du couple et de
Fr Forces hydrodynamiques générées par un sillage non déflexion ont été inclus dans leur calcul. Il a été constaté que la déflexion
uniforme de la lame augmente à mesure que l’épaisseur de la lame diminue. Young
Fv Forces hydrodynamiques générées par l'interaction fluide- [12, 13] a présenté une analyse d'interaction fluide-structure fortement
structure (vibrant dans un écoulement uniforme) entraînant couplée pour prédire les comportements hydroélastiques des hélices
la masse ajoutée et les matrices d'amortissement du fluide marines métalliques et composites dans des écoulements sous-cavitaires
̋ Vitesse de rotation de l'hélice { } et cavitaires. Un ensemble de conditions aux limites simplifiées a été
s La matrice de zone, s = diag s1, s2, . . ., sN introduit sur l’interaction de la pale et du fluide. Le fluide est modélisé à
k La fréquence des forces hydrodynamiques générées par un l'aide d'une méthode d'éléments limites 3D (BEM), qui a été combinée à
sillage non uniforme un modèle d'éléments finis de la pale. La modification des coefficients
je = √ d'influence due à la déformation de la pale a été examinée dans le cas de
je −1
U U = C−1k,vG, matrice des coefficients du panel l'hélice marine en composite flexible. Afin d'économiser les ressources de
calcul requises par les équations du système, les matrices asymétriques
E Matrice de coefficient d'interaction fluide-structure ajoutées dues au fluide ont été transformées en matrices de symétrie en
causée par (( ı · ∇ ) V 0 ) · n utilisant une technique de regroupement de type HRZ. Il a été constaté
UXL Matrices de coefficients de panel
que pour l’analyse des comportements hydrodynamiques d’hélices
V xl0 Matrice de vitesse constituée de V0xl (Ri) V0xl (Ri) composites flexibles, l’interaction fluide-structure fortement couplée doit
Composante de V0 dans la direction xi sur le point Ri Z1, être prise en compte. Neugebauer et coll. [14] ont étudié les problèmes
Z2 Matrices de fonctions de forme d’interaction fluide-structure des hélices en utilisant la méthode de
dynamique des fluides computationnelle. Les déplacements nodaux, les
vitesses et les accélérations à proximité de l'hélice ont été inclus dans leur
discussion de ces problèmes. Dans les méthodes fortement couplées calcul, et la différence entre les pressions hydrodynamiques obtenues par
(également appelées approches monolithiques), les équations du fluide et la méthode couplée et non couplée a été examinée. Lui et coll. [15] ont
des structures sont résolues simultanément de manière entièrement développé un algorithme de couplage FEM/CFD 3D pour étudier les
couplée. Le principal avantage est que les méthodes monolithiques ont comportements hydroélastiques d'une hélice modèle composite
tendance à être plus robustes que les méthodes faiblement couplées. hautement asymétrique. Les fréquences naturelles humides, les formes
modales, les coefficients de charge hydrodynamiques et les déplacements
Pour l'interaction fluide-structure unidirectionnelle, Taylor [3] a analysé des nœuds de pointe, etc., ont été examinés. Lee et coll. [16] ont introduit
la résistance d'une pale d'hélice en considérant la pale comme une poutre un algorithme basé sur BEM-FEM pour prédire les coefficients de
en porte-à-faux. Les résultats des réponses structurelles de la pale étaient poussée et de couple instables pour les hélices flexibles. Maljaars et
d'applications limitées en raison de la géométrie simple de la pale. Pour Kaminski [2] ont présenté une revue des recherches antérieures sur
améliorer les prédictions des contraintes de la pale, une coque mince

https://sider.ai/pdf-translate/E0G5UEKEG8A?fid=0 2/17
10/22/24, 7:18 PM PDF translate: -----A numerical method for predicting the hydroelastic response ofmarine propeller.pdf | Sider

190 J. Li et coll. / Recherche océanique appliquée 74 (2018) 188-204

les analyses hydroélastiques des hélices flexibles. Ils ont conclu que la
majorité des méthodes développées se limitaient aux analyses des
performances constantes des hélices flexibles. Concernant les
problèmes FSI des hélices sous-marines, Young [12, 13] a utilisé une
approche faiblement couplée pour résoudre les équations structurelles
et fluides couplées dans les deux sens, et le couplage bidirectionnel a
été implémenté en transférant les données des solveurs fluides et
structurels sur le surface des lames. Des efforts de recherche très
limités ont été consacrés aux méthodes fortement couplées pour
analyser les comportements hydrodynamiques des hélices sous-
marines. Kuo et Vorus [9] ont utilisé une telle méthode pour prédire la
contrainte dynamique des pales d'hélice. Les deux ensembles
d'équations fluides et structurelles ainsi que les termes de couplage dus
au FSI bidirectionnel ont été déterminés et combinés en un seul
ensemble d'équations, conduisant à une stratégie numérique plus
robuste et efficace. Une telle méthode fortement couplée sera utilisée
dans ce travail pour analyser les problèmes FSI de l'hélice. De plus, il
convient de noter que les méthodes couplées CFD-FEM sont
raisonnablement précises pour étudier les problèmes FSI des hélices
flexibles, mais l'effort de calcul requis par de telles méthodes est
considérable. Au contraire, des solutions précises peuvent être
obtenues en utilisant les méthodes couplées BEM et FEM avec un
effort de calcul
L'objectif relativement
de ce travail faible. une méthode numérique efficace
est de développer
pour prédire les réponses hydroélastiques des hélices marines dans le
sillage des sous-marins. Un ensemble de conditions aux limites Fig. 1. Le système de coordonnées fixes de l'hélice qui tourne avec l'hélice est illustré.
cinématiques raisonnablement précises sur l'interface fluide-structure a été
dérivé en utilisant les conditions de non-pénétration. Une méthode de
panneau dépendant du temps a été appliquée pour calculer les forces où ̈u, ̇u et u sont les vecteurs contenant respectivement les composantes
instables des hélices rigides, et la méthode des éléments finis combinée à nodales d'accélération, de vitesse et de déplacement. M et K sont les
une méthode de panneau dépendant de la fréquence a été utilisée pour matrices globales de masse et de rigidité assemblées par∫
formuler le modèle de structure de l'hélice et des matrices d'amortissement les matrices de chaque élément, données par : Me = Ve sNTNdV et
et de masse ajoutées en raison de le fluide. Une technique de réduction
modale combinée à la méthode Wilson a été utilisée pour déterminer les Ke=
réponses structurelles de l'hélice, ce qui permet de surmonter la faible ∫ Ve BTDBdV , où N, B = ˆDN et D sont la fonction de forme
efficacité numérique due aux matrices asymétriques ajoutées. Les effets
des fréquences d'excitation sur le problème FSI peuvent être explicitement matrice de tion, matrice de déformation-déplacement et matrice
pris en compte dans la méthode, qui fournit une grande compréhension constitutive du matériau, respectivement. Le modèle d'amortissement de
physique des interactions fluide-structure des hélices flexibles. Les Rayleigh est utilisé pour prendre en compte l'effet d'amortissement des
interactions fluide-structure unidirectionnelles et bidirectionnelles ont été pales, conduisant à la matrice d'amortissement globale C = M + ˇK, et et
prises en compte dans la présente analyse. Les résultats de la présente ˇ sont deux coefficients d'amortissement. C ̋ ̇u, K ̋u et F ̋ représentent les
méthode ont été validés par comparaison avec les solutions disponibles forces centrifuges et de Coriolis, qui sont négligées dans la présente
dans la littérature et celles obtenues à partir de logiciels commerciaux. analyse car la vitesse de rotation de l'hélice n'est pas grande. Fw est le
vecteur contenant les forces hydrodynamiques générées par l’afflux
spatial non uniforme après les coques.

2. Modèle mathématique 2.2. Gouverner les équations du fluide

2.1. Équations régissant les pales Une méthode de panel est utilisée pour déterminer le vecteur force
hydrodynamique Fw. L'écoulement est supposé incompressible, non
Un système de coordonnées cartésiennes o-xyz est établi pour analyser visqueux, et la vitesse totale Vtotal est exprimée par la somme du sillage
le mouvement de la pale, qui tourne avec l'hélice, comme le montre la d'afflux non uniforme V (obtenu par CFD comme le montre la Fig. 17) et
figure 1. L'axe y positif coïncide avec la ligne génératrice de la pale du débit potentiel de perturbation :
d'indexation de l'hélice, et l'axe x- L'axe est parallèle à l'axe central de V total = V +∇ (2)
l'arbre, positif ensuite. Le moyeu de l'hélice n'est pas pris en compte
dans l'analyse et les pales sont supposées fixées à l'emplanture. Les où est le potentiel de vitesse de perturbation correspondant au champ
pales de l'hélice sont constituées d'un matériau élastique isotrope, d'écoulement induit par l'hélice. Il est à noter que le champ
homogène et linéaire, et les propriétés du matériau sont données sous d'écoulement induit par l'hélice est supposé incompressible, non
la forme : coefficient de Poisson, module de Young E et densité s. visqueux et irrotationnel. Selon la théorie de l’écoulement potentiel et
La méthode des éléments finis est utilisée pour construire le modèle les conditions aux limites cinématiques à la surface de la pale, le
structurel des pales de l’hélice, qui sont spatialement discrétisées par problème des valeurs limites pour résoudre le potentiel de perturbation
des éléments isoparamétriques linéaires 3D. Chaque élément comporte peut s’écrire comme suit :
huit nœuds avec trois degrés de liberté de translation à chaque nœud. ∇2 = 0 (dans le domaine fluide ̋)
(3)
Basées sur la méthode des éléments finis, les équations discrétisées du
mouvement d'une lame élastique peuvent être obtenues en utilisant − ∂ ∂n = V · n − ∂ı∂ t · n + ((ı · ∇) V 0) · n (sur la
l'équation de Lagrange, donnée sous la forme : surface de la lame)
(4)
w (Rwake , t) = ( Rre , t − t′ )
M ̈u + (C ̋ + C) ̇u + (K ̋ + K )u = F ̋ + Fw (1) (état Morino Kutta) (5)

https://sider.ai/pdf-translate/E0G5UEKEG8A?fid=0 3/17
10/22/24, 7:18 PM PDF translate: -----A numerical method for predicting the hydroelastic response ofmarine propeller.pdf | Sider

J. Li et coll. / Recherche océanique appliquée 74 (2018) 188-204 191

où n est le vecteur normal unitaire extérieur. ı est le vecteur

déplacement des points nodaux sur la surface de la pale. V0 est la partie
uniforme du flux entrant. w (Rwake , t) représente le saut potentiel à
travers les feuilles de sillage, et (Rre , t − t′) désigne le saut potentiel à
travers la surface de la pale au niveau du bord de fuite, qui est égal au
potentiel du côté supérieur (côté aspiration ) moins celui du côté
inférieur (côté pression). t' est le temps nécessaire au fluide pour
parcourir la surface du sillage depuis le bord de fuite de la pale Rre
jusqu'au point de sillage Rwake. La condition limite cinématique (4)
(voir en annexe A) exige que l'écoulement soit tangent à la surface de la
pale oscillante. Cela conduit à une condition aux limites de type
Neumann. Les lecteurs peuvent se référer à la Réf. [17] pour plus de
détails sur la condition Kutta de Morino définie par l'équation. (5).
La distribution du potentiel de perturbation induite par l'hélice peut être
calculée en considérant des équations aux limites données. Par la suite,
les vitesses du champ d'écoulement peuvent être obtenues par l'équation.
(2). La pression du fluide peut être obtenue en utilisant l'équation de
Bernoulli [12,13], donnée par :

(1 ) Fig. 2. Discrétisation par méthode panel des surfaces de pale et de sillage.

2 V 2 − 12 V 2total
p = p0 + (6)
−∂∂t
2.3. Carter d'hélice rigide

où p0 est la pression hydrostatique absolue au point de calcul. est la densité du Les équations régissant le potentiel de perturbation dû à l'hélice rigide
fluide.
tournant dans un sillage non uniforme peuvent s'écrire comme suit :
Dans la présente analyse, la pression hydrostatique absolue n'est pas
important. Lorsque la vitesse de perturbation induite par l'hélice ∇ et ∇2 r = 0 (dans le domaine fluide ̋)
l'amplitude des termes non uniformes dans le sillage de la coque sont de (9)
petits termes d'ordre supérieur dans l'équation. (6) peut être négligé.
Dans de tels cas, l'équation. (6) peut être réécrit comme suit : − ∂ r∂n = V · n (sur la surface de la
(10)
pale)
( ) ( Rre , t − )
∂∂ t + V · ∇ + 12 (∇ )2 w,r (Réveil, t) = r t′ (état Morino Kutta) (11)
p=− (7)
où w,r (Rwake , t) représente le saut potentiel à travers les feuilles de
( ) sillage, et r (Rre , t − t′) désigne le saut potentiel à travers la surface de la
∂∂t+V0· pale au niveau du bord de fuite. t' est le temps nécessaire au fluide pour
p=−
∇
(8) se déplacer le long de la surface du sillage depuis le bord de fuite de la
pale Rre jusqu'au point de sillage Rwake. Le potentiel de vitesse de
perturbation r est régi par l'équation de Laplace. (9), et selon l'identité de
Notez l'équation. (7) et l'équation. (8) sont respectivement la formule de Green, il peut être calculé par l'expression intégrale suivante :
Bernoulli non linéaire complète et la formule de Bernoulli linéaire. ∫ ) ∂ G (RP , RQ)
r ( RQ , t
La fluctuation de pression p peut être obtenue soit par l'équation. (7) ou 2 r (RP , t) = ds
l'équation. (8), et le modèle hydrodynamique peut être combiné avec le ∂nQ
modèle structurel décrit dans l’équation. (1). Le vecteur force ∫ ∂ r ( RQ , t ) SB

hydrodynamique Fw peut être assemblé par les forces hydrodynamiques (RP, )

sur chaque élément, données comme suit : −
∂nQG RQ ds

SB
Peu = ZT ∫ ) ∂ G (RP , RQ)
w,r ( RQ , t
∫ pds,
+ ds (12)
∂nQ
Logiciel
où Z est la matrice d'interpolation des forces. La
perturbation induite par l'hélice où SB et Sw désignent la surface de la pale et de la feuille de sillage, respectivement
potentiel peut être
exprimé par : = r + v (voir Annexe A), et donc Fw peut être représenté tivement. Le modèle de sillage prescrit est considéré ici. L'angle de
par : Fr + Fv. Fr contient les forces hydrodynamiques dues à pas du sillage est supposé égal à l'angle de pas de la pale et la
l’écoulement non uniforme, qui est la force motrice du système de contraction radiale n'est pas prise en compte. L'indice Q représente le
couplage fluide-structure. Fv représente le vecteur de force point variable dans l'intégration, et l'indice P correspond au point de
hydrodynamique dû à l'interaction du fluide et de l'hélice, qui aboutit contrôle sur le panneau. RP et RQ sont les vecteurs de position des
aux matrices de masse ajoutée et d'amortissement induites par points sur la surface de la pale et de la feuille de sillage,
l'écoulement. Dans ce qui suit, la formulation détaillée et la mise en respectivement. G ( RP , RQ ) est la fonction de Green du Laplace
œuvre des deux parties des forces fluides sont présentées. Afin équation dans un domaine fluide 3D illimité, qui est donné par G ( RP ,
d'améliorer l'efficacité du calcul, une méthode de panneau dépendant du RQ ) = 1/|RP(RQ,
− RQt |. nQ est le vecteur normal unitaire extérieur,
) représente le saut potentiel à travers le sillage
temps est appliquée pour calculer la pression hydrodynamique Fr due à et r,w
la rotation des pales rigides, tandis qu'une méthode de panneau feuilles. La pale discrétisée et la surface de sillage sont illustrées à la Fig. 2.
dépendant de la fréquence est utilisée pour déterminer la pression L'intervalle radial du moyeu à la pointe est divisé en Nr sous-intervalles, et
hydrodynamique Fv due à la structure. déformation des lames l'intervalle dans le sens de la corde du bord d'attaque au bord d'attaque.
élastiques.

https://sider.ai/pdf-translate/E0G5UEKEG8A?fid=0 4/17
10/22/24, 7:18 PM PDF translate: -----A numerical method for predicting the hydroelastic response ofmarine propeller.pdf | Sider

192 J. Li et coll. / Recherche océanique appliquée 74 (2018) 188-204

problème, {r }
le bord de fuite est décomposé en Nc sous-intervalles. Les forces des |t est déterminé par l’équation. (13), et par la suite, le
dipôles et des sources sont supposées uniformes sur le panneau. les forces des dipôles et des sources à tout instant t peuvent être obtenues. Il
L'intervalle de sillage hélicoïdal depuis le bord de fuite jusqu'à une existe deux manières de calculer la pression instable, qui
distance suffisante en aval est divisé en Nw sous-intervalles. est régi par la formule de Bernoulli non linéaire, à savoir l'équation. (7).
Remplacement des équations. (10) et (11) dans l’équation. (12), le système L’une consiste à utiliser une opération différentielle sur la surface de la lame
d’équations discrétisé suivant peut être obtenu : pour calculer ∇ [18], et l’autre consiste à utiliser une méthode d’intégration
directe qui nécessite des traitements numériques robustes de six intégrales
Cr {r } |t = G { V · } |t + Cm {r }
singulières dans l’équation. (13); voir Réf. [19-23]. La méthode d'intégration
n Nw−1∑ |t−mt (13)
directe est utilisée ici. Ce faisant, ∇ r peut être calculé en prenant la lième
m=1
composante de la dérivée spatiale de l'équation. (12) comme suit :
où les éléments de la ième ligne et de la jième colonne dans les matrices ∫ ( ∂ G (RP, ))
de coefficients de panneau Cr, G et Cm sont les suivants : 2 ∂ r (RP , t) = r ( RQ , t )∂ RQ
ds
⎛ ⎞
∂ G ( Ri , ) ∂xl ∂xl ∂nQ
∫ ∫ ∂ G ( Ri , )
Rj ⎜⎜⎝∫ Rj1 ⎟⎟⎠ ∫ ∂ r ( RQ , t )SB (RP, )
Cr (je, j) =2 Je (je, j) −
∂njds − ∫s ds ∂nQ ∂G RQ
∂nj1 −
∂ xl ds
sj j1 1 SB

⎛ ⎞ ( ∂ G (RP, ))
∂ G ( Ri , ) ∫
w,r ( RQ , t
⎜⎜⎝∫Rj1 )∂ RQ
⎟⎟⎠ +
∂xl ∂nQ
ds (14)
+
∫s ∂nj1
ds
Logiciel
j1 2
Remplacement des équations. (10) et (11) dans l’équation. (14), l’équation
⎛ ⎞ discrétisée peut s’écrire :
∂ G ( Ri , ) { }
⎜⎜⎝∫Rjm+1 ⎟⎟⎠ ∂ r
|t =
{ r} |t + Gxl{V · n} |t Cxlm{ r} |t−mt
(15)
Cm (je, j) = ∂ xl
∫s ∂njm+1
ds
Cxlr +Nw −1∑ m=1
jm+1
⎞
1 où les éléments de la ième ligne et de la jième colonne dans les matrices
⎛ ∂ G ( Ri , ) de coefficients de panneau Cxlr , Gxl et Cxlm sont donnés comme suit :
⎜⎜⎝∫Rjm+1 ⎛
⎟⎟⎠ ∫∫ ( ∂ G ( Ri , ) )
−
∫s ∂njm+1
ds
Cxlr (je, j) = 1 2 ⎜⎝ ∂ Rj
ds
jm+1 2
∂xl ∂nj
sj
∫∫ ⎛ ⎞
1 ⎛ ∂ G ( Ri , ) ⎞
G (je, j) = |Ri − Rj| ds; je, j = 1, · · ·, N ⎜⎜⎝∫∂
+ ⎝ Rj1 ⎠ds ⎟⎟⎠
sj ∫ s ∂xl ∂nj1
j1
et { r ⎞ ⎞
1
⎛ ⎛ ∂ G ( Ri , ) ⎞
} |t = [ r (R1 ) , r (R2 ) , · · ·, r (RN )] T|t
⎜⎜⎝∫∂
− ⎝ Rj1 ⎠ds ⎟⎟⎠
⎟⎟⎠
{ Vn} ∫ s ∂xl ∂nj1
|t = [V (R1 ) · n (R1 ) , V (R2 ) · n (R2 ) · · ·, V (RN ) · n (RN )] T|t
j1 2

t = ̋ , N = 2KNr Nc
⎛⎛ ) ⎞ ⎞
⎛ ∂ G ( Ri ,
⎜⎜⎝∫∂
Cxlm (je, j) = 12 ⎜⎜⎝ ⎝ Rjm+1 ⎠ds ⎟⎟⎠
où I est une matrice d'identité N × N ; sj est l'aire du jème panneau ; ∫ s ∂xl ∂njm+1
s jj′ est l'aire du j'ième panneau sur la jième bande de vortex. nj est jm+1 1
le vecteur normal unitaire vers l'extérieur du jème panneau sur la ⎞ ⎞
⎛
pale, et njj′ est le vecteur normal unitaire du j'ème panneau sur la ⎛ ∂ G ( Ri , ) ⎞
jème bande de vortex. Ri représente le vecteur position d'un point ⎜⎜⎝∫∂ ⎝ Rjm+1 ⎟⎟⎠
⎟⎟⎠
⎠ds
de contrôle sur le ième panneau, et Rjj′ désigne le vecteur position −
∫ s ∂xl ∂njm+1
d'un point de contrôle sur le j'ième panneau. jm+1
panneau de la jième bande de vortex. L'indice 1 indique que l'équation 2

est calculée uniquement pour les panneaux situés sur le bord de fuite de ∫∫ ∂ G ( Ri , )
la surface d'aspiration, tandis que l'indice 2 indique que le calcul est Gxl (je, j) = 12 Rj
effectué uniquement pour les panneaux situés sur le bord de fuite de la ∂ xlds ; je, j = 1, · · ·, N
surface de pression. ̋ est la vitesse de rotation de l’hélice ; est l'angle sj
correspondant de chaque élément du panneau sur le sillage de l'hélice. N
est le nombre total de panneaux. K est le nombre de pales. et { } [ (RN )
Le processus de résolution est résumé comme suit : premièrement, la ∂ r ∂ r (R1 ) , ∂ (R2 ) , · · ·, ∂
résistance des dipôles et des sources est résolue en considérant le problème ∂ xl | = ∂ xl r∂xl r∂xl
]T |
constant des hélices ; deuxièmement, afin de résoudre l'instabilité t t

https://sider.ai/pdf-translate/E0G5UEKEG8A?fid=0 5/17
10/22/24, 7:18 PM PDF translate: -----A numerical method for predicting the hydroelastic response ofmarine propeller.pdf | Sider

J. Li et coll. / Recherche océanique appliquée 74 (2018) 188-204 193

Re (ik (RQ ) e−ict)

Il est à noter que ∂ r∂ t peut être obtenu en utilisant une différence convient, et v (RP , t) peut être exprimé comme suit :
v (RP , t) = Re ( k,v (RP ) e−ikt ) , où Re est la partie réelle du
finie entre des pas de temps adjacents. En appliquant l'équation. (7),
pr (t) peut être calculé à tout instant t et la force dynamique Fr (t) est expression. L'expression discrétisée de l'équation. (21)
calculée par :
{ Pr(t)} peut
k,v } s'écrire
= Vous {: je
{ + E { })
Fr (t) = Z1T s (16) (je sais } jek (22)

dans lequel
Par souci de commodité, le maillage de la méthode des éléments finis
coïncide avec celui de la méthode des panneaux sur la surface de la U = C−1k,vG
pale. Z1 est la matrice d'interpolation de force qui distribue la force de
pression du ième panneau aux nœuds structurels du modèle d'éléments ⎛ ⎞
finis. s est une matrice diagonale donnée par s = diag { s1, s2, . . ., sn } ∫∫ ∂ G ( Ri , ) ∂ G ( Ri )
∫∫
et si est la ième zone du panneau. En appliquant la transformée de Rj ⎜⎜⎝Nw∑
eikt , Rjj′
ds
⎟⎟⎠
Fourier pour Fr (t), la force dynamique Fr (k) dans le domaine fréquentiel
Ck,v (i, j) = 2 I (i, j) −
∂njds − jj'

peut être calculée. k est la fréquence.

j′ =1 ∂njj′
sj s j j′
1
⎛ ⎞
2.4. Carter d'hélice élastique ∂ G ( Ri )
∫∫
+
⎜⎜⎝Nw∑
eikt , Rjj′
jj' ds
⎟⎟⎠
Les équations régissant le potentiel de perturbation dû aux hélices j′ =1 ∂njj′
flexibles vibrant dans un écoulement uniforme peuvent s'écrire comme s j j′

suit : 2

⎡ ( ) ( ) ( ·n) ⎤
∂V 0 ∂ ∣∣∣∣ ∂V 0 ∂ ∣∣∣∣ ∂V
⎢⎢⎢⎢⎢⎢⎣
x·n 1 y·n 1 0 ∂z ∣∣∣∣ 1 ⎥⎥⎥⎥⎥⎥⎦
E= ...
( ) ( ) ( ·n)
∂V 0 ∂ ∣∣∣∣N ∂V 0 ∂ ∣∣∣∣N ∂V
x·n y·n 0 ∂z ∣∣∣∣N

{ je } = [ ık (R1 ) , ık (R2 ) , · · ·, ık ]T
(RN )
∇2 v = 0 (dans fluide, domaine ̋)
(17) { je = [ ık (R1 ) · n (R1 ) , ık (R2 ) · n (R2 ) , · · ·, ık (RN ) · ]T
sais } n (RN )
− ∂ v∂n = − ∂ı∂ t · n + ((ı · ∇) V 0) · n (sur la
surface de la lame) (18) { k, v } = [ k,v (R1 ) , k,v (R2 ) , · · ·, k,v (RN ]T
( Rre , t − ) )
w,v (Réveil, t) = v t′ (état Morino Kutta) (19)
où t jj 'représente le temps requis pour que le vortex se déplace
où w,v (Rwake , t) représente le saut potentiel à travers les feuilles de du jème panneau sur le bord de fuite jusqu'au j'ème panneau sur
sillage, et v (Rre , t − t′) désigne le saut potentiel à travers la surface de la
la jème bande de vortex.
pale au niveau du bord de fuite.
Le potentiel de vitesse de perturbation v est régi par le coefficient de Laplace La formule linéaire de Bernoulli Eq. (8) est appliqué pour
équation. En utilisant l’identité de Green, v peut être calculé par l’expression déterminer la pression due aux petites déformations des
intégrale suivante : pales. La
} même formule pour calculer ∇ v, donnée par : {
∂k,v = Ux l(ik {ıkn} + E {ık})
∫ ∫ ∂xl (23)
∂ v (RQ, t)
2 v (RP , t) =
v (RQ , t) ∂ G (RP , ds- ∂nQ G (RP , RQ )ds
RQ ) ∂nQ
où
SB
SB

∫ Uxl = Cxlk,vU + Gxl

+
w,v (RQ , t) ∂ G (RP , ds
RQ ) ∂nQ (20) ⎛
∫∫ ( ∂ G ( Ri , ) )
Logiciel
12 ∂ Rj
où SB, Sw, RP, RQ, G ( RP , RQ ) Cxlk,v (je, j) = ⎜⎝ ds
et nQ suivent la(RQ,
mêmet moyenne- ∂xl ∂nj
)
sj
opérations telles que définies dans la section précédente. w,v représente
⎛ ⎞
le saut potentiel à travers les feuilles de réveil. En
∫∫ ⎛ ∂ G ( Ri , ) ⎞
remplaçant l'équation. (18) dans l’équation. (20), on obtient
∫ ∫ +
⎜⎜⎝Nw∑eikt jj' ∂
⎝ Rjj′ ⎠ds ⎟⎟⎠
2 v (RP, t) −
v (RQ , t) ∂ G (RP , ds-
w,v (RQ , t) ∂ G (RP , ds
j′ =1 ∂xl
∂njj′
s jj
RQ ) ∂nQ RQ ) ∂nQ
1
∫ SB Sw
(21) ⎛ ⎞ ⎞
( · nQ + ı (RQ , t) · ∇ ) ) · nQ ⎛ ( Ri, ⎞
− ∂ı (RQ , (( V0 ) ∫∫ ⎝ ∂G Rjj′)
=
t) ∂t
G (RP, RQ) ds
⎜⎜⎝Nw∑ eikt jj' ∂
⎠ds ⎟⎟⎠
⎟⎟⎠
−
SB j′ =1 ∂xl
∂njj′
s jj
La fréquence de la force d’excitation pour le système dynamique linéaire 2
est obtenue dans la section 2.3. Nous supposons que ı ( RQ , t ) =

https://sider.ai/pdf-translate/E0G5UEKEG8A?fid=0 6/17
10/22/24, 7:18 PM PDF translate: -----A numerical method for predicting the hydroelastic response ofmarine propeller.pdf | Sider

194 J. Li et coll. / Recherche océanique appliquée 74 (2018) 188-204

⎧ 3∑ ⎫
{ } [ ⎪⎪⎪⎪⎨
k2Im(U)Z1 + l=1 Vxl0 Im(Uxl ⎪⎪⎪⎪⎬
(RN ) CW (k) = k Z1T s )EZ2+
∂k,v ∂k,v (R1 ) , · · ·, ∂ ]T [
⎪⎪⎪⎪⎩ ] ⎪⎪⎪⎪⎭
(28)
∂xl = ∂xl 3∑
k,v ∂ xl l=1 Vxl0
k −Re(U)EZ2 +
Re(Uxl )Z1
Remplacement des équations. (22) et (23) dans l’équation. (8), on obtient la
pression dynamique comme : où Im désigne la partie imaginaire d’une expression.

⎛ ⎛ ( ) ⎞⎞
−UE { je } + 3∑
l=1 V xl0 3. Algorithme numérique pour les équations gouvernantes

{ pv(k)} ⎜⎜⎜⎜⎝−⎜⎜⎜⎜⎝ik Uxl{ıkn} ⎟⎟⎟⎟⎠

⎟⎟⎟⎟⎠
Pour résoudre l’équation. (26), les méthodes d'intégration temporelle
= Ré
e−ict } + 3∑
+k2U { l=1 V xl0 directe, telles que les méthodes de Wilson et de Newmark, peuvent être
jekn utilisées. Cependant, les matrices de l'équation. (26) sont de grandes
Uxl E {ık} dimensions, et la ressource de calcul requise par ces méthodes pour
(24) résoudre l'équation. (26) directement est considérable. Il est à noter que la
méthode de superposition modale est une méthode efficace pour les
analyses de dynamique des structures. Malheureusement, les matrices de
où V0xl (Ri) sont les composantes de V0 dans la direction xl sur le point masse ajoutée et d'amortissement dans l'équation. (26) sont des matrices
Ri, données par V xl0 = diag { V0xl (R1 ) , · · ·, V0xl (RN ) } .La pression asymétriques. Par conséquent, l’équation. (26) ne peut pas être calculé
directement par la méthode de superposition modale. Une technique de
vecteur dans l'équation. (24) est exprimé en termes de vecteurs de regroupement de type HRZ [12] peut être utilisée, qui remplace MW (k) et
déplacement inconnus { ık } et { ıkn } , qui sont calculés aux centroïdes de CW (k) par les matrices diagonales ajoutées. Cependant, cette méthode
les panneaux de lames. Il faut exprimer le vecteur pression en fonction des peut s’avérer imprécise pour prédire les modes de vibration des pales de
déplacements nodaux
En supposant { } des pales. { je l’hélice. Ceci sera discuté dans les sections suivantes. Ainsi, une technique
ık = Z2uk et sais } = Z1u k, où Z1 et Z2 sont modale réduite est utilisée dans l’étude. Les premiers vecteurs mode
matrices de fonctions de forme et uk est le vecteur de déplacement nodal, la humide Nmode de l’hélice élastique peuvent être combinés en une matrice
force dynamique peut être calculée comme suit : modale, donnée par : k. Le vecteur déplacement peut être exprimé en
termes
des vecteurs modaux ik et coordonnées généralisées qik (t), étant donné
⎛ ⎛ ( ) ⎞⎞
3∑ comme : uk =
je
-UEZ2uk + l=1 V xl0 Uxl Nmode∑ je=1qik k. Ensuite, l'équation. (26) peut s’écrire :
⎜⎜⎜⎜⎝ ⎜⎜⎜⎜⎝ik ⎟⎟⎟⎟⎠
⎟⎟⎟⎟⎠ (t)
Déclenchement Z1uk
Fv (k) = Ré
du Z1T 3∑ ˆMk ̈qk (t) + ˆCk ̇qk (t) + ˆK kqk (t) = ˆF (29)
+k2UZ1uk + l=1 V xl0 Uxl
EZ2uk où
(25) k={ }
1 2k , · · ·, Nmode
k,
k

Si nous remplaçons la condition aux limites cinématique définie dans l’Eq. qk (t) = [ q1k (t) , q2k (t) , · · ·, (t) ]T
(18) avec celui de l’équation. (10), la méthode du panneau dépendant de la qNmodek
fréquence peut être utilisée pour calculer la pression hydrodynamique Fr
due à la rotation de la pale rigide.
ˆMk = kT [M + MW (k)] k

2.5. Matrices de masse ajoutée et d'amortissement ˆCk = kT [C + CW (k)] k

Afin de calculer le potentiel de perturbation vibratoire k,v et la force ˆK k = kTK k

hydroélastique Fv (k) via la méthode des panneaux, le vecteur de
déplacement nodal tridimensionnel uk est nécessaire à partir du modèle
d'éléments finis. À son tour, l’analyse par éléments finis des réponses ˆF = kTF

structurelles (déplacement uk, vitesse ̇uk et accélération ̈uk) de l’hélice

nécessite les vecteurs de force hydrodynamiques Fr (k) et Fv (k). Dans Maintenant, le nombre total de degrés de liberté impliqués dans
la présente analyse, le problème d’interaction fluide-structure est résolu l’équation. (26) est réduit au mode N, et l'équation. (29) peut être
en plaçant la force hydroélastique sur le côté gauche de l’équation du résolu très efficacement en utilisant la méthode Wilson. Une fois
mouvement de la pale, donnée par : qk(t) déterminé, uk peut être calculé
je
par : uk =concernant
Nmode∑
je=1qik k. Un organigramme le com-
(t)
[M + MW (k)] ̈uk + [C + CW (k)] ̇uk + K uk = Fr (k) (26) La présentation de la présente méthodologie est donnée sur la figure 3.

où MW (k) et CW (k) sont respectivement la matrice de masse ajoutée et la 4. Résultats et discussions

matrice d'amortissement ajouté, données sous la forme
4.1. Validation
⎧ 3∑ ⎫
⎪⎪⎪⎪⎨
k2Re(U)Z1 + l=1 Vxl0 ⎪⎪⎪⎪⎬Trois types d'hélices, à savoir l'hélice 4119, l'hélice 4381 et l'hélice 4382,
MW (k) = k2 Z1T s [ Re(Uxl )EZ2 ] sont considérés ici pour tester l'exactitude de la présente méthode de
⎪⎪⎪⎪⎩ 3∑
⎪⎪⎪⎪⎭(27)
prévision des performances en eau libre des hélices. Une technique de
−k −Im(U)EZ2 + l=1 Vxl0 progression géométrique est utilisée dans le sens de l'envergure, de la
Im(Uxl )Z1 racine à la pointe, et dans le sens de la longueur de l'arc.

https://sider.ai/pdf-translate/E0G5UEKEG8A?fid=0 7/17
10/22/24, 7:18 PM PDF translate: -----A numerical method for predicting the hydroelastic response ofmarine propeller.pdf | Sider

J. Li et coll. / Recherche océanique appliquée 74 (2018) 188-204 195

Fig. 3. Organigramme de la méthodologie d'analyse des forces de réaction instables des hélices rigides et des hélices élastiques en considérant une interaction fluide-structure unidirectionnelle et
bidirectionnelle.

Fig. 4. Discrétisation par panneau des surfaces de pale et de sillage. (a) : hélice 4119 ; (b) : hélice 4381 ; (c) : hélice 4382.

du bord d'attaque au bord de fuite, tandis qu'une segmentation semi- coefficient de poussée K héliceT = Thélice / n2D4 et coefficient de couple
cosinus au sillage proche combinée à une segmentation uniforme au sillage K propellerQ = Qpropeller / n2D5 par rapport au coefficient d'avance J =
lointain est appliquée dans la direction de la nappe tourbillonnaire de fuite ; V/nD pour les trois hélices sont représentés sur la Fig. 5. Thélice et
voir Fig. 4. Les variations des valeurs prédites et mesurées

https://sider.ai/pdf-translate/E0G5UEKEG8A?fid=0 8/17
10/22/24, 7:18 PM PDF translate: -----A numerical method for predicting the hydroelastic response ofmarine propeller.pdf | Sider

196 J. Li et coll. / Recherche océanique appliquée 74 (2018) 188-204

Fig. 5. Coefficients de poussée et de couple prévus et mesurés en fonction du rapport d'avance. (a)
: hélice 4119 ; (b) : hélice 4381 ; (c) : hélice 4382.

Fig. 6. Sillage de l'hélice HSP dans la direction (a) : x. (b) : direction y. (c) : direction z.

Qhélice représente respectivement la poussée et le couple des hélices, n

désigne la vitesse de rotation de l'hélice, V est la vitesse d'avancement et D
est le diamètre de l'hélice. Sur la figure 5, les courbes « exp » représentent les résultats mentaux sont plutôt bons. Cela indique que la méthode
les résultats expérimentaux obtenus par Kerwin et Lee [24] et Greeley et différentielle et la méthode d’intégration directe conviennent toutes
Kerwin [25] ; « diff » désigne les résultats obtenus par la méthode deux pour prédire les forces fluides constantes. De plus, si l'influence de
différentielle ; 'inte' représente les solutions prédites par la méthode la viscosité due à la couche limite et à la séparation des écoulements
d'intégration directe. Lorsque l'effet visqueux dû à la force visqueuse fs =
près du bord de fuite est prise en compte, les coefficients de poussée
12 CD | V total | V total si (CD = 0,0042 est le coefficient de frottement)
seront diminués et les coefficients de couple seront augmentés. Ces
est pris en compte, les solutions correspondantes sont étiquetées par un
résultats devraient être cohérents avec la théorie du panneau de l’hélice.
symbole « V », comme le montre la figure 5. On observe sur la figure 5
que l'accord entre les présentes solutions et l'expérience
Les coefficients de force obtenus par la méthode d'intégration directe
sont plus petits que ceux obtenus par la méthode différentielle.

https://sider.ai/pdf-translate/E0G5UEKEG8A?fid=0 9/17
10/22/24, 7:18 PM PDF translate: -----A numerical method for predicting the hydroelastic response ofmarine propeller.pdf | Sider

J. Li et coll. / Recherche océanique appliquée 74 (2018) 188-204 197

Fig. 8. Comparaisons des coefficients de force instable pour la lame de clé prédits par le
modèle actuel et par les résultats publiés dans le domaine fréquentiel. (a) : coefficients de
poussée ; (b) : coefficients de couple.
Fig. 7. Comparaisons des coefficients de force instable pour la lame de clé prédits par le
modèle actuel et par les résultats publiés dans le domaine temporel. (a) : coefficients de
poussée ; (b) : coefficients de couple. problème (voir équations (9) – (11)), la méthode différentielle [18] et la
méthode d'intégration directe sont utilisées pour calculer ∇ . Cela conduit
à un total de quatre méthodes de résolution, étiquetées comme suit :
Les coefficients de poussée instable et de couple prédits obtenus à partir
TIME DIFF, TIME INTE, FRE DIFF et FRE INTE. Dans la présente
de la présente méthode sont comparés aux résultats de Hoshino pour la
analyse, les méthodes TIME DIFF, TIME INTE et FRE INTE sont
pale de clé (voir sur la figure 6) de l'hélice HSP (hélice hautement
utilisées pour le calcul. De plus, afin d'obtenir les solutions basées sur
asymétrique). Le sillage non uniforme au niveau du disque de l'hélice est
différentes formules de Bernoulli, la méthode TIME INTE est
basé sur le modèle du 22e Comité de propulsion de l'ITTC, et la
décomposée en trois types de méthodes, données comme suit : TIME
géométrie de l'hélice est donnée sur la figure 6. Les coefficients de vitesse
INTE 1, TIME INTE 2 et TIME INTE 3. Dans nos cinq codes,
de sillage sur la figure 6 sont définis par ( 1 − Vx /V, Vy /V, Vz /V ) .
l'intégralité La formule de Bernoulli non linéaire (7) est appliquée dans
La segmentation du maillage suit la même configuration que celle utilisée TIME DIFF et TIME INTE 1, et la formule de Bernoulli p = − ( ∂ ∂ t + V
dans le cas de performances stables, sauf qu'une segmentation uniforme est 0 · ∇ + 12 |∇ |2 ) est utilisée dans TIME INTE 2
utilisée le long de la feuille de sillage. L’hélice est ici considérée comme
et TIME INTE 3. FRE INTE est basé sur la formule linéaire de Bernoulli (8).
rigide. Lorsque la ligne génératrice d'hélice de la lame de clé tourne dans le
sens des aiguilles d'une montre dans le sillage, le coefficient de poussée
instable K keybladeT = Tkeyblade / n2D4 et le coefficient de couple Concernant les résultats de TIME INTE 2 et TIME INTE 1, on observe que
la vitesse V dans la formule de Bernoulli non linéaire complète (7) peut
K keybladeQ = Qkeyblade / n2D5 pour la lame de clé changera être remplacée par V0 . De plus, comme le montrent les résultats de TIME
périodiquement. Par conséquent, les fréquences d’oscillation des INTE 3 et TIME INTE 2, la formule linéaire de Bernoulli (8) peut être
coefficients de force instationnaire incluront la fréquence de passage axial utilisée pour obtenir des solutions raisonnablement précises. Comparé aux
(APF Hz) et les harmoniques d’ordre élevé. Les coefficients de force dans résultats de TIME INTE 3 et FRE INTE, on constate que le résultat calculé
les domaines temporel et fréquentiel sont présentés sur les Fig. 7 et Fig. 8. par FRE INTE est la meilleure approximation pour la prédiction du
On observe sur la Fig. 7 que le changement des coefficients de force coefficient de poussée instationnaire et du coefficient de couple. Cette
instable se fait en un tour et que l'angle est entre la ligne du générateur de méthode sera utilisée dans la suite pour calculer les matrices de masse
l'hélice. de la lame de clé et la direction positive de l'axe y. Les résultats de ajoutée et d’amortissement. Les résultats de la méthode TIME DIFF
l'Institut de recherche et de conception pour la construction navale confirment en outre la validité de la méthode FRE INTE et des différentes
(ICEPRONAV) et du Samsung Ship Model Basin (SSMB) correspondent méthodes TIME INTE. Une grande erreur des coefficients instationnaires
respectivement aux limites supérieure et inférieure des résultats dans la plage de 250◦ à 300◦ est trouvée dans le cas de charges lourdes. La
numériques de HSP [26-28]. Pour calculer les forces instationnaires, la raison réside dans le fait que la condition de Kutta de Morino est utilisée
méthode du domaine temporel donnée dans la section 2.3 et la méthode du dans le calcul plutôt que la condition de Kutta, et qu'un modèle de sillage
domaine fréquentiel présentée dans la section 2.4 sont utilisées. Une fois prescrit est implicite ici ainsi que l'angle de tangage de
résolu à partir de la valeur limite

https://sider.ai/pdf-translate/E0G5UEKEG8A?fid=0 10/17
10/22/24, 7:18 PM PDF translate: -----A numerical method for predicting the hydroelastic response ofmarine propeller.pdf | Sider

198 J. Li et coll. / Recherche océanique appliquée 74 (2018) 188-204

Fig. 9. Discrétisation d'une pale de l'hélice 4119 et son modèle de calcul dans ANSYS.

Fig. 10. Comparaisons des résultats numériques prédits dans le vide par le modèle actuel et ANSYS
avec 40 éléments (dans le sens de l'envergure) × 40 (dans le sens de la corde) × 2 (épaisseur). (a) :
fréquences modales ; (b) : amplitudes de réaction ; (c) : amplitudes de déplacement.

le sillage est réglé égal à l'angle de pas de la pale plutôt qu'à l'angle de
pas hydrodynamique. Cependant, la figure 8 montre que presque tous
les résultats calculés par nos méthodes se situent entre les limites Le domaine fréquentiel peut être utilisé pour obtenir des solutions
supérieure et inférieure des résultats numériques du HSP. De plus, les raisonnablement précises pour les hélices élastiques.
coefficients de poussée et de couple instables prévus des cinq modèles Pour valider le modèle structurel, les résultats actuels calculés pour une
sont bien comparés aux résultats publiés dans la littérature. Cela pale de l'hélice 4119 sont comparés aux solutions numériques obtenues
confirme la validité des méthodes d'intégration directe dans les à partir du logiciel d'éléments finis ANSYS. La pale est fixée au pied et
domaines temporel et fréquentiel et des méthodes différentielles dans le la configuration géométrique de la pale est représentée sur la figure 9.
domaine temporel pour prédire les forces instationnaires. De plus, la L'effet d'amortissement structurel de la pale est négligé. La
formule de Bernoulli linéaire avec méthode d'intégration directe dans le segmentation du maillage sur la surface de la pale

https://sider.ai/pdf-translate/E0G5UEKEG8A?fid=0 11/17
10/22/24, 7:18 PM PDF translate: -----A numerical method for predicting the hydroelastic response ofmarine propeller.pdf | Sider

J. Li et coll. / Recherche océanique appliquée 74 (2018) 188-204 199

Fig. 11. Comparaisons des cinq premières formes modales dans le vide prédites par le modèle actuel et ANSYS avec 40 éléments (dans le sens de l'envergure) × 40 (dans le sens de la corde) × 2 (épaisseur).

Fig. 13. Fréquences fondamentales de deux hélices modèles à 7 pales à forte inclinaison.

les performances des hélices ainsi que les données expérimentales

rapportées par Tian et al. [29, 30] sont comparés sur la figure 13. Il ressort
de la figure 13 que les résultats numériques obtenus par la présente
méthode sont en très bon accord avec les données expérimentales.
Cette section est consacrée à la validation de l'analyse de l'interaction
Fig. 12. Maillage des éléments de frontière de l'hélice modèle à 7 pales à forte inclinaison. fluide-structure. Les matrices de masse ajoutée et d’amortissement ajouté
doivent être validées. Il convient de noter qu’une vérification directe des
éléments de ces deux matrices de couplage fluide-structure peut ne pas
est la même que celle utilisée dans le modèle de panneau, et une être simple. Cependant, la matrice de masse ajoutée peut être
segmentation uniforme est utilisée dans le sens de l’épaisseur. On observe indirectement validée en comparant les fréquences propres, les formes
à partir des Fig. 10 et 11 que les résultats calculés des fréquences propres modales et les réponses dynamiques de l'hélice dans un écoulement
et des formes modales de la pale sont bien comparés aux solutions de stationnaire.
Les fréquences propres, les formes modales et les réponses dynamiques
référence. sont calculées pour la même pale de l'hélice 4119 considérée dans la
Afin d’étudier les comportements dynamiques de la pale, une force validation FEM. Cependant, la lame est désormais immergée dans une eau
concentrée est appliquée à une extrémité de la pale dans la direction x, stationnaire. On peut déduire de l'Annexe B que les fréquences naturelles
comme le montre la Fig. 9. La force est définie comme : f = sin (2 f0t), où et les formes modales de la pale dérivées de l'analyse d'interaction fluide-
structure pour un écoulement stationnaire sont les mêmes que celles
la fréquence d'excitation est prise comme f0 = 1000 Hz. Les amplitudes
obtenues à partir d'une analyse d'interaction acoustique-solide, lorsque le
de déplacement dans trois directions en un même point sont données sur
nombre d'onde du son est défini. à zéro. Les résultats actuels calculés pour
la figure 10, et les forces de réaction mesurées au pied de la pale dans
l'interaction fluide-structure de la pale sont comparés aux solutions
trois directions sont également présentées. On observe que les résultats
obtenues à partir de l'analyse de l'interaction acoustique-structure en
présents sont en bon accord avec ceux obtenus par ANSYS.
utilisant ANSYS et Virtual. Acoustique du laboratoire. Comme le montrent
Les fréquences fondamentales de deux hélices modèles à 7 pales à forte les fig. 14 et 15, les présents résultats sont bien comparés à ces solutions de
inclinaison ont été calculées pour valider la précision de la présente référence. De plus, les fréquences propres déterminées par la méthode de
méthode. Les deux hélices en bronze et en matières plastiques ont les symétrisation de type HRZ sont également comparées aux présentes
mêmes configurations géométriques, comme le montre la Fig. 12. Le solutions. La technique de regroupement de type HRZ est obtenue en
diamètre des deux hélices est de 250 mm. Pour plus de détails sur les mettant à l'échelle les entrées diagonales des matrices cohérentes de masse
paramètres géométriques et matériels des hélices, le lecteur pourra se ajoutée et d'amortissement avec le rapport de la somme des entrées.
référer à Tian et al. [29,30]. La fréquence fondamentale calculée

https://sider.ai/pdf-translate/E0G5UEKEG8A?fid=0 12/17
10/22/24, 7:18 PM PDF translate: -----A numerical method for predicting the hydroelastic response ofmarine propeller.pdf | Sider

200 J. Li et coll. / Recherche océanique appliquée 74 (2018) 188-204

Tableau 1 Comparaisons des premiers 5èmes modes humides calculés

par ANSYS-Virtual.Lab, modèle actuel et modèle de type HRZ (Hz).

numéro de mode ANSYS-VIRTUALLAB modèle actuel De type HRZ

1 914.3 935.9 996.2

2 2135.3 2085.3 1663.3
3 2641.9 2760.2 2171.4
4 4066.9 4146.7 2672.1
5 4168.3 4217.4 2803.6

en utilisant ANSYS Fluent. Les équations de Navier-Stokes moyennées par

Reynolds sont résolues sur un système de grilles structurées en blocs avec un
domaine cylindrique circulaire. Les distances entre l'entrée et la proue, la
sortie vers la poupe et l'axe central de la coque jusqu'à la limite latérale sont
prises comme suit : L, 2L et 10Dm, où L et Dm désignent la longueur totale
et le diamètre maximal de la coque, comme indiqué. sur la figure 16,
respectivement. Les données géométriques de SUBOFF sont tirées de la
littérature [31]. La figure 16 montre les grilles sur la surface de la coque.
L'entrée et la limite latérale sont définies comme vitesse d'entrée, tandis
qu'une condition limite de pression de sortie est appliquée à la sortie. Le
modèle K standard et les fonctions de paroi standard sont utilisés pour la
fermeture des turbulences. Le y+ moyen des cellules adjacentes aux murs
sur la surface de la coque est d'environ 50.
La méthode combinée du panneau du domaine temps-fréquence et la
méthode FEM ont été utilisées pour calculer les réponses dynamiques de
l'hélice 4382 tournant dans le sillage d'un sous-marin SUBOFF, comme le
montre la figure 17. L'unité de la vitesse de sillage est m/s. Le diamètre
de l'hélice est de 250 mm et la vitesse angulaire de l'hélice est de 2746
tr/min. Les paramètres matériels de l'hélice sont donnés comme suit :
module d'Young E = 210 GPa, densité s = 7 800 kg/m3 et rapport de
Poisson = 0,3. La densité du fluide est = 1000 kg/m3 . L'amortissement
structurel de l'hélice est négligé. Pour l'hélice immergée dans le vide et
l'eau, les dix premières fréquences propres de l'hélice prédites par la
présente méthode sont répertoriées dans le tableau 2. Les forces de
réaction instables des hélices rigides et des hélices élastiques sont
déterminées en utilisant un fluide unidirectionnel et bidirectionnel. -
Fig. 14. Comparaisons des résultats numériques dans l'eau prédits par le modèle actuel et
analyses d'interactions de structures. Il est à noter que seule la force
ANSYS-Virtual.Lab avec 24 éléments (dans le sens de l'envergure) × 24 (dans le sens de la
corde) × 2 (épaisseur). (a) : fréquences du modèle ; (b) : amplitudes de déplacement. motrice Fr (k) du système est considérée dans la modélisation
unidirectionnelle, et les matrices ajoutées ne sont pas prises en compte.
Concernant la modélisation bidirectionnelle, les matrices ajoutées et la
qui contribuent au mouvement dans la même direction à la somme des force motrice Fr (k) sont prises en compte. Les résultats numériques pour
entrées diagonales qui contribuent au mouvement dans cette direction ; voir la comparaison de la lame rigide et de la lame élastique en considérant
Young [12]. Il ressort du tableau 1 que la méthode de type HRZ n'est pas l'interaction de structure fluide unidirectionnelle et bidirectionnelle sont
précise pour prédire les modes de vibration de la pale. présentés sur la figure 18. Il est à noter que BPF signifie la fréquence de
passage de la lame. De la figure 18 et du tableau 2, les conclusions
4.2. Analyse des forces instationnaires
suivantes peuvent être tirées :

Afin d’obtenir le sillage de SUBOFF au niveau du disque d’hélice sans l’hélice, 1) Les fréquences propres fondamentales de l'hélice 4382 dans le vide et dans
les équations régissant le champ d’écoulement sont résolues l'eau sont respectivement de 1138 Hz et 666 Hz, ce qui

Fig. 15. Comparaisons des formes modales dans l'eau prédites par le modèle actuel et ANSYS-Virtual.Lab avec 24 éléments (dans le sens de l'envergure) × 24 (dans le sens de la corde) × 2 (épaisseur).

https://sider.ai/pdf-translate/E0G5UEKEG8A?fid=0 13/17
10/22/24, 7:18 PM PDF translate: -----A numerical method for predicting the hydroelastic response ofmarine propeller.pdf | Sider

J. Li et coll. / Recherche océanique appliquée 74 (2018) 188-204 201

Fig. 16. Maillages de SUBOFF en CFD.

Fig. 17. Sillage du sous-marin SUBOFF pour hélice 4382.

Tableau 2 Comparaisons des premiers 10èmes modes calculés par l'interaction fréquences, ce qui peut augmenter l’ampleur des forces hydrodynamiques.
fluide-structure unidirectionnelle (modes secs) et l'interaction fluide-structure
bidirectionnelle pour une pale de l'hélice 4382 (modes humides).

numéro de mode sous vide (modes secs) dans l'eau (modes humides) Déclin(%)
Il est à noter qu'une fois que les matrices ajoutées MW (k) et CW (k)
sont calculées par la méthode du panel dépendant de la fréquence et
1 1137.7 666.3 0,41
2 2310.8 1575.3 0,32
que la force entraînée Fr (k) est obtenue par la méthode du panel
3 3394.6 2491.7 0,27 dépendant du temps, les équations discrétisées du le système entier
4 4897,5 3513.5 0,28 avec un nombre considérable de degrés de liberté sera transformé en
5 6220.2 4492.0 0,28 un ensemble d’équations discrétisées réduites en employant la
6 8229.1 5698.9 0,31
technique modale réduite. Les équations discrétisées réduites
7 10087.1 7100.1 0,30
8 10233.9 7521.4 0,27
comportent très peu de degrés de liberté, ce qui est efficace sur le plan
9 12349.8 8803.5 0,29
informatique. Par exemple, dans le cas de Nmode = 2000 pris dans la
10 14575.8 11008.4 0,24 présente analyse qui est suffisant pour obtenir des solutions
convergées, le temps de calcul de la présente formulation implémentée
dans les scripts MATLAB sur un PC Pentium IV-2,4 GHz est d'environ
10 s.
correspondent à la 5ème fréquence de passage des pales (1144 Hz) et à
la 3ème fréquence de passage des pales (686 Hz). Par conséquent, 5. Conclusions
quelles que soient les interactions fluide-structure unidirectionnelles ou
bidirectionnelles, la fréquence des forces hydrodynamiques étant Une méthode de panneau basée sur le potentiel 3D couplée à un FEM 3D
proche des fréquences naturelles de l'hydroptère (mode sec ou mode est développée pour l'analyse hydroélastique des hélices marines
humide), les amplitudes de réactions atteindront à des valeurs élastiques. La méthode du panneau dépendant du temps est utilisée pour
significatives. Ceci est conforme à la théorie classique des vibrations. Il déterminer les forces hydrodynamiques, et la méthode du panneau
est à noter que comme 1144 Hz est plus proche de 1138 Hz que de 686 dépendant de la fréquence est utilisée pour calculer la masse ajoutée et les
Hz à 666 Hz, les amplitudes de vibration de la lame élastique obtenues matrices d'amortissement. La méthode FEM est appliquée pour évaluer les
par analyse d'interaction fluide-structure unidirectionnelle sont vibrations structurelles et les forces de réaction. Une technique de
supérieures à celles déterminées par analyse fluide-structure réduction modale combinée à la méthode Wilson est utilisée pour
bidirectionnelle. méthode d'interaction avec la structure, comme le surmonter la faible efficacité numérique causée par les matrices fluides
montre la figure 18. asymétriques ajoutées. L'hélice est considérée comme rigide avec une
2) Les fréquences propres des hélices prédites par l'analyse rigidité structurelle infinie ou élastique avec une rigidité finie. Les
bidirectionnelle de la structure du fluide sont près de 29 % inférieures à interactions fluide-structure unidirectionnelles et bidirectionnelles sont
celles obtenues par une analyse unidirectionnelle. Cela indique qu'une prises en compte pour l'analyse hydroélastique des hélices élastiques. Les
analyse d'interaction fluide-structure bidirectionnelle fortement couplée matrices de masse et d’amortissement ajoutées sont négligées dans
doit être utilisée pour prédire les performances instables des hélices l’analyse unidirectionnelle de l’interaction fluide-structure, mais ces
élastiques, même si la déformation de l'hélice est faible. Il est à noter matrices sont superposées aux matrices structurelles de masse et
qu'une diminution de la fréquence fondamentale de l'hélice peut d’amortissement dans l’analyse bidirectionnelle de l’interaction fluide-
entraîner des résonances des hélices à des températures inférieures. structure. La validité de la présente méthode est confirmée en comparant
les chiffres prédits.

https://sider.ai/pdf-translate/E0G5UEKEG8A?fid=0 14/17
10/22/24, 7:18 PM PDF translate: -----A numerical method for predicting the hydroelastic response ofmarine propeller.pdf | Sider

202 J. Li et coll. / Recherche océanique appliquée 74 (2018) 188-204

la formule de Bernoulli non linéaire complète pour déterminer directement

la force motrice du système de couplage fluide-structure. L'application de
la méthode de panel purement dépendante de la fréquence présentée par
Kuo et Vorus [9] n'est pas simple pour traiter les termes non linéaires de la
formule de Bernoulli non linéaire. De plus, la présente méthode surmonte
les inconvénients du panel purement dépendant du temps présenté par
Young [10], qui fournit explicitement de grandes informations physiques
sur l'effet des matrices ajoutées et le coût de calcul requis par le modèle
actuel est considérablement plus faible. que celle de la méthode du panel
dépendant du temps. De plus, la présente méthode de panneau dépendant
de la fréquence est appliquée à la surface de la pale plutôt qu'à la surface de
cambrure comme indiqué dans la réf. [7]. Par conséquent, l’effet
d’épaisseur de la lame peut être directement pris en compte dans le présent
modèle. De plus, une technique de réduction modale est utilisée pour
surmonter l'énorme ressource de calcul requise par la méthode de Wilson et
les modes imprécis déterminés par la technique de regroupement de type
HRZ telle qu'employée dans la réf. [10].

À l’avenir, cette méthode sera étendue pour analyser les caractéristiques

des matrices de masse ajoutée et d’amortissement pour les hélices
élastiques dans différentes conditions d’écoulement et avec différents
paramètres géométriques et matériels.

Reconnaissance

Ce travail a été soutenu par la Fondation nationale des sciences naturelles de

Chine (numéro de subvention 11602138).

Annexe A. Conditions aux limites cinématiques sur la surface de la

pale

Dans le système de coordonnées cartésiennes illustré à la figure 1, un point

matériel arbitraire (x, y, z) sur la surface de la pale a un vecteur de
déplacement décrit comme ( ıx (x, y, z, t) , ıy (x , y, z, t) , ız (x, y, z, t) ) . Le
les coordonnées du point deviennent ( ̄x, ̄y, ̄z), qui s'exprime comme suit :
( ıx (x, y, z, t) , ıy (x, y, z, t) , ız (x, y, z, t) )
( ̄x, ̄y, ̄z) = (x, y, z) + (A.1)

Laissern (x, y, z) = 0 définir le géométrique fonction lequel

décrit la surface non déformée, puis l'équation de la surface déformée peut
s'écrire comme suit :
N ( ̄x, ̄y, ̄z) = n ( ̄x − ıx , ̄y − ıy , ̄z − ız )
=0 (A.2)

Supposons que le vecteur normal de la surface de la pale déformée au

point ( ̄x, ̄y, ̄z) soit ( N ̄x , N ̄y , N ̄z ) , alors :

N¯x = nx(1 − ıxx x¯x − ıxy y¯x − ıxz z¯x) + ny(−ıyx x¯x
− ıyy y¯x − ıyz z¯x) +nz (−ızx x¯x − ızy y¯zz¯x − ı )

N¯y = nx(−ıxx x¯y − ıxy y¯y − ıxz z¯y) + ny(1 − ıyx x¯y
− ıyy y¯y − ıyz z¯y) +nz (−ızx x¯y − ızy y¯z¯y − y¯ ı ) (A.3)

N¯z = nx(−ıxx x¯z − ıxy y¯z − ıxz z¯z ) + ny(−ıyx x¯z

− ıyy y¯z − ıyz z¯z ) +nz (1 − ızx x¯z − ızy y¯z − z¯z )
Fig. 18. Comparaisons de la pale rigide et de la pale élastique en considérant les harmoniques
de force instable de l'interaction fluide-structure unidirectionnelle et bidirectionnelle pour
l'hélice 4382. (a) : direction x ; (b) : direction y ; (c) : direction z.
où ıxx représente la dérivée spatiale de ıx par rapport à x, et x ̄x
représente la dérivée spatiale par rapport à ̄xdirection. D'autres
résultats calques avec ceux obtenus à partir de logiciels commerciaux et des symboles sont définis comme similaires à la définition de ıxx et x ̄x.
littératures. On constate que l'interaction fluide-structure bidirectionnelle En prenant la dérivée spatiale pour les deux côtés de l’équation. (A.1), on
doit être utilisée pour prédire les performances instables des hélices obtient l’équation suivante :
élastiques. ⎡ ⎤ ⎡ ⎤
L'algorithme d'interaction fluide-structure considéré ici est similaire à ceux ⎢⎣ 1 + ixy ixz ⎢⎣ xy xz
⎥⎦ = 1 0 0
⎤

présentés par Kuo et Vorus [9] et Young [10]. Cependant, le modèle actuel ıxx ıyx 1 + ıyy
iyz ⎥⎦x ̄x
y ̄x oui oui ⎡ ⎣ 0 1 0
⎦ (A.4)
utilise un panneau dépendant du temps combiné à une méthode de panneau
1 + izz 0 0 1
dépendant de la fréquence, qui peut résoudre izx ızy zx z ̄y zz

https://sider.ai/pdf-translate/E0G5UEKEG8A?fid=0 15/17
10/22/24, 7:18 PM PDF translate: -----A numerical method for predicting the hydroelastic response ofmarine propeller.pdf | Sider

J. Li et coll. / Recherche océanique appliquée 74 (2018) 188-204 203

En supposant que le vecteur déplacement inconnu est petit et en substituant Nous supposons ∼p (R, t) = p (R) eiωt , et l'équation de Helmholtz peut
l’équation. (A.4) dans l’équation. (A.3), on obtient : être obtenu sous la forme :

N¯x = nx − nx ıxx − ny ıyx ∇2p (R) + k2p (R) = 0 (B.2)

− nz ızx N¯y = ny − nx ıxy
(A.5)
− ny ıyy − nz ızy N¯z = nz où k est le nombre d'onde, donné par : k = ω/c 0 , et c0 est la vitesse du son.
− nx ıxz − ny ıyz − nz ızz
Afin de calculer l'amplitude de pression acoustique p (R) dans l'équation. (B.2),
La composante de vitesse normale aux frontières solides de la lame déformée l’expression intégrale suivante est employée :
doit être nulle, ce qui conduit à l'expression suivante :
∫ p ( RQ ) ∂ G ( ∫ ∂p( )
RQ (RP, )
( )
V ( ̄x, ̄y, ̄z, t) + ∇ ( ̄x, ̄y, ̄z, t) − ∂ı∂ t ( ̄x, ̄y, ̄z,
2 p (RP ) =
RP , RQ ) ∂nQ ds- ∂nQG RQ ds

t) SB SB
(B.3)
· ( N ̄xi + N ̄yj + N )
̄z k =0 (A.6)

où SB désigne la surface de la lame. L'indice Q correspond au point

où V ( ̄x, ̄y, ̄z, t) représente le vecteur vitesse de variable dans l'intégration, et l'indice p correspond au point de
couler. ∇ ( ̄x, ̄y, ̄z, t) est la vitesse induite générée par la déformation contrôle sur le panneau. RP et RQ sont les vecteurs de position des
lame. ∂ı∂ t ( ̄x, ̄y, ̄z, t) représente le vecteur vitesse du point ( ̄x, ̄y, ̄z). points sur la lame. G ( RP , RQ ) est le vert
fonction pour l'équation de Helmholtz dans un domaine fluide 3D illimité,
qui est donnée par G ( RP , RQ ) = e−ik|RP −RQ | |RP−RQ | . nQ est l'out-
En appliquant le développement de Taylor pour V ( ̄x, ̄y, ̄z,
t), ∇ ( ̄x, ̄y, ̄z, t) et ∂ı ∂ t ( ̄x, ̄y, ̄z, t), on obtient : vecteur normal de l'unité de service. Si et seulement si k = 0, la fonction de
( je · ) Green dégénérera en solution fondamentale de l'équation de Laplace.
V ( ̄x, ̄y, ̄z, t) = V (x, y, z, t) + ∇ V (x, y, z, t)
La condition limite cinématique sur la surface de la pale pour le système de
( je · ) couplage acoustique-solide nécessite que la vitesse et l'accélération
∇ ( ̄x, ̄y, ̄z, t) = ∇ (x, y, z, t) + ∇ ∇ (x, y, z, t)
(A.7)
normales d'un point matériel sur la surface de la pale soient égales à celles
de la particule fluide au même endroit, exprimées comme suit :
∂ı ∂ t ( ̄x, ̄y, ̄z, t) = ∂ı∂ t (x, y, z, t) + (ı · ∇) ∂ı∂
t (x, y, z, t) ∇ ∼p (R, t) = − = − ∂2∼ı
(sur la surface de la lame) (B.4)
Supposons que ∇ (x, y, z, t) et que les parties non uniformes du flux
∂∼v (R, t) ∂ t (R, t) ∂ t2
entrant soient petites. En remplaçant l'équation. (A.5) et l’équation. (A.7) Éq. (B.4) peut être réécrit comme suit :
dans l’équation. (A.6), on obtient la condition aux limites cinématique
sur la surface de la pale, exprimée par : ∂ ̃p ∂n = − ∂2 ̃ın
(B.5)
(R, t) ∂ t2
∂ı ∂ t · n = ∂ ∂n + V · n + ((ı · ∇) V 0 − (V 0 ·
(A.8) En supposant que ̃ın (R, t) = ın (R) eiωt
∇) ı) · n conviendrait avec
̃p (R, t) = p (R) eiωt , Éq. (B.5) peut être exprimé comme suit :
où n est le unité normale vecteur de le où-
indiquer (x, y, z), ∂p
formé lame surface
( nx à
, ny , ).V0 et n est exprimé
√ 1 ∂n = ω2ın (B.6)
comme: nz représente le uniforme parties de
n2x +n2y +n2z

le flux entrant. La lame est supposée présenter une faible déformation élastique, En remplaçant l'équation. (B.6) dans l’équation. (B.3), la forme discrétisée de
et l’Eq. (B.3) peut s’écrire ainsi :
on obtient donc l'équation suivante : ( [ ])−1
∂G
− ∂ ∂n∼= V · n − ∂ı∂ t · n+ ((ı · {p} = − ω2 2 Je −
∂n Gin (B.7)
(A.9)
∇) V 0) · n
et on obtient la force nodale suivante :
Par souci de commodité, les conditions aux limites décrites par l’équation ci-
( [ ])−1
dessus peuvent être décomposées en deux conditions aux limites différentes, ∂G
données comme suit : f = ZT s {p} = − ω2ZT s 2 Je −
∂n Gzu (B.8)

− ∂ r∂n = V ·
(A.10) L'expression détaillée de la fonction de Green G est donnée par
n cos (k|RP − RQ | ) − je pèche ( k|RP )
− ∂ v∂n = − ∂ı∂ t · n+ ((ı · ∇) G (RP, RQ ) = e−ik|RP −RQ | − RQ |
(A.11) |RP-RQ | =
V 0) · n |RP-RQ |

(B.9)
où l'équation. (A.10) et l’équation. (A.11) sont les conditions cinématiques aux
limites pour la lame rigide et la lame élastique, respectivement.
La force nodale pourrait être réécrite comme
Annexe B. Comparaisons entre l'interaction fluide-structure et
f = −(Continu + iCau) (B.10)
l'interaction acoustique-solide
où maman et Ca sont la masse ajoutée et l'amortissement ajouté
Le champ de pression acoustique ̃p (R, t) est régi par une équation d’onde 3D matrice.
dans le domaine temporel, donnée par
Après une dérivation rigoureuse, l’énoncé suivant a pu être obtenu :
∇2 ̃p (R, t) − ∂2 ̃p (R, =0 (B.1)
1 c20 t) ∂ t2 Ma = Mw |V0=0, Ca = Cw | V 0 = 0 = 0, pour k = 0 (B.11)

https://sider.ai/pdf-translate/E0G5UEKEG8A?fid=0 16/17
10/22/24, 7:18 PM PDF translate: -----A numerical method for predicting the hydroelastic response ofmarine propeller.pdf | Sider

204 J. Li et coll. / Recherche océanique appliquée 74 (2018) 188-204

En résumé, la matrice de masse ajoutée déterminée par l'analyse [15] XD He, Y. Hong, RG Wang, Optimisation hydroélastique d'une hélice marine composite dans
d'interaction fluide-structure pour le cas d'un écoulement stationnaire est un sillage non uniforme, Ocean Eng. 39 (2012) 14-23, http://dx.doi.
org/10.1016/j.oceaneng.2011.10.007.
égale à celle obtenue à partir d'une analyse d'interaction acoustique-solide
[16] H. Lee, M. Song, J. Suh, B. Chang, Analyse hydroélastique des hélices marines basée sur
lorsque le nombre d'onde du fluide est nul. Les matrices d’amortissement un algorithme FSI couplé BEM-FEM, Int. J. Nav. Archit. Océan Ing. 6 (2014) 562-577,
ajouté disparaissent toutes deux dans un tel cas. Ainsi, les équations http://dx.doi.org/10.2478/IJNAOE-2013-0198.
gouvernantes de l’interaction fluide-structure auront les mêmes solutions [17] L. Morino, C.-C. Kuo, Potentiel aérodynamique subsonique pour les
configurations complexes : une théorie générale, AIAA J. 12 (1974) 191-197.
que l’équation gouvernante de l’interaction acoustique-solide.
[18] T. Hoshino, Analyse hydrodynamique d'hélices en écoulement constant à l'aide d'une
méthode de panneau de surface, J. Soc. Nav. Archit. Japon 165 (1989).
[19] X.-W. Gao, Évaluation des intégrales de domaine régulières et singulières avec
Références
discrétisation aux limites uniquement – ​théorie et code Fortran, J. Comput. Appl.
Mathématiques. 175 (2005) 265-290.
[1] R. Kamakoti, W. Shyy, Interaction fluide-structure pour les applications aéroélastiques, Prog. [20] MS Tarafder, GK Saha, ST Mehedi, Analyse de l'écoulement potentiel autour
Aérosp. Sci. 40 (2004) 535-558. d'hydroptères tridimensionnels par la méthode de panneau basée sur une source combinée
[2] PJ Maljaars, ML Kaminski, Analyse hydroélastique des hélices flexibles : un aperçu, dans : et un dipôle, J. Mar. Sci. Technologie. 18 (2010) 376-384.
Fourth Int Symp. Mar. Propulsors, Austin, Texas, États-Unis, 2015. [21] JL Hess, AMO Smith, Calcul du flux potentiel de non-soulèvement sur des corps
[3] D. Taylor, La vitesse et la puissance des navires, Washington DC, 1933. [4] J. Conolly, tridimensionnels arbitraires, J. Ship Res. 8 (1964) 22-44.
Force des hélices, Trans. R. Inst. Nav. Archit. 103 (1961) [22] JL Hess, AMO Smith, Calcul du flux potentiel sur les corps arbitraires, Prog. Aéronaute. Sci.
139-204. [5] P. Genalis, Résistance élastique des hélices — Analyse par méthodes matricielles, 8 (1967) 1-138.
[23] Hassan Ghassemi, AR Kohansal, Évaluation numérique de différents niveaux
1970. [6] P. Atkinson, La prédiction de la distorsion et des contraintes des hélices marines à l'aide d'intégrales singulières, découlant du BEM et son application dans l'analyse des hydroptères,
d'un Appl. Mathématiques. Calculer. 213 (2009) 277-289.
modèle d'éléments finis à coque épaisse super paramétrique, Trans. R. Inst. Nav. Archit. 115 [24] JE Kerwin, C.-S. Lee, Prédiction des performances stables et instables des hélices
(1973) 359-375. marines par la théorie numérique de la surface de levage, SNAME Trans. 86 (1978) 218-253.
[7] J. Ma, Contraintes dans les hélices marines, J. Ship Res. 18 (1974) 252. [8] T. Sontvedt,
Contraintes des pales d'hélice, application des méthodes d'éléments finis, [25] DS Greeley, Justin E. Kerwin, Méthodes numériques pour la conception et l'analyse d'hélices
Calculer. Structure. 4 (1974) 193-204. [9] J. Kuo, W. Vorus, Contrainte dynamique des pales en débit constant, SNAME Trans. 90 (1982) 415-453.
d'hélice, dans : Tenth Sh. Technologie. Rés. [26] Y. Ukon, T. Kudo, Y. Kurobe, Mesure de la distribution de pression sur des hélices à
Symp., Norfolk, 1985, pp. [10] P. Atkinson, J. Glover, Effets hydroélastiques de l'hélice. grande échelle, dans : Proceeding Propellers/Shafting'91, SNAME, Virginia Beach, 1991.
transactions, Trans. Soc.
Nav. Archit. Mar. Ing. 21 (1988). [11] H. Lin, J. Lin, Comportement hydroélastique non [27] Y. Ukon, T. Kudo, et al., Mesure de la répartition de la pression sur une hélice à grande
linéaire des hélices utilisant un échelle, mesure sur une hélice très asymétrique, J. Soc. Nav. Archit. Japon 170 (1990).
méthode des éléments et théorie des surfaces de levage, J. Mar. Sci. Technologie. 1
(1996) 114-124. [28] Y. Ukon, H. Yuasa, Répartition de la pression et contrainte des pales sur une hélice fortement
[12] YL Young, Analyse hydroélastique dépendant du temps des propulseurs cavitants, J. asymétrique, dans : 19e Symp Nav. Hydrodyn., Séoul, Corée, 1992.
Fluids Struct. 23 (2007) 269-295, http://dx.doi.org/10.1016/j.jfluidstructs. 2006.09.003. [29] J. Tian, ​P. Croaker, J. Li, H. Hua, Études expérimentales et numériques sur la vibration
induite par fx1flow des pales d'hélice sous un afflux non uniforme, Proc. Inst. Mécanique. Ing.
[13] YL Young, Analyse d'interaction fluide-structure d'hélices marines composites flexibles, Partie M J. Eng. Marit. Environ. 231 (2016), http://dx.doi.org/10. 1177/1475090216654306.
J. Fluids Struct. 24 (2008) 799-818, http://dx.doi.org/10.1016/j. jfluidstructs.2007.12.010.
[30] J. Tian, ​Z. Zhang, Z. Ni, H. Hua, Analyse des vibrations induites par l'écoulement
[14] J. Neugebauer, M. Abdel-maksoud, M. Braun, Interaction fluide-structure des hélices, 2008. d'hélices élastiques dans un afflux cyclique : une étude expérimentale et numérique, Appl.
Océan Rés. 65 (2017) 47-59, http://dx.doi.org/10.1016/j.apor.2017.03.014.
[31] NC Groves, TT Huang, MS Chang, Caractéristiques géométriques des modèles DARPA
SUBOFF, 1989.

https://sider.ai/pdf-translate/E0G5UEKEG8A?fid=0 17/17