CAD Galileo 973 Zienkievictz and Campbell 2) 994 3) 4) 7) 3 8) W R n, n 2, 3 T (X): X W x W Q T D (X) T (X) DV O( D ) () V T (X) Eule

Similar documents
応力とひずみ.ppt


untitled

1

( š ) œ 525, , , , ,000 85, , ,810 70,294 4,542,050 18,804,052 () 178,710 1,385, , ,792 72,547 80,366

untitled

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


untitled

( ) œ ,475, ,037 4,230,000 4,224,310 4,230,000 4,230,000 3,362,580 2,300, , , , , , ,730 64,250 74

„¤‰ƒ‰IŠv‚æ‡S−ª†{“Å‘IB5-97

第8章 位相最適化問題

<4D F736F F D2088CF88F589EF8E9197BF F690EC816A2E646F63>

š š o š» p š î å ³å š š n š š š» š» š ½Ò š ˆ l ˆ š p î å ³å š î å» ³ ì š š î å š o š š ½ ñ š å š š n n å š» š m ³ n š

untitled

untitled

ロシア語便覧 1

官報(号外第196号)

Ë,, ÌÓ ÏÓÈ ÂÈ? ÚÓÚ, ÚÓÚ


untitled

ú r(ú) t n [;t] [;t=n]; (t=n; 2t=n]; (2t=n; 3t=n];:::; ((nä 1)t=n;t] n t 1 (nä1)t=n e Är(t)=n (nä 2)t=n e Är(t)=n e Är((nÄ1)t=n)=n t e Är(t)=n e Är((n

<4D F736F F D2092B28DB882C982C282A282C42E646F63>

( ) 1,771,139 54, , ,185, , , , ,000, , , , , ,000 1,000, , , ,000

Ÿ ( Ÿ ) Ÿ šœš 100,000 10,000,000 10,000,000 3,250,000 1,000,000 24,350,000 5,000,000 2,500,000 1,200,000 1,000,000 2,960,000 7,000,000 1,500,000 2,200

š ( š ) ,400, , ,893, ,743, ,893, ,893, ,658,475 35,884,803 2,167,037 7,189,060 55,417,575 81,08

š ( š ) (6) 11,310, (3) 34,146, (2) 3,284, (1) 1,583, (1) 6,924, (1) 1,549, (3) 15,2

73

.....Z...^.[ \..

<4D F736F F D F8DE98BCA8CA797A78FAC8E9988E397C3835A E815B82CC8A E646F63>

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

Ÿ Ÿ ( ) Ÿ , , , , , , ,000 39,120 31,050 30,000 1,050 52,649, ,932,131 16,182,115 94,75


untitled

エジプト、アブ・シール南丘陵頂部・石造建造物のロータス柱の建造方法

Ÿ ( ) Ÿ ,195,027 9,195,027 9,195, ,000 25, ,000 30,000 9,000,000 9,000, ,789, ,000 2,039,145 3,850,511 2,405,371

ID 3) 9 4) 5) ID 2 ID 2 ID 2 Bluetooth ID 2 SRCid1 DSTid2 2 id1 id2 ID SRC DST SRC 2 2 ID 2 2 QR 6) 8) 6) QR QR QR QR

Ÿ ( ) ,166,466 18,586,390 85,580,076 88,457,360 (31) 1,750,000 83,830,000 5,000,000 78,830, ,388,808 24,568, ,480 6,507,1


Proposal of Driving Torque Control Method for Electric Vehicle with In-Wheel Motors Masataka Yoshimura (Yokohama National University) Hiroshi Fujimoto

“LŁñ‡¤‡½‡Ã1„”“ƒ‡¨‡Ł‡è

Microsoft Word - 99

<4D F736F F D F8DE98BCA8CA797A78FAC8E9988E397C3835A E815B82CC8A E646F63>

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

<4D F736F F D2088CF88F589EF8E9197BF81698CA28E9490E78DCE816A2D312E646F63>

<4D F736F F D EC08E7B8FF38BB BD90AC E A837A815B B83578C668DDA97702E646F63>


…_…C…L…fi…J…o†[fiü“ePDF/−mflF™ƒ

WINS クラブ ニュース

untitled

untitled

p1_5.pmd

untitled

6kg 1.1m 1.m.1m.1 l λ ϵ λ l + λ l l l dl dl + dλ ϵ dλ dl dl + dλ dl dl 3 1. JIS 1 6kg 1% 66kg 1 13 σ a1 σ m σ a1 σ m σ m σ a1 f f σ a1 σ a1 σ m f 4



320_…X…e†Q“õ‹øfiÁ’F

<4D F736F F D2088CF88F589EF8E9197BF816991E596EC927C A2E646F63>

<4D F736F F D F8DE98BCA8CA797A78FAC8E9988E397C3835A E815B82CC8A E646F63>

( ) g 900,000 2,000,000 5,000,000 2,200,000 1,000,000 1,500, ,000 2,500,000 1,000, , , , , , ,000 2,000,000

Microsoft Word - ’V‘é−gŁš.doc


„¤‰ƒ‰IŠv‚æ‡S−ª†{“Å‘IB5-97

< F31332D8B638E FDA8DD E F1292E6A>

V(x) m e V 0 cos x π x π V(x) = x < π, x > π V 0 (i) x = 0 (V(x) V 0 (1 x 2 /2)) n n d 2 f dξ 2ξ d f 2 dξ + 2n f = 0 H n (ξ) (ii) H

Microsoft Word - −C−…−gŁš.doc

<4D F736F F D BB388E78CA48B B E6328AAA D655F92B290AE82B382E782C E646F63>

2012専門分科会_new_4.pptx

ii

<4D F736F F D C8E86817A91E682528E9F97A28FAF92AC926A8F978BA493AF8E5189E68AEE967B8C7689E6>

Microsoft Word - GrCadSymp1999.doc

2 G(k) e ikx = (ik) n x n n! n=0 (k ) ( ) X n = ( i) n n k n G(k) k=0 F (k) ln G(k) = ln e ikx n κ n F (k) = F (k) (ik) n n= n! κ n κ n = ( i) n n k n

50 2 I SI MKSA r q r q F F = 1 qq 4πε 0 r r 2 r r r r (2.2 ε 0 = 1 c 2 µ 0 c = m/s q 2.1 r q' F r = 0 µ 0 = 4π 10 7 N/A 2 k = 1/(4πε 0 qq

ÍÂ~ÊÂ ÊÊ ÇÍ ÌÉÊÊÌÊÇÍÂÈÍ Ê ÊÌÊÊÍÉÉÉÆÉÉÍÆÂsÊÂ ÌÉÊ~ÊsÊÆÇ ÉÉÊsÆÍÆÊÉ~ÇÈÉÇÉÉÊsÉÆÆjÇÆÇÉÉÉÆÉÉÍ ÆÂ ÊÊÍÉÂÇÍÌÉÊsÊÊÇÉÂÊÍÍÉwÊÊÂÌÉ t ÊwÎÔ ÑÊÔÖÏÑ Ö Ñ ÑÒÔÇ ÈÍÍÇÉÊÊÍÂÇ

( š ) š 13,448 1,243,000 1,249,050 1,243,000 1,243,000 1,249,050 1,249, , , ,885

0304_ふじみ野地福_本編_01


š ( š ) ,148,770 3,147,082 1, ,260 1,688 1,688 10,850 10, , ,

& Vol.5 No (Oct. 2015) TV 1,2,a) , Augmented TV TV AR Augmented Reality 3DCG TV Estimation of TV Screen Position and Ro

<4D F736F F D BB388E78CA48B B E6338AAA2B92B290AE2B E646F63>

untitled

情報理論 第5回 情報量とエントロピー

PPP_‚Ü‚Æ‚ß.pdf

š ( š ) 7,930,123,759 7,783,750, ,887, ,887 3,800,369 2,504,646,039 i 200,000,000 1,697,600, ,316.63fl 306,200,

Microsoft Word - p2-11堀川先生_紀要原稿_ final.doc

鉄鋼協会プレゼン

5 1.2, 2, d a V a = M (1.2.1), M, a,,,,, Ω, V a V, V a = V + Ω r. (1.2.2), r i 1, i 2, i 3, i 1, i 2, i 3, A 2, A = 3 A n i n = n=1 da = 3 = n=1 3 n=1

ÊÈÌÊ fêôöôï Ö É É ~ Œ ~ Œ ÈÍÉÆÍ s Ê É Â Ê ÉÉÆÍÇÉ Ê Ê É Ê ÈÍv ÈÍ É ÈÍ Â ÇÍ vèé Ê Ê É ÈÉËÈÆ ÊÌÉ Ê~Æ Ê Ê ÈÍfÆ Ê ÊÉÆÉÊ Ê Ê ÈÍ Ê ÈÉËÈÆ

2

untitled

1


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

力学的性質

T rank A max{rank Q[R Q, J] t-rank T [R T, C \ J] J C} 2 ([1, p.138, Theorem 4.2.5]) A = ( ) Q rank A = min{ρ(j) γ(j) J J C} C, (5) ρ(j) = rank Q[R Q,

車両開発における構造・機構のCAE

fm

untitled

<4D F736F F D F8DE98BCA8CA797A78FAC8E9988E397C3835A E815B82CC8A E646F63>

Ò ÑÔÏÓ ÐÎ ÆÉ z uññòõ w g ÌÊÉÇÍ ˆ ˆ Ð Ö Ò z Ò ÑÔÏÓ Ð ÓÑÐÒÒ ÎÔÖÏÖ ÎÖÐÖÑÕ uôöðöõ Î~ËÍÂÌÉÂ ÑÑÒÕÊ ÉÊÍ ÌÆÇÇ Î Ê ÈÂÊÈÇÊÓÑÐÒÒ ÇÂ z uêèéæíçî ÍÇÊÈÍÂ t Ê Ç ÈÍÂ Â

Transcription:

EOIRS OF SHONAN INSTITUTE OF TECHNOLOGY Vol. 37, No., 2003 * A Geneal Numeical Solution to Ditibuted Shape Optimization Poblem in Stuctual Deign aatohi SHIODA* In thi pape, a numeical hape optimization method of continua i peented fo typical tength, igidity and vibation poblem in tuctual deign. A the tength poblem, the minimization poblem of maximum te and the hape detemination poblem that achieve a given deied te ditibution ae fomulated. The igidity poblem involve the minimization poblem of extenal wok and the hape detemination poblem that achieve a given deied diplacement ditibution. Alo, the vibation poblem involve maximization of eigen fequency with mode tacking. Each poblem i fomulated and enitivity function ae deived uing the Lagangian multiplie method and the mateial deivative method. The taction method, which i a hape optimization method, i employed to find the optimal domain vaiation that educe the objective functional. The popoed numeical analyi method make it poible to deign optimal tuctue efficiently. Example of computed eult ae peented to how the validity and pactical utility of the popoed method. GUI (Gaphical Ue Inteface) * 4 0 CAD ) 7

37 2 2. CAD 2.2 600 Galileo 973 Zienkievictz and Campbell 2) 994 3) 4) 7) 3 8) W R n, n 2, 3 T (X): X W x W Q T D (X) T (X) DV O( D ) () V T (X) Eule V C Q {V C (W: R n ) V 0 in Q} (2) O( D ) lim O( D )/ D 0 D Æ0 J f J W f dx Eule J () (3) 8

J f dx f n dg n W G (4) È ( x ) max Í, (8) x ŒW Î a n n n i V i n G W Eintein ( ),i ( )/ x i ( ) Eule ( ) Lagange 8) J f a KS 9) G J f dg (5) W ( x) KS( ( x)) ÏÔ Ê ˆ Ô ln Ìexp dx Á Ë ÓÔ a Ô (9) Eule J J { (, ini ) } d f f f k n G G k 2 3 4 4. (6) KS (x) 5 200 KS l p 4..2 Fig. W G W V W G W G deign f f i, i=,, n P P i W G (9) Given W, f in W, P on G, o R (0) find W (o V ) () {( xx yy) ( yy zz) ( zz xx) 2 2 2 2 2 2 2 2 6( )} xy yz xz (7) 4.. Fig.. Domain vaiation of continuum. 9

37 that minimize ln W ÏÔ Ê exp ( x ) ˆ Ô Ì Á dx Ë ÓÔ a Ô (2) A a W Ê expá Ë a ˆ dx (9) ubject to 0 0 (3) a(v,w) l(w) fo all w U (4) a(v,w) l(w) () (2) v v i w w i U, 0 R W a( vw, ) eijklvk, lwi, jdw l( w) f w dw Pw dg i i i i W G (5) (6) e e ijkl L w(x), L ÏÔ Ê ( x ) ˆ Ô L( W, vw,, L) ln Ìexp dx Á W Ë ÓÔ a Ô l( w) a( v, w) L( ) 0 (7) n L v,w L a (v,w ) l (w ) fo all w U (20) Ê ˆ a( v, w) exp ijdx A Á Ë A vdx fo all v ŒU k W W Ê expá Ë (2) L ( 0 ) 0, 0 0, L 0. (22)(23)(24) (20) v (2) w (22) (24) Kuhn Tucke L v, w L L l G (V ) (25) V l G (V ) a a ˆ ij ij ij k (n l V l 0 o G ) (f 0) L V L l( w ) a( v, w ) A W Ê ˆ exp Á ijdx a (, ) Ë v w a ÏÔ Ê ˆ Ô Ì eijklvk, lwi, j fiwi exp A Á L G ÓÔ Ë a Ô nvd l l G L ( 0 ), V ŒCQ (8) ij G i i G l ( V ) GVdG (26) Ï ÔÊ Ê ˆ ˆ Ô G ÌÁ eijklvk, lwi, j fiwi exp n A Á L a ÓÔ Ë Ë Ô on G deign G \G fix (27) G v (20) L (22) (24) V 20

4.2 (x) (x) (28) W C W C (2) ij a( v, w ) 2 { ( x) ( x) } vdx k v Wc { } W ( x) ( x) dx C fo all v U (29) (26),(27) V l G (V ) G, G 2 (30) (32) W G 2 l ( V ) G VdG G VdG G i i i i c 2 G { ( x) ( x) } n G 2 ( e v w f w ) n ijkl k, l i, j i i L (30) (3) (32) 5 5. v (6) l(v) l(v) V l G (V ) 2 ij k G (33) (34) (35) a( v, w) l( v ) fo all v ŒU G ( 2 fv e v v L) n on G G \ G i i ijkl k, l i, j deign fix (33) (34) (35) (33) (20) v w G 2 5.2 4.2 G D G i i G l ( V ) GVdG {( v ( x) v ( x)}{( v ( x) v ( x)} dg i i i i (36) G D ṽ i (37) a( v, w ) 2 { v ( x) v ( x)} v dg i i i WC (37) V l G (V ) G G i i G l ( V ) GVdG G ( e v w f w ) n ijkl k, l i, j i i L (38) (39) 2

37 6 v () l () l () (4) a(v (),w) l () b(v (),w), fo all w U (40) b(v (),w) ( ) i i W b( v, w ) v w dx (4) w Lagange L v () w L a(v (),w ) l () b(v (),w ), fo all w U (42) a(v (),w) l () b(v (),w), fo all v () U (43) b(v (),w) (44) L ( 0 ) 0, 0 0, L 0. (45)(46)(47) (42) (44) v () (43) w v () w (45) (47) (22) (24) Kuhn Tucke v () w L Lagange L V G G i i G l ( V ) GVdG G ( e v v l v L) n ijkl ( ) k, l ( ) i, j ( ) ( ) i (48) (49) 7 2 4 6 7. V (50) a(v,w) l(w) fo all w C Q (50) (50) G V (50) V L L DL l G (DV ) O( D ) (5) (50) (5) e ijkl a(v,w) $ a 0 : a( x, x) a x fo all x ŒU (52) D DL a(dv,dv ) 0 (53) (50) V L 0) 7.2 FE 2 22

FE FE (50) D G G deign Fig. 2 7.3 v v () w (2) (25) (54) [B] T [D] w 2{ (x) (x)} ij T [ B] [ D][ B] dx{ x} W W Ï T Ô Ê ˆ [ B] [ D] Ì exp dx A Á ÓÔ Ë a Ô ij Ô Fig. 2. Schematic flow chat of hape optimization ytem. 2{v i (x)-ṽ i (x)} w 7.4 PID PID ) L G L 8 23

37 Fig. 3. Solid am poblem. Fig. 5. Iteation hitoie. Fig. 6. Notched plate poblem. Fig. 4. Initial & optimal hape and te ditibution (Pa). Fig. 7. Taget hape. 4 8 8. 8.. Fig. 3 Fig. 3(a) P Fig. 3(b) Fig. 4 Fig. 5 30% 8..2 Fig. 6 W C (x) Fig. 7 Fig. 6(a) 24

Fig. 8. Shape identification eult. Fig. 0. Iteation hitoie. Fig.. Toion am poblem. Fig. 9. Compaion of te ditibution along A B. Fig. 6(b) Fig. 8 Fig. 9 A B Fig. 0 Fig. 7 8.2 8.2. Fig. Fig. (a) Fig. (b) Fig. 2 Fig. 3 2 Fig. 4 Fig. 5 8.2.2 Fig. 6 G D y ṽ 2 Fig. 6(a) Fig. 6(b) Fig. 7 Fig. 8 25

37 Fig. 4. Iteation hitoie. Fig. 2. Initial hape. Fig. 5. Optimal hape fo min-max te poblem. Fig. 3. Optimal hape. Fig. 9 26

Fig. 6. Table uppot poblem. Fig. 9. Iteation hitoie. Fig. 7. Initial hape and defomation mode. 8.3 Fig. 20(a) (Fig. 20(b)) PC (ulti Point Containt) Fig. 2(a) (b) 20% 9 Fig. 8. Shape identification eult. 27

37 Fig. 20. Engine mount backet poblem. Fig. 2. Initial & optimal hape and vibation mode. AI ) 60-574, A(994), 479. 2) Zienkiewicz, O. C. and Campbell, J., in Optimal Stuctual Deign, Gallaghe, R. H. and Zienkiewicz, O. C. (Ed.), John Wiley, New Yok (973). 3) Cea, J., in Optimization of Ditibuted Paamete Stuctue, Haug, E. J. and Cea, J., Vol. 2, Sijthoff & Noodhoff, Alphen aan den Rijn (98). 4), 62-604, A (996), 283 2837 5), 63-607, A (997), 60 67. 6), 6-587, C (995), 49 54.. 7) Ihaa, H., Azegami, H. and Shimoda,., in Compute Aided Optimum Deign of Stuctue VI, WIT Pe, 28

Southampton, (999), 87 95. 8) Sokolowki, J. and Zoleio, J. P., in Intoduction to Shape Optimization, Shape Senitivity Analyi, Spinge-Velag, New Yok, (99). 9) Keielmeie, G. and Steinhaue, R., in Sytematic Contol Deign by Optimizing a Vecto Pefomance Index, Intenational Fedeation of Active Contol Sympoium on Compute-Aided Deign of Contol Sytem, Zuich, Switzeland, Augut 29 3, (979), 3 7. 0) Vol., 996 5 ). ) Shimoda,., Azegami, H., Sakuai, T., in Poceeding of 9 th. Intenational Confeence on Vehicle Stuctual echanic and CAE, Toy, Apil 4 6, SAE (995), 223. 29