O metodě konečných prvkůLect_7.ppt

M. Okrouhlík

Ústav termomechaniky, AV ČR, PrahaPlzeň, 2010

Řešení statických úloh

kruk =[2 1 3 2 4 3 5 4]

okr = [1 2 9 10]

% 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]

E:\edu_mkp_liberec_2\mtl_prog\prut_statika.m

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

vsesily =

1.0e+003 *

0.2887 0.5000 0 0 0.0000 -1.0000 0.0000 0 -0.2887 0.5000

0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2-0.1

0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

Podrobné odvození:

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

…………………….

% 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

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

2 4 6 8 100

0.2

0.4

0.6

0.8

1

1.2x 10

-3 pruhyby na hornim povrchu [m]

poradova cisla uzlu na povrchu

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

Four-node bilinear element, plane stress, full integration

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

Jak na to?