3D stress精选.ppt

3D stress Mostly from Ragan and Schultz, unpublished Tractions from stress Traction normal and shear components Traction from principal stresses Full stress tensor rotation example Example application to South Mountains faults * Note sign convention change Tractions from stress T=S*N (MATLAB) sigma=dot(T,N) (MATLAB) Traction normal and shear components S = [-50 -20 10; -20 -30 -15; 10 -15 -120] [l,m,n] =plunge_trend_to_dir_cosines(30,290); N=[l;m;n]; T=S*N; %equation 13.11 T_mag = sqrt(sum(T.^2)); %normalize components of T to get its direction cosines lt=T(1)./T_mag; mt = T(2)./T_mag; nt = T(3)./T_mag; %plot traction vector [plunge, trend] = dir_cosines_to_plunge_trend2(lt, mt, nt); %we know the orientation of the normal traction, but what is its magnitude? sigma = dot(T,N);%equation 13.13 %Now for the shear traction; use the McKenzie construction B = cross(T,N); %vector normal to the plane containing T and N B_mag = sqrt(B(1)^2 + B(2)^2 + B(3)^2); lb = B(1)./B_mag; mb = B(2)./B_mag; nb = B(3)./B_mag; Ts = cross(N,B); %shear traction direction Ts_mag = sqrt(Ts(1)^2 + Ts(2)^2 + Ts(3)^2); Ts(1) = Ts(1)./Ts_mag; Ts(2) = Ts(2)./Ts_mag; Ts(3) = Ts(3)./Ts_mag; %lets check that the normal and shear are components of the traction testmag = sqrt(sum(sigma.^2 + Ts_mag.^2)); traction vector components are 6.4660 10.9900 -44.8311 traction magnitude 46.6091 traction vector direction cosines 0.1387 0.2358 -0.9619 traction plunge = 74.1 trend = 239.5 normal traction mag -29.44 shear traction mag 36.13 check that components make same length as traction: 46.6091 =?= 46.6091 N, sigma tau T B Traction from principal stresses alpha1 = 45; alpha2 = 60; S = [120 0 0; 0 80 0; 0 0 20] l = cosd(alpha1); m = cosd(alpha2); n = sqrt(1 - l.^2 - m.^2); N =[l;m;n]; [plunge, trend] = dir_cosines_to_plunge_trend(l, m, n);%pole to plane ld = -l; md = -m; nd = cosd(asind(n)); [dip, dipdir] = dir_cosines_to_plunge_trend(ld, md, nd);%plane itself T=S*N; %equation 13.11 T