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
33
34
35
36
|
import matplotlib.pyplot as plt
import numpy as np
import math
from matplotlib import rcParams
rcParams['font.family'] = 'sans-serif'
rcParams['font.sans-serif'] = ['Yu Gothic', 'Meirio']
k = [0.048, 0.024] # outer, inner
k_label = ["outer", "inner"]
dt = 10
thickness = [18, 35, 70] # [um]
width = np.arange(0.1, 5, 0.1)
current_outer = []
current_inner = []
for i in thickness:
current_outer.append(k[0]*dt**0.44*np.power(i*width*1.55, 0.725))
# 1.55 = 1u*1m*39370.1^2
for i in thickness:
current_inner.append(k[1]*dt**0.44*np.power(i*width*1.55, 0.725))
fig = plt.figure()
ax = fig.add_subplot(111)
ax.set_title("current (delta T = 10degC) vs line width")
for i, c in enumerate(current_outer):
ax.plot(width, c, label="{}_{}um".format(k_label[0], thickness[i]))
for i, c in enumerate(current_inner):
ax.plot(width, c, linestyle = "dashed", label="{}_{}um".format(k_label[1], thickness[i]))
ax.set_xlim(0, 5)
ax.set_ylim(0, 10)
ax.set_xticks(np.arange(0, 5, 0.2), minor=True)
ax.set_yticks(np.arange(0, 10, 0.5), minor=True)
ax.set_xlabel("line width [mm]")
ax.set_ylabel("Current [A]")
ax.legend()
ax.grid(which="both")
|
