function [sxx,sxy,syx,syy] = stress_trans(srr,srt,str,stt,axr,axt,ayr,ayt)
% 2-D stress transformation function
sxx = axr.*axr.*srr + axr.*axt.*srt + axt.*axr.*str + axt.*axt.*stt;
sxy = axr.*ayr.*srr + axr.*ayt.*srt + axt.*ayr.*str + axt.*ayt.*stt;
syx = ayr.*axr.*srr + ayr.*axt.*srt + ayt.*axr.*str + ayt.*axt.*stt;
syy = ayr.*ayr.*srr + ayr.*ayt.*srt + ayt.*ayr.*str + ayt.*ayt.*stt;