Explication formule fonction gaussien

[weld]

Bonjour à tous,

Je sollicite l'aide de la communauté digischool au sujet de la programmation d'une source de chaleur modélisant une opération de soudage sur un tube (trajectoire de la source circulaire sur un plan) . En gros, j'aimerais avoir un éclaircissement sur l'opération mathématique réalisée à partir de la ligne 8; une explication sur le repère X3,Y3 ainsi que X4,Y4 . Je tiens à préciser que le programme est un extrait de l'ensemble que l'on calcul avec python afin de modéliser la propagation thermique sur le tube. Le programme fonctionne ainsi je suis chargé d'expliquer ce que l'on fait mathématiquement.

La seule incohérence que l'on a c'est que la température minimum calculé est de -373 °C..... Tous scientifiques sait que c'est impossible !!

1# def gaussien(X,Y,Z,t,FY,T0):
2# X1 = (X-R1*cos(T0))
3# Y1 = (Y-R1*sin(T0))
4# Z1 = Z-3.006
5# X2 = cos(T0)*X1+sin(T0)*Y1
6# Y2 = cos(T0)*Y1-sin(T0)*X1
7# Z2 = Z1
8# X3 = (X-R1*cos(T0+vit*(t-tmaintien)))
9# Y3 = (Y-R1*sin(T0+vit*(t-tmaintien)))
10# X4 = cos(T0+vit*(t-tmaintien))*X3+sin(T0+vit*(t-tmaintien))*Y3
11# Y4 = cos(T0+vit*(t-tmaintien))*Y3-sin(T0+vit*(t-tmaintien))*X3
12# Q=((6*sqrt(3)*ETA*U*I)/((ah*bh*ch)*3.1*sqrt(3.1)))*exp(-((3*X4**2/ch**2)+(3*Y4**2/ah**2)+(3*(Z1)**2
13# /bh**2)))
14# if(t< tmaintien):
15# Q=((6*sqrt(3)*ETA*U*I)/((ah*bh*ch)*3.1*sqrt(3.1)))*exp(-((3*X2**2/ch**2)+(3*Y2**2/ah**2)+(3*Z2**2
16# /bh**2)))
17# if(t > tfin):
18# Q=0.
19# return Q


Jamais le soleil ne voit l'ombre. - Léonard De Vinci

[pascal16]

un petit dessin du repère serait le bienvenu.
ce que sont X, Y..et le Q retourné.
quelques explications sur TO qui ressemble à un angle, puis à une température.
U, I ETA ah,bh,ch sont définies nulle part
tu passes FY en variable, et elle ne sert pas.

Vu de loin, Tmaintient est une température de maintient, donc ta température finale est la température de maintient.
X1 Y1et Z1, un calcul de changement d'origine en coordonnées cylindrique (normal pour faire des cercles)