Диссертация (1141568), страница 21
Текст из файла (страница 21)
Подготовка к испытаниюРисунок Б.8 – Балка БС-II-0 после испытанияРисунок Б.9 – Балка БС-II-1. Подготовка к испытанию184Рисунок Б.10 – Балка БС-II-1 после испытанияРисунок Б.11 – Балка БС-II-2. Подготовка к испытаниюРисунок Б.12 – Балка БС-II-2 после испытания185Рисунок Б.13 – Балка БС-III-0. Подготовка к испытаниюРисунок Б.14 – Балка БС-III-0 после испытанияРисунок Б.15 – Балка БС-III-1.
Подготовка к испытанию186Рисунок Б.16 – Балка БС-III-1 после испытанияРисунок Б.17 – Балка БС-III-2. Подготовка к испытаниюРисунок Б.18 – Балка БС-III-2 после испытания187Приложение ВПрограмма «Нелинейный расчет железобетонных элементов с косвеннымармированием». Язык: встроенный в ПК Matlabclearclc%1.ПОСТРОЕНИЕ ДИАГРАММ%1.1.Исходные данные%1.1.1.Параметры бетонаdisp('ПОСТРОЕНИЕ ДИАГРАММ:')disp('ПАРАМЕТРЫ БЕТОНА:')disp('Прочность бетона на сжатие (МПа):')Rb=input('Rb=');disp('Модуль упругости бетона (МПа):')Eb0=input('Eb0=');disp('Относительные деформации в вершине диаграммы сжатия бетона (по умолчаниюepsilon_b0=0.002):')epsilon_b0=input('epsilon_b0=');disp('Предельные относительные деформации сжатия бетона (по умолчаниюepsilon_bu=0.0035):')epsilon_bu=input('epsilon_bu=');disp('Начальный коэффициент упругости бетона (по умолчанию vb0=1):')vb0=input('vb0=');disp('Начальный коэффициент упругости бетона (по умолчанию ub0=0.2):')ub0=input('ub0=');%1.1.2.Параметры косвенного армированияdisp('ПАРАМЕТРЫ КОСВЕННОГО АРМИРОВАНИЯ:')disp('Прочность стали сеток косвенного армирования на растяжение (МПа):')Rsxy=input('Rsxy=');disp('Модуль упругости стали сеток косвенного армирования (МПа):')Esxy=input('Esxy=');disp('Коэффициент косвенного армирования:')mu_xy=input('mu_xy=');%1.2.Построение диаграммы сжатия бетона (Карпенко)vbu=Rb/Eb0/epsilon_b0;omega_1=2-2.5*vbu;omega_2=1-omega_1;for i=1:30sigma_b=Rb/30*i;ita=sigma_b/Rb;vb=vbu+(vb0-vbu)*sqrt(1-omega_1*ita-omega_2*ita^2);epsilon_b=sigma_b/vb/Eb0;s(1,i)=sigma_b;e(1,i)=epsilon_b;end188s=[0,s];e=[0,e];%1.2.1.Аппроксимация МНК в виде ?=p1*?b^2+q2*?b%Коэффициенты многочлена:p2=(sum(s)*e(end)^2-s(end)*sum(e.^2))/(sum(e)*e(end)^2-e(end)*sum(e.^2));p1=(s(end)-p2*e(end))/e(end)^2;%1.2.2.Построение графикаplot(e,s);hold on;grid on;plot([epsilon_b0,epsilon_bu],[Rb,Rb]);%1.3.Построение диаграммы сжатия для бетона, усиленного сетками:ro_xy=0.375*mu_xy*Rsxy/Rb; %Конструктивный коэффициентRb3=Rb*((1-ro_xy)/2+sqrt(((1-ro_xy)/2)^2+9*ro_xy)) %Прочность на сжатие%бетона, усиленного сетками (МПа)q1=0.9*log(Rb)-2.61;q2=0.9*(Rb3-Rb)/Rb3+0.1*log(Rb)-1.29;q3=0.1*(Rb3-Rb)/Rb3;m=(-q2-(q2^2-4*q1*q3)^0.5)/(2*q1);k=1/(0.1+0.9*m)+2.9-log(Rb);sigma_bxyu=m*Rb3;%Усилие поперечного обжатия, создаваемое сетками%Коэффициент увеличения деформативности:en=exp((2.9224-0.00408*Rb)*(0.9*sigma_bxyu/Rb)^(0.3124+0.0022*Rb));epsilon_b03=epsilon_b0*en %Относительные деформации в вершине диаграммы%сжатия ограниченного бетонаepsilon_bu3=epsilon_bu*en %Предельные относительные деформации в вершине%диаграммы сжатия ограниченного бетона%1.3.1.Коэффициенты многочлена:epsilon_b13=epsilon_b03/1000;epsilon_b23=0.23*epsilon_b03;gorizont=((-0.38*Rb+3.31)*mu_xy+exp(0.039*Rb-1));if gorizont>=1gorizont=1;endepsilon_b33=gorizont*epsilon_b03sig1=Eb0;sig2=0.62*Rb3;sig3=0;sig4=Rb3;t1=4*epsilon_b13^3;t2=3*epsilon_b13^2;t3=2*epsilon_b13;t4=1;189t5=epsilon_b23^4;t6=epsilon_b23^3;t7=epsilon_b23^2;t8=epsilon_b23;t9=4*epsilon_b33^3;t10=3*epsilon_b33^2;t11=2*epsilon_b33;t12=1;t13=epsilon_b33^4;t14=epsilon_b33^3;t15=epsilon_b33^2;t16=epsilon_b33;epsilonmntt=[t1 t2 t3 t4; t5 t6 t7 t8; t9 t10 t11 t12; t13 t14 t15 t16];sigmamntt=[sig1; sig2; sig3; sig4];tt=inv(epsilonmntt)*sigmamntt;p3=tt(1,1);p4=tt(2,1);p5=tt(3,1);p6=tt(4,1);p7=-0.2*Rb3/(epsilon_bu3-epsilon_b03);p8=Rb3-p7*epsilon_b03;%1.3.2.Построение графикаfor i=1:30epsilon_b3=epsilon_b33/30*i;e3(1,i)=epsilon_b3;s3(1,i)=p3*epsilon_b3^4+p4*epsilon_b3^3+p5*epsilon_b3^2+p6*epsilon_b3;ends3=[0,s3];e3=[0,e3];plot(e3,s3);plot([epsilon_b33,epsilon_b03],[Rb3,Rb3]);plot([epsilon_b03,epsilon_bu3],[Rb3,0.8*Rb3]);%2.РАСЧЕТ НОРМАЛЬНОГО СЕЧЕНИЯ ПО ПРОЧНОСТИ%2.1.Исходные данные%2.1.1.Параметры продольного армирования:disp('РАСЧЕТ НОРМАЛЬНОГО СЕЧЕНИЯ ПО ПРОЧНОСТИ:')disp('ПАРАМЕТРЫ ПРОДОЛЬНОГО АРМИРОВАНИЯ:')disp('Прочность стали продольного армирования на растяжение (МПа):')Rs=input('Rs=');disp('Модуль упругости стали продольного армирования (МПа)')Es=input('Es=');disp('Площадь продольного армирования (см2)')190As=input('As=');disp('Предельные относительные деформации арматуры')epsilon_s2=input('epsilon_s2=');%2.1.2.Параметры сечения:disp('ПАРАМЕТРЫ СЕЧЕНИЯ:')disp('Момент трещинообразования:')Mcrc=input('Mcrc=');disp('Ширина сечения (см):')b=input('b=');disp('Высота сечения (см):')h=input('h=');disp('Суммарная толщина бокового защитного слоя косвенного армирования (см):')c=input('c=');disp('Толщина верхнего защитного слоя для косвенного армирования (см):')c3=input('c3=');disp('Привязка центра тяжести продольного армирования относительно растянутой грани(см):')a=input('a=');%2.2.Расчет сечения:h0=h-a; %Рабочая высота сечения (см)h1=h0-c3; %Рабочая высота сечения после разрушения защитного слоя (см)epsilon_s0=Rs/Es; %Относительные деформации продольной арматуры при%пределе текучести%2.2.1.Задаем массивы для дальнейшего заполнения:%Вводим величины, которые будут пересчитываться в циклеeb1=0; %Деформации крайнего сжатого волокнаeb2=0; %Деформакии крайнего сжатого волокна армированного ядраes=0;%Деформации арматурыx=0;%Высота сжатой зоныx1=0;%Высота сжатой зоны армированного ядраkriviznacr=0; %Кривизна сечения с трещинойM=0;%Момент, воспринимаемый сечением%Создаются пустые однострочные массивы для последующего заполнения%элементами:etap_m=[];eb1_m=[];eb2_m=[];es_m=[];x_m=[];x1_m=[];kriviznacr_m=[];M_m=[];%2.2.2.Вычисление массивов данныхshag=epsilon_bu3/1000; %шаг, с которым будут увеличиваться деформации%бетонаpe=1000;%коэффициент для перевода единиц измерения момента в [кН*м]191while eb1<epsilon_bu3 & es<epsilon_s2eb1=eb1+shag;if eb1<epsilon_b0 & eb2<=epsilon_b0etap=1;%этап 1%Решенме уравненияq1=b*(p1/3*eb1+p2/2)+(b-c)*(p3/5*eb1^3+p4/4*eb1^2+p5/3*eb1+p6/2p1/3*eb1-p2/2);q2=b*c3*(p1*eb1+p2)+Es*As;q3=b*c3^2*(p1*eb1+p2/2)-Es*As*h1;q4=b*p1/3*c3^3*eb1;xxx=[q1 q2 q3 q4];korni=roots(xxx);korni=korni(abs(imag(korni))<1e-5);korni=real(korni);korni=korni(korni>0);x1=min(korni);%Нахождение данных для текущего шагаx=x1+c3;es=(h1-x1)/x1*eb1;eb2=x/(h0-x)*es;kriviznacr=(eb1+es)/h1;if es<=epsilon_s0; %Усилие в растянутой арматуреsig=es*Es;else sig=Rs;Es=Rs/es;endM=((b-c)*(p3/6*eb1^4+p4/5*eb1^3+p5/4*eb1^2+p6/3*eb1-p1/4*eb1^2p2/3*eb1)*x1^2+b*(p1/4*eb2^2+p2/3*eb2)*x^2+sig*As*(h0-x))/pe;%Заполнение массивов данныхeb1_m=[eb1_m,eb1];eb2_m=[eb2_m,eb2];es_m=[es_m,es];x1_m=[x1_m,x1];x_m=[x_m,x];kriviznacr_m=[kriviznacr_m,kriviznacr];etap_m=[etap_m,etap];M_m=[M_m,M];elseif eb1<=epsilon_b0 & eb2>epsilon_b0etap=2;%этап 2%Решенме уравненияq1=b*(p1*epsilon_b0^3/3/eb1+p2*epsilon_b0^2/2/eb1)+(bc)*(p3/5*eb1^4+p4/4*eb1^3+p5/3*eb1^2+p6/2*eb1-p1/3*eb1^2-p2/2*eb1)+b*Rb*(1epsilon_b0/eb1);q2=b*c3*Rb+eb1*Es*As;q3=-h1*eb1*Es*As;xx=[q1 q2 q3];korni=roots(xx);korni=korni(abs(imag(korni))<1e-5);192korni=real(korni);korni=korni(korni>0);x1=min(korni);%Нахождение данных для текущего шагаx=x1+c3;es=(h1-x1)/x1*eb1;eb2=x/(h0-x)*es;kriviznacr=(eb1+es)/h1;if es<=epsilon_s0; %Усилие в растянутой арматуреsig=es*Es;else sig=Rs;Es=Rs/es;endx0=epsilon_b0/eb1*x1;M=((b-c)*(p3/6*eb1^4+p4/5*eb1^3+p5/4*eb1^2+p6/3*eb1-p1/4*eb1^2p2/3*eb1)*x1^2+b*(p1*epsilon_b0^4/4/eb1^2+p2*epsilon_b0^3/3/eb1^2)*x1^2+b*Rb*(xx0)*(x+x0)/2+sig*As*(h0-x))/pe;%Заполнение массивов данныхeb1_m=[eb1_m,eb1];eb2_m=[eb2_m,eb2];es_m=[es_m,es];x1_m=[x1_m,x1];x_m=[x_m,x];kriviznacr_m=[kriviznacr_m,kriviznacr];etap_m=[etap_m,etap];M_m=[M_m,M];elseif eb1>epsilon_b0 & eb2<=epsilon_buetap=3;%этап 3%Решенме уравненияq1=c*(p1*epsilon_b0^3/3/eb1+p2*epsilon_b0^2/2/eb1)+(bc)*(p3/5*eb1^4+p4/4*eb1^3+p5/3*eb1^2+p6/2*eb1)+Rb*c*(1-epsilon_b0/eb1);q2=b*c3*Rb+eb1*Es*As;q3=-h1*eb1*Es*As;xx=[q1 q2 q3];korni=roots(xx);korni=korni(abs(imag(korni))<1e-5);korni=real(korni);korni=korni(korni>0);x1=min(korni);%Нахождение данных для текущего шагаx=x1+c3;es=(h1-x1)/x1*eb1;eb2=x/(h0-x)*es;kriviznacr=(eb1+es)/h1;if es<=epsilon_s0; %Усилие в растянутой арматуреsig=es*Es;else sig=Rs;Es=Rs/es;endx0=epsilon_b0/eb1*x1;193M=((bc)*(p3/6*eb1^4+p4/5*eb1^3+p5/4*eb1^2+p6/3*eb1)*x1^2+c*(p1*epsilon_b0^4/4/eb1^2+p2*epsilon_b0^3/3/eb1^2)*x1^2+c*Rb*(x1-x0)*(x1+x0)/2+b*c3*Rb*(x1+c3/2)+sig*As*(h0-x))/pe;%Заполнение массивов данныхeb1_m=[eb1_m,eb1];eb2_m=[eb2_m,eb2];es_m=[es_m,es];x1_m=[x1_m,x1];x_m=[x_m,x];kriviznacr_m=[kriviznacr_m,kriviznacr];etap_m=[etap_m,etap];M_m=[M_m,M];elseif eb1<=epsilon_b33etap=4;%этап 4%Решенме уравненияq1=(b-c)*(p3/5*eb1^4+p4/4*eb1^3+p5/3*eb1^2+p6/2*eb1);q2=eb1*Es*As;q3=-h1*eb1*Es*As;xx=[q1 q2 q3];korni=roots(xx);korni=korni(abs(imag(korni))<1e-5);korni=real(korni);korni=korni(korni>0);x1=min(korni);%Нахождение данных для текущего шагаx=x1;es=(h1-x1)/x1*eb1;kriviznacr=(eb1+es)/h1;if es<=epsilon_s0; %Усилие в растянутой арматуреsig=es*Es;else sig=Rs;Es=Rs/es;endM=((b-c)*(p3/6*eb1^4+p4/5*eb1^3+p5/4*eb1^2+p6/3*eb1)*x1^2+sig*As*(h1x1))/pe;%Заполнение массивов данныхeb1_m=[eb1_m,eb1];es_m=[es_m,es];x1_m=[x1_m,x1];x_m=[x_m,x];kriviznacr_m=[kriviznacr_m,kriviznacr];etap_m=[etap_m,etap];M_m=[M_m,M];elseif eb1>epsilon_b33 & eb1<=epsilon_b03etap=5;%Решенме уравнения%этап 5194q1=(bc)*(p3*epsilon_b33^5/5/eb1+p4*epsilon_b33^4/4/eb1+p5*epsilon_b33^3/3/eb1+p6*epsilon_b33^2/2/eb1+Rb3*(1-epsilon_b33/eb1));q2=eb1*Es*As;q3=-h1*eb1*Es*As;xx=[q1 q2 q3];korni=roots(xx);korni=korni(abs(imag(korni))<1e-5);korni=real(korni);korni=korni(korni>0);x1=min(korni);%Нахождение данных для текущего шагаx=x1;es=(h1-x1)/x1*eb1;kriviznacr=(eb1+es)/h1;if es<=epsilon_s0; %Усилие в растянутой арматуреsig=es*Es;else sig=Rs;Es=Rs/es;endx33=epsilon_b33/eb1*x1;M=((bc)*(p3*epsilon_b33^6/6/eb1^2+p4*epsilon_b33^5/5/eb1^2+p5*epsilon_b33^4/4/eb1^2+p6*epsilon_b33^3/3/eb1^2)*x1^2+(b-c)*Rb3*(x1-x33)*(x1+x33)/2+sig*As*(h1-x1))/pe;%Заполнение массивов данныхeb1_m=[eb1_m,eb1];es_m=[es_m,es];x1_m=[x1_m,x1];x_m=[x_m,x];kriviznacr_m=[kriviznacr_m,kriviznacr];etap_m=[etap_m,etap];M_m=[M_m,M];elseif eb1>epsilon_b03%Решенме уравненияq1=(bc)*(p3*epsilon_b33^5/5/eb1+p4*epsilon_b33^4/4/eb1+p5*epsilon_b33^3/3/eb1+p6*epsilon_b33^2/2/eb1+p7/2*(eb1-epsilon_b03^2/eb1)+p8*(1epsilon_b03/eb1)+Rb3/eb1*(epsilon_b03-epsilon_b33));q2=eb1*Es*As;q3=-h1*eb1*Es*As;xx=[q1 q2 q3];korni=roots(xx);korni=korni(abs(imag(korni))<1e-5);korni=real(korni);korni=korni(korni>0);x1=min(korni);%Нахождение данных для текущего шагаx=x1;es=(h1-x1)/x1*eb1;195kriviznacr=(eb1+es)/h1;if es<=epsilon_s0; %Усилие в растянутой арматуреsig=es*Es;else sig=Rs;Es=Rs/es;endx33=epsilon_b33/eb1*x1;x03=epsilon_b03/eb1*x1;M=((bc)*(p3*epsilon_b33^6/6/eb1^2+p4*epsilon_b33^5/5/eb1^2+p5*epsilon_b33^4/4/eb1^2+p6*epsilon_b33^3/3/eb1^2+p7/3*(eb1-epsilon_b03^3/eb1^2)+p8/2*(1epsilon_b03^2/eb1^2))*x1^2+(b-c)*Rb3*(x03-x33)*(x03+x33)/2+sig*As*(h1-x1))/pe;%Заполнение массивов данныхeb1_m=[eb1_m,eb1];es_m=[es_m,es];x1_m=[x1_m,x1];x_m=[x_m,x];kriviznacr_m=[kriviznacr_m,kriviznacr];etap_m=[etap_m,etap];M_m=[M_m,M];endendEs=200000;%Проредим массивы (удаление начала этапа 4) для обеспечения плавного%деформирования арматуры:i=2;while i<=numel(es_m)if es_m(i)<es_m(i-1)es_m(i)=[];eb1_m(i)=[];x1_m(i)=[];x_m(i)=[];kriviznacr_m(i)=[];etap_m(i)=[];M_m(i)=[];elsei=i+1;endend%Массив момента для построения графика деформаций крайнего волокна%защитного слоя до его разрушения:krivizna_m2=[]; %Укороченный массив кривизны до начала этапа 4for i=1:numel(eb2_m)krivizna_m2=[krivizna_m2,kriviznacr_m(i)];end%2.2.3.Вычислим итоговых значений моментов%Вычислим максимальный момент до разрушения защитного слояMrz_m=[];%Часть массива момента для этапа 3for i=1:numel(etap_m)196if etap_m(i)<4Mrz_m=[Mrz_m,M_m(i)];endendMrz=real(max(Mrz_m))%Момент разрушения защитного слоя%Вычислим остаточную несущую способность:%Собираем массив моментов для этапов 4-5Most_m=[];for i=1:numel(M_m)if etap_m(i)>=4Most_m=[Most_m,M_m(i)];endend%Проверяем работает ли балка после разрушения защитного слоя (не пуст%ли массив Most):proverka=isempty(Most_m);if proverka==1Most=0;elseMost=max(Most_m);end%Вычислим разрушающий момент:Mult=max(M_m)%2.4.Запись в файлfilename='results.xlsx';zagolovki1={'Rb3=';'Mult=';'Mrz=';'Most=';'Mcrc=';'-';'Mn=';'psi=';'1/R='};zagolovki2={'Etap','M','eb1','eb2','es','x1','x','1/Rcr'};xlswrite(filename,zagolovki1,'Лист1','A2');xlswrite(filename,zagolovki2,'Лист1','D1');xlswrite(filename,Rb3,'Лист1','B2');xlswrite(filename,Mult,'Лист1','B3');xlswrite(filename,Mrz,'Лист1','B4');xlswrite(filename,Most,'Лист1','B5');xlswrite(filename,Mcrc,'Лист1','B6');xlswrite(filename,etap_m','Лист1','D2');xlswrite(filename,M_m','Лист1','E2');xlswrite(filename,eb1_m','Лист1','F2');xlswrite(filename,eb2_m','Лист1','G2');xlswrite(filename,es_m','Лист1','H2');xlswrite(filename,x1_m','Лист1','I2');xlswrite(filename,x_m','Лист1','J2');xlswrite(filename,kriviznacr_m','Лист1','K2');%3.РАСЧЕЧ СРЕДНЕЙ КРИВИЗНЫ НОРМАЛЬНОГО СЕЧЕНИЯ%3.1.Исходные данныеdisp('РАСЧЕЧ СРЕДНЕЙ КРИВИЗНЫ НОРМАЛЬНОГО СЕЧЕНИЯ:')disp('ИСХОДНЫЕ ДАННЫЕ:')197disp('Изгибающий момент, при котором требуется определить кривизну (кН*м):')Mn=input('Mn=');disp('%Коэффициент, отвечающий за профиль арматуры:')fi_sl=input('fi_sl=');%3.2.Вычисляем коэффициент psi%Определяем усилие в арматуре в трещине непосредственно после образованияfor i=1:numel(M_m)if M_m(i)>=Mcrcepsilon_scrc=(es_m(i)-es_m(i-1))/(M_m(i)-M_m(i-1))*(Mcrc-M_m(i-1))+es_m(i-1)breakendend%Определяем усилие в арматуре в трещине при заданном моментеfor i=1:numel(M_m)if M_m(i)>=Mnepsilon_sn=(es_m(i)-es_m(i-1))/(M_m(i)-M_m(i-1))*(Mn-M_m(i-1))+es_m(i-1)breakendend%Коэффициент psi:fi_m=epsilon_scrc/epsilon_sn;psi=1.25-1*fi_sl*fi_mif psi>1psi=1end%3.3.Задаем приведенный модуль упругости арматуры (актуален только для%заданного момента)Es=Es/psi;%3.4.Вычисляем массивы данных (фиктивные)%3.4.1.Задаем массивы для дальнейшего заполнения:%Вводим величины, которые будут пересчитываться в циклеeb1m=0; %Деформации крайнего сжатого волокнаeb2m=0; %Деформакии крайнего сжатого волокна армированного ядраesm=0;%Деформации арматурыxm=0;%Высота сжатой зоныx1m=0;%Высота сжатой зоны армированного ядраMm=0;%Момент, воспринимаемый сечением%Создаются пустые однострочные массивы для последующего заполнения%элементами:etapm_m=[];eb1m_m=[];eb2m_m=[];esm_m=[];xm_m=[];x1m_m=[];Mm_m=[];%3.4.2.Вычисление массивов данных198while Mm<=Mneb1m=eb1m+shag;if eb1m<epsilon_b0 & eb2m<=epsilon_b0etapm=1;%этап 1%Решенме уравненияq1=b*(p1/3*eb1m+p2/2)+(b-c)*(p3/5*eb1m^3+p4/4*eb1m^2+p5/3*eb1m+p6/2p1/3*eb1m-p2/2);q2=b*c3*(p1*eb1m+p2)+Es*As;q3=b*c3^2*(p1*eb1m+p2/2)-Es*As*h1;q4=b*p1/3*c3^3*eb1m;xxx=[q1 q2 q3 q4];korni=roots(xxx);korni=korni(abs(imag(korni))<1e-5);korni=real(korni);korni=korni(korni>0);x1m=min(korni);%Нахождение данных для текущего шагаxm=x1m+c3;esm=(h1-x1m)/x1m*eb1m;eb2m=xm/(h0-xm)*esm;sig=esm*Es;Mm=((b-c)*(p3/6*eb1m^4+p4/5*eb1m^3+p5/4*eb1m^2+p6/3*eb1m-p1/4*eb1m^2p2/3*eb1m)*x1m^2+b*(p1/4*eb2m^2+p2/3*eb2m)*xm^2+sig*As*(h0-xm))/pe;%Заполнение массивов данныхeb1m_m=[eb1m_m,eb1m];eb2m_m=[eb2m_m,eb2m];esm_m=[esm_m,esm];x1m_m=[x1m_m,x1m];xm_m=[xm_m,xm];etapm_m=[etapm_m,etapm];Mm_m=[Mm_m,Mm];elseif eb1m<=epsilon_b0 & eb2m>epsilon_b0etapm=2;%этап 2%Решенме уравненияq1=b*(p1*epsilon_b0^3/3/eb1m+p2*epsilon_b0^2/2/eb1m)+(bc)*(p3/5*eb1m^4+p4/4*eb1m^3+p5/3*eb1m^2+p6/2*eb1m-p1/3*eb1m^2-p2/2*eb1m)+b*Rb*(1epsilon_b0/eb1m);q2=b*c3*Rb+eb1m*Es*As;q3=-h1*eb1m*Es*As;xx=[q1 q2 q3];korni=roots(xx);korni=korni(abs(imag(korni))<1e-5);korni=real(korni);korni=korni(korni>0);x1m=min(korni);%Нахождение данных для текущего шага199xm=x1m+c3;esm=(h1-x1m)/x1m*eb1m;eb2m=xm/(h0-xm)*esm;sig=esm*Es;x0=epsilon_b0/eb1m*x1m;Mm=((b-c)*(p3/6*eb1m^4+p4/5*eb1m^3+p5/4*eb1m^2+p6/3*eb1m-p1/4*eb1m^2p2/3*eb1m)*x1m^2+b*(p1*epsilon_b0^4/4/eb1m^2+p2*epsilon_b0^3/3/eb1m^2)*x1m^2+b*Rb*(xm-x0)*(xm+x0)/2+sig*As*(h0-xm))/pe;%Заполнение массивов данныхeb1m_m=[eb1m_m,eb1m];eb2m_m=[eb2m_m,eb2m];esm_m=[esm_m,esm];x1m_m=[x1m_m,x1m];xm_m=[xm_m,xm];etapm_m=[etapm_m,etapm];Mm_m=[Mm_m,Mm];elseif eb1m>epsilon_b0 & eb2m<=epsilon_buetapm=3;%этап 3%Решенме уравненияq1=c*(p1*epsilon_b0^3/3/eb1m+p2*epsilon_b0^2/2/eb1m)+(bc)*(p3/5*eb1m^4+p4/4*eb1m^3+p5/3*eb1m^2+p6/2*eb1m)+Rb*c*(1-epsilon_b0/eb1m);q2=b*c3*Rb+eb1m*Es*As;q3=-h1*eb1m*Es*As;xx=[q1 q2 q3];korni=roots(xx);korni=korni(abs(imag(korni))<1e-5);korni=real(korni);korni=korni(korni>0);x1m=min(korni);%Нахождение данных для текущего шагаxm=x1m+c3;esm=(h1-x1m)/x1m*eb1m;eb2m=xm/(h0-xm)*esm;sig=esm*Es;x0=epsilon_b0/eb1m*x1m;Mm=((bc)*(p3/6*eb1m^4+p4/5*eb1m^3+p5/4*eb1m^2+p6/3*eb1m)*x1m^2+c*(p1*epsilon_b0^4/4/eb1m^2+p2*epsilon_b0^3/3/eb1m^2)*x1m^2+c*Rb*(x1mx0)*(x1m+x0)/2+b*c3*Rb*(x1m+c3/2)+sig*As*(h0-xm))/pe;%Заполнение массивов данныхeb1m_m=[eb1m_m,eb1m];eb2m_m=[eb2m_m,eb2m];esm_m=[esm_m,esm];x1m_m=[x1m_m,x1m];xm_m=[xm_m,xm];etapm_m=[etapm_m,etapm];Mm_m=[Mm_m,Mm];200endend%3.5.Вычисляем искомую осредненную кривизнуepsilon_s=(esm_m(end)-esm_m(end-1))/(Mm_m(end)-Mm_m(end-1))*(Mn-Mm_m(end1))+esm_m(end-1);epsilon_b1=(eb1m_m(end)-eb1m_m(end-1))/(Mm_m(end)-Mm_m(end-1))*(Mn-Mm_m(end1))+eb1m_m(end-1);epsilon_b2=(eb2m_m(end)-eb2m_m(end-1))/(Mm_m(end)-Mm_m(end-1))*(Mn-Mm_m(end1))+eb2m_m(end-1);krivizna=(epsilon_s+epsilon_b1)/h1epsilon_s=epsilon_s*100000epsilon_b2=epsilon_b2*100000%3.6.Запись в файлxlswrite(filename,Mn,'Лист1','B8');xlswrite(filename,psi,'Лист1','B9');xlswrite(filename,krivizna,'Лист1','B10');201Приложение Г202Приложение Г.