WTS (formerly Leonov) + GUI: dp=dplas(sig,tau0,A0,H0,T) - File Exchange

Code covered by the BSD License

Highlights from
WTS (formerly Leonov) + GUI

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.

dp=dplas(sig,tau0,A0,H0,T)

function dp=dplas(sig,tau0,A0,H0,T)

%Calculates the plastic strain rate as a function of the stress tensor sig
if taueq(sig)==0
    dp=dev(sig)/ (2 * A(A0,H0,T) * tau0);
else
    dp=dev(sig)* sinh(taueq(sig)/tau0) / (2 * A(A0,H0,T) * taueq(sig));
end

Highlights from WTS (formerly Leonov) + GUI

Highlights from
WTS (formerly Leonov) + GUI