WTS (formerly Leonov) + GUI: s=sigma(J,Be,K,G) - 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.

s=sigma(J,Be,K,G)
function s=sigma(J,Be,K,G) %Calculates the Cauchy stress tensor s=K(J-1)eye(3)+G*dev(Be);