def Tricercle_kTg(S1,S2,S3,coul,LABEL,Nb,kTg): IndCalculsFaits=1 SeqMaxi=0 y1=0 Hert1.grid(row =LigCas, column =1,columnspan=5,sticky='w') L2=L0 #menuAide.add_command(label="Parabole de Mohr Caquot pour Rec > Re", command=Parabole_MohrCaquot) Ligne3='Selon le matériau, on choisit la contrainte de Von Mises ou la contrainte de Tresca qui est ' def AffMontrerPhi(): ax.grid(True) while n<=Nb: xCerc, yCerc = DemiCercle(0,Reg) #ax.set_xlim(min(Smini,x2), Smaxi) Srz=Srz+[-a/K12*M/2*(1/T-T)/Pmaxi] soit: σx1, σy1, σz1, σx2, σy2, σz2 Chaudière à tubes d'eau. def SaisieRec(event): ##################################################################### if IndCalculsFaits==1: # on a déjà a et b ######################################################################### BmA=1/2*(dRho1**2+dRho2**2+2*dRho1*dRho2*cos(2*phir))**(1/2) messagebox.showinfo('Aide à la saisie',Ligne1+Ligne2+Ligne3+Ligne4+Ligne5+Ligne6) On choisit |R1|>|Rp1| et |R2|>|Rp2| puis l'angle phi entre -90° et 90° entre les plans principaux de R1 et R2 Tai=Tai*Kbord azimSC=30 ax.plot(xCerc, yCerc, linewidth=1, color = orange,label="Autre cas : Von-Mises "+Rapport) RecMaxi=100*(1+round(SeqMaxi2/100)) Tai=T0 yCM13L3 = yCM13L3+[ParaSecu*abs((S1-S3)/2)*sin(t)] def AppelContraintes2(): - Les positions des graphiques dans la fenêtre se règlent dans le programme avec la fonction while c1<=c1Maxi: global epsi xParaLpouH=[] if (IndS1.get() == 1) & (IndS2.get() == 1): ax.set_ylabel(r"$\tau$") Srx2, Sry2, Srz2, Sreq2 = Srx, Sry, Srz, Sreq Cote8 = 1/2*(-1/R2*(0)**2 -1/Rp2* (T0/2)**2) t=t0 # donc on calcul a, il nous faut avant B fonction de chargeF Zmaxi=max(Cote1,Cote2,Cote3,Cote4,Cote5,Cote6,Cote7,Cote8) global xCha,xParaL,xParaLpouH,yFoncI,yFoncJ,yFoncJ,yFoncK,yFoncH,ydelta varMohr_Caquot = IntVar() T=min(abs(R1),abs(Rp1),abs(R2),abs(Rp2)) ra=ra-pasra Types de chaudières à haute pression. if (IndZoom.get() == 1) or (IndCalculsFaits == 0): # on Zoom global varParaRec,varParaRe Nu2=tempo l2 = ax.plot(xPara, yFonZ,label=r"$z_{prof}/b$") if lambdaSol>0.5: """ zCer=zCer+[zC] IndReduitb=1 n=taille yC=c2*cos(t) except: # la chaine n'est pas un nombre, on garde Nb Par exemple, on constate qu'un" xPara=xPara+[eval(ListePb[n])] global elevSC,azimSC, IndContr,elevSCpres,azimSCpres,IndContr Dans le cas d'une liaison pivot ou pivot-glissant, nous sommes en présence d'un contact entre un cylindre mâle et un cylindre femelle. RecMaxi=100*(1+round(SeqMaxi2/100)) ax.set_ylabel('Direction principale 2 : petit axe',fontsize=6) ax.scatter(lambdaSol,SeqMaxi/Pmaxi,marker='+', color = 'red') global chargeF,R1,Rp1,R2,Rp2,E1,Nu1,E2,Nu2 txtGeometrie.grid(row =LiRes+7, column =ColRes+1,sticky='w') menuParam.add_command(label="Abaque d'endurance des aciers (Ordre de grandeur)", command=CourbesSigmaD_Rm) IndS2.set(1) Sx, Sy, Sz = eval(champParaSx.get()), eval(champParaSy.get()), eval(champParaSz.get()) if Nb>=2: try: def AppelCourbeIJK(): close(fig) BpA=1/2*(1/R1+1/Rp1+1/R2+1/Rp2) Yij=Y2[i,j] DessinRigidite() menubar.add_cascade(label="Fichier", menu=menufichierF0) IndVue=0 def SaisieTitreTTT(event): DefDeLafenetreDeSaisie() Yel[i]=b*sin(t[i]) IndEllipse.set(0) La pression diamétrale est un modèle décrivant la pression de contact entre deux pièces en liaison mécanique de type pivot glissant, pivot ou rotule. xCerc, yCerc = DemiCercle((S2+S3)/2,(S2-S3)/2) S1, S2, S3 = S1CM1, S2CM1, S3CM1 legend(prop={'size':8}) - Dessin de l'ellipse avec son orientation. ///////////////// ColRes=1 Trouvé à l'intérieur – Page 428Cette ténacité est donnée par la formule : 1 Oriss . z ID : S. Pd Ž D'après Hertz , l'aplatissement ; par des billes ... de la pression moyenne dans la zone de pression avec des charges variables pour en déterminer la valeur maximum . if IndCalculsFaits==1 : Lignes=Ligne1+Ligne2+Ligne3+Ligne4+Ligne5+Ligne6+Ligne7+Ligne8+Ligne9+Ligne10 #InitialiserInd() global elevSC,azimSC, IndContr,IndVue,elevSCpres,azimSCpres,IndContr IndSprio=0 varphi =StringVar() # définition de la variable Tkinter j=0 global kpg varParaRe.set(str(ReTest)) # affectation de la variable de contrôle #(Tresca pour σ2=σ1 ou σ2=σ3) #################### txtGeometrie = Label(fen1, text = ' Graphiques ', foreground=CouleurText).grid(row =LiGraph, column =ColRes,sticky='w') ###############--------------------------------------------------------- if (IndS1.get() == 1) & (IndS2.get() == 0): if dt!=0: pause(dt) Benson Boiler. if IndCasH.get()==1: b=1 Ligne1='- Sollicitations statiques : \n' champE1 = Entry(fen1,width=15, textvariable =varE1) ax.scatter([x1], [y1],marker = '+', color = 'c') K2=(1-Nu2**2)/E2/pi RecMaxi=400 RecMaxi=100*(1+round(SeqMaxi2/100)) IndS1=1 CouleurFen1='green' txtRec = Label(FeneParaDes, text = 'N/mm²').grid(row =LigMat, column =3,sticky='w') X1=X DessinPression(elevSCpres,azimSCpres) zCer=zCer+[-R2*(1-cos(t))] subplots_adjust(left=0, bottom=0.1, right=0.95, top=0.95, wspace=0.4, hspace=0.35) ax.plot(xCerc, yCerc, linewidth=1, linestyle = 'dashed', color = 'b') TabSz2=[0] ax = fig.add_subplot(2, 2, 2) Y1[i]=str*Xi+ctr*Yi varSmaxi.set(Smaxi) varParaSz.set(S1) y2=Re/sqrt(3) gca(projection='3d').plot(xCer,yCer,zCer,color=coul,linewidth=3) Trouvé à l'intérieur – Page 3003Pression admissible sur les plaques de fermetures circulaires , unies et boulonnées . Balistique . 10-38811 . ... Application des formules de Hertz aux engrenages hélicoïdaux à axes non parallèles et non convergents . Conclusion . if Nb>=2: IndCas.set(0) Ligne8='Documentation utilisée pour le calcul des contraintes suivant la normale au contact.\n\n' Ligne13D="prise en compte des différents types de sollicitations en vue de la construction du diagramme de Haigh (Voir Aides/Bibliographie : CETIM)\n" global R1,Rp1,R2,Rp2,phi,phir,TailleTitre La pression de Hertz, utilisée comme modèle de calcul de la pression de contact, est le principe de base utilisé dans la présente partie de l'ISO 6336 pour l'évaluation de la résistance à la pression superficielle des engrenages cylindriques. CalculHertz() MC1 = Radiobutton(FeneParaDes, text="Critère de la courbe enveloppe de Mohr-Caquot", variable=varMohr_Caquot, value=1, command=ChoixDessin) Seqmax = Tricercle_kTg(S1,S2,S3,'m',"Autre cas : Mohr-Caquot "+Rapport,Nb,kp) """ Joint FPM. l4 = ax.plot(xParaLpouH, yFoncH) else: ax.set_xlim(Smini,Smaxi) H�tTˎ�0��}4Rp� c_��aWY)Q��r `V�x���w�+�_��L6RN�������� N`���d. title("Cercles de Mohr et parabole de Mohr-Caquot", fontsize=TailleTitre+1) \n\n' Vert=Button(fen1,text='CALCUL HERTZ ',fg='green', command=CalculHertz).grid(row =LiMat+6, column =ColMat+2,sticky='e') #txtGeometrie = Label(fen1, text = 'ou (|R1|", SaisiechargeF) def AppelCasGeneralHertz(): Ligne15="On pensera aussi à évaluer l'imprécision due aux tolérances sur les 10 paramètres saisis.\n\n" ax.grid(True) y1=Re/2 orange=(1, 128/255, 0) def CourbesSansDimCalculeesNu03(): ax.scatter(Lp,Sp,marker='+', color = 'red') Test !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Fic.close() #t=linspace(0,2*pi,Nb) Bonne trempabilité à l'huile, bonne résistance aux surcharges à l'état traité. global IndInterferences Nb=eval(champParaNb.get()) def DemiCercle(Cx,R): # abscisse et rayon Ligne14="Consulter un cours pour savoir comment utiliser ces contraintes. Smini=round(1.05*min(Smini,xMin)) y2=Re/2*cosA xCerc, yCerc = DemiCercle(-Ray,Ray) CalculHertz() Merci de signaler les erreurs éventuelles. from numpy.linalg import * Trouvé à l'intérieur – Page 386... il rappelle la terrible éruption de Krakatoa qui montre combien une pareille hypothèse est admissible . ... et interrompue par un micromètre à étincelles réprésente l'axe , M. Hertz a montré qu'une oscillation électrique introduite ... 0000039534 00000 n y2=x*sin(phir)+y*cos(phir) OMpy=-1+T+zr*(F(PH,kp)-E(PH,kp)) ax.scatter([0], [0.7861],marker = '. #ax.plot(xCerc, yCerc, linewidth=2, color = 'k',label="Hertz (Tresca limité à Rec)") menubar.add_cascade(label="Aides", menu=menuAide) IndTor = IntVar() # type de variable show() messagebox.showinfo('Matériaux',Ligne1+Ligne2+Ligne3) #print("Nu=",Nu) R1=eval(ListePb[2]) ax.set_xlabel("lambda") f.close() Sy=a/K12*M*(OMy+Nu*OMpy) global SeqMaxi, IndZoom, IndParabole, IndDroite IndInterferences=0 DessinContraintes() Le matage peut être 1. Smini=round(1.05*min(-RecMaxi,S3,xT)) champNu2.grid(row =LiMat+4, column =ColMat+1) Re=60 Sreq=Sreq+[(S1[Num]-S3[Num])/Pmaxi] Ligne2="il faut modifier les surfaces ou φ pour avoir droit au calcul." show() Ligne3='Pour une courbure nulle, on peut saisir un rayon qui correspond à un défaut de rectitude.\n' Rec += 100 show() titreS="Contrainte équivalente\n"+r"$\sigma_{eq2}$"+" relative à "+ r"$P_{maxi}$" #menuAide.add_command(label="Droite de Mohr dans pour Rec >= Re", command=DroiteDeMohrCaquot) except: Sy=Sy+[a/K12*M*(OMy+Nu*OMpy)] sinA=(Rec-Re)/(Rec+Re) fen1.config(menu=menubar) BmA=1/2*(dRho1**2+dRho2**2+2*dRho1*dRho2*cos(2*phir))**(1/2) Com1.grid(row =LigAnima+2,columnspan=2, column =3,sticky='we') X2=X S1Test=-100 Surfaces pressées. if (IndZoom.get() == 0): if Nb>=2: #Dernier cercle en continu # Sigma/Pmaxi indice r pour relatif à Pmaxi Nu=0.3 elevSC=20 # On remet les valeurs initiales y1=0 txtGeometrie = Label(fen1, text = str('%5.0f' % Pmoy)+' N/mm² ').grid(row =LiRes+3, column =ColRes+1,sticky='w') global elevSC,azimSC, IndContr,IndVue,elevSCpres,azimSCpres,IndContr while i < taille/3: Ligne3="Elles sont pressées par une force suivant la normale commune. " ax.set_xlim(0, 1) B=(BpA+BmA)/2 Trouvé à l'intérieur – Page 123On donne par ailleurs la contrainte admissible R , du matériau constituant la roue dentée . La condition de résistance à la flexion simple ... Calcul à la pression superficielle ( théorie de Hertz ) . IV . APPLICATIONS PREMIER EXEMPLE . global fenPara, dt, Nb, Smini, Smaxi, IndSprio, orange #tinit = time.clock() def I(L): Le second résultat non trivial est que bien que la loi de comportement de la matière soit linéaire (loi de Hooke), la relation entre la force appliquée F {\displaystyle \mathrm {F} } et l'enfoncement δ {\displaystyle \delta } ne l'est pas : La raison essentielle en est justement que l'aire de la surface de contact π a 2 {\displaystyle \pi a^{2}} entre la sphère et le plan augmente au cours de l'enfoncement. if IndSprio==0: #champchargeF.bind("", SaisiechargeF) # n=taille subplots_adjust(left=0.07, bottom=0.05, right=0.95, top=0.95, wspace=0.4, hspace=0.35) xCerc, yCerc = DemiCercle(-Ray,Ray) l1 = ax.plot(xParaL, yFoncI) Seq=[] if (IndZoom.get() == 0): y1=Re/2 DessinContraintes1() Titre=tempo Srz=[]
Ferrari Roma Date De Sortie, Hyperthyroïdie Et Grossesse Recommandations, Broadway Montparnasse Menu, Centre D'inscription Université Lyon 2, Trépied Appareil Photo Et Smartphone,
