Tải bản đầy đủ (.pdf) (43 trang)

chương 2 đại số ma trận mô hình máy tính trong thiết kế máy bay

Bạn đang xem bản rút gọn của tài liệu. Xem và tải ngay bản đầy đủ của tài liệu tại đây (4.05 MB, 43 trang )

<span class="text_page_counter">Trang 2</span><div class="page_container" data-page="2">

<b>Chương2.</b>



<b>ĐẠISỐMATRẬN</b>

<b>MơhìnhMáytínhtrongThiết ếMáybayk</b>

Đểthiếtkế ếth hệmáybaythương mại vàquân sựkếtiếp,các kỹsưt i    ạ Boing’sPhantom Work sử dụng mơhình 3D và tính tốn động lực chấ  ỏng (CFD). ọ nghiên cứu dịngt l H khơngkhíquanhmộ máybayảođểtrả ờt l icho nh ngcâuh ivềthiếtkếquantr ngtrước ữ ỏ ọ  khit oạ nhữngmơhìnhvậ lý.Điềunàyđãt giảmđáng ểsốchutrìnhthiế  ếvàgiáthành–k t k vàđạisốtuyế tínhđóngn vai trị quan tr ng trong tiếntrình này. ọ    

Máybayảolúcđầunhưmộtmơ hình tốn“khung dây” mà chỉt nt i     ồ ạ trong bộnhớmáy  tínhvàthiếtbịđầucuố hiể thịđồ ọ (Mộ mơhìnhBoingi n h a. t 777đãchỉranhưvậy.)Mơ hìnhtốnhọcnàysắ  ế vàtácđộng p x p tới mỗibướcthiết k ếvàchế ạ máybaycảtrong ẫ t o l n ngồi.PhântíchCFDquantâm ớ t i bềmặ tbênngồi.

Mặcdùbềngồilớ vỏcủ máybaykhip a hồnthànhtrơngcóvẻtrơ tru,tuynhiêncấun trúcbề mặtlàkháphứctạp.Bên ạnhcánhvàthânmáyc bay, một máy bay cịn cóvỏ, bộ      phậngiữthăngbằng,bộ ậnđiph ều ướng,…Cáchh màkhơngkhí ịchchuyể xungd n quanh cácbộphậnnàyquyếtđịnhcách mà máy bay di chuyển trên bầ     utr i. Các phương trình mơờ     tảsự ịchchuyể này làd n   ấ phức ạp,và ọphả tínhtốnđếnđầuvào vàr t t h i  đầuracủ độnga cơ,vành ngữ vệtcịnl iạ docánhcủamáybay.Đểnghiêncứu luồng khí, các kỹsưcầ    nmột mô ảt cả itiến pcấ caocủ bềmặ máybay.a t

M tộ máytính ạ  mộ mơhìnhbềmặ bằngcáchthêm ướ bachiềucủ các ộ trongt o ra t t l i a h p mơhình ngun th y.ủ Cáchộptrong ướinàyhoặchồntồn ằmtronghoặcl n ngồimáy bay,hoặcchúnggiaovớibềmặtcủamáybay.Máytính ự chọ cáckhố giao ớibềm tl a n i v ặ và chianhỏchúng,chỉgiữl iạ những ộ  ỏh nh pnh ơ màvẫ cịngiaovớ bềmặt.Qtrìnhn i chia nhỏnàyđượcl p l iặ  ạ chođến khilướ i là cựcm n. M tị  ộ l iđiểnướ hình có thểchứa trên 400,000    hộp.

Qtrìnhtìmluồngkhíquanhmáybay bao  ồmviệcgiả  ặ đi l p lg i l p  ặ  ạicác ệphươngh

<b>trìnhtuyếntính x</b><i>A</i> =<b>b</b>mà có thểliên quanđến hai<b>triệuphươ trìnhvàbiế số Véctơb</b>ng n . thayđổimỗ  ần,và ự vào ữliệucủ  ướ vànghiệmcủ cácphươngtrìnhtrướcđó.Sửi l d a d a l i a dụngmộ máytínhnhanhnhấ ngồithịtrường,mộtnhómcủ FantomWorkscóth dànht t a ể  từmộ vàigiờđế mộ vàingàyđểthiếtt n t lập và giải quyết   mộ ấtv nđềvềlu ng khí.Sau khiồ   nhómnàyphântíchcácgiảipháp,họcó ểcóth mộtsốđiều ỉnh nhch  ỏtớibềmặtcủamáy baytồnbộqtrìnhmộ  ầ t l n nữa.Cóthểcầnđến hàngngàn lần chạ  y CFD.

Chương này trìnhbàyhaikháiniệmquantrọ hng ỗtrợviệctìmnghiệm ủ  ữ hệc anh ng phươngtrìnhđồsộ ưvậy:nh

 Khốimatr n:M thệCFDđi nậ ộ ể hìnhcómatrậ  ệsố“thư vớ  ầu ế cácphầ  ửn h a” i h h t n t khơng.Nhómcácbiế mộ cáchchínhxác ẫ  ớ mộ matrậnkhố vớ nhiềukhốin t d n t i t i i không.Mục giớ thiệunhữngmatrậ nhưvậ vàmô ảcácứng ụngcủ chúng.2.4 i n y t d a  Nhântửhóamatrận:Th mchíkhiđượcviếtậ vớ nhữngmatrậnkhối, ệcácphươngi h

trìnhnàylàrấtphứctạ Đểđơnp. giả  ơn ữ cáctínhtốn,phầ m mn h n a n ề CFDcủ a Boingsửdụngcáig iọ lànhân ửhóaLUcủ ma ật a trnhệs .Mục2.5 ảoluậnvềLUố th và vi ệcnhântửhóacácmatrậ hữun íchkhác.Chitiếthơnvềnhântửhóaxuấ hiệ ởt n mộtsố ụcsauđó.m

</div><span class="text_page_counter">Trang 3</span><div class="page_container" data-page="3">

3 Đểphântích nghiệm của một hệluồng    khí, các kỹsưmuốnhình dung ra luồng khí      trên bề mặt máy bay. ọsử dụng đồH họa máytính và  đại sốtuyế  tính để cung cấpn phươngtiệnchođồhọa.Mẫukhungdâycủabềmặtmáybayđược ưutrữnhưlàl dữ ệ li u trongnhiềumatrận.Mộtkhimộ hìnhảnhđượcxuấtramànt hìnhmáytính,cáckỹsưcó thểthayđổiquymơcủanó, phóng to ơcjthunhh ỏcácvùng ỏ,vàxoaynhữnghìnhnh ảnhđểxemcácphần mà có thểbịchekhuất. Mỗi    thao tác nàyđược  hồn thành bởiviệ c nhâncácmatrậnthíchhợp.Mục2.7giả ithíchcácýtưởngcơbả n.

Khả năngphân tích và giải các phương trình của chúng tơi sẽ đượ  nângc cao khi chúngtôithựchiệncácbiế đổiđạin sốcácmatrận. ơ  ữa,cácđịnhnghĩavàđịnhlýH n n trongchương này cungcấpmộtsốcông ụcơbảnđểc giảiquyếtnhi u ngdụngcủađạiề ứ sốtuyếntínhbaogồmhai hay nhiềumatrận.V iớ nh ngữ matrậnvng,Định lý Ma trận nghịchđảotrong mục 2.3 kết nối hầuhết các khái niệm ãđượctrìnhbàytrướcđó.Mụ         đ c 2.4và2.5kiểmtracácmatrận ốikh vàviệcnhân ửhóamat trận,màxuấ hiệ tronghầut n h tế cáccáchdunghiệ đạicủađạin sốtuyế tính.Mục2.6n và2.7mơ ảt haiứngdụ củang đạisốmatrận, cho tới kinhtếhọc và    đồhọamáy tính.

Nếu <i>A</i>làmộtmatrậncỡ<i>m</i>´<i>n</i>,ngh làĩa mộtmatrận<i>vớ mhàngvàncột,khiđóphần tửở</i>i

<i>hàng ivàc t cộ  j ủaAđượckýhiệubởia<small>ij</small></i>vàđượckýhiệu là phần tử( , )   <i>i jcủaA.</i>Chẳng hạn, phần tử(3,2)làsố<i>a</i><sub>32</sub>ởhàngba, cộ thai.Mỗi c t của ộ <i>Alàmột dãy msốthực màcó thểđồng</i>    nhấtvớ mộ véctơtrongi t <i><small>m</small></i>.Thơngthường,nhữngcộtnày<b>đượckýhiệubởia a</b><sub>1</sub>, <sub>2</sub>¼,<b>a</b><i><sub>n</sub></i>và

<i>củamatrậnA .Một</i>matrậ đườngchéolàmộtn matrậ vngmàn các phần tửngồiđường  chéođềubằng0.Mộtmatrậ màcácphầ tửđề bằngkhơngđượcn n u g iọ là ma  trận không và  viếtlà0.Cỡcủamatrậ khôngthườngđượcchorõràngtùyvào ừngn t ngữcảnh.

Sốhọcvềvéctơđãđượcmơt tả ừtrướ cósựmởrộngc tựnhiênchocácmatrậ Tanóin. hai ma trận là bằng nhau  ếu chúng có cùng cỡ (nghĩa làn  có cùng sốhàng và sốcộ  các cộtt) tươngứngcủachúnglàbằ nhau, ặcng ho nóimộtcáchkháccácphầ  ử ủachúnglàbằ n t c ng nhau.<i>NếuA Blàcácmatrậncỡ</i>và <i>m</i>´<i>n,thìtổngA B</i>+ làmatrận cỡ <i>m</i>´<i>n</i>màcáccộtc aủ  nólàtổngcủacác ộ  ươ ứngc t t ng của <i>A</i>và<i>B.Từ ệ  ộ cácvéctơcộtđược</i>vi c cng thự hiệntheoc

</div><span class="text_page_counter">Trang 4</span><div class="page_container" data-page="4">

4 t<i>ừngphầ  ử,mỗ phầ  ửcủaA B</i>n t i n t + làtổng tươngứng củacácphần t<i> ử ủaAvàB.Tổ </i>c ng

<i>A</i>+<i>B</i>đượcđịnhnghĩa<i>chỉkhi A   </i> <i>và B có cùng cỡ.</i> 

 Khiđó

<i>nhưngA C</i>+ khơngthểxácđịnhđược vì chúng khơng có cùng cỡ.     

<i>Nếurlàm tộ vôhướngvàAlàm tộ matrận,thìphépnhânvơhướngrAlàmatrậncó các</i> cột là tích của <i>r</i>với các ột tương ứngc <i> của A. Cũ  nhưvéctơ, A</i>ng - <i>được dùng cho( 1)A</i>- , và

</div><span class="text_page_counter">Trang 5</span><div class="page_container" data-page="5">

5 Rõ ràng rằng sốcácc t<i>ộ  của A phả  ươ  thích</i>i t ng <i> vớ  sốcác hàng củ  B đểcho một t h</i>i a  ổ ợp tuyếntínhnhư<i>Ab</i><sub>1</sub>t nồ <i>tại. Cũngnhưvậy,định nghĩacủa ABchỉra rằng AB cócùng sốhàng</i>      

<i>vớiAvàcócùngsốcộtvớ B.</i>i

baonhiêunếuchúng ồ  ại?t n t

<b>GIẢI</b><i>TừA </i>có 5c tộ <i>vàBcó 5hàng, tích của ABlà t ntạivàlàmộtmatrậncỡ3 2</i>    ồ ´ :

Tích<i>BA</i>khơngt n t iồ  ạ vìhaicộ của t <i>B</i>khơng ươngthíchvớ 3t i <i>hàngcủ A.</i>a

 Địnhnghĩacủa<i>AB</i>làquatrọ chotínhtốnlýthuyếtvàứngdụng,tuynhiênquyng tắcsauđâychomộ phươngt thức ệuquảchoviệctínhtốncácphầ tửriêng ẻcủahi n l <i>AB</i>khi

Xemxétcácphần tửliênquansẽchota th ấyvìsaohaiphươngpháptínhtíchmatrậncho cùngmộtkếtquả.

<b>GIẢI</b>Đểtìm<i>phầntửởhàng 1 cột 3 củama trận AB, ta chỉra hàng1 của A và cột 3 của B.</i>                Nhântươngứngcácphầ  ửvàcộngcáckế quả ạ nhưđãthấydướ đn t t l i i ây:

</div><span class="text_page_counter">Trang 6</span><div class="page_container" data-page="6">

ĐỊNHLÝ2<i>GiảsửAlàmột matrậncỡ</i> <i>m</i>´<i>n,BvàClàcácmatrậncócỡsaochocácphép</i>

tốncộngvànhânđượ chỉralàthự hiệnđược.c c

 Tínhtrái‐phảitrongphépnhânlàbắ buộcvìt <i>AB</i>và<i>BA</i>khơnglngiốngnhau.Điềunày khơngcógì llà <i>ạvìchúngtabiế rằngcáccộ củaAB các</i>t t là t hổ ợptuyế tínhcủan cácc cột ủa

<i>A,trongkhiđó,cácc</i>ột ủac <i>BA</i>lạiđượct oạ thànhtừcác ột cc <i> ủaB. ịtrí ủa</i>V c cácthừ sốa

<i>trongtíchABđượcnhấnmạnhbở tagọ Alàđượcnhânphải</i>i i <i>bở B,cịn </i>i <i>B</i>làđượcnhântrái

<i>bởiA.NếuAB</i>=<i>BA</i>thì <i>ta nóiA và B là giao hốn với nhau.</i>      

<i>x</i>Ỵ  thì<i>A x<small>k</small></i> làkết quả của phép nhân  trái <i>x</i>bởi<i>A ặp ạ  </i>l l i l<i>k</i> ần,Nếu<i>k</i>=0thì <small>0</small>

<i>A x</i>là<i>chính x.Dođó</i> <small>0</small>

<i>A</i>đượcthểhiện nhưma trận 

</div><span class="text_page_counter">Trang 7</span><div class="page_container" data-page="7">

1.Từviệcmỗivéctơtrong<i><small>n</small></i> có thểxemnhưmột ma trận cỡ   <i>n</i>´1,cáctínhchấtvềsự chuyểnvịtrongĐịnhlý3cũngđượ ápdụngchocácc véctơ.Cho

<b>Tính x</b>(<i>A</i>) ,<i><small>T</small></i> <b>x</b><i><small>T</small>A<small>T</small></i>,<b>xx</b><i><small>T</small></i><b>vàx x</b><i><small>T</small></i> .Cót nồ t i<b>ạ  x</b><i>A<small>TT</small></i>không?

2.Cho <i>A</i>làmộtmatrận 4 4´    <b>và x là một véctơtrong</b> <small>4</small><b>.Cáchnàonhanhnhấtđểtính x</b><i>A</i><small>2</small> ? Hãytínhsốphépnhân.

</div><span class="text_page_counter">Trang 8</span><div class="page_container" data-page="8">

<b>2.2.NGHỊCHĐẢOCỦAMATRẬN </b>

Một ma trận cỡ<i>n</i>´<i>nA được gọi<b> là khảnghịch ếu ồ  ạ  ma trận C sao choCA I</b></i>n t n t i = và

<i>AC</i> =<i>I</i>với<i>I</i>=<i>I<sub>n</sub></i>làmatrậnđơn<i><b>vịcấpn.Trongtrường ợpnàyClàmộ nghịchđảocủ </b></i>h t a

<i>A</i>. Thực tế,  được xác định<i>C</i>  duynhất<i> bởi A, vì ếu B là mộ  nghịch đảo khác của A thì</i>n t

 <i>Giátrịad bc</i>- đượcgọi là nh đị thức<i>của ma trận A,và taviết det( )</i>    <i>A</i>=<i>ad</i>-<i>bc</i>.Định lý4chỉrarằng <i>A</i>khảnghịchnếuvà ỉnếu det( )ch <i>A</i>¹0.

độngtạicácđiểm 1, 2, 3 nhưhình    <b>vẽ. Gọi ftrong</b>  <small>3</small>làtập tất cảcác lực tạinhữngđiểm này,  

<b>vàgọiytrong</b> <small>3</small>

làtổng sốcácbiếndạng (nghĩalà cácchuyểnđộng)củathanhtạibađiểm   đó.Theođịnh luật Hookecủa vật lý, có thểchỉra rằng      

với <i>Dlàmatrận u n</i> ố (flexibilitymatrix).Nghịchđảocủanóđượcgọilàmatrậncứng(stiffness matrix).Hãymơtảcáchiệ  ượngvậ lýcn t t <i>ủacác ộ củaDvà</i>ct <small>1</small>

<i>D</i><small>-</small>.

</div><span class="text_page_counter">Trang 9</span><div class="page_container" data-page="9">

Bi<b>ểu ễnvéctơe1</b>di =(1, 0, 0)nhưlàmộtđơn vịlực tácđộng  xuố dưới t ing  ạ điểm1trênthanh (với l ựcbằngkhơng ạ t i hai<b>điểmkia).Khiđó e1</b><i>D</i> ,c<i>ộtthứnhất ủaD,làtập</i>c hợ cácđộuố p n củathanhdomộtđơnvị ựctácđộngl t iạ 1.Tươ tng ựtacócácmơ ảchocộ thứt t haivàthứ

<b>Bâygiờ,tabiểudiễne</b><sub>1</sub>nhưlàmộtvéctơuố Khiđón. <small>1</small>

<i>D</i><small>-</small>baogồm cáclực sinh rasựuốn.   Ngh là,ĩa cộtđầutiêncủa <small>1</small>

<i>D</i><small>-</small>baogồmcác lựcđược tácđộng   tại bađiểmđểsinhramộtđơn vịuốn t i ạ điểm1vàđộuốnbằng không tạicác iểm  đ khác. Tương tự ư  nh vậy,cột 2 và3của 

<i>D</i><small>-</small>baogồm những lực  cầnđếnđểtạo ra mộtđơn  vịuố tn ươngứng ạ đ ểm2và3.Tạ mỗ t i i i i cột,mộthoặ hailự phảitráidấu(hướnglêntrên)đểc c tạo ra mộtđơnvịuốn tại iểmtương đ ứngvàđộuốnbằng khơng tạihai i  đểmcịn ại.Nếuđộl uốn ođược, chđ  ẳnghạ bằninchen, trênpound,thìcácphần t ửcủamatrậ cn ứngđượ đobằngpoundc trêninche.

<i>2. NếuA, B</i>làcácmatrậncỡ<i>n</i>´<i>n</i>khảnghịchthì<i>AB</i>cũng khảnghịch, vànghịchđảo  của<i>AB</i>là tích ủa nghịch đảoc <i> của A và B</i> theo thứ ự ngượct lại, nghĩ a là

Mộtmatrậncơsởlàmộtmatrậnnhậnđượcbằngcáchthự hiệnc mộ phépbiế đổisơcấpt n hàngcơbảntrênmatrậ đơ vị.Vídụsauđâyminh ọaba ạin n h lo matrậ cơsở.n

Tíchcácmatrậ nkhảnghịchlàkhảnghịch,vànghịchđảolàtíchcác nghịchđảotheochiềungược ạl i.

</div><span class="text_page_counter">Trang 10</span><div class="page_container" data-page="10">

10  Hãytính<i>E A E A E Avà</i><sub>1</sub> , <sub>2</sub> , <sub>3</sub> mơ tảlàm thếnào mà các tích này lại cóthểnhậnđượcbằngcác        phépbiếnđổ sơci <i>ấphàng ơbả trênA.</i>c n

<b>GIẢI</b>Tathấy rằng 

 Cộng 4- lầnhàng <i>1 củaA vàohàng 3</i>  choma trận <i>E A(đây </i><sub>1</sub> là một toán tửthaythếhàng).  Đổichỗhàng 1 và hàng 2 cho ta      <i>E A</i><sub>2</sub> <i>,vànhânhàng3củaAvới5tađượcE A.</i><sub>3</sub>

Phépnhântráibởimatrậ n<i>E</i><small>1</small>trong <i>Ví dụ5 có tác dụng tương tựtrên bất kỳma trận 3 n</i>        ´  nào.Nócộng 4- lần hàng 1vào hàng3 của ma trậnđó.Nóiriêng,từ    <i>E I</i><sub>1</sub> =<i>E</i><sub>1</sub>,tacó ểth thấyrằng<i>E</i><sub>1</sub>tựnóđượcchobởiphépbiế đổisơcấphàngtn ương ự ừmatrậ đơnt t n vị.Bở i

</div><span class="text_page_counter">Trang 11</span><div class="page_container" data-page="11">

Trongthựchành, hiếmkhiđượctính,trừkhi cầta nđến cácphầntử của .Đểtính ảc và mấ  ấ t g p ba ầ l n sốphéptốngiải  bằngphươngphápbiếnđổisơcấphàng.

<b>THUẬTTỐNTÌM</b><i>A</i> 

Biến đổi hàng ma trậ  mở rộng n [<i>A I</i>]. Nếu <i>Atương đương hàng với I thì</i>[<i>A I</i> ]tương đươnghàngvới [<i>I A</i><small>-1</small>].Ngượclại, <i>A</i>khơngcónghịchđảo.

<i>A</i><small>-</small>chínhxáclà các nghiệm củacác hệtrong (2). Quan sát nàylà có íchvì          mộtsốvấnđềứngdụ cóng thể ỉucầutìmmột ặcch ho hai ột ủac c <i>A</i><small>-1</small>.Trongtrườnghợp đó,chỉnhững hệtươngứng trong (2)   cần phải giải mà thôi.   

   

</div>

×