暖阳三杯
2008-01-26, 22:37
这是函数中的命令:
function PlotStress(iStress)
%显示应力云图
%iStress----应力分量指示,可以为:
%1----x方向正应力
%2----y方向正应力
%3----剪应力
%4----第一主应力
%5----第二主应力
%6----最大剪应力
global gNode gElement gNodeStress
switch iStress
case 1
title='x方向正应力';
case 2
title='y方向正应力';
case 3
title='剪应力';
case 4
title='最大主应力';
case 5
title='最小主应力';
case 6
title='最大剪应力';
end
figure;
axis equal;
axis off;
set(gcf,'NumberTitle','off');
set(gcf,'Name',title);
Minstress=gNodeStress(:,iStress)
Maxstress=max( gNodeStress(:,iStress));
caxis([Minstress,Maxstress]);
colormap('gray');
当运行PlotStress(4)
则出现下列错误提示:
??? Attempted to access gNodeStress(:,4); index out of bounds because size(gNodeStress)=[0,0].
Error in ==> PlotStress at 32
Minstress=gNodeStress(:,iStress)
不明白哪儿出问题了,请高手帮忙看看,谢谢啊!!
function PlotStress(iStress)
%显示应力云图
%iStress----应力分量指示,可以为:
%1----x方向正应力
%2----y方向正应力
%3----剪应力
%4----第一主应力
%5----第二主应力
%6----最大剪应力
global gNode gElement gNodeStress
switch iStress
case 1
title='x方向正应力';
case 2
title='y方向正应力';
case 3
title='剪应力';
case 4
title='最大主应力';
case 5
title='最小主应力';
case 6
title='最大剪应力';
end
figure;
axis equal;
axis off;
set(gcf,'NumberTitle','off');
set(gcf,'Name',title);
Minstress=gNodeStress(:,iStress)
Maxstress=max( gNodeStress(:,iStress));
caxis([Minstress,Maxstress]);
colormap('gray');
当运行PlotStress(4)
则出现下列错误提示:
??? Attempted to access gNodeStress(:,4); index out of bounds because size(gNodeStress)=[0,0].
Error in ==> PlotStress at 32
Minstress=gNodeStress(:,iStress)
不明白哪儿出问题了,请高手帮忙看看,谢谢啊!!