1 / 19

O metod ě konečných prvků Lect_ 7

O metod ě konečných prvků Lect_ 7 .ppt. Řešení statických úloh. M. Okrouhlík Ústav termomechaniky, AV ČR , Praha Plze ň , 2010. % kodova cisla d=[1 2 3 4; 1 2 5 6; 3 4 5 6; 3 4 7 8; 5 6 7 8; 5 6 9 10; 7 8 9 10]. kruk =[2 1 3 2 4 3 5 4] okr = [1 2 9 10].

galen
Download Presentation

O metod ě konečných prvků Lect_ 7

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. O metodě konečných prvkůLect_7.ppt Řešení statických úloh M. Okrouhlík Ústav termomechaniky, AV ČR, Praha Plzeň, 2010

  2. % kodova cisla d=[1 2 3 4; 1 2 5 6; 3 4 5 6; 3 4 7 8; 5 6 7 8; 5 6 9 10; 7 8 9 10] kruk =[2 1 3 2 4 3 5 4] okr = [1 2 9 10]

  3. E:\edu_mkp_liberec_2\mtl_prog\prut_statika.m vsesily = 1.0e+003 * 0.2887 0.5000 0 0 0.0000 -1.0000 0.0000 0 -0.2887 0.5000 pos = 1.0e-005 * 0.0700 -0.2021 0.0000 -0.4042 -0.0700 -0.2021 sila = 288.6751 0.0000 -288.6751 -577.3503 -288.6751 -0.0000 288.6751

  4. Podrobné odvození: Stejskal, V., Okrouhlík, M.: Kmitání s Matlabem, Vydavatelství ČVUT, Praha 2002, ISBN 80-01-02435-0

  5. residuals and external forces resi = Columns 1 through 5 -3.6380e-012 1.8190e-011 0 1.0914e-011 -2.5466e-011 Columns 6 through 10 -7.2760e-012 -1.0914e-011 7.2760e-012 7.2760e-012 3.6380e-012 Columns 11 through 15 5.6389e-011 0 4.3656e-011 7.2760e-011 1.7462e-010 Columns 16 through 20 -7.2760e-012 -1.4552e-010 -2.9104e-011 -1.8190e-012 -2.1828e-011 Columns 21 through 25 -8.7311e-011 4.3656e-011 9.0949e-012 -4.3656e-011 2.9104e-010 Columns 26 through 30 3.6380e-011 2.6193e-010 7.2760e-011 2.7649e-010 -1.4552e-010 Columns 31 through 35 -2.0009e-010 -1.8917e-010 1.0914e-010 1.4552e-011 -1.0550e-010 Columns 36 through 40 5.8208e-011 5.0000e+002 -5.4570e-011 5.0000e+002 1.4552e-010 ……………………. % static solution dis = xk\p; % check clc r=xk*dis; % infate displacements dist=[0 0 0 0 dis']; disp('inflated displacements') disf=dist point=[1 7:4:2*lmax-1]; displo=disf(point); dismax=max(disf) disp('residuals and external forces') resi = r' plot(displo); grid title('pruhyby na hornim povrchu [m]') xlabel('poradova cisla uzlu na povrchu') % pause % end of beam1.m

  6. E:\edu_mkp_liberec_2\mtl_prog\beam_statika\beam1.m

  7. Four-node bilinear element, plane stress, full integration

  8. Jak na to? www.it.cas.cz/cs/elektronicka-kniha-numerical-methods-computation-mechanics E:\edu_mkp_liberec_2\pdf_jpg_my_old_texts\skripta_jaderna\aplik_mechanika_kontinua_1989.pdf

More Related