Parent Directory
|
Revision Log
|
Revision Graph
initial checkin of Timour's MatLAB scripts
1 | edhill | 1.1 | Heat=1.0; |
2 | Nx=10; | ||
3 | Ny=20; | ||
4 | for i=1:Nx | ||
5 | for j=1:Ny | ||
6 | heating(i,j)=-Heat*(sin(pi*(2*j/Ny)));; | ||
7 | end | ||
8 | end | ||
9 | Stress=1.0; | ||
10 | Nx=10; | ||
11 | Ny=20; | ||
12 | for i=1:Nx | ||
13 | for j=1:Ny | ||
14 | stress(i,j)=-Stress*(cos(pi*(2*j/Ny))); | ||
15 | end | ||
16 | end | ||
17 | for i=1:Nx | ||
18 | for j=1:Ny | ||
19 | stressy(i,j)=-Stress*(cos(pi*(2*j/Ny))); | ||
20 | end | ||
21 | end | ||
22 | for i=1:Nx | ||
23 | for j=1:Ny | ||
24 | stressx(i,j)=-Stress*(cos(pi*(2*i/Nx)));; | ||
25 | end | ||
26 | end | ||
27 | stress=stressx+stressy; |
ViewVC Help | |
Powered by ViewVC 1.1.22 |