Create a circular 3D surf plot

i have the following code:
%% Assign constant
v=0.33;% poisson ratio
E=73.1*10^9; % young's modulus [N]
ti=0.1185;% thickness of bulkhead intact [m]
td=0.1185;% thickness of bulkhead Damage[m]
d=3.11; % diameter of bulkhead[m]
ai=d/2; % radius of bulkhead[m]
r= linspace(-ai,ai) ;
z=linspace(0,ti);
q0=5; % pressure distribution [N/m^2]
%% compute stress for Bulkhead intact state
Mr=0.0625*q0*(3+v)*(ai^2-r.^2);
Sigmari = (Mr.*z)/(ti^3/12);
the output " Sigmari" reflects the stress at some point on a circular plate as a function of thr radius r and the thickness z. I want to create a plot that would the plate and the shear distribution on the plate with a color legend.
Thank you in advance.

回答(1 个)

%% Assign constant
v=0.33;% poisson ratio
E=73.1*10^9; % young's modulus [N]
ti=0.1185;% thickness of bulkhead intact [m]
td=0.1185;% thickness of bulkhead Damage[m]
d=3.11; % diameter of bulkhead[m]
ai=d/2; % radius of bulkhead[m]
r = linspace(-ai,ai).' ;
z=linspace(0,ti);
q0=5; % pressure distribution [N/m^2]
%% compute stress for Bulkhead intact state
Mr=0.0625*q0*(3+v)*(ai^2-r.^2);
Sigmari = (Mr.*z)/(ti^3/12);
[x, y] = pol2cart(r, z);
pcolor(x, y, Sigmari);
axis equal

类别

帮助中心File Exchange 中查找有关 Animation 的更多信息

产品

版本

R2021b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by