1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
|
isomat 0 E 1000 nu 0.0 { density = 1234;};
Default flxelement::thickness = 0.1;
Default flxelement::planestr::gausspoints = 4;
Default flxelement::planestr::sri::gausspoints = 1;
const F = 0.01;
const ex = 400;
const ey = 10;
const dx = 10/ex;
const dy = 0.2/ey;
timer define tim1;
timer start tim1;
for ( i = 0; ey+1-i; i+1) {
node i [0,dy*i,0] fix wx fix wy fix wz;
};
for ( j = 1; ex+1-j; j+1) {
for ( i = 0; ey+1-i; i+1) {
node j*(ey+1)+i [dx*j,dy*i,0] fix wz;
};
};
for ( j = 0; ex-j; j+1) {
for ( i = 0; ey-i; i+1) {
elplanestr j*ey+i j*(ey+1)+i (j+1)*(ey+1)+i (j+1)*(ey+1)+i+1 j*(ey+1)+i+1;
};
};
nload ex*(ey+1) wy 1*F/(ey*2);
for ( i = 1; ey-i; i+1) {
nload ex*(ey+1)+i wy 1*F/(ey);
};
nload ex*(ey+1)+ey wy 1*F/(ey*2);
fem assdof;
fem assstf;
fem assforces;
fem solve {
printsol=false;
eps=1e-5;
iter=100000;
pcn=1;
};
timer stop tim1;
timer print tim1;
timer delete tim1;
#for ( i = 0; ey+1-i; i+1) {
# calc ex*(ey+1)+i;
# calc nodeu(ex*(ey+1)+i,wx);
# calc nodeu(ex*(ey+1)+i,wy);
#};
filestream res ("res.txt");
elplot x,y,z,E,n,G {stream=res;};
# reference displacement: 50
|