0001 function plot_drivingstress(md,options_structure,width,i);
0002
0003
0004
0005
0006
0007
0008
0009
0010 [sx sy s]=drivingstress(md);
0011
0012
0013 [x y z elements is2d]=processmesh(md,options_structure);
0014 [dstress isongrid]=processdata(md,s,options_structure);
0015 dstress=dstress/1000;
0016
0017
0018 subplot(width,width,i);
0019
0020
0021 if ~isnan(options_structure.edgecolor),
0022 edgecolor=options_structure.edgecolor;
0023 else
0024 edgecolor='none';
0025 end
0026
0027
0028 if ~isongrid
0029 if is2d
0030 A=elements(:,1); B=elements(:,2); C=elements(:,3);
0031 patch( 'Faces', [A B C], 'Vertices', [x y z],'CData', dstress,'FaceColor','flat','EdgeColor',edgecolor);
0032 else
0033 A=elements(:,1); B=elements(:,2); C=elements(:,3); D=elements(:,4); E=elements(:,5); F=elements(:,6);
0034 patch( 'Faces', [A B C], 'Vertices', [x y z],'CData', dstress,'FaceColor','flat','EdgeColor',edgecolor);
0035 patch( 'Faces', [D E F], 'Vertices', [x y z],'CData', dstress,'FaceColor','flat','EdgeColor',edgecolor);
0036 patch( 'Faces', [A B E D], 'Vertices', [x y z],'CData', dstress,'FaceColor','flat','EdgeColor',edgecolor);
0037 patch( 'Faces', [B E F C ], 'Vertices', [x y z],'CData', dstress,'FaceColor','flat','EdgeColor',edgecolor);
0038 patch( 'Faces', [C A D F ], 'Vertices', [x y z],'CData', dstress,'FaceColor','flat','EdgeColor',edgecolor);
0039 end
0040
0041 elseif isongrid
0042 if is2d
0043 A=elements(:,1); B=elements(:,2); C=elements(:,3);
0044 patch( 'Faces', [A B C], 'Vertices', [x y z],'FaceVertexCData', dstress,'FaceColor','interp','EdgeColor',edgecolor);
0045 else
0046 if options_structure.layer>=1,
0047 A=elements(:,1); B=elements(:,2); C=elements(:,3);
0048 patch( 'Faces', [A B C], 'Vertices', [x y z],'FaceVertexCData', dstress,'FaceColor','interp','EdgeColor',edgecolor);
0049 else
0050 A=elements(:,1); B=elements(:,2); C=elements(:,3); D=elements(:,4); E=elements(:,5); F=elements(:,6);
0051 patch( 'Faces', [A B C], 'Vertices', [x y z],'FaceVertexCData', dstress,'FaceColor','interp','EdgeColor',edgecolor);
0052 patch( 'Faces', [D E F], 'Vertices', [x y z],'FaceVertexCData', dstress,'FaceColor','interp','EdgeColor',edgecolor);
0053 patch( 'Faces', [A B E D], 'Vertices', [x y z],'FaceVertexCData', dstress,'FaceColor','interp','EdgeColor',edgecolor);
0054 patch( 'Faces', [B E F C ], 'Vertices', [x y z],'FaceVertexCData', dstress,'FaceColor','interp','EdgeColor',edgecolor);
0055 patch( 'Faces', [C A D F ], 'Vertices', [x y z],'FaceVertexCData', dstress,'FaceColor','interp','EdgeColor',edgecolor);
0056 end
0057 end
0058 end
0059
0060
0061 if isnan(options_structure.title)
0062 options_structure.title='Driving stress [kPa]';
0063 end
0064 applyoptions(md,dstress,options_structure);