WTS (formerly Leonov) + GUI: f=Vact(a,b,truestrain) - 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.

f=Vact(a,b,truestrain)

function 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.
%Units: Vact[nm3]

f=(a - b *  (1/3) * (exp(-2 * truestrain) + 2 * exp(truestrain)).^(0.5) .* (exp(2 * truestrain) + 2 * exp(-truestrain)).^(0.5) );


%check if Vact gets negative. If YES show error message !
if f<=0
    error('NEGATIVE activation volume Vact ! Please reedit the depedence of Vact on strain such that Vact stays positive')
end

Highlights from WTS (formerly Leonov) + GUI

Highlights from
WTS (formerly Leonov) + GUI