Approche de modélisation :
Je modélise sur Matlab-Simulink une structure flexible très fine. Tous les points du modèle sont reliés entre eux par des ressorts et des amortisseurs de cette façon (sans les attaches au centre):
Nom : Description du maillage.png
Affichages : 161
Taille : 104,6 Ko

L'équation générale de mon modèle appliquée à chaque point du maillage est la suivante:
Nom : Formule d'un système dynamique mass-ressort-ammortisseur.PNG
Affichages : 174
Taille : 1,2 Ko
Avec k la rigidité des ressorts, et c l'amortissement des amortisseurs.

Pour adapter les propriétés physiques du matériau que je veux modéliser, la rigidité du ressort a été fixée à une valeur très élevée, autour de k = 5000. Cela signifie que mes liens à ressort sont très réactifs à toute déformation.

Problème :
Ceci mène à mon problème : les liens à haute rigidité induisent un déplacement à haute fréquence que je peux considérer comme du bruit dans la simulation.

  • La simulation est beaucoup plus lente car le pas de temps variable que j'utilise doit être très faible.
  • Ces déplacements à haute fréquence (autour de 160 Hz, qui est la fréquence de résonance des ressorts) restent tout au long des simulations.
    Voici une simulation de ma structure tournant à une vitesse angulaire constante :


Nom : Evolution temporelle des coordonnées sphériques d'un des points du maillage.png
Affichages : 164
Taille : 27,8 Ko

On peut voir que R vibre à une fréquence très élevée. Cependant, l'amplitude du déplacement est clairement négligeable.
Pour accélérer la simulation, je veux supprimer ces vibrations !

Investigation :
Pour les supprimer, j'étudie les techniques de filtrage du signal, principalement le filtrage passe-bas. A chaque boucle de notre simulation, et ce qui doit entrer dans notre filtre sont les données de tous mes points dans tous mes axes.

Nom : Block simulink d'un filtre passe-bas.png
Affichages : 167
Taille : 2,9 Ko

La version continue du filtre passe-bas de la bibliothèque Simulink a été testée sur l'accélération, la vitesse et la position, avec plusieurs fréquences de coupure de 100 Hz à 500 Hz.
Par exemple, pour une fréquence de coupure de 200Hz et en filtrant la position à t=0.6 sec j'ai :

Nom : Evolution temporelle filtrée des coordonnées sphériques d'un des points du maillage.png
Affichages : 154
Taille : 19,3 Ko
C'est un mouvement dans le plan donc je n'ai pas d'angle d'élévation, mais l'angle d'azimut et la distance du point au centre sont complètement divergents.
Slon moi, le problème peut provenir de :
  • Le fait que je sois dans un système en boucle fermée.
  • Le fait que pour mon maillage, le filtre reçoit 81 vecteurs de 3*1 à chaque pas de temps et peut-être que le bloc filtre n'est pas fait pour fonctionner avec cela.


Question principale :
  • Existe-t-il des techniques de filtrage pour les systèmes en boucle fermée et à entrées multiples qui pourraient résoudre mon problème ?


Merci d'avance pour votre aide !