0001 function pressure=PressureElemCompute(m,inputs,type);
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011 global cluster gridset
0012
0013
0014 elements=m.elements;
0015 grids=m.grids;
0016 materials=m.materials;
0017 loads=m.loads;
0018 gridset=m.gridset;
0019
0020
0021 [n1,n2]=GetNumberOfActiveElements(elements);
0022
0023
0024 pressure=zeros((n2-n1)+1,1);
0025
0026
0027 if strcmpi(type,'2d')
0028 disp('pressure not computed for 2d meshes')
0029 return
0030 end
0031
0032
0033 for n=n1:n2,
0034 if ~isempty(elements(n).element),
0035 stress_tensor=Stress(elements(n).element,grids,materials,inputs);
0036 trace=stress_tensor(1)+stress_tensor(2)+stress_tensor(3);
0037 pressure(n)=-1/3*trace;
0038 end
0039 end