Video on Hill48 yield function parameter optimization:
• Hill48 anisotropic yie...
MATLAB Code for Hill48 yield function parameter optimization:
www.mathworks.com/matlabcentr...
MATLAB code for Yield surface plot in this video:
%% Von Mises
syms x y % x=sigma1; y=sigma2
YF=sqrt((x^2+y^2+(x-y)^2)/2);
figure (1)
Mises = fcontour(YF,'k')
Mises.LevelList = 10 % Initial Yield stress
Mises.YRange = [-20,20];
Mises.XRange = [-20,20];
%% Hill48
F=0.8404; G=0.7324; H=0.2677;
YF=(G+H)*x^2+(F+H)*y^2-2*H*x*y;
hold on
set(gca,'fontsize', 18)
Hill = fcontour(YF,'b')
Hill.LevelList = 100
Hill.YRange = [-20,20];
Hill.XRange = [-20,20];
%% Experimental
grid on
plot([0,10],[9.5,0],'*r')
legend ('von Moses','Hill 48','Experimental')