IdentifiantMot de passe
Loading...
Mot de passe oublié ?Je m'inscris ! (gratuit)
Navigation

Inscrivez-vous gratuitement
pour pouvoir participer, suivre les réponses en temps réel, voter pour les messages, poser vos propres questions et recevoir la newsletter

MATLAB Discussion :

Création d'une sphère ouverte des 2 cotés


Sujet :

MATLAB

  1. #1
    Membre régulier Avatar de Décembre
    Inscrit en
    Avril 2010
    Messages
    277
    Détails du profil
    Informations forums :
    Inscription : Avril 2010
    Messages : 277
    Points : 110
    Points
    110
    Par défaut Création d'une sphère ouverte des 2 cotés
    Salut,

    Afin de créer une sphère centrée en (x,y,z) de rayon r=1.5mm alors j'ai fait ce qui suit:
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
     
     
    [xo,yo,zo] = sphere(30);
     
     
    r=1.5;
    phi=linspace(0,pi,30);
    theta=linspace(0,2*pi,40);
    [phi,theta]=meshgrid(phi,theta);
    xo=r*sin(phi).*cos(theta);
    yo=r*sin(phi).*sin(theta);
    zo=r*cos(phi);
    surf(xo+x, yo+y, zo+z)
    Maintenat je voudrais faire des ouvertures des deux coté de la sphère. comment pourrais-je le faire?
    merci

  2. #2
    Rédacteur/Modérateur

    Avatar de Jerome Briot
    Homme Profil pro
    Freelance mécatronique - Conseil, conception et formation
    Inscrit en
    Novembre 2006
    Messages
    20 302
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Haute Garonne (Midi Pyrénées)

    Informations professionnelles :
    Activité : Freelance mécatronique - Conseil, conception et formation

    Informations forums :
    Inscription : Novembre 2006
    Messages : 20 302
    Points : 53 160
    Points
    53 160
    Par défaut
    Tu ne nous dis rien sur la forme de ces ouvertures...

    Essaie ceci :

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
    14
    15
    16
    17
    18
    r = 1.5;
    
    phi = linspace(pi*.1,pi*.9,30);
    theta = linspace(0,2*pi,40);
    
    [phi,theta] = meshgrid(phi,theta);
    
    x = r*sin(phi).*cos(theta);
    y = r*sin(phi).*sin(theta);
    z = r*cos(phi);
    
    xo = rand;
    yo = rand;
    zo = rand;
    
    figure
    surf(xo+x, yo+y, zo+z)
    axis equal vis3d
    Ingénieur indépendant en mécatronique - Conseil, conception et formation
    • Conception mécanique (Autodesk Fusion 360)
    • Impression 3D (Ultimaker)
    • Développement informatique (Python, MATLAB, C)
    • Programmation de microcontrôleur (Microchip PIC, ESP32, Raspberry Pi, Arduino…)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

  3. #3
    Membre régulier Avatar de Décembre
    Inscrit en
    Avril 2010
    Messages
    277
    Détails du profil
    Informations forums :
    Inscription : Avril 2010
    Messages : 277
    Points : 110
    Points
    110
    Par défaut
    Merci beaucoup Dut

    Voila, grâce à vous j'ai obtenu le ressort que je voulais avoir:


    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
    14
    15
    16
    17
    18
    19
    20
    21
    22
    23
    24
    t = 0:pi/30:3*pi;
    xs=sin(t);
    ys=cos(t);
    t1(:,1)=xs;
    t1(:,2)=ys;
    t1(:,3)=t;
    for i=1: length(t1)
        r = 1;
        phi = linspace(pi/2,pi/2,30);
        theta = linspace(0,2*pi,40);
     
        [phi,theta] = meshgrid(phi,theta);
     
        xo = r*sin(phi).*cos(theta);
        yo = r*sin(phi).*sin(theta);
        zo = r*cos(phi);
     
        %xo = rand;
        %yo = rand;
        %zo = rand;
     
     
        mesh(xo+t1(i,1), yo+t1(i,2), zo+t1(i,3))
    end
    Maintenant j'essaie de transformer ce mesh en isosurface mais je n'y arrive pas voila ce que j'ai ajouté mais ça ne donne rien:
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    patch(isosurface(xo+t1(i,1), yo+t1(i,2), zo+t1(i,3),-3))
    Pouvez-vous m'aidez s'il vous plait?

  4. #4
    Rédacteur/Modérateur

    Avatar de Jerome Briot
    Homme Profil pro
    Freelance mécatronique - Conseil, conception et formation
    Inscrit en
    Novembre 2006
    Messages
    20 302
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Haute Garonne (Midi Pyrénées)

    Informations professionnelles :
    Activité : Freelance mécatronique - Conseil, conception et formation

    Informations forums :
    Inscription : Novembre 2006
    Messages : 20 302
    Points : 53 160
    Points
    53 160
    Par défaut
    Voici ton code optimisé :

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
    14
    15
    16
    17
    18
    19
    20
    21
    22
    23
    24
    25
    26
    27
    28
    29
    30
    31
    32
    figure
    hold on
     
    t = 0:pi/30:3*pi;
    xs = sin(t);
    ys = cos(t);
     
    r = 1;
    phi = linspace(pi/2,pi/2,30);
    theta = linspace(0,2*pi,40);
     
    [phi,theta] = meshgrid(phi,theta);
     
    xo = r*sin(phi).*cos(theta);
    yo = r*sin(phi).*sin(theta);
    zo = r*cos(phi);
     
    n = numel(xs);
     
    X = zeros(40,n);
    Y = zeros(40,n);
    Z = zeros(40,n);
     
    for i = 1:n    
        X(:,i) = xo(:,1)+xs(i);
        Y(:,i) = yo(:,1)+ys(i);
        Z(:,i) = zo(:,1)+t(i);
    end
     
    mesh(X,Y,Z)
     
    axis equal vis3d
    Que cherches-tu à faire ensuite ?
    Ingénieur indépendant en mécatronique - Conseil, conception et formation
    • Conception mécanique (Autodesk Fusion 360)
    • Impression 3D (Ultimaker)
    • Développement informatique (Python, MATLAB, C)
    • Programmation de microcontrôleur (Microchip PIC, ESP32, Raspberry Pi, Arduino…)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

  5. #5
    Membre régulier Avatar de Décembre
    Inscrit en
    Avril 2010
    Messages
    277
    Détails du profil
    Informations forums :
    Inscription : Avril 2010
    Messages : 277
    Points : 110
    Points
    110
    Par défaut
    Merci Dut pour l'optimisation

    En fait je cherche à exporter ce ressort en format *.STL, pour le lire dans un logiciel de modélisation 3D tel que blender par exemple.


    Pour cela il faudrait que je crée l'isosurface pour avoir la liste des vertex et des faces, j'ai lu sur le help de matlab que c'est faisable avec la fonction isosurface mais l'exemple montre le cas d'une matrice 3D or dans mon cas je n'ai qu'un mesh plotté

  6. #6
    Rédacteur/Modérateur

    Avatar de Jerome Briot
    Homme Profil pro
    Freelance mécatronique - Conseil, conception et formation
    Inscrit en
    Novembre 2006
    Messages
    20 302
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Haute Garonne (Midi Pyrénées)

    Informations professionnelles :
    Activité : Freelance mécatronique - Conseil, conception et formation

    Informations forums :
    Inscription : Novembre 2006
    Messages : 20 302
    Points : 53 160
    Points
    53 160
    Par défaut
    Citation Envoyé par Décembre Voir le message
    Pour cela il faudrait que je crée l'isosurface pour avoir la liste des vertex et des faces,
    Utilise la fonction SURF2PATCH :

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    f = surf2patch(X,Y,Z,'triangles');
    Ingénieur indépendant en mécatronique - Conseil, conception et formation
    • Conception mécanique (Autodesk Fusion 360)
    • Impression 3D (Ultimaker)
    • Développement informatique (Python, MATLAB, C)
    • Programmation de microcontrôleur (Microchip PIC, ESP32, Raspberry Pi, Arduino…)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

  7. #7
    Membre régulier Avatar de Décembre
    Inscrit en
    Avril 2010
    Messages
    277
    Détails du profil
    Informations forums :
    Inscription : Avril 2010
    Messages : 277
    Points : 110
    Points
    110
    Par défaut
    Je ne sais pas comment te remercier Dut toutes les formules de remerciement du monde ne peuvent pas exprimer ma gratitude

    Merci

  8. #8
    Membre régulier Avatar de Décembre
    Inscrit en
    Avril 2010
    Messages
    277
    Détails du profil
    Informations forums :
    Inscription : Avril 2010
    Messages : 277
    Points : 110
    Points
    110
    Par défaut
    Bonsoir,

    Je n'arrive pas à arrêter de poser des questions

    Grâce à l'aide de dut j'ai réussi à modéliser un ressort à partir d'une sphère, en le comparant avec les ressorts que j'obtiens avec les logiciels de modélisation 3D j'ai constater que le modéliser avec des sphères était une erreur alors j'ai penser à refaire la même chose avec un cylindre:

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
     
    t = 0:0.01:pi;
    xs = sin(pi*t);
    ys = cos(pi*t);
     
    [x,y,z] = cylinder(0.5)
     
    for i=1:length(t)
        hold on
        surf(xs(i)+x, ys(i)+y, t(i)+z)
    end
    Si vous exécuter le programme vous allez comparer que les cylindres sont mis dans le mauvais sens, alors j'aimerais bien vous demander s'il y a moyen de les inverser

+ Répondre à la discussion
Cette discussion est résolue.

Discussions similaires

  1. [AC-97] Création d'une table à partir des résultats d'une requete analyse croisée
    Par docjo dans le forum Requêtes et SQL.
    Réponses: 2
    Dernier message: 11/11/2009, 22h46
  2. [Data] Création d'une session avec des données du runtime
    Par NeedToHelp dans le forum Spring
    Réponses: 2
    Dernier message: 17/07/2009, 12h13
  3. Réponses: 4
    Dernier message: 22/11/2007, 20h23
  4. une caméra, une sphère et des collisions
    Par Kurisu dans le forum Développement 2D, 3D et Jeux
    Réponses: 1
    Dernier message: 22/06/2007, 16h28
  5. prb création d'une fonction utilisant des fonctions de postgis
    Par Benjamin_es dans le forum PostgreSQL
    Réponses: 2
    Dernier message: 13/03/2007, 12h34

Partager

Partager
  • Envoyer la discussion sur Viadeo
  • Envoyer la discussion sur Twitter
  • Envoyer la discussion sur Google
  • Envoyer la discussion sur Facebook
  • Envoyer la discussion sur Digg
  • Envoyer la discussion sur Delicious
  • Envoyer la discussion sur MySpace
  • Envoyer la discussion sur Yahoo