%%*************************************************************************
%StVenant_Kirchhoff
[Dtan,stress_gp,stress_vmisses_gp,strain_gp,strain_vmisses_gp,stress_trial]=...
    StVenant_Kirchhoff_response(StVenant_Kirchhoff,material_type_e,F,...
    stress_gp,stress_vmisses_gp,strain_gp,strain_vmisses_gp,I,L,L1,IdyadI);

%%*************************************************************************
%Neo_Hookean_Compressible_response
[Dtan,stress_gp,stress_vmisses_gp,strain_gp,strain_vmisses_gp,stress_trial]=...
    NeoHookean_Compressible_response(NeoHookean_Compressible,material_type_e,F,...
    stress_gp,stress_vmisses_gp,strain_gp,strain_vmisses_gp);

%%*************************************************************************
%Neo_Hookean_Incompressible_respons
[Dtan,stress_gp,stress_vmisses_gp,strain_gp,strain_vmisses_gp,stress_trial] =...
    NeoHookean_Incompressible_response(NeoHookean_Incompressible,material_type_e,F,...
    stress_gp,stress_vmisses_gp,strain_gp,strain_vmisses_gp);

%%*************************************************************************
%Ferrara-Pandolfi_Incompressible_response
vecNormal=zeros(1,3);
vecN1=zeros(1,3);
Gpnt = N' * node(sctr,:); % global GP
vecNormal(1,1) = Gpnt(1); vecNormal(1,2) = Gpnt(2); vecNormal(1,3) = 0; % Normal vector of Gauss point
vecN1(1,:)=[0 0 1];
vecMbar=cross(vecN1(1,:),vecNormal(1,:)); % local circumferential axis
vecNbar1=cross(vecNormal(1,:),cross(vecN1(1,:),vecNormal(1,:))); % Axial geometry vector
Mbar=vecMbar/norm(vecMbar);
Nbar1=vecNbar1/norm(vecNbar1);

[Dtan,stress_gp,stress_vmisses_gp,strain_gp,strain_vmisses_gp,J3,stress_trial] =...
    FerraraPandolfi_Incomp_response(Ferrara_Pandolfi_Macro,material_type_e,...
    F,stress_gp,stress_vmisses_gp,strain_gp,strain_vmisses_gp,Mbar,Nbar1);

%%*************************************************************************
%RVE_response (Homogenization)
[Dtan,stress_gp,stress_vmisses_gp,strain_gp,strain_vmisses_gp,RVEu_gp,stress_trial]=...
    RVE_response(RVEVol,RVEnumelem,RVEelement,RVEnode,RVEu_gp,RVETDOF,RVEact,RVE_material_type,RVEIndexer,...
    stress_gp,F,strain_gp,stress_vmisses_gp,strain_vmisses_gp,...
    RVEW,RVEQ,RVEelemType,...
    StVenant_Kirchhoff,NeoHookean_Compressible,NeoHookean_Incompressible,Ferrara_Pandolfi,...
    I,L,Lz,L1,IdyadI,max_iter,global_toll);