begin process at 2012 02 05 02:13:38
  Trouver un code source :
 
dans
 
Accueil > 

Code

 > 

Math & Algorithmes

 > COLLISIONS EN 3D AVEC VISUAL PYTHON

COLLISIONS EN 3D AVEC VISUAL PYTHON


 Information sur la source

Note :
9 / 10 - par 2 personnes
9,00 / 10

  • 1

  • 2

  • 3

  • 4

  • 5

  • 6

  • 7

  • 8

  • 9

  • 10
Catégorie :Math & Algorithmes Classé sous :collisions, billard Niveau :Débutant Date de création :24/01/2006 Vu / téléchargé :4 963 / 370

Auteur : HCD

Ecrire un message privé
Commentaire sur cette source (5)
Ajouter un commentaire et/ou une note

 Description

Ce code fait suite à "COLLISIONS en cascades" déposé fin 2005 sur le site.
Il met en oeuvre l'extension VISUAL de PYTHON, très intéressante pour les représentation graphiques en 3D.

Source

  • #! /usr/bin/env python
  • # -*- coding: Latin-1 -*-
  • "COLLISIONS_VPYTHON"
  • "Ce script est démarqué de l'exemple fourni dans 'GAS' de l'extension VISUAL de PYTHON"
  • from visual import *
  • from time import clock
  • from visual.graph import *
  • from random import random
  • from visual.controls import*
  • #################################################################################
  • win=500
  • # Valeurs typiques
  • L=1.0
  • gray=(0.7,0.7,0.7) # couleur des bordures
  • Raxes=0.01 # rayon des bordures
  • scene=display(title="COLLISIONS_VPYTHON",width=win,height=win,background=(0,1,1),x=0,y=0,range=L/4,center=(L/2.,L/2.,L/2.))
  • scene.lights=[vector(1,1,1)]
  • scene.ambient=0
  • xaxis=curve(pos=[(0,0,0),(L,0,0)],color=gray,radius=Raxes)
  • yaxis=curve(pos=[(0,0,0),(0,L,0)],color=gray,radius=Raxes)
  • xaxis2=curve(pos=[(0,L,0),(L,L,0)],color=gray,radius=Raxes)
  • yaxis2=curve(pos=[(L,0,0),(L,L,0)],color=gray,radius=Raxes)
  • dt=0.002
  • #################################################################################
  • def boules():
  • ""
  • global pos,v,radius,boule,N,R,Nsteps,t,dt
  • N =10 # valeur à modifier pour augmenter ou diminuer le nombre de boules
  • R=0.05
  • boule=[]
  • colors=[color.red,color.green,color.blue,color.yellow,color.cyan,color.magenta,color.white]
  • poslist,vlist,rlist=[],[],[]
  • for i in range(N):
  • Lmin=R
  • Lmax=L-Lmin
  • x,y,z,r=R+(Lmax-Lmin)*random(),R+(Lmax-Lmin)*random(),0,R
  • boule=boule+[sphere(pos=(x,y,z),radius=r,color=colors[i % 7])]
  • theta,phi=pi*random(),2*pi*random()
  • vx,vy,vz=sin(theta)*cos(phi),sin(theta)*sin(phi),0
  • poslist.append((x,y,z))
  • vlist.append((vx,vy,vz))
  • rlist.append(r)
  • pos=array(poslist)
  • v=array(vlist)
  • radius=array(rlist)
  • t=0.0
  • Nsteps=0
  • pos=pos+v*(dt/2.)
  • time=clock()
  • ################################################################################
  • def move_boules():
  • ""
  • global pos,v,radius,boule,Nsteps,t,dt
  • while 1:
  • rate(100)
  • pos =pos+ v*dt
  • r = pos-pos[:,NewAxis]
  • rmag = sqrt(add.reduce(r*r,-1))
  • hit = less_equal(rmag,radius+radius[:,NewAxis])-identity(N)
  • hitlist = sort(nonzero(hit.flat)).tolist()
  • # si collisions:
  • for ij in hitlist:
  • i,j= divmod(ij,N)
  • hitlist.remove(j*N+i)
  • vtot=v[i]+v[j]
  • vi,vj=v[i],v[j]
  • ri,rj=boule[i].radius,boule[j].radius
  • a=mag(vj-vi)**2
  • if a==0: continue
  • b=2*dot(pos[i]-pos[j],vj-vi)
  • c=mag(pos[i]-pos[j])**2-(ri+rj)**2
  • d=b**2-4.*a*c
  • if d<0: continue
  • deltat=(-b+sqrt(d))/(2.*a) # t-deltat lorsque il y a contact
  • pos[i],pos[j]=pos[i]-v[i]*deltat,pos[j]-v[j]*deltat # retour à la configuration de contact
  • vcmi,vcmj=v[i]-vtot/2,v[j]-vtot/2
  • rrel=norm(pos[j]-pos[i])
  • vcmi,vcmj=vcmi-2*dot(vcmi,rrel)*rrel,vcmj-2*dot(vcmj,rrel)*rrel
  • v[i],v[j]=vcmi+vtot/2,vcmj+vtot/2
  • pos[i],pos[j]=pos[i]+v[i]*deltat,pos[j]+v[j]*deltat
  • # rebonds sur les bordures
  • outside=less_equal(pos,R+Raxes)
  • v1=v*outside
  • v=v-v1+abs(v1)
  • outside=greater_equal(pos,L-R-Raxes)
  • v1=v*outside
  • v=v-v1-abs(v1)
  • # mise à jour des positions
  • for i in range(N):
  • boule[i].pos=pos[i]
  • Nsteps=Nsteps+1
  • t=t+dt
  • def lancer():
  • ""
  • boules()
  • move_boules()
  • lancer()
#! /usr/bin/env python 
# -*- coding: Latin-1 -*- 

"COLLISIONS_VPYTHON"
"Ce script est démarqué de l'exemple fourni dans 'GAS' de l'extension VISUAL de PYTHON"

from visual import *
from time import clock
from visual.graph import *
from random import random
from visual.controls import*

#################################################################################

win=500
# Valeurs typiques
L=1.0 
gray=(0.7,0.7,0.7) # couleur des bordures
Raxes=0.01 # rayon des bordures
scene=display(title="COLLISIONS_VPYTHON",width=win,height=win,background=(0,1,1),x=0,y=0,range=L/4,center=(L/2.,L/2.,L/2.))
scene.lights=[vector(1,1,1)]
scene.ambient=0
xaxis=curve(pos=[(0,0,0),(L,0,0)],color=gray,radius=Raxes)
yaxis=curve(pos=[(0,0,0),(0,L,0)],color=gray,radius=Raxes)
xaxis2=curve(pos=[(0,L,0),(L,L,0)],color=gray,radius=Raxes)
yaxis2=curve(pos=[(L,0,0),(L,L,0)],color=gray,radius=Raxes)
dt=0.002

#################################################################################

def boules():
    ""
    global pos,v,radius,boule,N,R,Nsteps,t,dt
    N =10 # valeur à modifier pour augmenter ou diminuer le nombre de boules
    R=0.05 
    boule=[]
    colors=[color.red,color.green,color.blue,color.yellow,color.cyan,color.magenta,color.white]
    poslist,vlist,rlist=[],[],[]
    for i in range(N):
        Lmin=R
        Lmax=L-Lmin
        x,y,z,r=R+(Lmax-Lmin)*random(),R+(Lmax-Lmin)*random(),0,R
        boule=boule+[sphere(pos=(x,y,z),radius=r,color=colors[i % 7])]
        theta,phi=pi*random(),2*pi*random()
        vx,vy,vz=sin(theta)*cos(phi),sin(theta)*sin(phi),0
        poslist.append((x,y,z))
        vlist.append((vx,vy,vz))
        rlist.append(r)
    pos=array(poslist)
    v=array(vlist)
    radius=array(rlist)
    t=0.0
    Nsteps=0
    pos=pos+v*(dt/2.) 
    time=clock()

################################################################################

def move_boules():
    ""
    global pos,v,radius,boule,Nsteps,t,dt
    while 1:
        rate(100)
        pos =pos+ v*dt
        r = pos-pos[:,NewAxis] 
        rmag = sqrt(add.reduce(r*r,-1))
        hit = less_equal(rmag,radius+radius[:,NewAxis])-identity(N)
        hitlist = sort(nonzero(hit.flat)).tolist()
       # si collisions:
        for ij in hitlist:
            i,j= divmod(ij,N)
            hitlist.remove(j*N+i)
            vtot=v[i]+v[j]
            vi,vj=v[i],v[j]
            ri,rj=boule[i].radius,boule[j].radius
            a=mag(vj-vi)**2
            if a==0: continue 
            b=2*dot(pos[i]-pos[j],vj-vi)
            c=mag(pos[i]-pos[j])**2-(ri+rj)**2
            d=b**2-4.*a*c
            if d<0: continue 
            deltat=(-b+sqrt(d))/(2.*a) # t-deltat lorsque il y a contact
            pos[i],pos[j]=pos[i]-v[i]*deltat,pos[j]-v[j]*deltat # retour à la configuration de contact
            vcmi,vcmj=v[i]-vtot/2,v[j]-vtot/2 
            rrel=norm(pos[j]-pos[i])
            vcmi,vcmj=vcmi-2*dot(vcmi,rrel)*rrel,vcmj-2*dot(vcmj,rrel)*rrel
            v[i],v[j]=vcmi+vtot/2,vcmj+vtot/2
            pos[i],pos[j]=pos[i]+v[i]*deltat,pos[j]+v[j]*deltat
        # rebonds sur les bordures
        outside=less_equal(pos,R+Raxes)
        v1=v*outside
        v=v-v1+abs(v1)
        outside=greater_equal(pos,L-R-Raxes)
        v1=v*outside
        v=v-v1-abs(v1)
        # mise à jour des positions
        for i in range(N):
            boule[i].pos=pos[i]
        Nsteps=Nsteps+1
        t=t+dt

def lancer():
    ""
    boules()
    move_boules()
lancer()



 Fichier Zip

Les Membres Club peuvent télécharger directement un fichier contenu dans le zip sans télécharger le zip en entier !

Télécharger le zip


 Sources du même auteur

Source avec Zip MASTERMIND (PAR JULIE ET HCD)
Source avec Zip Source avec une capture RÉSOLUTION DES ÉQUATIONS DU 4ÈME DEGRÉ
Source avec Zip Source avec une capture LE CENTRE DU CERCLE RETROUVÉ PAR NAPOLÉON !
Source avec Zip Source avec une capture CERCLE ET DROITE D' EULER ANIMÉS DANS UN TRIANGLE
Source avec Zip Source avec une capture UN SERPENT PYTHON ... PAS COMME LES AUTRES !

 Sources de la même categorie

Source avec Zip BASE64 ENCRYPT/DECRYPT PYTHON BY MAXOU56800 par Maxou56800
Source avec Zip Source avec une capture TRIANGULATION par mecrosoft
Source avec Zip Source avec une capture COURBE DE BEZIER par mecrosoft
Source avec Zip Source avec une capture CALCUL D'AIRE D'UN TRIANGLE [INTERFACE GRAPHIQUE] par SeventhSon
Source avec Zip Source avec une capture SUITE DE FIBONACCI [INTERFACE GRAPHIQUE] par SeventhSon

 Sources en rapport avec celle ci

Source avec Zip Source avec une capture UN SERPENT PYTHON ... PAS COMME LES AUTRES ! par HCD
Source avec Zip Source avec une capture COLLISIONS EN CASCADES : PLUS D'UNE CENTAINE DE BOULES S'ENT... par HCD
Source avec Zip BILLARD FRANÇAIS par HCD

Commentaires et avis

Commentaire de tony1239456 le 17/04/2006 22:00:14

Salut, je voulais dire que ton code est pas mal, mais connais-tu un bon tuto sur le web permettant d'apprendre a se servir de visual python en français ? Merci d'avance.

Commentaire de HCD le 17/04/2006 23:16:21

TONY 1239456

Voici la réponse à ta question :

http://guigui.developpez.com/Tutoriel/Python/frvpython/webdoc/

Commentaire de tony1239456 le 18/04/2006 12:56:59

Merci pour le tuto, mais comme je voudrais créer un mini-jeu en 3D, j'aimerais savoir comment importer ses images au lieu de mettre des couleurs rouge, bleu... Merci d'avance

Commentaire de zorg724 le 01/11/2007 20:39:20

On pourra mettre des image dans la version 2008 de Vpython (ajout de texture et de transparence)

Commentaire de geek79soa le 08/01/2011 11:40:51

Bonjour!!
J'ai un projet à faire sur la collision de deux boules: choc élastique. A vrai dire je dois programmer cette collision sur Python mais je sais vraiment pas quoi faire.
Quelqu'un pourrait m'aider s'il vous plait .
Merci

 Ajouter un commentaire




Nos sponsors


Sondage...

Comparez les prix

CalendriCode

Février 2012
LMMJVSD
  12345
6789101112
13141516171819
20212223242526
272829    

Consulter la suite du CalendriCode

 
Développement réalisé par Nicolas SOREL (Nix) avec l'aide de : Cyril DURAND et Emmanuel (EBArtSoft), Merci à Vincent pour ses précieux conseils.
CodeS-SourceS.com© Toute reproduction même partielle est interdite sauf accord écrit du Webmaster
CodeS-SourceS.com© est une marque déposée tous droits réservés

Google Coop CodeS-SourceS Google Coop CodeS-SourceS
Temps d'éxécution de la page : 1,825 sec (3)

Nous contacter | Annoncer sur CodeS-SourceS | Mentions légales