cambios menores

Comparar_Resultados/M800-50A.mdb

@@ -124,129 +124,3 @@ alphaValue_Bertotti=3.447
 alphaValue_Bertotti_Maxwell=2
 betavalue_Steinmetz=0
 LaminationThickness=0.5
-[M1000-65A]
-Type=Fixed_Solid
-Solid Type=Steel
-Thermal Conductivity=30
-Specific Heat=460
-Density=7650
-Notes=
-ElectricalResistivity=1.8E-7
-TempCoefElectricalResistivity=0
-PoissonsRatio=0.3
-YoungsCoefficient=210000
-YieldStress=295
-BValue[0]=0
-HValue[0]=0
-BValue[1]=0.1
-HValue[1]=83.3
-BValue[2]=0.2
-HValue[2]=107
-BValue[3]=0.3
-HValue[3]=119
-BValue[4]=0.4
-HValue[4]=130
-BValue[5]=0.5
-HValue[5]=140
-BValue[6]=0.6
-HValue[6]=150
-BValue[7]=0.7
-HValue[7]=160
-BValue[8]=0.8
-HValue[8]=172
-BValue[9]=0.9
-HValue[9]=185
-BValue[10]=1
-HValue[10]=200
-BValue[11]=1.1
-HValue[11]=217
-BValue[12]=1.2
-HValue[12]=237
-BValue[13]=1.3
-HValue[13]=275
-BValue[14]=1.4
-HValue[14]=368
-BValue[15]=1.5
-HValue[15]=604
-BValue[16]=1.6
-HValue[16]=1360
-BValue[17]=1.7
-HValue[17]=3370
-BValue[18]=1.8
-HValue[18]=7010
-BValue[19]=1.9
-HValue[19]=16457.33
-BValue[20]=2
-HValue[20]=30342.8247
-BValue[21]=2.1
-HValue[21]=50000
-BValue[22]=2.2
-HValue[22]=75000
-BValue[23]=2.3
-HValue[23]=120000
-Frequency[0]=50
-LossDensity[0]=0.06
-FluxDensity[0]=0.1
-Frequency[1]=50
-LossDensity[1]=0.24
-FluxDensity[1]=0.2
-Frequency[2]=50
-LossDensity[2]=0.5
-FluxDensity[2]=0.3
-Frequency[3]=50
-LossDensity[3]=0.81
-FluxDensity[3]=0.4
-Frequency[4]=50
-LossDensity[4]=1.16
-FluxDensity[4]=0.5
-Frequency[5]=50
-LossDensity[5]=1.56
-FluxDensity[5]=0.6
-Frequency[6]=50
-LossDensity[6]=1.99
-FluxDensity[6]=0.7
-Frequency[7]=50
-LossDensity[7]=2.47
-FluxDensity[7]=0.8
-Frequency[8]=50
-LossDensity[8]=3.01
-FluxDensity[8]=0.9
-Frequency[9]=50
-LossDensity[9]=3.6
-FluxDensity[9]=1
-Frequency[10]=50
-LossDensity[10]=4.26
-FluxDensity[10]=1.1
-Frequency[11]=50
-LossDensity[11]=4.99
-FluxDensity[11]=1.2
-Frequency[12]=50
-LossDensity[12]=5.87
-FluxDensity[12]=1.3
-Frequency[13]=50
-LossDensity[13]=6.82
-FluxDensity[13]=1.4
-Frequency[14]=50
-LossDensity[14]=7.9
-FluxDensity[14]=1.5
-Frequency[15]=50
-LossDensity[15]=9.05
-FluxDensity[15]=1.6
-Frequency[16]=50
-LossDensity[16]=10
-FluxDensity[16]=1.7
-Frequency[17]=50
-LossDensity[17]=10.8
-FluxDensity[17]=1.8
-KcValue=0.65
-KhValue_Steinmetz=0.06257
-KhValue_Bertotti=0.01137
-KhValue_Bertotti_Maxwell=0.02027
-KeddyValue_Steinmetz=1E-0005
-KexcValue=0.008469
-KexcValue_Maxwell=0.00403
-alphaValue_Steinmetz=1.477
-alphaValue_Bertotti=3.245
-alphaValue_Bertotti_Maxwell=2
-betavalue_Steinmetz=0.2514
-LaminationThickness=0.65

graficar.py

@@ -1,2 +1,40 @@
 import numpy as np
 import matplotlib.pyplot as plt
+import json
+
+
+
+
+rutaPLA = 'Comparar_Resultados/PLAcond.mdb'
+rutaMetal = 'Comparar_Resultados/M800-50A.mdb'
+
+with open("h_values/h_values.json", "r") as f : hPLA = np.array(json.load(f))
+
+muVacio = 4 * np.pi * 1e-7
+bPLA = muVacio * hPLA 
+
+bPLAcond = []
+hPLAcond = []
+
+with open(rutaPLA, 'r') as f:
+    for linea in f.readlines():
+        if linea.startswith("BValue["): bPLAcond.append(float(linea.split("=")[1]))
+        if linea.startswith("HValue["): hPLAcond.append(float(linea.split("=")[1]))
+
+bM800 = []
+hM800 = []
+
+with open(rutaMetal, 'r') as f:
+    for linea in f.readlines():
+        if linea.startswith("BValue["): bM800.append(float(linea.split("=")[1]))
+        if linea.startswith("HValue["): hM800.append(float(linea.split("=")[1]))
+
+print(bM800)
+
+plt.figure()
+plt.plot(hPLA, bPLA, label = 'PLA')
+plt.plot(hPLAcond, bPLAcond, label = 'PLACond')
+plt.plot(hM800, bM800, label = 'M800-50A')
+plt.title('Comparacion entre Curvas para Motor CAD')
+plt.legend()
+plt.show()