WTS (formerly Leonov) + GUI: d=dev(a) - 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.

d=dev(a)
function d=dev(a) %This function determines the deviatoric part of a d=a- (1/3) trace(a)eye(3);