Last change
on this file since 4 was
4,
checked in by zhangj, 11 years ago
|
Initial import--MML version from SOLEIL@2013
|
File size:
444 bytes
|
Line | |
---|
1 | function pass = circplate |
---|
2 | % Deflection of a circular plate |
---|
3 | |
---|
4 | tol = chebfunpref('eps'); |
---|
5 | |
---|
6 | d = domain(0,1); |
---|
7 | f = chebfun('exp(-50*r)',d); % loading function |
---|
8 | |
---|
9 | r = chebfun('r',d); |
---|
10 | D = diff(d); I = eye(d); |
---|
11 | B = diag(r.^3)*D^4 + diag(2*r.^2)*D^3 - diag(r)*D^2 + D; % r^3*(biharm) |
---|
12 | |
---|
13 | B.lbc(1)=D; B.lbc(2)=D^3; % symmetry at origin |
---|
14 | B.rbc(1)=I; B.rbc(2)=D; % clamped at boundary |
---|
15 | |
---|
16 | u1 = B\(r.^3.*f); |
---|
17 | load circplate |
---|
18 | |
---|
19 | pass = norm(u-u1) < 1e-10*(tol/eps); |
---|
Note: See
TracBrowser
for help on using the repository browser.