WTS (formerly Leonov) + GUI: t=taueq(sigma) - File Exchange

Code covered by the BSD License

WTS_GUI(varargin) End initialization code - DO NOT EDIT
A=A(A0,H0,T)
[T33,truestrain,J,S33,R33]=WT... This function determines the total true axial stress in a glassy polymer in uniaxial deformation according to the WTS model
d=dev(a) This function determines the deviatoric part of a
d=symm(l) Calculates the symmetric part of l
dp=dplas(sig,tau0,A0,H0,T) Calculates the plastic strain rate as a function of the stress tensor sig
dy=evolution(time,y,strainrat...
f=B(epsilon,poisson); isochoric part of the chauchy green tensor
f=F(epsilon,poisson); deformation gradient tensor F
f=Vact(a,b,truestrain) In this model, the activation volume is a function of invariants of the strain tensor. A Poisson's ratio of 0.5 is assumed.
s=sigma(J,Be,K,G) Calculates the Cauchy stress tensor
t=taueq(sigma) This function calculates the equivalent stress
w=asym(l) Calculates the anti-symmetric part of l
example.m
View all files

from WTS (formerly Leonov) + GUI by Michael Wendlandt
Calculates stress vs. strain curves of plastically deformed glassy polymers. GUI available.

t=taueq(sigma)
function t=taueq(sigma) %This function calculates the equivalent stress t=sqrt(0.5trace(dev(sigma)dev(sigma)));