質問編集履歴
1
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -100,9 +100,9 @@
|
|
|
100
100
|
|
|
101
101
|
FL = -((ρ*U*A)/2*M)*dCL
|
|
102
102
|
|
|
103
|
-
ML =
|
|
103
|
+
ML = FL/2
|
|
104
|
-
|
|
104
|
+
|
|
105
|
-
eq = sp.Eq( sp.diff(y, z ,2)+(ML/(E*I))*sp.diff(y, z, 1), 0 )
|
|
105
|
+
eq = sp.Eq( sp.diff(y, z ,2)+(ML/(E*I))*(y^2)*sp.diff(y, z, 1), 0 )
|
|
106
106
|
|
|
107
107
|
yrez = sp.dsolve(eq)
|
|
108
108
|
|
|
@@ -219,3 +219,45 @@
|
|
|
219
219
|
IndexError: At least one sheet must be visible
|
|
220
220
|
|
|
221
221
|
```
|
|
222
|
+
|
|
223
|
+
|
|
224
|
+
|
|
225
|
+
追記
|
|
226
|
+
|
|
227
|
+
手計算で解く手法を記載します。
|
|
228
|
+
|
|
229
|
+
初期値とかは適当なので、上記プログラムとは一致してません。
|
|
230
|
+
|
|
231
|
+
|
|
232
|
+
|
|
233
|
+
d²y/dx²+ay²dy/dx=0
|
|
234
|
+
|
|
235
|
+
|
|
236
|
+
|
|
237
|
+
dy/dx=pとおくと
|
|
238
|
+
|
|
239
|
+
d²y/dx²
|
|
240
|
+
|
|
241
|
+
=dp/dx
|
|
242
|
+
|
|
243
|
+
=pdp/dy
|
|
244
|
+
|
|
245
|
+
|
|
246
|
+
|
|
247
|
+
pdp/dy+ay²p=0
|
|
248
|
+
|
|
249
|
+
p(dp/dy+ay²)=0
|
|
250
|
+
|
|
251
|
+
|
|
252
|
+
|
|
253
|
+
p=0のときy=C (Cは任意定数)
|
|
254
|
+
|
|
255
|
+
|
|
256
|
+
|
|
257
|
+
p≠0のとき
|
|
258
|
+
|
|
259
|
+
dp/dy+ay²=0
|
|
260
|
+
|
|
261
|
+
|
|
262
|
+
|
|
263
|
+
p=-(a/3)y³+A
|