๐Ÿฅ‡๊ณต๋Œ€์ด๊ฑฐ์ €๊ฒƒ(๋ง‰ํ•™๊ธฐ)/๐ŸŒCAE

๐Ÿฅ‡๊ณต๋Œ€์ด๊ฑฐ์ €๊ฒƒ(๋ง‰ํ•™๊ธฐ)/๐ŸŒCAE

    ๋ฏธ๋ถ„๋ฐฉ์ •์‹ ์ˆ˜์น˜ํ•ด๋ฒ•-์˜ค์ผ๋Ÿฌ,midpoint, heun method, RK4

    kookyungmin.github.io/language/2017/09/16/matlab13/ [matlab] 13.๋ฏธ๋ถ„๋ฐฉ์ •์‹์˜์ˆ˜์น˜ํ•ด๋ฒ•(1) ๋ฏธ๋ถ„๋ฐฉ์ •์‹์˜ ์ˆ˜์น˜ํ•ด๋ฒ• - ์ดˆ๊ธฐ ๊ฐ’ ๋ฌธ์ œ kookyungmin.github.io ์ข‹์€ ์ž๋ฃŒ๊ฐ€ ์žˆ์–ด์„œ ๊ณต์œ ํ•˜๊ณ ์ž ํ•œ๋‹ค.

    ๋‰ดํ„ด์˜ ์—ญ์ œ๊ณฑ์˜ ๋ฒ•์น™- ๋งŒ์œ ์ธ๋ ฅ

    ๋‰ดํ„ด์˜ ์—ญ์ œ๊ณฑ์˜ ๋ฒ•์น™- ๋งŒ์œ ์ธ๋ ฅ

    ์—ญ์ œ๊ณฑ์˜ ๋ฒ•์น™ : ์–ด๋–ค ๋ฌผ๋ฆฌ๋Ÿ‰์ด ๊ฑฐ๋ฆฌ์˜ ์ œ๊ณฑ์— ๋ฐ˜๋น„๋ก€ ๋‹ฌ ํ‘œ๋ฉด์—์„œ์˜ ์ค‘๋ ฅ์€ ์ง€๊ตฌ์˜ 1/6๋ฐ–์— ๋˜์ง€ ์•Š๋Š”๋‹ค. ์ด ๋•Œ๋ฌธ์— ๋‹ฌ์— ๊ฐ„ ์šฐ์ฃผ์ธ๋“ค์˜ ๋™์ž‘์ด ์ž์—ฐ์Šค๋Ÿฝ์ง€ ๋ชปํ•˜๋‹ค. ๋งŒ์œ ์ธ๋ ฅ์˜ ๊ฐ€์žฅ ํฐ ํŠน์ง•์€ ์ค‘๋ ฅ์ด ๋‘ ๋ฌผ์ฒด ์‚ฌ์ด ๊ฑฐ๋ฆฌ์˜ ์ œ๊ณฑ์— ๋ฐ˜๋น„๋ก€ํ•œ๋‹ค๋Š” ์ ์ด๋‹ค. ๋ฌผ๋ฆฌํ•™์—์„œ๋Š” ์ด์ฒ˜๋Ÿผ ์–ด๋–ค ๋ฌผ๋ฆฌ๋Ÿ‰์ด ๊ฑฐ๋ฆฌ์˜ ์ œ๊ณฑ์ด ๋ฐ˜๋น„๋ก€ํ•˜๋Š” ๊ฒฝ์šฐ๊ฐ€ ๊ฐ„ํ˜น ์žˆ๋Š”๋ฐ, ์ด๋ฅผ ์—ญ์ œ๊ณฑ์˜ ๋ฒ•์น™(inverse square law)์ด๋ผ๊ณ  ๋ถ€๋ฅธ๋‹ค. ๋‰ดํ„ด์€ [ํ”„๋ฆฐํ‚คํ”ผ์•„] 1๊ถŒ์—์„œ ํ–‰์„ฑ์ด ์—ญ์ œ๊ณฑ์˜ ํž˜์„ ๋ฐ›๋Š”๋‹ค๋Š” ๊ฐ€์ •ํ•˜์— ์ผ€ํ”Œ๋Ÿฌ์˜ ์„ธ ๊ฐ€์ง€ ๋ฒ•์น™์„ ์œ ๋„ํ–ˆ๋‹ค. ๋˜ํ•œ ํž˜์ด ์ •ํ™•ํ•˜๊ฒŒ ๊ฑฐ๋ฆฌ์˜ ์ œ๊ณฑ์— ๋ฐ˜๋น„๋ก€ํ•˜๋ฉด ๊ทธ ๊ถค๋„๋Š” ๋‹ซํžŒ ๊ถค๋„์ž„์„ ์‰ฝ๊ฒŒ ๋ณด์ผ ์ˆ˜ ์žˆ๋‹ค. ๊ทธ๋ ‡๋‹ค๋ฉด ์ค‘๋ ฅ์€ ์™œ ๊ฑฐ๋ฆฌ์˜ ์ œ๊ณฑ์— ๋ฐ˜๋น„๋ก€ํ• ๊นŒ? ์ค‘๋ ฅ์„ ์ง๊ด€์ ์œผ๋กœ ์ดํ•ดํ•˜๊ธฐ ์œ„ํ•ด์„œ ์šฐ์„  ํผ์งํ•˜๊ณ  ๋‘ฅ๊ทผ ์‚ฌ๊ณผ๋ฅผ ํ•˜๋‚˜ ์ค€๋น„..

    Plot ์ฐธ๊ณ 

    Plot ์ฐธ๊ณ 

    1. Specifiers 2.์„ ๋‘๊ป˜ plot(x,y1,x,y2,'LineWidth',2.0) 3.legend ํŠน์ • ๊ทธ๋ž˜ํ”„๋ฅผ ์ง€์ •ํ•˜์—ฌ legend๋ฅผ ๋ถ™์ด๋Š” ๋ฐฉ์‹ plot3(x,y,z,'line specifiers' , 'property Name' , 'property value') ์ด ๋ช…๋ น์–ด๋ฅผ ์‚ฌ์šฉํ•˜๋Š” ์กฐ๊ฑด ๋˜ํ•œ ์ฐจ์ด์ ์€ ์—†์Šต๋‹ˆ๋‹ค. ๋ช…๋ น์–ด๋ฅผ ์‚ฌ์šฉํ•˜๋ฉด ์š”๋Ÿฐ ์‹์œผ๋กœ ๊ทธ๋ž˜ํ”„๊ฐ€ ์ถœ๋ ฅ์„ํ•ฉ๋‹ˆ๋‹ค. ์‰ฝ์ฃ ? ๋งŒ์ผ '์‚ฌ์šฉ์ž ๋””์ž์ธ' ์ž…๋ ฅ์— ๋Œ€ํ•œ ๋‚ด์šฉ์ด ๊ธฐ์–ต์ด ์•ˆ๋‚˜์‹œ๋ฉด '2์ฐจ์› ๊ทธ๋ž˜ํ”„' ํฌ์ŠคํŒ…์—์„œ ์„ค๋ช…์„ ๋ณด๊ณ  ์˜ค์‹œ๋ฉด ๋ฉ๋‹ˆ๋‹ค. ๊ทธ๋ž˜ํ”„๋ฅผ ์ถœ๋ ฅํ–ˆ์ง€๋งŒ ๋’ท๋ฐฐ๊ฒฝ์— ์„ ์ด ๊ตฌ๋ถ„์„ ์ด ์—†์œผ์‹œ๋‹ค๋ฉด 'grid on' ๋ช…๋ น์–ด๋ฅผ ์ž…๋ ฅํ•˜์‹œ๋ฉด ๋ฉ๋‹ˆ๋‹ค. (๊ตฌ๋ถ„์„ ์„ ์ง€์šฐ๊ณ  ์‹ถ๋‹ค๋ฉด 'grid off') ํ–‰๋ ฌ์˜ ๋ฐ์ดํ„ฐ๊ฐ€ ์•„๋‹Œ ๋ฒ”์œ„ ํ˜•์‹์œผ๋กœ ์ฃผ์–ด์ง„ ๋ฐ์ดํ„ฐ๊ฐ€ ์žˆ๋‹ค..

    ๋ณด๊ฐ„๋ฒ•

    ๋ณด๊ฐ„๋ฒ•

    1. ๋ณด๊ฐ„๋ฒ• ๋ณด๊ฐ„๋ฒ•์ด๋ž€? ์‹ค๋ณ€์ˆ˜ x์˜ ํ•จ์ˆ˜ f(x)์˜ ๋ชจ์–‘์€ ๋ฏธ์ง€์ด๋‚˜, ์–ด๋–ค ๊ฐ„๊ฒฉ์„ ๊ฐ€์ง€๋Š” 2๊ฐœ ์ด์ƒ์ธ ๋ณ€์ˆ˜์˜ ๊ฐ„ x์— ๋Œ€ํ•œ ํ•จ์ˆ˜๊ฐ’์ด ์•Œ๋ ค์ ธ ์žˆ์„ ๊ฒฝ์šฐ, ๊ทธ ์‚ฌ์ด์˜ ์ž„์˜์˜ x์— ๋Œ€ํ•œ ํ•จ์ˆ˜ ๊ฐ’์„ ์ถ”์ •ํ•˜๋Š” ๊ฒƒ -๋„ค์ด๋ฒ„ ๋ฐฑ๊ณผ์ถœ์ฒ˜ ์˜ˆ๋ฅผ ๋“ค๋ฉด ์‹œ๊ฐ„ 30์ดˆ 1๋ถ„ 1๋ถ„30์ดˆ 2๋ถ„30์ดˆ 3๋ถ„ 4๋ถ„ ์ „๋ ฅ 110 270 350 260 200 150 ์ด๋Ÿฐ ์ž๋ฃŒ๊ฐ€ ์žˆ์„ ๋•Œ ์‹œ๊ฐ„์— ๋”ฐ๋ฅธ ์ „๋ จ๋Ÿ‰์— ๋Œ€ํ•œ ๊ด€๊ณ„๋ฅผ ๋ชจ๋ฅด๊ธฐ ๋•Œ๋ฌธ์— 2๋ถ„์ผ ๋•Œ ์ „๋ ฅ์ด๋‚˜ 3๋ถ„ 30์ดˆ์ผ๋•Œ ์ „๋ ฅ์„ ์•Œ ์ˆ˜๊ฐ€ ์—†๋‹ค. ์ด๋•Œ ์‚ฌ์šฉํ•˜๋Š” ๊ฒƒ์ด ๋ณด๊ฐ„๋ฒ•์ด๋‹ค. ๋งคํŠธ๋žฉ์—์„œ๋Š” ์„ ํ˜• ๋ณด๊ฐ„๋ฒ•์œผ๋กœ interp ํ•จ์ˆ˜๋ฅผ ์‚ฌ์šฉํ•œ๋‹ค. interp1 : ๋ณ€์ˆ˜๊ฐ€ ํ•œ ๊ฐœ์ธ ํ•จ์ˆ˜๋ฅผ ๋‚ด์‚ฝํ• ๋•Œ interp2: ๋ณ€์ˆ˜๊ฐ€ ๋‘ ๊ฐœ์ธ ํ•จ์ˆ˜๋ฅผ ๋‚ด์‚ฝํ• ๋•Œ interp1(x,y,x_int)๊ฐ€ ๊ธฐ๋ณธํ˜•์ด๊ณ  x๋Š” ๋…๋ฆฝ๋ณ€์ˆ˜ ..

    ODE + MATLAB

    ODE + MATLAB

    1 ์˜ค์ผ๋Ÿฌ : first-order 2 Heun's method 3 Midpoint method 4 Runge-Kutta methods(Seconds, Third, Fourth, Higher) 1. F = ma ๋‰ดํ„ด์˜ ์ œ 2๋ฒ•์น™ (ํž˜์˜ ํ•ฉ๋ ฅ 0์ด ์•„๋‹๋•Œ) 2. ํ‘ธ๋ฆฌ์— ์—ด ๋ฐฉ์ •์‹ : ์—ด์ด ์šด๋™ ํŠน์„ฑ, ์—ด๋ฅ˜ ๋ฐ€๋„ q ๊ตฌํ• ๋•Œ, T๋Š” ์˜จ๋„์˜ ๋ฏธ๋ถ„: ๋†’์€๊ณณ์—์„œ ๋‚ฎ์€๊ณณ ๊ธฐ์šธ๊ธฐ -๊ฐ€ ๋œ๋‹ค. ๋‹จ์œ„ K. ์—ฌ๊ธฐ์„œ qยฏ๋Š” ์—ด๋ฅ˜ ๋ฐ€๋„(heat flow density)[๋‹จ์œ„: W/ใŽก], ฮบ๋Š” ์—ด ์ „๋„๋„(heat conductivity: ๋ฌผ์งˆ์˜ ๊ณ ์œ ํ•œ ํŠน์„ฑ)[๋‹จ์œ„: W/K/m], T๋Š” ์˜จ๋„(temperature)[๋‹จ์œ„: K]์ด๋‹ค. ์‹ (1)์„ ์ž˜ ์ดํ•ดํ•˜๋ ค๋ฉด ๊ตฌ๋ฐฐ ์—ฐ์‚ฐ์ž(gradient operator) โˆ‡ยฏ๋ฅผ ๋ด์•ผ ํ•œ๋‹ค. ๊ตฌ๋ฐฐ๋Š”..

    [Matlab]getframe

    take a snapshot of scenario www.youtube.com/watch?v=3I1_5M7Okqo&ab_channel=ChristopherLum

    ๋ช…๋ น๋ฌธ ์ตํžˆ๊ธฐ

    CAE ๋ฅผ ํ•˜๋‹ค๋ณด๋ฉด ๋‹ค์–‘ํ•œ ๋ช…๋ น๋ฌธ์„ ์‚ฌ์šฉํ•˜๊ณ , ์–ธ์–ด์˜ ๊ธฐ๋ณธ ํ‹€์ธ ๋ฌธ๋ฒ•์„ ๋ฐฐ์šฐ๊ฒŒ ๋˜๋Š”๋ฐ ์ •๋ฆฌ๊ฐ€ ์ž˜ ๋˜์–ด ์žˆ๋Š” ๋ธ”๋กœ๊ทธ๋“ค์ด ์ฐธ ๋งŽ์€ ๊ฒƒ ๊ฐ™๋‹ค. m.blog.naver.com/cj3024/221111764415

    1์ฃผ์ฐจ-์ž์œ ๋‚™ํ•˜ ์šด๋™,

    1์ฃผ์ฐจ-์ž์œ ๋‚™ํ•˜ ์šด๋™,

    CAE์˜ ์ˆ˜์—…์ด ์‹œ์ž‘๋˜์—ˆ๋‹ค. ๋ฌธ์ œ๋ฅผ ํ’€๊ธฐ์ „ ๊ธฐ๋ณธ ์ ์ธ ๊ฐœ๋…์„ ์ •๋ฆฌํ•˜๊ณ  ๊ฐˆ ํ•„์š”๊ฐ€ ์žˆ์–ด์„œ, ์˜›๋‚  ๊ธฐ์–ต๋„ ๋˜์‚ด๋ฆด๊ฒธ ์ •๋ฆฌํ•ด๋ณด์•˜๋‹ค. 1. ์งˆ๋Ÿ‰๊ณผ ๋ฌด๊ฒŒ, SI unit, MKS ๋‹จ์œ„๊ณ„ ์งˆ๋Ÿ‰ (mass) ๋ฌผ์ฒด์˜ ๊ณ ์œ ํ•œ ์–‘; kg ํ˜น์€ g ์–‘ํŒ” ์ €์šธ์ด๋‚˜ ์œ—์ ‘์‹œ ์ €์šธ์„ ์ด์šฉํ•˜์—ฌ ์งˆ๋Ÿ‰์„ ์ธก์ •ํ•œ๋‹ค. ๋ฌผ์ฒด์˜ ๊ณ ์œ ํ•œ ์–‘์ด๋ฏ€๋กœ ์žฅ์†Œ์— ๋”ฐ๋ผ ๋ณ€ํ•˜์ง€ ์•Š๋Š”๋‹ค. ๋ฌด๊ฒŒ (weight) ์ค‘๋ ฅ์œผ๋กœ ์ธํ•ด ๋ฌผ์ฒด๋ฅผ ๋Œ์–ด๋‹น๊ธฐ๋Š” ํž˜(force): N. ํ˜น์€ kgf, gf ์šฉ์ˆ˜์ฒ  ์ €์šธ์ด๋‚˜ ์•‰์€๋ฑ…์ด ์ €์šธ์„ ์ด์šฉํ•˜์—ฌ ์ค‘๋ ฅ์˜ ํž˜์˜ ํฌ๊ธฐ๋ฅผ ์•Œ์•„๋‚ธ๋‹ค. ์ค‘๋ ฅ์€ ๋ฌผ์ฒด๋ฅผ ๋Œ์–ด๋‹น๊ธฐ๋Š” ํž˜์ด๋ฏ€๋กœ ์ค‘๋ ฅ์ด ๋‹ค๋ฅธ ์žฅ์†Œ์—๋Š” ๋ฌด๊ฒŒ ๋˜ํ•œ ๋ณ€ํ•œ๋‹ค. ์˜ˆ) ์ง€๊ตฌ์—์„œ์˜ ๋ฌด๊ฒŒ = 60kg x 9.8m/s2 = 588N = 60kgf ๋‹ฌ์—์„œ์˜ ๋ฌด๊ฒŒ = 60kg x 1.622m/s2 = 97..