質問編集履歴
5
補足情報の更新
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -1,4 +1,4 @@
|
|
|
1
|
-
|
|
1
|
+
実現したいこと
|
|
2
2
|
|
|
3
3
|
|
|
4
4
|
|
4
補足情報の更新
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -1,4 +1,4 @@
|
|
|
1
|
-
##
|
|
1
|
+
##現したいこと
|
|
2
2
|
|
|
3
3
|
|
|
4
4
|
|
|
@@ -177,3 +177,7 @@
|
|
|
177
177
|
コード中の式は、工学系論文の式を引用しています。
|
|
178
178
|
|
|
179
179
|
引用した論文では、最終的にC_ttの値は、上記の変数の時におよそ1.124E-12の結果となるようです。
|
|
180
|
+
|
|
181
|
+
|
|
182
|
+
|
|
183
|
+
|
3
補足情報の更新
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -172,4 +172,8 @@
|
|
|
172
172
|
|
|
173
173
|
コードを全文掲載しました。
|
|
174
174
|
|
|
175
|
+
|
|
176
|
+
|
|
177
|
+
コード中の式は、工学系論文の式を引用しています。
|
|
178
|
+
|
|
175
|
-
最終的に
|
|
179
|
+
引用した論文では、最終的にC_ttの値は、上記の変数の時におよそ1.124E-12の結果となるようです。
|
2
コード全文を掲載
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -44,17 +44,113 @@
|
|
|
44
44
|
|
|
45
45
|
### 該当のソースコード
|
|
46
46
|
|
|
47
|
+
```Python
|
|
47
48
|
|
|
49
|
+
import numpy as np
|
|
50
|
+
|
|
51
|
+
from scipy import integrate
|
|
52
|
+
|
|
53
|
+
|
|
54
|
+
|
|
55
|
+
# 変数を定義
|
|
56
|
+
|
|
57
|
+
delta_in = 0.15E-3
|
|
58
|
+
|
|
59
|
+
delta_out = 1.00E-3
|
|
60
|
+
|
|
61
|
+
delta_tc = 0.6E-3
|
|
62
|
+
|
|
63
|
+
Do = 1.03E-3
|
|
64
|
+
|
|
65
|
+
Di = 1.00E-3
|
|
66
|
+
|
|
67
|
+
eps_0 = 8.854187E-12
|
|
68
|
+
|
|
69
|
+
eps_r = 4.0
|
|
70
|
+
|
|
71
|
+
h_core = 12.3E-3
|
|
72
|
+
|
|
73
|
+
dia_in_core = 17.6E-3
|
|
74
|
+
|
|
75
|
+
dia_out_core = 32.8E-3
|
|
76
|
+
|
|
77
|
+
|
|
78
|
+
|
|
79
|
+
r_in = (dia_in_core/2) - delta_tc - Do/2
|
|
80
|
+
|
|
81
|
+
r_out = (dia_out_core/2) + delta_tc + Do/2
|
|
82
|
+
|
|
83
|
+
|
|
84
|
+
|
|
85
|
+
|
|
86
|
+
|
|
87
|
+
# 積分１
|
|
48
88
|
|
|
49
89
|
y_in = lambda phi: np.sin(phi)/(phi*(delta_in+Do*(1-np.cos(phi))))
|
|
50
90
|
|
|
51
91
|
y_out = lambda phi: np.sin(phi)/(phi*(delta_out+Do*(1-np.cos(phi))))
|
|
52
92
|
|
|
53
|
-
|
|
54
|
-
|
|
55
93
|
integ_in = integrate.quad(y_in, (-1/2)*np.pi, (1/2)*np.pi)
|
|
56
94
|
|
|
57
95
|
integ_out = integrate.quad(y_out, (-1/2)*np.pi, (1/2)*np.pi)
|
|
96
|
+
|
|
97
|
+
|
|
98
|
+
|
|
99
|
+
L_T_in = h_core + delta_tc + Do
|
|
100
|
+
|
|
101
|
+
L_T_out = L_T_in
|
|
102
|
+
|
|
103
|
+
L_T = h_core*2+ (dia_out_core-dia_in_core) + Do
|
|
104
|
+
|
|
105
|
+
|
|
106
|
+
|
|
107
|
+
|
|
108
|
+
|
|
109
|
+
|
|
110
|
+
|
|
111
|
+
C_tt_air_in = (Do/2.0)*eps_0*L_T_in * integ_in[0]
|
|
112
|
+
|
|
113
|
+
C_tt_air_out = (Do/2.0)*eps_0*L_T_out * integ_out[0]
|
|
114
|
+
|
|
115
|
+
|
|
116
|
+
|
|
117
|
+
|
|
118
|
+
|
|
119
|
+
|
|
120
|
+
|
|
121
|
+
# 積分２
|
|
122
|
+
|
|
123
|
+
f = lambda phi, r: np.sin(phi)/(phi*( ((delta_out-delta_in)/(r_out-r_in))*(r-r_in + delta_in + Do*(1-np.cos(phi)))))
|
|
124
|
+
|
|
125
|
+
val, err = integrate.dblquad(f, r_in, r_out, lambda phi: (-1/2)*np.pi, lambda phi: (1/2)*np.pi)
|
|
126
|
+
|
|
127
|
+
|
|
128
|
+
|
|
129
|
+
C_tt_air_l = (Do/2.0)*eps_0 * val
|
|
130
|
+
|
|
131
|
+
print("C_tt_air_l", C_tt_air_l)
|
|
132
|
+
|
|
133
|
+
|
|
134
|
+
|
|
135
|
+
C_tt_ins = (np.pi*eps_0*eps_r*L_T)/(2*np.log(Do/Di))
|
|
136
|
+
|
|
137
|
+
print("C_tt_ins", C_tt_ins)
|
|
138
|
+
|
|
139
|
+
|
|
140
|
+
|
|
141
|
+
C_tt_in = (C_tt_ins*C_tt_air_in)/(C_tt_ins+C_tt_air_in)
|
|
142
|
+
|
|
143
|
+
C_tt_out = (C_tt_ins*C_tt_air_out)/(C_tt_ins + C_tt_air_out)
|
|
144
|
+
|
|
145
|
+
C_tt_l = (C_tt_ins*C_tt_air_l)/(C_tt_ins + C_tt_air_l)
|
|
146
|
+
|
|
147
|
+
|
|
148
|
+
|
|
149
|
+
C_tt = C_tt_in + C_tt_out + 2*C_tt_l
|
|
150
|
+
|
|
151
|
+
```
|
|
152
|
+
|
|
153
|
+
|
|
58
154
|
|
|
59
155
|
|
|
60
156
|
|
|
@@ -74,4 +170,6 @@
|
|
|
74
170
|
|
|
75
171
|
|
|
76
172
|
|
|
77
|
-
|
|
173
|
+
コードを全文掲載しました。
|
|
174
|
+
|
|
175
|
+
最終的に、C_ttの値は、上記の変数の時にはおよそ1.124E-12の結果にしたいです。
|
1
試したことを更新
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -64,7 +64,9 @@
|
|
|
64
64
|
|
|
65
65
|
|
|
66
66
|
|
|
67
|
+
積分範囲を、0からnp.pi、(-1/2)*np.piから0などにすると積分できましたが、
|
|
68
|
+
|
|
67
|
-
|
|
69
|
+
(-1/2)*np.piから(1/2)*np.piに関しては上記のようなエラーが出ます。
|
|
68
70
|
|
|
69
71
|
|
|
70
72
|
|