%Solving not linear task to destroy contact at using finite element method.
%In program uses the rectangular multiplex element.

%
%Initial data (input data)

    %Parameters material contactive parts
        %сталь 20л/ steel 430A
        mu=0.25;        %(0,24..0,28) Poisson's coeff.
        sigma_cr=216e6; %sigma_0.2
        E=200e9;        %Young's modulus (for this task elastic modulus)

    %Body geometry
        L=1;    %length
        H;    %height parts body

    %Grid options
        N=11;   %points at length
        M=11;   %points at height part body

%
%Elements for the program

    element_number = zeros((M-1)*2+1 ,N-1);
    point_global_number_first_part = zeros(N*M ,4);
    point_global_number_second_part = zeros(N*M ,4);
    point_global_number_contact_part = zeros(2*N ,4);

    coordinate_first_part = zeros(N*M ,4);
    coordinate_second_part = zeros(N*M ,4);
    coordinate_contact_part = zeros(2*N ,4);

    sigma_first_part = zeros(N*M ,4);
    sigma_second_part = zeros(N*M ,4);
    sigma_contact_part = zeros(2*N ,4);

%
%Program body

    %