Changeset 9423 for issm/trunk/src/m/model/mechanicalproperties.m
- Timestamp:
- 08/25/11 15:11:19 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/m/model/mechanicalproperties.m
r8298 r9423 4 4 % this routine computes the components of the stress tensor 5 5 % strain rate tensor and their respective principal directions. 6 % the results are in the model md: md. stress and md.strainrate6 % the results are in the model md: md.results 7 7 % 8 8 % Usage: … … 94 94 stress.principalaxis1=directionsstress(:,3:4); 95 95 stress.effectivevalue=1/sqrt(2)*sqrt(stress.xx.^2+stress.yy.^2+2*stress.xy.^2); 96 md. stress=stress;96 md.results.stress=stress; 97 97 98 98 strainrate=struct('xx',[],'yy',[],'xy',[],'principalvalue1',[],'principalaxis1',[],'principalvalue2',[],'principalaxis2',[],'effectivevalue',[]); … … 105 105 strainrate.principalaxis1=directionsstrain(:,3:4); 106 106 strainrate.effectivevalue=1/sqrt(2)*sqrt(strainrate.xx.^2+strainrate.yy.^2+2*strainrate.xy.^2); 107 md. strainrate=strainrate;107 md.results.strainrate=strainrate; 108 108 109 109 deviatoricstress=struct('xx',[],'yy',[],'xy',[],'principalvalue1',[],'principalaxis1',[],'principalvalue2',[],'principalaxis2',[],'effectivevalue',[]); … … 116 116 deviatoricstress.principalaxis1=directionsstress(:,3:4); 117 117 deviatoricstress.effectivevalue=1/sqrt(2)*sqrt(stress.xx.^2+stress.yy.^2+2*stress.xy.^2); 118 md. deviatoricstress=deviatoricstress;118 md.results.deviatoricstress=deviatoricstress;
Note:
See TracChangeset
for help on using the changeset viewer.