数値計算:有限要素法

Size: px
Start display at page:

Download "数値計算:有限要素法"

Transcription

1 ( ) 1 / 61

2 ( ) 2 / 61

3 ( ) 3 / 61

4 P(0) P(x) u(x) P(L) f P(0) P(x) P(L) ( ) 4 / 61

5 L P(x) E(x) A(x) x P(x) P(x) u(x) P(x) u(x) (0 x L) ( ) 5 / 61

6 u(x) 0 L x ( ) 6 / 61

7 P(0) P(L) f d dx ( EA du dx ) = 0 u(0) = 0 E(L)A(L) du dx (L) = f (boundary value problem) ( ) 7 / 61

8 (finite element method; FEM) ( ) ( ) 8 / 61

9 U U = f W L 0 ( ) 2 1 du 2 EA dx, dx W = f u(l) min I = U W, subject to u(0) = 0 ( ) 9 / 61

10 [0, L] 6 h = L/6 Pi Pj 0 xi xj L x (nodal point) x 0 = 0, x 1 = h, x 2 = 2h,, x 6 = L u i = u(xi ) u(x) 7 u 0, u 1, u 2,, u 6 ( ) 10 / 61

11 [x i, x j ] u(x) u(x) = u i N i,j (x) + u j N j,i (x), (x i x x j ) u i, u j x i, x j xi xj L N i,j (x) = x j x { h 1 (x = xi ) = 0 (x = x j ) x 0 0 xi xj L 1 N j,i (x) = x x i { h 1 (x = xj ) = 0 (x = x i ) x ( ) 11 / 61

12 u(x) = u(x) (0 x L) u 0 N 0,1 (x) + u 1 N 1,0 (x) (x 0 x x 1 ) u 1 N 1,2 (x) + u 2 N 2,1 (x) (x 1 x x 2 ) u 2 N 2,3 (x) + u 3 N 3,2 (x) (x 2 x x 3 ). u 5 N 5,6 (x) + u 6 N 6,5 (x) (x 5 x x 6 ) u(x) 7 u 0, u 1,, u 6 ( ) 12 / 61

13 u(x) 0 h 2h 3h 4h 5h L P0 P1 P2 P3 P4 P5 P6 x ( ) 13 / 61

14 u(x) u5 u6 u2 u4 u0 u1 u3 0 h 2h 3h 4h 5h L P0 P1 P2 P3 P4 P5 P6 x ( ) 13 / 61

15 [x i, x j ] N i,j (x) = x j x, N j,i (x) = x x i h h N i,j(x) = 1 h, N j,i(x) = 1 h u(x) = u i N i,j (x) + u j N j,i (x) du = u i N dx i,j(x) + u j N j,i(x) 1 = u i h + u 1 j h = u i + u j h ( ) 14 / 61

16 = 7 u 0, u 1,, u 6 L 0 = u N = u 0 u 1 u 2. u 6 x1 x 0 x2 x 1 x3 + + x x6 x 5 ( ) 15 / 61

17 E xj ( ) 2 1 du xj ( ) 2 x i 2 EA 1 dx = dx x i 2 EA ui + u j dx h = 1 E xj 2 h ( u 2 i + u j ) 2 A dx x i = 1 [ [ ] E V ui u i,j j 2 h 2 V i,j V i,j = xj x i A dx = [x i, x j ] V i,j V i,j ] [ ui u j ], ( ) 16 / 61

18 L ( ) 2 1 du U = 0 2 EA dx dx = 1 [ ] [ [ ] E V u0 u 0,1 V 0,1 u0 1 2 h 2 V 0,1 V 0,1 u 1 [ ] [ 1 [ ] E V u1 u 1,2 V 1,2 u1 2 2 h 2 V 1,2 V 1,2 u 2 [ ] [ 1 [ ] E V u2 u 2,3 V 2,3 u2 3 2 h 2 V 2,3 V 2,3 u 3 [ ] [ 1 [ ] E V u5 u 5,6 V 5,6 u5 6 2 h 2 V 5,6 V 5,6 u 6 ] + ] + ] + + ] ( ) 17 / 61

19 K = E h 2 U = 1 2 u N T Ku N V 0,1 V 0,1 V 0,1 V 0,1 + V 1,2 V 1,2 V 1,2 V 1,2 + V 2,3 V 2,3 V 2, V 4,5 + V 5,6 V 5,6 V 5,6 V 5,6 ( ) 18 / 61

20 K = E h 2 U = 1 2 u N T Ku N V 0,1 V 0,1 V 0,1 V 0,1 + V 1,2 V 1,2 V 1,2 V 1,2 + V 2,3 V 2,3 V 2, V 4,5 + V 5,6 V 5,6 V 5,6 V 5,6 ( ) 18 / 61

21 K = E h 2 U = 1 2 u N T Ku N V 0,1 V 0,1 V 0,1 V 0,1 + V 1,2 V 1,2 V 1,2 V 1,2 + V 2,3 V 2,3 V 2, V 4,5 + V 5,6 V 5,6 V 5,6 V 5,6 ( ) 18 / 61

22 K = E h 2 U = 1 2 u N T Ku N V 0,1 V 0,1 V 0,1 V 0,1 + V 1,2 V 1,2 V 1,2 V 1,2 + V 2,3 V 2,3 V 2, V 4,5 + V 5,6 V 5,6 V 5,6 V 5,6 ( ) 18 / 61

23 A V i,j = Ah K = EA h ( ) 19 / 61

24 A V i,j = Ah K = EA h ( ) 20 / 61

25 A V i,j = Ah K = EA h ( ) 20 / 61

26 f W = f u(l) = f u 6 W = f T u N f = f ( ) 21 / 61

27 u(0) = u 0 = 0 a T u N = 0 a = ( ) 22 / 61

28 min I = U W, subject to u(0) = 0 min I (u N ) = 1 2 ut N K u N f T u N, subject to a T u N = 0 ( ) 23 / 61

29 ( ) min I (u N ) = 1 2 ut N K u N f T u N, subject to a T u N = 0 min J(u N, λ) = I (u N ) λa T u N = 1 2 ut N K u N f T u N λa T u N ( ) 24 / 61

30 ( ) J u N = Ku N f λa = 0 J λ = at u N = 0 [ K a a T 0 ] [ un λ ] = [ f 0 ] ( ) 25 / 61

31 ( ) J u N = Ku N f λa = 0 J λ = at u N = 0 [ K a a T 0 ] [ un λ ] = [ f 0 ] ( ) 25 / 61

32 ( ) J u N = Ku N f λa = 0 J λ = at u N = 0 [ K a a T 0 ] [ un λ ] = [ f 0 ] = u N λ ( ) 25 / 61

33 A(x) = a 2bx (a, b a 2bL > 0 ) E V i,j = xj x i (a 2bx) dx = [ ax bx 2] xj x i = {a b(x i + x j )} h V 0,1 = (a bh)h, V 1,2 = (a 3bh)h, V 2,3 = (a 5bh)h, [0, L] 4 a bh K = E h a + bh a + bh 2a 4bh a + 3bh a + 3bh 2a 8bh a + 5bh a + 5bh 2a 12bh a + 7bh a + 7bh a 7bh ( ) 26 / 61

34 P(x) r(x) = 3(a bx) (a, b a bl > 0 ) A(x) = πr 2 = 3π(a bx) 2 xj ] x=xj V i,j = 3π(a bx) 2 (a bx)3 dx = [π x i b x=x i = π { 3a 2 3ab(x i + x j ) + b 2 (xi 2 + x i x j + xj 2 ) } h [0, L] 4 (x i = ih h = L/4) V 0,1 = π(3a 2 3abh + b 2 h 2 )h V 1,2 = π(3a 2 9abh + 7b 2 h 2 )h V 2,3 = π(3a 2 15abh + 19b 2 h 2 )h V 3,4 = π(3a 2 21abh + 37b 2 h 2 )h ( ) 27 / 61

35 ( ) 28 / 61

36 P(0) P(x) u(x,t) 0 P(L) f(t) P(0) P(x) t P(L) ( ) 29 / 61

37 T = L 0 ( ) 2 1 u L 2 ρa 1 dx = t 0 2 ρa u2 dx 7 u 0, u 1,, u 6 [x i, x j ] u(x, t) u(x, t) = u i (t) N i,j (x) + u j (t) N j,i (x), (x i x x j ) u(x, t) = u i (t) N i,j (x) + u j (t) N j,i (x), (x i x x j ) ( ) 30 / 61

38 xj x i 1 2 ρa u2 dx = 1 2 m i,j; i,j = xj = 1 2 m i,j; j,i = m j,i; i,j = xj ρa { u i (t) N i,j (x) + u j (t) N j,i (x)} 2 dx x i [ ] [ ] [ ] m ui u i,j; i,j m i,j; j,i ui j m j,i; i,j m j,i; j,i u j x i ρa (N i,j ) 2 dx, m j,i; j,i = xj x i ρa N i,j N j,i dx xj x i ρa (N j,i ) 2 dx, ( ) 31 / 61

39 L T = 1 ρa u 2 dx 2 0 = 1 [ ] [ ] [ m u0 u 0,1; 0,1 m 0,1; 1,0 u0 1 2 m 1,0; 0,1 m 1,0; 1,0 u 1 1 [ ] [ ] [ m u1 u 1,2; 1,2 m 1,2; 2,1 u1 2 2 m 2,1; 1,2 m 2,1; 2,1 u 2 1 [ ] [ ] [ m u2 u 2,3; 2,3 m 2,3; 3,2 u2 3 2 m 3,2; 2,3 m 3,2; 3,2 u 3 1 [ ] [ ] [ m u5 u 5,6; 5,6 m 5,6; 6,5 u5 6 2 m 6,5; 5,6 m 6,5; 6,5 u 6 ] + ] + ] + + ] = 1 2 u N T M u N ( ) 32 / 61

40 M = m 0,1; 0,1 m 0,1; 1,0 m 1,0; 0,1 m 1,0; 1,0 + m 1,2; 1,2 m 1,2; 2,1. m 2,1; 1,2.. m4,5; 5,4 m 5,4; 4,5 m 5,4; 5,4 + m 5,6; 5,6 m 5,6; 6,5 m 6,5; 5,6 m 6,5; 6,5 ( ) 33 / 61

41 A ρ xj x i (N i,j ) 2 dx = xj x i N i,j N j,i dx = M = ρah 6 xj x i (N j,i ) 2 dx = h 3 xj x i N j,i N i,j dx = h ( ) 34 / 61

42 A ρ M = ρah ( ) 35 / 61

43 A ρ M = ρah ( ) 36 / 61

44 A ρ M = ρah ( ) 36 / 61

45 A(x) = a 2bx (a, b a 2bL > 0 ) ρ xj x (N i,j ) 2 dx = h x i 12 (3x i + x j ), xj x (N j,i ) 2 dx = h x i 12 (x i + 3x j ), xj x i x N i,j N j,i dx = xj x i x N j,i N i,j dx = h 12 (x i + x j ) m i,j; i,j = ρh 6 {2a b(3x i + x j )}, m j,i; j,i = ρh 6 {2a b(x i + 3x j )}, m i,j; j,i = m j,i; i,j = ρh 6 {a b(x i + x j )} ( ) 37 / 61

46 [0, L] 4 2a bh M = ρh 6 a bh a bh 4a 8bh a 3bh a 3bh 4a 16bh a 5bh a 5bh 4a 24bh a 7bh a 7bh 2a 15bh ( ) 38 / 61

47 L(u N, u N ) = T U + W + λr = 1 2 u N T M u N 1 2 u N T Ku N + f T u N + λa T u N L d L = Ku N + f + λa Mü N = 0 u N dt u N ( ) 39 / 61

48 R = a T u N Ṙ = a T u N R = a T ü N a T ü N + 2νa T u N + ν 2 a T u N = 0 ( ) 40 / 61

49 Mü N λa = Ku N + f a T ü N = a T (2ν u N + ν 2 u N ) u N = v N [ M a a T ] [ vn λ ] [ = Ku N + f a T (2νv N + ν 2 u N ) ] ( ) 41 / 61

50 λ, µ ( ) 42 / 61

51 P i P i P i Pk Pi Pj ( ) 43 / 61

52 P i P i P i Pk Pi Pj ( ) 43 / 61

53 P i P i P i Pk Pi uk ui Pj uj P i P i P i u i = [ ui v i ] [ uj, u j = v j ] [ uk, u k = v k ] ( ) 43 / 61

54 U i,j,k = 1 [ u T 2 i uj T uk T K i,j,k (6 6) ] Ki,j,k u i u j u k K i,j,k = λj i,j,k λ J i,j,k λ, Jµ i,j,k (6 6) + µj i,j,k µ ( ) 44 / 61

55 1 P2 P2 1 1 P3 P0 1 P1 P1 ( ) 45 / 61

56 J 0,1,2 λ = J 3,2,1 λ = 1 2 J 0,1,2 µ = J 3,2,1 µ = ( ) 46 / 61

57 P2 P3 P2 P2 P3 1 = P0 1 P1 P0 1 P1 P1 ( ) 47 / 61

58 U = U 0,1,2 + U 3,2,1 = 1 2 [ u T 0 u T 1 u T 2 u T 3 ] K u 0 u 1 u 2 u 3 K (8 8) J λ, J µ (8 8) K = λj λ + µj µ ( ) 48 / 61

59 J λ (0, 0) J 0,1,2 λ = ( ) 49 / 61

60 J λ (0, 1) J 0,1,2 λ = ( ) 49 / 61

61 J λ (0, 2) J 0,1,2 λ = ( ) 49 / 61

62 J λ (1, 0) J 0,1,2 λ = ( ) 49 / 61

63 J λ (1, 1) J 0,1,2 λ = ( ) 49 / 61

64 J λ (1, 2) J 0,1,2 λ = ( ) 49 / 61

65 J λ (2, 0) J 0,1,2 λ = ( ) 49 / 61

66 J λ (2, 1) J 0,1,2 λ = ( ) 49 / 61

67 J λ (2, 2) J 0,1,2 λ = ( ) 49 / 61

68 J 0,1,2 λ J λ J λ = ( ) 50 / 61

69 J λ (3, 3) J 3,2,1 λ = ( ) 51 / 61

70 J λ (3, 2) J 3,2,1 λ = ( ) 51 / 61

71 J λ (3, 1) J 3,2,1 λ = ( ) 51 / 61

72 J λ (2, 3) J 3,2,1 λ = ( ) 51 / 61

73 J λ (2, 2) J 3,2,1 λ = ( ) 51 / 61

74 J λ (2, 1) J 3,2,1 λ = ( ) 51 / 61

75 J λ (1, 3) J 3,2,1 λ = ( ) 51 / 61

76 J λ (1, 2) J 3,2,1 λ = ( ) 51 / 61

77 J λ (1, 1) J 3,2,1 λ = ( ) 51 / 61

78 J 3,2,1 λ J λ J λ = ( ) 52 / 61

79 J 0,1,2 λ J 3,2,1 λ J λ = ( ) 53 / 61

80 J 0,1,2 µ J µ J µ = ( ) 54 / 61

81 J 3,2,1 µ J µ J µ = ( ) 55 / 61

82 J 0,1,2 µ J 1,3,2 µ J µ = ( ) 56 / 61

83 U = 1 2 u N T Ku N, T = 1 2 u N T M u N = = ( ) 57 / 61

84 x = x 1 x 2 x 3 y x y x y 1 = y x x 2 y x 3 ( ) 58 / 61

85 b = b 1 b 2 b 3 y = b T x = b 1 x 1 + b 2 x 2 + b 3 x 3 y x 1 = b 1, y x 2 = b 2, y x 3 = b 3 ( b T x ) = b x ( ) 59 / 61

86 A = a 11 a 12 a 13 a 12 a 22 a 23 a 13 a 23 a 33 y = 1 2 x T Ax = 1 2 { a11 x a 22 x a 33 x a 12 x 1 x 2 + 2a 13 x 1 x 3 + 2a 23 x 2 x 3 } y x 1 = 1 2 {a 11 2x 1 + 2a 12 x 2 + 2a 13 x 3 } = a 11 x 1 + a 12 x 2 + a 13 x 3 ( ) 60 / 61

87 y x 1 = a 11 x 1 + a 12 x 2 + a 13 x 3 y x 2 = a 12 x 1 + a 22 x 2 + a 23 x 3 y x 3 = a 13 x 1 + a 23 x 2 + a 33 x 3 ( ) 1 x 2 x T Ax = Ax ( ) 61 / 61

9. 05 L x P(x) P(0) P(x) u(x) u(x) (0 < = x < = L) P(x) E(x) A(x) P(L) f ( d EA du ) = 0 (9.) dx dx u(0) = 0 (9.2) E(L)A(L) du (L) = f (9.3) dx (9.) P

9. 05 L x P(x) P(0) P(x) u(x) u(x) (0 < = x < = L) P(x) E(x) A(x) P(L) f ( d EA du ) = 0 (9.) dx dx u(0) = 0 (9.2) E(L)A(L) du (L) = f (9.3) dx (9.) P 9 (Finite Element Method; FEM) 9. 9. P(0) P(x) u(x) (a) P(L) f P(0) P(x) (b) 9. P(L) 9. 05 L x P(x) P(0) P(x) u(x) u(x) (0 < = x < = L) P(x) E(x) A(x) P(L) f ( d EA du ) = 0 (9.) dx dx u(0) = 0 (9.2) E(L)A(L)

More information

1 4 1 ( ) ( ) ( ) ( ) () 1 4 2

1 4 1 ( ) ( ) ( ) ( ) () 1 4 2 7 1995, 2017 7 21 1 2 2 3 3 4 4 6 (1).................................... 6 (2)..................................... 6 (3) t................. 9 5 11 (1)......................................... 11 (2)

More information

73

73 73 74 ( u w + bw) d = Ɣ t tw dɣ u = N u + N u + N 3 u 3 + N 4 u 4 + [K ] {u = {F 75 u δu L σ (L) σ dx σ + dσ x δu b δu + d(δu) ALW W = L b δu dv + Aσ (L)δu(L) δu = (= ) W = A L b δu dx + Aσ (L)δu(L) Aσ

More information

Part () () Γ Part ,

Part () () Γ Part , Contents a 6 6 6 6 6 6 6 7 7. 8.. 8.. 8.3. 8 Part. 9. 9.. 9.. 3. 3.. 3.. 3 4. 5 4.. 5 4.. 9 4.3. 3 Part. 6 5. () 6 5.. () 7 5.. 9 5.3. Γ 3 6. 3 6.. 3 6.. 3 6.3. 33 Part 3. 34 7. 34 7.. 34 7.. 34 8. 35

More information

all.dvi

all.dvi 5,, Euclid.,..,... Euclid,.,.,, e i (i =,, ). 6 x a x e e e x.:,,. a,,. a a = a e + a e + a e = {e, e, e } a (.) = a i e i = a i e i (.) i= {a,a,a } T ( T ),.,,,,. (.),.,...,,. a 0 0 a = a 0 + a + a 0

More information

A A = a 41 a 42 a 43 a 44 A (7) 1 (3) A = M 12 = = a 41 (8) a 41 a 43 a 44 (3) n n A, B a i AB = A B ii aa

A A = a 41 a 42 a 43 a 44 A (7) 1 (3) A = M 12 = = a 41 (8) a 41 a 43 a 44 (3) n n A, B a i AB = A B ii aa 1 2 21 2 2 [ ] a 11 a 12 A = a 21 a 22 (1) A = a 11 a 22 a 12 a 21 (2) 3 3 n n A A = n ( 1) i+j a ij M ij i =1 n (3) j=1 M ij A i j (n 1) (n 1) 2-1 3 3 A A = a 11 a 12 a 13 a 21 a 22 a 23 a 31 a 32 a 33

More information

W u = u(x, t) u tt = a 2 u xx, a > 0 (1) D := {(x, t) : 0 x l, t 0} u (0, t) = 0, u (l, t) = 0, t 0 (2)

W u = u(x, t) u tt = a 2 u xx, a > 0 (1) D := {(x, t) : 0 x l, t 0} u (0, t) = 0, u (l, t) = 0, t 0 (2) 3 215 4 27 1 1 u u(x, t) u tt a 2 u xx, a > (1) D : {(x, t) : x, t } u (, t), u (, t), t (2) u(x, ) f(x), u(x, ) t 2, x (3) u(x, t) X(x)T (t) u (1) 1 T (t) a 2 T (t) X (x) X(x) α (2) T (t) αa 2 T (t) (4)

More information

AHPを用いた大相撲の新しい番付編成

AHPを用いた大相撲の新しい番付編成 5304050 2008/2/15 1 2008/2/15 2 42 2008/2/15 3 2008/2/15 4 195 2008/2/15 5 2008/2/15 6 i j ij >1 ij ij1/>1 i j i 1 ji 1/ j ij 2008/2/15 7 1 =2.01/=0.5 =1.51/=0.67 2008/2/15 8 1 2008/2/15 9 () u ) i i i

More information

gr09.dvi

gr09.dvi .1, θ, ϕ d = A, t dt + B, t dtd + C, t d + D, t dθ +in θdϕ.1.1 t { = f1,t t = f,t { D, t = B, t =.1. t A, tdt e φ,t dt, C, td e λ,t d.1.3,t, t d = e φ,t dt + e λ,t d + dθ +in θdϕ.1.4 { = f1,t t = f,t {

More information

1 u t = au (finite difference) u t = au Von Neumann

1 u t = au (finite difference) u t = au Von Neumann 1 u t = au 3 1.1 (finite difference)............................. 3 1.2 u t = au.................................. 3 1.3 Von Neumann............... 5 1.4 Von Neumann............... 6 1.5............................

More information

5. [1 ] 1 [], u(x, t) t c u(x, t) x (5.3) ξ x + ct, η x ct (5.4),u(x, t) ξ, η u(ξ, η), ξ t,, ( u(ξ,η) ξ η u(x, t) t ) u(x, t) { ( u(ξ, η) c t ξ ξ { (

5. [1 ] 1 [], u(x, t) t c u(x, t) x (5.3) ξ x + ct, η x ct (5.4),u(x, t) ξ, η u(ξ, η), ξ t,, ( u(ξ,η) ξ η u(x, t) t ) u(x, t) { ( u(ξ, η) c t ξ ξ { ( 5 5.1 [ ] ) d f(t) + a d f(t) + bf(t) : f(t) 1 dt dt ) u(x, t) c u(x, t) : u(x, t) t x : ( ) ) 1 : y + ay, : y + ay + by : ( ) 1 ) : y + ay, : yy + ay 3 ( ): ( ) ) : y + ay, : y + ay b [],,, [ ] au xx

More information

n ξ n,i, i = 1,, n S n ξ n,i n 0 R 1,.. σ 1 σ i .10.14.15 0 1 0 1 1 3.14 3.18 3.19 3.14 3.14,. ii 1 1 1.1..................................... 1 1............................... 3 1.3.........................

More information

A

A A 2563 15 4 21 1 3 1.1................................................ 3 1.2............................................. 3 2 3 2.1......................................... 3 2.2............................................

More information

6.1 (P (P (P (P (P (P (, P (, P.

6.1 (P (P (P (P (P (P (, P (, P. (011 30 7 0 ( ( 3 ( 010 1 (P.3 1 1.1 (P.4.................. 1 1. (P.4............... 1 (P.15.1 (P.16................. (P.0............3 (P.18 3.4 (P.3............... 4 3 (P.9 4 3.1 (P.30........... 4 3.

More information

4 2016 3 8 2.,. 2. Arakawa Jacobin., 2 Adams-Bashforth. Re = 80, 90, 100.. h l, h/l, Kármán, h/l 0.28,, h/l.., (2010), 46.2., t = 100 t = 2000 46.2 < Re 46.5. 1 1 4 2 6 2.1............................

More information

II A A441 : October 02, 2014 Version : Kawahira, Tomoki TA (Kondo, Hirotaka )

II A A441 : October 02, 2014 Version : Kawahira, Tomoki TA (Kondo, Hirotaka ) II 214-1 : October 2, 214 Version : 1.1 Kawahira, Tomoki TA (Kondo, Hirotaka ) http://www.math.nagoya-u.ac.jp/~kawahira/courses/14w-biseki.html pdf 1 2 1 9 1 16 1 23 1 3 11 6 11 13 11 2 11 27 12 4 12 11

More information

ct_root.dvi

ct_root.dvi 7 1 1 3 ( ) (i) ( ) (ii) ( ) (iii) 1.1 2 1.2 1.3 1.4 1.3 1.4 1.1 1.1.1 ( ) ( ) 1) 1 x x X R X =[0, x] x b(x) b( ) 2 b(0) = 0 x 0 b (0) = + b (x) =0 x 0 b (x) < 0 x 2 ( ) Θ Θ ={θ 0,θ 1 } 0

More information

all.dvi

all.dvi 29 4 Green-Lagrange,,.,,,,,,.,,,,,,,,,, E, σ, ε σ = Eε,,.. 4.1? l, l 1 (l 1 l) ε ε = l 1 l l (4.1) F l l 1 F 30 4 Green-Lagrange Δz Δδ γ = Δδ (4.2) Δz π/2 φ γ = π 2 φ (4.3) γ tan γ γ,sin γ γ ( π ) γ tan

More information

all.dvi

all.dvi 72 9 Hooke,,,. Hooke. 9.1 Hooke 1 Hooke. 1, 1 Hooke. σ, ε, Young. σ ε (9.1), Young. τ γ G τ Gγ (9.2) X 1, X 2. Poisson, Poisson ν. ν ε 22 (9.) ε 11 F F X 2 X 1 9.1: Poisson 9.1. Hooke 7 Young Poisson G

More information

nakata/nakata.html p.1/20

nakata/nakata.html p.1/20 http://www.me.titech.ac.jp/ nakata/nakata.html p.1/20 1-(a). Faybusovich(1997) Linear systems in Jordan algebras and primal-dual interior-point algorithms,, Euclid Jordan p.2/20 Euclid Jordan V Euclid

More information

ばらつき抑制のための確率最適制御

ばらつき抑制のための確率最適制御 ( ) http://wwwhayanuemnagoya-uacjp/ fujimoto/ 2011 3 9 11 ( ) 2011/03/09-11 1 / 46 Outline 1 2 3 4 5 ( ) 2011/03/09-11 2 / 46 Outline 1 2 3 4 5 ( ) 2011/03/09-11 3 / 46 (1/2) r + Controller - u Plant y

More information

III 1 (X, d) d U d X (X, d). 1. (X, d).. (i) d(x, y) d(z, y) d(x, z) (ii) d(x, y) d(z, w) d(x, z) + d(y, w) 2. (X, d). F X.. (1), X F, (2) F 1, F 2 F

III 1 (X, d) d U d X (X, d). 1. (X, d).. (i) d(x, y) d(z, y) d(x, z) (ii) d(x, y) d(z, w) d(x, z) + d(y, w) 2. (X, d). F X.. (1), X F, (2) F 1, F 2 F III 1 (X, d) d U d X (X, d). 1. (X, d).. (i) d(x, y) d(z, y) d(x, z) (ii) d(x, y) d(z, w) d(x, z) + d(y, w) 2. (X, d). F X.. (1), X F, (2) F 1, F 2 F F 1 F 2 F, (3) F λ F λ F λ F. 3., A λ λ A λ. B λ λ

More information

kou05.dvi

kou05.dvi 2 C () 25 1 3 1.1........................................ 3 1.2..................................... 4 1.3..................................... 7 1.3.1................................ 7 1.3.2.................................

More information

変 位 変位とは 物体中のある点が変形後に 別の点に異動したときの位置の変化で あり ベクトル量である 変位には 物体の変形の他に剛体運動 剛体変位 が含まれている 剛体変位 P(x, y, z) 平行移動と回転 P! (x + u, y + v, z + w) Q(x + d x, y + dy,

変 位 変位とは 物体中のある点が変形後に 別の点に異動したときの位置の変化で あり ベクトル量である 変位には 物体の変形の他に剛体運動 剛体変位 が含まれている 剛体変位 P(x, y, z) 平行移動と回転 P! (x + u, y + v, z + w) Q(x + d x, y + dy, 変 位 変位とは 物体中のある点が変形後に 別の点に異動したときの位置の変化で あり ベクトル量である 変位には 物体の変形の他に剛体運動 剛体変位 が含まれている 剛体変位 P(x, y, z) 平行移動と回転 P! (x + u, y + v, z + w) Q(x + d x, y + dy, z + dz) Q! (x + d x + u + du, y + dy + v + dv, z +

More information

直交座標系の回転

直交座標系の回転 b T.Koama x l x, Lx i ij j j xi i i i, x L T L L, L ± x L T xax axx, ( a a ) i, j ij i j ij ji λ λ + λ + + λ i i i x L T T T x ( L) L T xax T ( T L T ) A( L) T ( LAL T ) T ( L AL) λ ii L AL Λ λi i axx

More information

6.1 (P (P (P (P (P (P (, P (, P.101

6.1 (P (P (P (P (P (P (, P (, P.101 (008 0 3 7 ( ( ( 00 1 (P.3 1 1.1 (P.3.................. 1 1. (P.4............... 1 (P.15.1 (P.15................. (P.18............3 (P.17......... 3.4 (P................ 4 3 (P.7 4 3.1 ( P.7...........

More information

1 1 u m (t) u m () exp [ (cπm + (πm κ)t (5). u m (), U(x, ) f(x) m,, (4) U(x, t) Re u k () u m () [ u k () exp(πkx), u k () exp(πkx). f(x) exp[ πmxdx

1 1 u m (t) u m () exp [ (cπm + (πm κ)t (5). u m (), U(x, ) f(x) m,, (4) U(x, t) Re u k () u m () [ u k () exp(πkx), u k () exp(πkx). f(x) exp[ πmxdx 1 1 1 1 1. U(x, t) U(x, t) + c t x c, κ. (1). κ U(x, t) x. (1) 1, f(x).. U(x, t) U(x, t) + c κ U(x, t), t x x : U(, t) U(1, t) ( x 1), () : U(x, ) f(x). (3) U(x, t). [ U(x, t) Re u k (t) exp(πkx). (4)

More information

106 4 4.1 1 25.1 25.4 20.4 17.9 21.2 23.1 26.2 1 24 12 14 18 36 42 24 10 5 15 120 30 15 20 10 25 35 20 18 30 12 4.1 7 min. z = 602.5x 1 + 305.0x 2 + 2

106 4 4.1 1 25.1 25.4 20.4 17.9 21.2 23.1 26.2 1 24 12 14 18 36 42 24 10 5 15 120 30 15 20 10 25 35 20 18 30 12 4.1 7 min. z = 602.5x 1 + 305.0x 2 + 2 105 4 0 1? 1 LP 0 1 4.1 4.1.1 (intger programming problem) 1 0.5 x 1 = 447.7 448 / / 2 1.1.2 1. 2. 1000 3. 40 4. 20 106 4 4.1 1 25.1 25.4 20.4 17.9 21.2 23.1 26.2 1 24 12 14 18 36 42 24 10 5 15 120 30

More information

3. :, c, ν. 4. Burgers : t + c x = ν 2 u x 2, (3), ν. 5. : t + u x = ν 2 u x 2, (4), c. 2 u t 2 = c2 2 u x 2, (5) (1) (4), (1 Navier Stokes,., ν. t +

3. :, c, ν. 4. Burgers : t + c x = ν 2 u x 2, (3), ν. 5. : t + u x = ν 2 u x 2, (4), c. 2 u t 2 = c2 2 u x 2, (5) (1) (4), (1 Navier Stokes,., ν. t + B: 2016 12 2, 9, 16, 2017 1 6 1,.,,,,.,.,,,., 1,. 1. :, ν. 2. : t = ν 2 u x 2, (1), c. t + c x = 0, (2). e-mail: iwayama@kobe-u.ac.jp,. 1 3. :, c, ν. 4. Burgers : t + c x = ν 2 u x 2, (3), ν. 5. : t +

More information

6 2 2 x y x y t P P = P t P = I P P P ( ) ( ) ,, ( ) ( ) cos θ sin θ cos θ sin θ, sin θ cos θ sin θ cos θ y x θ x θ P

6 2 2 x y x y t P P = P t P = I P P P ( ) ( ) ,, ( ) ( ) cos θ sin θ cos θ sin θ, sin θ cos θ sin θ cos θ y x θ x θ P 6 x x 6.1 t P P = P t P = I P P P 1 0 1 0,, 0 1 0 1 cos θ sin θ cos θ sin θ, sin θ cos θ sin θ cos θ x θ x θ P x P x, P ) = t P x)p ) = t x t P P ) = t x = x, ) 6.1) x = Figure 6.1 Px = x, P=, θ = θ P

More information

< 1 > (1) f 0 (a) =6a ; g 0 (a) =6a 2 (2) y = f(x) x = 1 f( 1) = 3 ( 1) 2 =3 ; f 0 ( 1) = 6 ( 1) = 6 ; ( 1; 3) 6 x =1 f(1) = 3 ; f 0 (1) = 6 ; (1; 3)

< 1 > (1) f 0 (a) =6a ; g 0 (a) =6a 2 (2) y = f(x) x = 1 f( 1) = 3 ( 1) 2 =3 ; f 0 ( 1) = 6 ( 1) = 6 ; ( 1; 3) 6 x =1 f(1) = 3 ; f 0 (1) = 6 ; (1; 3) < 1 > (1) f 0 (a) =6a ; g 0 (a) =6a 2 (2) y = f(x) x = 1 f( 1) = 3 ( 1) 2 =3 ; f 0 ( 1) = 6 ( 1) = 6 ; ( 1; 3) 6 x =1 f(1) = 3 ; f 0 (1) = 6 ; (1; 3) 6 y = g(x) x = 1 g( 1) = 2 ( 1) 3 = 2 ; g 0 ( 1) =

More information

1 Ricci V, V i, W f : V W f f(v ) = Imf W ( ) f : V 1 V k W 1

1 Ricci V, V i, W f : V W f f(v ) = Imf W ( ) f : V 1 V k W 1 1 Ricci V, V i, W f : V W f f(v = Imf W ( f : V 1 V k W 1 {f(v 1,, v k v i V i } W < Imf > < > f W V, V i, W f : U V L(U; V f : V 1 V r W L(V 1,, V r ; W L(V 1,, V r ; W (f + g(v 1,, v r = f(v 1,, v r

More information

,. Black-Scholes u t t, x c u 0 t, x x u t t, x c u t, x x u t t, x + σ x u t, x + rx ut, x rux, t 0 x x,,.,. Step 3, 7,,, Step 6., Step 4,. Step 5,,.

,. Black-Scholes u t t, x c u 0 t, x x u t t, x c u t, x x u t t, x + σ x u t, x + rx ut, x rux, t 0 x x,,.,. Step 3, 7,,, Step 6., Step 4,. Step 5,,. 9 α ν β Ξ ξ Γ γ o δ Π π ε ρ ζ Σ σ η τ Θ θ Υ υ ι Φ φ κ χ Λ λ Ψ ψ µ Ω ω Def, Prop, Th, Lem, Note, Remark, Ex,, Proof, R, N, Q, C [a, b {x R : a x b} : a, b {x R : a < x < b} : [a, b {x R : a x < b} : a,

More information

( )/2 hara/lectures/lectures-j.html 2, {H} {T } S = {H, T } {(H, H), (H, T )} {(H, T ), (T, T )} {(H, H), (T, T )} {1

( )/2   hara/lectures/lectures-j.html 2, {H} {T } S = {H, T } {(H, H), (H, T )} {(H, T ), (T, T )} {(H, H), (T, T )} {1 ( )/2 http://www2.math.kyushu-u.ac.jp/ hara/lectures/lectures-j.html 1 2011 ( )/2 2 2011 4 1 2 1.1 1 2 1 2 3 4 5 1.1.1 sample space S S = {H, T } H T T H S = {(H, H), (H, T ), (T, H), (T, T )} (T, H) S

More information

simx simxdx, cosxdx, sixdx 6.3 px m m + pxfxdx = pxf x p xf xdx = pxf x p xf x + p xf xdx 7.4 a m.5 fx simxdx 8 fx fx simxdx = πb m 9 a fxdx = πa a =

simx simxdx, cosxdx, sixdx 6.3 px m m + pxfxdx = pxf x p xf xdx = pxf x p xf x + p xf xdx 7.4 a m.5 fx simxdx 8 fx fx simxdx = πb m 9 a fxdx = πa a = II 6 ishimori@phys.titech.ac.jp 6.. 5.4.. f Rx = f Lx = fx fx + lim = lim x x + x x f c = f x + x < c < x x x + lim x x fx fx x x = lim x x f c = f x x < c < x cosmx cosxdx = {cosm x + cosm + x} dx = [

More information

1 1.1 / Fik Γ= D n x / Newton Γ= µ vx y / Fouie Q = κ T x 1. fx, tdx t x x + dx f t = D f x 1 fx, t = 1 exp x 4πDt 4Dt lim fx, t =δx 3 t + dxfx, t = 1

1 1.1 / Fik Γ= D n x / Newton Γ= µ vx y / Fouie Q = κ T x 1. fx, tdx t x x + dx f t = D f x 1 fx, t = 1 exp x 4πDt 4Dt lim fx, t =δx 3 t + dxfx, t = 1 1 1.1......... 1............. 1.3... 1.4......... 1.5.............. 1.6................ Bownian Motion.1.......... Einstein.............. 3.3 Einstein........ 3.4..... 3.5 Langevin Eq.... 3.6................

More information

ランダムウォークの境界条件・偏微分方程式の数値計算

ランダムウォークの境界条件・偏微分方程式の数値計算 B L06(2018-05-22 Tue) : Time-stamp: 2018-05-22 Tue 21:53 JST hig,, 2, multiply transf http://hig3.net L06 B(2018) 1 / 38 L05-Q1 Quiz : 1 M λ 1 = 1 u 1 ( ). M u 1 = u 1, u 1 = ( 3 4 ) s (s 0)., u 1 = 1

More information

18 ( ) I II III A B C(100 ) 1, 2, 3, 5 I II A B (100 ) 1, 2, 3 I II A B (80 ) 6 8 I II III A B C(80 ) 1 n (1 + x) n (1) n C 1 + n C

18 ( ) I II III A B C(100 ) 1, 2, 3, 5 I II A B (100 ) 1, 2, 3 I II A B (80 ) 6 8 I II III A B C(80 ) 1 n (1 + x) n (1) n C 1 + n C 8 ( ) 8 5 4 I II III A B C( ),,, 5 I II A B ( ),, I II A B (8 ) 6 8 I II III A B C(8 ) n ( + x) n () n C + n C + + n C n = 7 n () 7 9 C : y = x x A(, 6) () A C () C P AP Q () () () 4 A(,, ) B(,, ) C(,,

More information

,, Mellor 1973),, Mellor and Yamada 1974) Mellor 1973), Mellor and Yamada 1974) 4 2 3, 2 4,

,, Mellor 1973),, Mellor and Yamada 1974) Mellor 1973), Mellor and Yamada 1974) 4 2 3, 2 4, Mellor and Yamada1974) The Turbulence Closure Model of Mellor and Yamada 1974) Kitamori Taichi 2004/01/30 ,, Mellor 1973),, Mellor and Yamada 1974) Mellor 1973), 4 1 4 Mellor and Yamada 1974) 4 2 3, 2

More information

f(x) = f(x ) + α(x)(x x ) α(x) x = x. x = f (y), x = f (y ) y = f f (y) = f f (y ) + α(f (y))(f (y) f (y )) f (y) = f (y ) + α(f (y)) (y y ) ( (2) ) f

f(x) = f(x ) + α(x)(x x ) α(x) x = x. x = f (y), x = f (y ) y = f f (y) = f f (y ) + α(f (y))(f (y) f (y )) f (y) = f (y ) + α(f (y)) (y y ) ( (2) ) f 22 A 3,4 No.3 () (2) (3) (4), (5) (6) (7) (8) () n x = (x,, x n ), = (,, n ), x = ( (x i i ) 2 ) /2 f(x) R n f(x) = f() + i α i (x ) i + o( x ) α,, α n g(x) = o( x )) lim x g(x) x = y = f() + i α i(x )

More information

第5章 偏微分方程式の境界値問題

第5章 偏微分方程式の境界値問題 October 5, 2018 1 / 113 4 ( ) 2 / 113 Poisson 5.1 Poisson ( A.7.1) Poisson Poisson 1 (A.6 ) Γ p p N u D Γ D b 5.1.1: = Γ D Γ N 3 / 113 Poisson 5.1.1 d {2, 3} Lipschitz (A.5 ) Γ D Γ N = \ Γ D Γ p Γ N Γ

More information

Sturm-Liouville Green KEN ZOU Hermite Legendre Laguerre L L [p(x) d2 dx 2 + q(x) d ] dx + r(x) u(x) = Lu(x) = 0 (1) L = p(x) d2 dx

Sturm-Liouville Green KEN ZOU Hermite Legendre Laguerre L L [p(x) d2 dx 2 + q(x) d ] dx + r(x) u(x) = Lu(x) = 0 (1) L = p(x) d2 dx Sturm-Liouville Green KEN ZOU 2006 4 23 1 Hermite Legendre Lguerre 1 1.1 2 L L p(x) d2 2 + q(x) d + r(x) u(x) = Lu(x) = 0 (1) L = p(x) d2 2 + q(x) d + r(x) (2) L = d2 2 p(x) d q(x) + r(x) (3) 2 (Self-Adjoint

More information

(3) (2),,. ( 20) ( s200103) 0.7 x C,, x 2 + y 2 + ax = 0 a.. D,. D, y C, C (x, y) (y 0) C m. (2) D y = y(x) (x ± y 0), (x, y) D, m, m = 1., D. (x 2 y

(3) (2),,. ( 20) ( s200103) 0.7 x C,, x 2 + y 2 + ax = 0 a.. D,. D, y C, C (x, y) (y 0) C m. (2) D y = y(x) (x ± y 0), (x, y) D, m, m = 1., D. (x 2 y [ ] 7 0.1 2 2 + y = t sin t IC ( 9) ( s090101) 0.2 y = d2 y 2, y = x 3 y + y 2 = 0 (2) y + 2y 3y = e 2x 0.3 1 ( y ) = f x C u = y x ( 15) ( s150102) [ ] y/x du x = Cexp f(u) u (2) x y = xey/x ( 16) ( s160101)

More information

D v D F v/d F v D F η v D (3.2) (a) F=0 (b) v=const. D F v Newtonian fluid σ ė σ = ηė (2.2) ė kl σ ij = D ijkl ė kl D ijkl (2.14) ė ij (3.3) µ η visco

D v D F v/d F v D F η v D (3.2) (a) F=0 (b) v=const. D F v Newtonian fluid σ ė σ = ηė (2.2) ė kl σ ij = D ijkl ė kl D ijkl (2.14) ė ij (3.3) µ η visco post glacial rebound 3.1 Viscosity and Newtonian fluid f i = kx i σ ij e kl ideal fluid (1.9) irreversible process e ij u k strain rate tensor (3.1) v i u i / t e ij v F 23 D v D F v/d F v D F η v D (3.2)

More information

Untitled

Untitled II 14 14-7-8 8/4 II (http://www.damp.tottori-u.ac.jp/~ooshida/edu/fluid/) [ (3.4)] Navier Stokes [ 6/ ] Navier Stokes 3 [ ] Reynolds [ (4.6), (45.8)] [ p.186] Navier Stokes I 1 balance law t (ρv i )+ j

More information

211 kotaro@math.titech.ac.jp 1 R *1 n n R n *2 R n = {(x 1,..., x n ) x 1,..., x n R}. R R 2 R 3 R n R n R n D D R n *3 ) (x 1,..., x n ) f(x 1,..., x n ) f D *4 n 2 n = 1 ( ) 1 f D R n f : D R 1.1. (x,

More information

TOP URL 1

TOP URL   1 TOP URL http://amonphys.web.fc2.com/ 1 30 3 30.1.............. 3 30.2........................... 4 30.3...................... 5 30.4........................ 6 30.5.................................. 8 30.6...............................

More information

4. ϵ(ν, T ) = c 4 u(ν, T ) ϵ(ν, T ) T ν π4 Planck dx = 0 e x 1 15 U(T ) x 3 U(T ) = σt 4 Stefan-Boltzmann σ 2π5 k 4 15c 2 h 3 = W m 2 K 4 5.

4. ϵ(ν, T ) = c 4 u(ν, T ) ϵ(ν, T ) T ν π4 Planck dx = 0 e x 1 15 U(T ) x 3 U(T ) = σt 4 Stefan-Boltzmann σ 2π5 k 4 15c 2 h 3 = W m 2 K 4 5. A 1. Boltzmann Planck u(ν, T )dν = 8πh ν 3 c 3 kt 1 dν h 6.63 10 34 J s Planck k 1.38 10 23 J K 1 Boltzmann u(ν, T ) T ν e hν c = 3 10 8 m s 1 2. Planck λ = c/ν Rayleigh-Jeans u(ν, T )dν = 8πν2 kt dν c

More information

untitled

untitled - k k k = y. k = ky. y du dx = ε ux ( ) ux ( ) = ax+ b x u() = ; u( ) = AE u() = b= u () = a= ; a= d x du ε x = = = dx dx N = σ da = E ε da = EA ε A x A x x - σ x σ x = Eε x N = EAε x = EA = N = EA k =

More information

II ( ) (7/31) II ( [ (3.4)] Navier Stokes [ (6/29)] Navier Stokes 3 [ (6/19)] Re

II ( ) (7/31) II (  [ (3.4)] Navier Stokes [ (6/29)] Navier Stokes 3 [ (6/19)] Re II 29 7 29-7-27 ( ) (7/31) II (http://www.damp.tottori-u.ac.jp/~ooshida/edu/fluid/) [ (3.4)] Navier Stokes [ (6/29)] Navier Stokes 3 [ (6/19)] Reynolds [ (4.6), (45.8)] [ p.186] Navier Stokes I Euler Navier

More information

Microsoft Word - 表紙.docx

Microsoft Word - 表紙.docx 黒住英司 [ 著 ] サピエンティア 計量経済学 訂正および練習問題解答 (206/2/2 版 ) 訂正 練習問題解答 3 .69, 3.8 4 (X i X)U i i i (X i μ x )U i ( X μx ) U i. i E [ ] (X i μ x )U i i E[(X i μ x )]E[U i ]0. i V [ ] (X i μ x )U i i 2 i j E [(X i

More information

1. 1 BASIC PC BASIC BASIC BASIC Fortran WS PC (1.3) 1 + x 1 x = x = (1.1) 1 + x = (1.2) 1 + x 1 = (1.

1. 1 BASIC PC BASIC BASIC BASIC Fortran WS PC (1.3) 1 + x 1 x = x = (1.1) 1 + x = (1.2) 1 + x 1 = (1. Section Title Pages Id 1 3 7239 2 4 7239 3 10 7239 4 8 7244 5 13 7276 6 14 7338 7 8 7338 8 7 7445 9 11 7580 10 10 7590 11 8 7580 12 6 7395 13 z 11 7746 14 13 7753 15 7 7859 16 8 7942 17 8 Id URL http://km.int.oyo.co.jp/showdocumentdetailspage.jsp?documentid=

More information

II 2 3.,, A(B + C) = AB + AC, (A + B)C = AC + BC. 4. m m A, m m B,, m m B, AB = BA, A,, I. 5. m m A, m n B, AB = B, A I E, 4 4 I, J, K

II 2 3.,, A(B + C) = AB + AC, (A + B)C = AC + BC. 4. m m A, m m B,, m m B, AB = BA, A,, I. 5. m m A, m n B, AB = B, A I E, 4 4 I, J, K II. () 7 F 7 = { 0,, 2, 3, 4, 5, 6 }., F 7 a, b F 7, a b, F 7,. (a) a, b,,. (b) 7., 4 5 = 20 = 2 7 + 6, 4 5 = 6 F 7., F 7,., 0 a F 7, ab = F 7 b F 7. (2) 7, 6 F 6 = { 0,, 2, 3, 4, 5 },,., F 6., 0 0 a F

More information

構造と連続体の力学基礎

構造と連続体の力学基礎 II 37 Wabash Avenue Bridge, Illinois 州 Winnipeg にある歩道橋 Esplanade Riel 橋6 6 斜張橋である必要は多分無いと思われる すぐ横に道路用桁橋有り しかも塔基部のレストランは 8 年には営業していなかった 9 9. 9.. () 97 [3] [5] k 9. m w(t) f (t) = f (t) + mg k w(t) Newton

More information

名称未設定

名称未設定 ! = ( u v w = u i u = u 1 u u 3 u = ( u 1 u = ( u v = u i! 11! 1! 13 % T = $! 1!! 3 ' =! ij #! 31! 3! 33 & 1 0 0%! ij = $ 0 1 0 ' # 0 0 1& # % 1! ijm = $ 1 & % 0 (i, j,m = (1,,3, (, 3,1, (3,1, (i, j,m

More information

: 1g99p038-8

: 1g99p038-8 16 17 : 1g99p038-8 1 3 1.1....................................... 4 1................................... 5 1.3.................................. 5 6.1..................................... 7....................................

More information

23 7 28 i i 1 1 1.1................................... 2 1.2............................... 3 1.2.1.................................... 3 1.2.2............................... 4 1.2.3 SI..............................

More information

微分方程式の解を見る

微分方程式の解を見る ( ) norikazu[at]ms.u-tokyo.ac.jp http://www.infsup.jp/saito/materials/829gairon.pdf 28 2 9 NS ( ) 28 2 9 / 5 Newton 2 3 Navier-Stokes 4 5 6 NS ( ) 28 2 9 2 / 5 2 x 2 + bx + c = x = b ± b 2 4c 2 3 x 3 +

More information

lecture

lecture 5 3 3. 9. 4. x, x. 4, f(x, ) :=x x + =4,x,.. 4 (, 3) (, 5) (3, 5), (4, 9) 95 9 (g) 4 6 8 (cm).9 3.8 6. 8. 9.9 Phsics 85 8 75 7 65 7 75 8 85 9 95 Mathematics = ax + b 6 3 (, 3) 3 ( a + b). f(a, b) ={3 (a

More information

, = = 7 6 = 42, =

, = = 7 6 = 42, = http://www.ss.u-tokai.ac.jp/~mahoro/2016autumn/alg_intro/ 1 1 2016.9.26, http://www.ss.u-tokai.ac.jp/~mahoro/2016autumn/alg_intro/ 1.1 1 214 132 = 28258 2 + 1 + 4 1 + 3 + 2 = 7 6 = 42, 4 + 2 = 6 2 + 8

More information

2013 25 9 i 1 1 1.1................................... 1 1.2........................... 2 1.3..................................... 3 1.4..................................... 4 2 6 2.1.................................

More information

Sample function Re random process Flutter, Galloping, etc. ensemble (mean value) N 1 µ = lim xk( t1) N k = 1 N autocorrelation function N 1 R( t1, t1

Sample function Re random process Flutter, Galloping, etc. ensemble (mean value) N 1 µ = lim xk( t1) N k = 1 N autocorrelation function N 1 R( t1, t1 Sample function Re random process Flutter, Galloping, etc. ensemble (mean value) µ = lim xk( k = autocorrelation function R( t, t + τ) = lim ( ) ( + τ) xk t xk t k = V p o o R p o, o V S M R realization

More information

第1章 微分方程式と近似解法

第1章 微分方程式と近似解法 April 12, 2018 1 / 52 1.1 ( ) 2 / 52 1.2 1.1 1.1: 3 / 52 1.3 Poisson Poisson Poisson 1 d {2, 3} 4 / 52 1 1.3.1 1 u,b b(t,x) u(t,x) x=0 1.1: 1 a x=l 1.1 1 (0, t T ) (0, l) 1 a b : (0, t T ) (0, l) R, u

More information

keisoku01.dvi

keisoku01.dvi 2.,, Mon, 2006, 401, SAGA, JAPAN Dept. of Mechanical Engineering, Saga Univ., JAPAN 4 Mon, 2006, 401, SAGA, JAPAN Dept. of Mechanical Engineering, Saga Univ., JAPAN 5 Mon, 2006, 401, SAGA, JAPAN Dept.

More information

n ( (

n ( ( 1 2 27 6 1 1 m-mat@mathscihiroshima-uacjp 2 http://wwwmathscihiroshima-uacjp/~m-mat/teach/teachhtml 2 1 3 11 3 111 3 112 4 113 n 4 114 5 115 5 12 7 121 7 122 9 123 11 124 11 125 12 126 2 2 13 127 15 128

More information

TOP URL 1

TOP URL   1 TOP URL http://amonphys.web.fc2.com/ 1 6 3 6.1................................ 3 6.2.............................. 4 6.3................................ 5 6.4.......................... 6 6.5......................

More information

B ver B

B ver B B ver. 2017.02.24 B Contents 1 11 1.1....................... 11 1.1.1............. 11 1.1.2.......................... 12 1.2............................. 14 1.2.1................ 14 1.2.2.......................

More information

20 9 19 1 3 11 1 3 111 3 112 1 4 12 6 121 6 122 7 13 7 131 8 132 10 133 10 134 12 14 13 141 13 142 13 143 15 144 16 145 17 15 19 151 1 19 152 20 2 21 21 21 211 21 212 1 23 213 1 23 214 25 215 31 22 33

More information

IV (2)

IV (2) COMPUTATIONAL FLUID DYNAMICS (CFD) IV (2) The Analysis of Numerical Schemes (2) 11. Iterative methods for algebraic systems Reima Iwatsu, e-mail : iwatsu@cck.dendai.ac.jp Winter Semester 2007, Graduate

More information

3. :, c, ν. 4. Burgers : u t + c u x = ν 2 u x 2, (3), ν. 5. : u t + u u x = ν 2 u x 2, (4), c. 2 u t 2 = c2 2 u x 2, (5) (1) (4), (1 Navier Stokes,.,

3. :, c, ν. 4. Burgers : u t + c u x = ν 2 u x 2, (3), ν. 5. : u t + u u x = ν 2 u x 2, (4), c. 2 u t 2 = c2 2 u x 2, (5) (1) (4), (1 Navier Stokes,., B:,, 2017 12 1, 8, 15, 22 1,.,,,,.,.,,,., 1,. 1. :, ν. 2. : u t = ν 2 u x 2, (1), c. u t + c u x = 0, (2), ( ). 1 3. :, c, ν. 4. Burgers : u t + c u x = ν 2 u x 2, (3), ν. 5. : u t + u u x = ν 2 u x 2,

More information

ii 3.,. 4. F. (), ,,. 8.,. 1. (75%) (25%) =7 20, =7 21 (. ). 1.,, (). 3.,. 1. ().,.,.,.,.,. () (12 )., (), 0. 2., 1., 0,.

ii 3.,. 4. F. (), ,,. 8.,. 1. (75%) (25%) =7 20, =7 21 (. ). 1.,, (). 3.,. 1. ().,.,.,.,.,. () (12 )., (), 0. 2., 1., 0,. 24(2012) (1 C106) 4 11 (2 C206) 4 12 http://www.math.is.tohoku.ac.jp/~obata,.,,,.. 1. 2. 3. 4. 5. 6. 7.,,. 1., 2007 (). 2. P. G. Hoel, 1995. 3... 1... 2.,,. ii 3.,. 4. F. (),.. 5... 6.. 7.,,. 8.,. 1. (75%)

More information

numb.dvi

numb.dvi 11 Poisson kanenko@mbkniftycom alexeikanenko@docomonejp http://wwwkanenkocom/ , u = f, ( u = u+f u t, u = f t ) 1 D R 2 L 2 (D) := {f(x,y) f(x,y) 2 dxdy < )} D D f,g L 2 (D) (f,g) := f(x,y)g(x,y)dxdy (L

More information

p = mv p x > h/4π λ = h p m v Ψ 2 Ψ

p = mv p x > h/4π λ = h p m v Ψ 2 Ψ II p = mv p x > h/4π λ = h p m v Ψ 2 Ψ Ψ Ψ 2 0 x P'(x) m d 2 x = mω 2 x = kx = F(x) dt 2 x = cos(ωt + φ) mω 2 = k ω = m k v = dx = -ωsin(ωt + φ) dt = d 2 x dt 2 0 y v θ P(x,y) θ = ωt + φ ν = ω [Hz] 2π

More information

n Y 1 (x),..., Y n (x) 1 W (Y 1 (x),..., Y n (x)) 0 W (Y 1 (x),..., Y n (x)) = Y 1 (x)... Y n (x) Y 1(x)... Y n(x) (x)... Y n (n 1) (x) Y (n 1)

n Y 1 (x),..., Y n (x) 1 W (Y 1 (x),..., Y n (x)) 0 W (Y 1 (x),..., Y n (x)) = Y 1 (x)... Y n (x) Y 1(x)... Y n(x) (x)... Y n (n 1) (x) Y (n 1) D d dx 1 1.1 n d n y a 0 dx n + a d n 1 y 1 dx n 1 +... + a dy n 1 dx + a ny = f(x)...(1) dk y dx k = y (k) a 0 y (n) + a 1 y (n 1) +... + a n 1 y + a n y = f(x)...(2) (2) (2) f(x) 0 a 0 y (n) + a 1 y

More information

untitled

untitled 40 4 4.3 I (1) I f (x) C [x 0 x 1 ] f (x) f 1 (x) (3.18) f (x) ; f 1 (x) = 1! f 00 ((x))(x ; x 0 )(x ; x 1 ) (x 0 (x) x 1 ) (4.8) (3.18) x (x) x 0 x 1 Z x1 f (x) dx ; Z x1 f 1 (x) dx = Z x1 = 1 1 Z x1

More information

Holton semigeostrophic semigeostrophic,.., Φ(x, y, z, t) = (p p 0 )/ρ 0, Θ = θ θ 0,,., p 0 (z), θ 0 (z).,,,, Du Dt fv + Φ x Dv Φ + fu +

Holton semigeostrophic semigeostrophic,.., Φ(x, y, z, t) = (p p 0 )/ρ 0, Θ = θ θ 0,,., p 0 (z), θ 0 (z).,,,, Du Dt fv + Φ x Dv Φ + fu + Holton 9.2.2 semigeostrophic 1 9.2.2 semigeostrophic,.., Φ(x, y, z, t) = (p p 0 )/ρ 0, Θ = θ θ 0,,., p 0 (z), θ 0 (z).,,,, Du Dt fv + Φ x Dv Φ + fu + Dt DΘ Dt + w dθ 0 dz = 0, (9.2) = 0, (9.3) = 0, (9.4)

More information

x,, z v = (, b, c) v v 2 + b 2 + c 2 x,, z 1 i = (1, 0, 0), j = (0, 1, 0), k = (0, 0, 1) v 1 = ( 1, b 1, c 1 ), v 2 = ( 2, b 2, c 2 ) v

x,, z v = (, b, c) v v 2 + b 2 + c 2 x,, z 1 i = (1, 0, 0), j = (0, 1, 0), k = (0, 0, 1) v 1 = ( 1, b 1, c 1 ), v 2 = ( 2, b 2, c 2 ) v 12 -- 1 4 2009 9 4-1 4-2 4-3 4-4 4-5 4-6 4-7 4-8 4-9 4-10 c 2011 1/(13) 4--1 2009 9 3 x,, z v = (, b, c) v v 2 + b 2 + c 2 x,, z 1 i = (1, 0, 0), j = (0, 1, 0), k = (0, 0, 1) v 1 = ( 1, b 1, c 1 ), v 2

More information

y π π O π x 9 s94.5 y dy dx. y = x + 3 y = x logx + 9 s9.6 z z x, z y. z = xy + y 3 z = sinx y 9 s x dx π x cos xdx 9 s93.8 a, fx = e x ax,. a =

y π π O π x 9 s94.5 y dy dx. y = x + 3 y = x logx + 9 s9.6 z z x, z y. z = xy + y 3 z = sinx y 9 s x dx π x cos xdx 9 s93.8 a, fx = e x ax,. a = [ ] 9 IC. dx = 3x 4y dt dy dt = x y u xt = expλt u yt λ u u t = u u u + u = xt yt 6 3. u = x, y, z = x + y + z u u 9 s9 grad u ux, y, z = c c : grad u = u x i + u y j + u k i, j, k z x, y, z grad u v =

More information

( ) sin 1 x, cos 1 x, tan 1 x sin x, cos x, tan x, arcsin x, arccos x, arctan x. π 2 sin 1 x π 2, 0 cos 1 x π, π 2 < tan 1 x < π 2 1 (1) (

( ) sin 1 x, cos 1 x, tan 1 x sin x, cos x, tan x, arcsin x, arccos x, arctan x. π 2 sin 1 x π 2, 0 cos 1 x π, π 2 < tan 1 x < π 2 1 (1) ( 6 20 ( ) sin, cos, tan sin, cos, tan, arcsin, arccos, arctan. π 2 sin π 2, 0 cos π, π 2 < tan < π 2 () ( 2 2 lim 2 ( 2 ) ) 2 = 3 sin (2) lim 5 0 = 2 2 0 0 2 2 3 3 4 5 5 2 5 6 3 5 7 4 5 8 4 9 3 4 a 3 b

More information

05-5.dvi

05-5.dvi 131 71 71 71 71 71 71 71 7 1 71 71 71 71 71 7 1 71 71 71 71 71 71 7 1 71 71 71 71 71 71 71 71 7 1 1 71 71 71 71 71 71 71 71 71 7 1 75(468) 1 71 71 7 517.95 1 7.1 7. 1 71 71 71 71 71 71 71 7, 1 7.1 7. 1

More information

(1.2) T D = 0 T = D = 30 kn 1.2 (1.4) 2F W = 0 F = W/2 = 300 kn/2 = 150 kn 1.3 (1.9) R = W 1 + W 2 = = 1100 N. (1.9) W 2 b W 1 a = 0

(1.2) T D = 0 T = D = 30 kn 1.2 (1.4) 2F W = 0 F = W/2 = 300 kn/2 = 150 kn 1.3 (1.9) R = W 1 + W 2 = = 1100 N. (1.9) W 2 b W 1 a = 0 1 1 1.1 1.) T D = T = D = kn 1. 1.4) F W = F = W/ = kn/ = 15 kn 1. 1.9) R = W 1 + W = 6 + 5 = 11 N. 1.9) W b W 1 a = a = W /W 1 )b = 5/6) = 5 cm 1.4 AB AC P 1, P x, y x, y y x 1.4.) P sin 6 + P 1 sin 45

More information

x () g(x) = f(t) dt f(x), F (x) 3x () g(x) g (x) f(x), F (x) (3) h(x) = x 3x tf(t) dt.9 = {(x, y) ; x, y, x + y } f(x, y) = xy( x y). h (x) f(x), F (x

x () g(x) = f(t) dt f(x), F (x) 3x () g(x) g (x) f(x), F (x) (3) h(x) = x 3x tf(t) dt.9 = {(x, y) ; x, y, x + y } f(x, y) = xy( x y). h (x) f(x), F (x [ ] IC. f(x) = e x () f(x) f (x) () lim f(x) lim f(x) x + x (3) lim f(x) lim f(x) x + x (4) y = f(x) ( ) ( s46). < a < () a () lim a log xdx a log xdx ( ) n (3) lim log k log n n n k=.3 z = log(x + y ),

More information

24 I ( ) 1. R 3 (i) C : x 2 + y 2 1 = 0 (ii) C : y = ± 1 x 2 ( 1 x 1) (iii) C : x = cos t, y = sin t (0 t 2π) 1.1. γ : [a, b] R n ; t γ(t) = (x

24 I ( ) 1. R 3 (i) C : x 2 + y 2 1 = 0 (ii) C : y = ± 1 x 2 ( 1 x 1) (iii) C : x = cos t, y = sin t (0 t 2π) 1.1. γ : [a, b] R n ; t γ(t) = (x 24 I 1.1.. ( ) 1. R 3 (i) C : x 2 + y 2 1 = 0 (ii) C : y = ± 1 x 2 ( 1 x 1) (iii) C : x = cos t, y = sin t (0 t 2π) 1.1. γ : [a, b] R n ; t γ(t) = (x 1 (t), x 2 (t),, x n (t)) ( ) ( ), γ : (i) x 1 (t),

More information

tokei01.dvi

tokei01.dvi 2. :,,,. :.... Apr. - Jul., 26FY Dept. of Mechanical Engineering, Saga Univ., JAPAN 4 3. (probability),, 1. : : n, α A, A a/n. :, p, p Apr. - Jul., 26FY Dept. of Mechanical Engineering, Saga Univ., JAPAN

More information

) ] [ h m x + y + + V x) φ = Eφ 1) z E = i h t 13) x << 1) N n n= = N N + 1) 14) N n n= = N N + 1)N + 1) 6 15) N n 3 n= = 1 4 N N + 1) 16) N n 4

) ] [ h m x + y + + V x) φ = Eφ 1) z E = i h t 13) x << 1) N n n= = N N + 1) 14) N n n= = N N + 1)N + 1) 6 15) N n 3 n= = 1 4 N N + 1) 16) N n 4 1. k λ ν ω T v p v g k = π λ ω = πν = π T v p = λν = ω k v g = dω dk 1) ) 3) 4). p = hk = h λ 5) E = hν = hω 6) h = h π 7) h =6.6618 1 34 J sec) hc=197.3 MeV fm = 197.3 kev pm= 197.3 ev nm = 1.97 1 3 ev

More information

IA hara@math.kyushu-u.ac.jp Last updated: January,......................................................................................................................................................................................

More information

Part y mx + n mt + n m 1 mt n + n t m 2 t + mn 0 t m 0 n 18 y n n a 7 3 ; x α α 1 7α +t t 3 4α + 3t t x α x α y mx + n

Part y mx + n mt + n m 1 mt n + n t m 2 t + mn 0 t m 0 n 18 y n n a 7 3 ; x α α 1 7α +t t 3 4α + 3t t x α x α y mx + n Part2 47 Example 161 93 1 T a a 2 M 1 a 1 T a 2 a Point 1 T L L L T T L L T L L L T T L L T detm a 1 aa 2 a 1 2 + 1 > 0 11 T T x x M λ 12 y y x y λ 2 a + 1λ + a 2 2a + 2 0 13 D D a + 1 2 4a 2 2a + 2 a

More information

() Remrk I = [0, ] [x i, x i ]. (x : ) f(x) = 0 (x : ) ξ i, (f) = f(ξ i )(x i x i ) = (x i x i ) = ξ i, (f) = f(ξ i )(x i x i ) = 0 (f) 0.

() Remrk I = [0, ] [x i, x i ]. (x : ) f(x) = 0 (x : ) ξ i, (f) = f(ξ i )(x i x i ) = (x i x i ) = ξ i, (f) = f(ξ i )(x i x i ) = 0 (f) 0. () 6 f(x) [, b] 6. Riemnn [, b] f(x) S f(x) [, b] (Riemnn) = x 0 < x < x < < x n = b. I = [, b] = {x,, x n } mx(x i x i ) =. i [x i, x i ] ξ i n (f) = f(ξ i )(x i x i ) i=. (ξ i ) (f) 0( ), ξ i, S, ε >

More information

: , 2.0, 3.0, 2.0, (%) ( 2.

: , 2.0, 3.0, 2.0, (%) ( 2. 2017 1 2 1.1...................................... 2 1.2......................................... 4 1.3........................................... 10 1.4................................. 14 1.5..........................................

More information

Microsoft Word - 11問題表紙(選択).docx

Microsoft Word - 11問題表紙(選択).docx A B A.70g/cm 3 B.74g/cm 3 B C 70at% %A C B at% 80at% %B 350 C γ δ y=00 x-y ρ l S ρ C p k C p ρ C p T ρ l t l S S ξ S t = ( k T ) ξ ( ) S = ( k T) ( ) t y ξ S ξ / t S v T T / t = v T / y 00 x v S dy dx

More information

2 1 1 α = a + bi(a, b R) α (conjugate) α = a bi α (absolute value) α = a 2 + b 2 α (norm) N(α) = a 2 + b 2 = αα = α 2 α (spure) (trace) 1 1. a R aα =

2 1 1 α = a + bi(a, b R) α (conjugate) α = a bi α (absolute value) α = a 2 + b 2 α (norm) N(α) = a 2 + b 2 = αα = α 2 α (spure) (trace) 1 1. a R aα = 1 1 α = a + bi(a, b R) α (conjugate) α = a bi α (absolute value) α = a + b α (norm) N(α) = a + b = αα = α α (spure) (trace) 1 1. a R aα = aα. α = α 3. α + β = α + β 4. αβ = αβ 5. β 0 6. α = α ( ) α = α

More information

difgeo1.dvi

difgeo1.dvi 1 http://matlab0.hwe.oita-u.ac.jp/ matsuo/difgeo.pdf ver.1 8//001 1 1.1 a A. O 1 e 1 ; e ; e e 1 ; e ; e x 1 ;x ;x e 1 ; e ; e X x x x 1 ;x ;x X (x 1 ;x ;x ) 1 1 x x X e e 1 O e x x 1 x x = x 1 e 1 + x

More information

m dv = mg + kv2 dt m dv dt = mg k v v m dv dt = mg + kv2 α = mg k v = α 1 e rt 1 + e rt m dv dt = mg + kv2 dv mg + kv 2 = dt m dv α 2 + v 2 = k m dt d

m dv = mg + kv2 dt m dv dt = mg k v v m dv dt = mg + kv2 α = mg k v = α 1 e rt 1 + e rt m dv dt = mg + kv2 dv mg + kv 2 = dt m dv α 2 + v 2 = k m dt d m v = mg + kv m v = mg k v v m v = mg + kv α = mg k v = α e rt + e rt m v = mg + kv v mg + kv = m v α + v = k m v (v α (v + α = k m ˆ ( v α ˆ αk v = m v + α ln v α v + α = αk m t + C v α v + α = e αk m

More information

f : R R f(x, y) = x + y axy f = 0, x + y axy = 0 y 直線 x+y+a=0 に漸近し 原点で交叉する美しい形をしている x +y axy=0 X+Y+a=0 o x t x = at 1 + t, y = at (a > 0) 1 + t f(x, y

f : R R f(x, y) = x + y axy f = 0, x + y axy = 0 y 直線 x+y+a=0 に漸近し 原点で交叉する美しい形をしている x +y axy=0 X+Y+a=0 o x t x = at 1 + t, y = at (a > 0) 1 + t f(x, y 017 8 10 f : R R f(x) = x n + x n 1 + 1, f(x) = sin 1, log x x n m :f : R n R m z = f(x, y) R R R R, R R R n R m R n R m R n R m f : R R f (x) = lim h 0 f(x + h) f(x) h f : R n R m m n M Jacobi( ) m n

More information

SO(2)

SO(2) TOP URL http://amonphys.web.fc2.com/ 1 12 3 12.1.................................. 3 12.2.......................... 4 12.3............................. 5 12.4 SO(2).................................. 6

More information

Introduction to Numerical Analysis of Differential Equations Naoya Enomoto (Kyoto.univ.Dept.Science(math))

Introduction to Numerical Analysis of Differential Equations Naoya Enomoto (Kyoto.univ.Dept.Science(math)) Introduction to Numerical Analysis of Differential Equations Naoya Enomoto (Kyoto.univ.Dept.Science(math)) 2001 1 e-mail:s00x0427@ip.media.kyoto-u.ac.jp 1 1 Van der Pol 1 1 2 2 Bergers 2 KdV 2 1 5 1.1........................................

More information

A

A A04-164 2008 2 13 1 4 1.1.......................................... 4 1.2..................................... 4 1.3..................................... 4 1.4..................................... 5 2

More information

() x + y + y + x dy dx = 0 () dy + xy = x dx y + x y ( 5) ( s55906) 0.7. (). 5 (). ( 6) ( s6590) 0.8 m n. 0.9 n n A. ( 6) ( s6590) f A (λ) = det(a λi)

() x + y + y + x dy dx = 0 () dy + xy = x dx y + x y ( 5) ( s55906) 0.7. (). 5 (). ( 6) ( s6590) 0.8 m n. 0.9 n n A. ( 6) ( s6590) f A (λ) = det(a λi) 0. A A = 4 IC () det A () A () x + y + z = x y z X Y Z = A x y z ( 5) ( s5590) 0. a + b + c b c () a a + b + c c a b a + b + c 0 a b c () a 0 c b b c 0 a c b a 0 0. A A = 7 5 4 5 0 ( 5) ( s5590) () A ()

More information