Файловый архив студентов. Файловый архив студентов.
1261 вуз, 4614 предмет.
/ Регистрация
FAQ Обратная связь Вопросы и предложения
Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Московский государственный технический университет им. H.Э.Баумана
Предмет:
[НЕСОРТИРОВАННОЕ]
Файл:
дз Зарубин.docx
Скачиваний:
23
Добавлен:
09.02.2015
Размер:
485.37 Кб
Скачать
☆
►Содержание►

7. Заключение

В данной работе рассматривались кристаллы с кубической и гексагональной решетками. По матрице коэффициентов упругости были найдены матрицы коэффициентов податливости, построены графики зависимостей линейной податливости от направлений единичного вектора для различных случаев расположения относительно плоскостей кристалла.

Для поликристалла найдены оценки механических характеристик в предположениях хаотической ориентации зерен и статистически усредненной сферической формы кристаллических зерен.

Аналогичные расчеты были проведены и для случаев пористого сплава-смеси из материалов с кубической и гексагональной кристаллическими решетками. Построены графики зависимостей характеристик смеси от объемных долей металлов в сплаве при разных значениях пористости. Показана некорректность нижних оценок при наличии пор.

8. Список литературы

  1. Зарубин В.С., Кувыркин Г.Н. Математические модели механики и электродинамики сплошной среды. – М.: Изд-во МГТУ им. Н.Э. Баумана, 2008. – 512 с.: ил.(Математическое моделирование в технике и в технологии).

  2. Зарубин В.С. Прикладные задачи термопрочности элементов конструкций. М.: Машиностроение, 1985. 296 с.

9. Приложение

n=6;

CCo=Table[If[(ij)&&(i n/2),c11Co,If[(j n/2-1)&&(i n/2-1),c12Co,If[(j n/2)&&(i n/2),c13Co,If[(ij)&&(i> n/2),c44Co,0]]]],{i,1,n},{j,1,n}];

CCo[[3,3]]=c33Co;

CCo[[6,6]]=1/2*(c11Co-c12Co);

c11Co=303.5;

c12Co=161.2;

c13Co=99.8;

c33Co=354.2;

c44Co=71.9;

MatrixForm[CCo]

SCo=Table[If[(ij)&&(i n/2),s11Co,If[(j n/2-1)&&(i n/2-1),s12Co,If[(j n/2)&&(i n/2),s13Co,If[(ij)&&(i> n/2),s44Co,0]]]],{i,1,n},{j,1,n}];

s11Co=(-c13Co c13Co+c11Co c33Co)/((c11Co-c12Co)*(-2 c13Co c13Co+(c11Co+c12Co)*c33Co));

s12Co=(c13Co c13Co-c33Co c12Co)/((c11Co-c12Co)*(-2 c13Co c13Co+(c11Co+c12Co)*c33Co));

s13Co=-c13Co/c33Co (2 s12Co+1/(c11Co-c12Co));

s33Co=(-s13Co (c11Co+c12Co))/c13Co;

SCo[[3,3]]=s33Co;

s44Co=1/c44Co;

SCo[[6,6]]=2*(s11Co-s12Co);

MatrixForm[SCo.CCo]

e3[_]:=Sin[];

SECo[_]:=s11Co (1-e3[]^2)2+s33Co e3[]4+(2 s13Co+s44Co)e3[]2 (1-e3[]2)

Plot[SECo[],{,0,6},PlotRangeAll,PlotStyle{Thick,Black},AxesStyle->{Directive[Black,13],Directive[Black,13]},AxesLabel{Style[", рад"],Style["n, ГПа-1"]}, AxesOrigin{0,0}]

cmmppCo=(2c11Co+c33Co)+2(2c13Co+c12Co) (* machine precision *)

smmppCo=(2s11Co+s33Co)+2(2s13Co+s12Co)

cmpmpCo=(2c11Co+c33Co)+2(2c44Co+CCo[[6,6]])

smpmpCo=(2s11Co+s33Co)+1/2 (2s44Co+SCo[[6,6]])

FCo=1/15 (2 cmmppCo-cmpmpCo);

FCo=1/30 (3 cmpmpCo-cmmppCo);

sRCo=1/30 (3 smpmpCo-smmppCo);

RCo=1/(4sRCo);

GRCo=1/smmppCo

GFCo=FCo+2/3 FCo

RCo=15/(2*(3*smpmpCo-smmppCo))

FCo

RCo=(3*GRCo-2*RCo)/(6*GRCo+2RCo)

FCo=(3*GFCo-2*FCo)/(6*GFCo+2FCo)

ERCo=2 RCo (1+RCo)

EFCo=2 FCo (1+FCo)

CLi=Table[If[(ij)&&(i n/2),c11Li,If[(j n/2)&&(i n/2),c12Li,If[(ij)&&(i> n/2),c44Li,0]] ],{i,1,n},{j,1,n}];

c11Li=13.5;

c12Li=11.4;

c44Li=9.6;

SLi=Table[If[(ij)&&(i n/2),s11Li,If[(j n/2)&&(i n/2),s12Li,If[(ij)&&(i> n/2),s44Li,0]] ],{i,1,n},{j,1,n}];

s11Li=(c11Li+c12Li)/((c11Li-c12Li)*(c11Li+2*c12Li));

s12Li=-c12Li/(c11Li2+c11Li*c12Li-2*c12Li2);

s44Li=1/c44Li;

MatrixForm[SLi]

s44Li-2*(s11Li-s12Li

MatrixForm[SLi.CLi]

e1[_]:=Cos[];

e2[_]:=Sin[];

e3[_]:=0;

SELi[_]:=s11Li-(2 s11Li-2 s12Li-s44Li)*(e1[]2 e2[]2+e1[]2 e3[]2+e2[]2 e3[]2)

Plot[SELi[],{,0,6},PlotRange{0,0.4},PlotStyle{Thick,Black},AxesStyle->{Directive[Black,13],Directive[Black,13]},AxesLabel{Style[", рад"],Style["n, ГПа-1"]}, AxesOrigin{0,0}]

e1[_]:=1/Sqrt[2] Cos[];

e2[_]:=1/ Sqrt[2] Cos[];

e3[_]:=Sin[];

SELi[_]:=s11Li-(2 s11Li-2 s12Li-s44Li)*(e1[]2 e2[]2+e1[]2 e3[]2+e2[]2 e3[]2)

Plot[SELi[],{,0,6},PlotRangeAll,PlotStyle{Thick, Black},AxesStyle->{Directive[Black,13],Directive[Black,13]},AxesLabel{Style[", рад"],Style["n, ГПа-1"]}, AxesOrigin{0,0}]

[_]:=/6-;

e1[_]:=1/ Sqrt[6] Cos[[]]+1/Sqrt[2]Sin[[]];

e2[_]:=1/ Sqrt[6] Cos[[]]-1/Sqrt[2]Sin[[]];

e3[_]:=Sqrt[2/3] Cos[[]];

SELi[_]:=s11Li-(2 s11Li-2 s12Li-s44Li)*(e1[]2 e2[]2+e1[]2 e3[]2+e2[]2 e3[]2)

Plot[SELi[],{,0,6},PlotRangeAll,PlotStyle{Thick, Black},AxesStyle->{Directive[Black,13],Directive[Black,13]},AxesLabel{Style[", рад"],Style["n, ГПа-1"]}, AxesOrigin{0,0}]

cmmppLi=3(c11Li+2c12Li)

smmppLi=3(s11Li+2s12Li)

cmpmpLi=3(c11Li+2c44Li)

smpmpLi=3(s11Li+1/2 s44Li)

FLi=1/15 (2 cmmppLi-cmpmpLi);

FLi=1/30 (3 cmpmpLi-cmmppLi);

sRLi=1/30 (3 smpmpLi-smmppLi);

RLi=1/(4sRLi)

GRLi=1/smmppLi

GFLi=FLi+2/3 FLi

RLi

RLi=15/(2*(3*smpmpLi-smmppLi))

FLi

RLi=(3*GRLi-2*RLi)/(6*GRLi+2RLi)

FLi=(3*GFLi-2*FLi)/(6*GFLi+2FLi)

ERLi=2 RLi(1+RLi)

EFLi=2 FLi (1+FLi)

CP=Table[0.0,{i,1,n},{j,1,n}];

c011[G_,_]:=G+4/3 ;

c012[G_,_]:=G-2/3 ;

c044[G_,_]:=;

C0=Table[If[(ij)&&(i n/2),c011[G,],If[(j n/2)&&(i n/2),c012[G,],If[(ij)&&(i> n/2),c044[G,],0]] ],{i,1,n},{j,1,n}];

MatrixForm[C0]

[G_,_]:=(3*G-2*)/(6*G+2)

I1=Table[If[(i 3)&& (j 3),1,0],{i,6},{j,6}];

I2=1/2*Table[If[(i 3)&& (j 3)&& (ji),2,If[ (ji),1,0]],{i,6},{j,6}];

w=Simplify[3/2*(1-[G,])/(4-5[G,]) ((3*(4-5[G,]))/(1+[G,])*I1+5/2*Inverse[I2])];

[CC_]:=Inverse[CC-C0+C0.w].(C0-CC);

kkmm[_]:=Sum[Sum[[[i,j]],{i,1,3}],{j,1,3}];

kmkm[_]:=Sum[[[i,i]],{i,1,3}]+2*Sum[[[i,i]],{i,4,6}];

{RLi,GRLi}

{,G}/.FindMinimum[kmkm[[CLi]]^2+kkmm[[CLi]]^2,{{,FLi},{G,GFLi}}][[2]]

{FLi,GFLi}

esLi={,G}[[1]]/.FindMinimum[kmkm[[CLi]]^2+kkmm[[CLi]]^2,{{,FLi},{G,GFLi}}][[2]]

GesLi={,G}[[2]]/.FindMinimum[kmkm[[CLi]]^2+kkmm[[CLi]]^2,{{,FLi},{G,GFLi}}][[2]]

esLi=(3*GesLi-2*esLi)/(6*GesLi+2*esLi)

EesLi=2 esLi(1+esLi)

{RCo,GRCo}

{,G}/.FindMinimum[{kmkm[[CCo]]2+kmkm[[CCo]]2},{{,FCo},{G,GFCo}}][[2]]

{FCo,GFCo}

esCo={,G}[[1]]/.FindMinimum[kmkm[[CCo]]^2+kkmm[[CCo]]^2,{{,FCo},{G,GFCo}}][[2]]

GesCo={,G}[[2]]/.FindMinimum[kmkm[[CCo]]^2+kkmm[[CCo]]^2,{{,FCo},{G,GFCo}}][[2]]

esCo=(3*GesCo-2*esCo)/(6*GesCo+2*esCo)

EesCo=2 esCo(1+esCo)

V3=0.0;

0=Table[{t-1,{,G}[[1]]/.FindMinimum[{(a*kmkm[[CLi]]+(1-a-V3)*kmkm[[CCo]]+V3*kmkm[[CP]])^2+(a*kkmm[[CLi]]+(1-a-V3)*kkmm[[CCo]]+V3*kkmm[[CP]])^2}/.a0.1*(t-1)*(1-V3) ,{{,oc /.{oca FLi+(1-a-V3)*FCo}/.a0.1*(t-1)*(1-V3)},{G,ocG /.{ocGa GFLi+(1-a-V3)*GFCo}/.a0.1*(t-1)*(1-V3) }},MaxIterations50] [[2]]},{t,1,11}];

m01=ListPlot[0,PlotRangeAll,PlotStyle{Black,PointSize[0.02]},AxesOrigin{0,0}];

m02=Plot[{a FLi+(1-a-V3)*FCo}/.a0.1*t(1-V3),{t,0,10},PlotRangeAll,PlotStyle{Black,Thick},AxesOrigin{0,0}];

m03=Plot[{1/(a/RLi +((1-a-V3)/RCo))}/.a0.1*t(1-V3),{t,0,10},PlotRangeAll,PlotStyle{Dashed,Black,Thick},AxesOrigin{0,0}];

Show[m02,m03,m01,PlotRangeAll,AxesOrigin {0,0},AxesStyle->{Directive[Black,13],Directive[Black,13]},AxesLabel{ StandardForm["V1/(V1+V2)"],", ГПа"}]

G0=Table[{t-1,{,G}[[2]]/.FindMinimum[{(a*kmkm[[CLi]]+(1-a-V3)*kmkm[[CCo]]+V3*kmkm[[CP]])^2+(a*kkmm[[CLi]]+(1-a-V3)*kkmm[[CCo]]+V3*kkmm[[CP]])^2}/.a0.1*(t-1)*(1-V3) ,{{,oc /.{oca FLi+(1-a-V3)*FCo}/.a0.1*(t-1)*(1-V3)},{G,ocG /.{ocGa GFLi+(1-a-V3)*GFCo}/.a0.1*(t-1)*(1-V3) }},MaxIterations50] [[2]]},{t,1,11}];

s01=ListPlot[G0,PlotRangeAll,PlotStyle{Black,PointSize[0.02]},AxesOrigin{0,0}];

s02=Plot[{a GFLi+(1-a-V3)*GFCo}/.a0.1*t*(1-V3),{t,0,10},PlotRangeAll,PlotStyle{Black,Thick},AxesOrigin{0,0}];

s03=Plot[{1/(a/GRLi +((1-a-V3)/GRCo))}/.a0.1*t*(1-V3),{t,0,10},PlotRangeAll,PlotStyle{Dashed,Black,Thick},AxesOrigin{0,0}];

Show[s01,s02,s03,PlotRangeAll,AxesOrigin {0,0},AxesStyle->{Directive[Black,13],Directive[Black,13]},AxesLabel{ StandardForm["V1/(V1+V2)"],"k, ГПа"}]

0=Table[{t-1,(3*G0[[t,2]]-2*0[[t,2]])/(6*G0[[t,2]]+2*0[[t,2]])},{t,1,11}];

01=ListPlot[0,PlotRangeAll,PlotStyle{Black,PointSize[0.02]},AxesOrigin{0,0}];

02=Plot[{a FLi+(1-a-V3)*FCo}/.a0.1*t(1-V3),{t,0,10},PlotRangeAll,PlotStyle{Black,Thick},AxesOrigin{0,0}];

03=Plot[{1/(a/RLi +((1-a-V3)/RCo))}/.a0.1*t(1-V3),{t,0,10},PlotRangeAll,PlotStyle{Dashed,Black,Thick},AxesOrigin{0,0}];

Show[01,02,03,PlotRangeAll,AxesOrigin {0,0},AxesStyle->{Directive[Black,13],Directive[Black,13]},AxesLabel{ StandardForm["V1/(V1+V2)"], ""}]

E0=Table[{t-1,2 0[[t,2]](1+0[[t,2]])},{t,1,11}];

E01=ListPlot[E0,PlotRangeAll,PlotStyle{Black,PointSize[0.02]},AxesOrigin{0,0}];

E02=Plot[{a EFLi+(1-a-V3)*EFCo}/.a0.1*t(1-V3),{t,0,10},PlotRangeAll,PlotStyle{Black,Thick},AxesOrigin{0,0}];

E03=Plot[{1/(a/ERLi +((1-a-V3)/ERCo))}/.a0.1*t(1-V3),{t,0,10},PlotRangeAll,PlotStyle{Dashed,Black,Thick},AxesOrigin{0,0}];

Show[E01,E02,E03,PlotRangeAll,AxesOrigin {0,0},AxesStyle->{Directive[Black,13],Directive[Black,13]},AxesLabel{ StandardForm["V1/(V1+V2)"], "E, ГПа"}]

V3=0.1;

1=Table[{t-1,{,G}[[1]]/.FindMinimum[{(a*kmkm[[CLi]]+(1-a-V3)*kmkm[[CCo]]+V3*kmkm[[CP]])^2+(a*kkmm[[CLi]]+(1-a-V3)*kkmm[[CCo]]+V3*kkmm[[CP]])^2}/.a0.1*(t-1)*(1-V3) ,{{,oc /.{oca FLi+(1-a-V3)*FCo}/.a0.1*(t-1)*(1-V3)},{G,ocG /.{ocGa GFLi+(1-a-V3)*GFCo}/.a0.1*(t-1)*(1-V3) }},MaxIterations50] [[2]]},{t,1,11}];

m11=ListPlot[1,PlotRangeAll,PlotStyle{Black,PointSize[0.02]},AxesOrigin{0,0}];

m12=Plot[{a FLi+(1-a-V3)*FCo}/.a0.1*t*(1-V3),{t,0,10},PlotRangeAll,PlotStyle{Black,Thick},AxesOrigin{0,0}];

(*m13=Plot[{1/(a/RLi +((1-a-V3)/RCo))}/.a0.1*t(1-V3),{t,0,10},PlotRangeAll,PlotStyle{Green,Thick},AxesOrigin{0,0}];*)

Show[m11,m12,PlotRangeAll,AxesOrigin {0,0},AxesStyle->{Directive[Black,13],Directive[Black,13]},AxesOrigin{0,0},AxesLabel{ StandardForm["V1/(V1+V2)"],",ГПа"}]

G1=Table[{t-1,{,G}[[2]]/.FindMinimum[{(a*kmkm[[CLi]]+(1-a-V3)*kmkm[[CCo]]+V3*kmkm[[CP]])^2+(a*kkmm[[CLi]]+(1-a-V3)*kkmm[[CCo]]+V3*kkmm[[CP]])^2}/.a0.1*(t-1)*(1-V3) ,{{,oc /.{oca FLi+(1-a-V3)*FCo}/.a0.1*(t-1)*(1-V3)},{G,ocG /.{ocGa GFLi+(1-a-V3)*GFCo}/.a0.1*(t-1)*(1-V3) }},MaxIterations50] [[2]]},{t,1,11}];

s11=ListPlot[G1,PlotRangeAll,PlotStyle{Black,PointSize[0.02]},AxesOrigin{0,0}];

s12=Plot[{a GFLi+(1-a-V3)*GFCo}/.a0.1*t(1-V3),{t,0,10},PlotRangeAll,PlotStyle{Black,Thick},AxesOrigin{0,0}];

(*s13=Plot[{1/(a/GRLi +((1-a-V3)/GRCo))}/.a0.1*t(1-V3),{t,0,10},PlotRangeAll,PlotStyle{Green,Thick},AxesOrigin{0,0}];*)

Show[s11,s12,PlotRangeAll,AxesOrigin {0,0},AxesStyle->{Directive[Black,13],Directive[Black,13]},AxesLabel{ StandardForm["V1/(V1+V2)"],"k, ГПа"}]

1=Table[{t-1,(3*G1[[t,2]]-2*1[[t,2]])/(6*G1[[t,2]]+2*1[[t,2]])},{t,1,11}];

11=ListPlot[1,PlotRangeAll,PlotStyle{Black,PointSize[0.02]},AxesOrigin {0,0}];

12=Plot[{a FLi+(1-a-V3)*FCo}/.a0.1*t*(1-V3),{t,0,10},PlotRangeAll,PlotStyle{Black,Thick},AxesOrigin {0,0}];

(*13=Plot[{1/(a/RLi +((1-a-V3)/RCo))}/.a0.1*t(1-V3),{t,0,10},PlotRangeAll,PlotStyle{Green,Thick},AxesOrigin {0,0}];*)

Show[11,12,PlotRangeAll,AxesOrigin {0,0},AxesStyle->{Directive[Black,13],Directive[Black,13]},AxesLabel{ StandardForm["V1/(V1+V2)"],""}]

E1=Table[{t-1,2 1[[t,2]](1+1[[t,2]])},{t,1,11}];

E11=ListPlot[E1,PlotRangeAll,PlotStyle{Black,PointSize[0.02]},AxesOrigin {0,0}];

E12=Plot[{a EFLi+(1-a-V3)*EFCo}/.a0.1*t*(1-V3),{t,0,10},PlotRangeAll,PlotStyle{Black,Thick},AxesOrigin {0,0}];

(*E13=Plot[{1/(a/ERLi +((1-a-V3)/ERCo))}/.a0.1*t(1-V3),{t,0,10},PlotRangeAll,PlotStyle{Green,Thick},AxesOrigin {0,0}];*)

Show[E11,E12,PlotRangeAll,AxesOrigin {0,0},AxesStyle->{Directive[Black,13],Directive[Black,13]},AxesLabel{ StandardForm["V1/(V1+V2)"], "E, ГПа"}]

V3=0.2;

2=Table[{t-1,{,G}[[1]]/.FindMinimum[{(a*kmkm[[CLi]]+(1-a-V3)*kmkm[[CCo]]+V3*kmkm[[CP]])^2+(a*kkmm[[CLi]]+(1-a-V3)*kkmm[[CCo]]+V3*kkmm[[CP]])^2}/.a0.1*(t-1)*(1-V3) ,{{,oc /.{oca FLi+(1-a-V3)*FCo}/.a0.1*(t-1)*(1-V3)},{G,ocG /.{ocGa GFLi+(1-a-V3)*GFCo}/.a0.1*(t-1)*(1-V3) }},MaxIterations50] [[2]]},{t,1,11}];

m21=ListPlot[2,PlotRangeAll,PlotStyle{Black,PointSize[0.02]},AxesOrigin {0,0}];

m22=Plot[{a FLi+(1-a-V3)*FCo}/.a0.1*t*(1-V3),{t,0,10},PlotRangeAll,PlotStyle{Black,Thick},AxesOrigin {0,0}];

(*m23=Plot[{1/(a/RLi +((1-a-V3)/RCo))}/.a0.1*t*(1-V3),{t,0,10},PlotRangeAll,PlotStyle{Green,Thick},AxesOrigin {0,0}];*)

Show[m21,m22,PlotRangeAll,AxesOrigin {0,0},AxesStyle->{Directive[Black,13],Directive[Black,13]},AxesLabel{ StandardForm["V1/(V1+V2)"],", ГПа"}]

G2=Table[{t-1,{,G}[[2]]/.FindMinimum[{(a*kmkm[[CLi]]+(1-a-V3)*kmkm[[CCo]]+V3*kmkm[[CP]])^2+(a*kkmm[[CLi]]+(1-a-V3)*kkmm[[CCo]]+V3*kkmm[[CP]])^2}/.a0.1*(t-1)*(1-V3) ,{{,oc /.{oca FLi+(1-a-V3)*FCo}/.a0.1*(t-1)*(1-V3)},{G,ocG /.{ocGa GFLi+(1-a-V3)*GFCo}/.a0.1*(t-1)*(1-V3) }},MaxIterations50] [[2]]},{t,1,11}];

s21=ListPlot[G2,PlotRangeAll,PlotStyle{Black,PointSize[0.02]},AxesOrigin {0,0}];

s22=Plot[{a GFLi+(1-a-V3)*GFCo}/.a0.1*t(1-V3),{t,0,10},PlotRangeAll,PlotStyle{Black,Thick},AxesOrigin {0,0}];

(*s23=Plot[{1/(a/GRLi +((1-a-V3)/GRCo))}/.a0.1*t*(1-V3),{t,0,10},PlotRangeAll,PlotStyle{Green,Thick},AxesOrigin {0,0}];*)

Show[s21,s22,PlotRangeAll,AxesOrigin {0,0},AxesStyle->{Directive[Black,13],Directive[Black,13]},AxesLabel{ StandardForm["V1/(V1+V2)"],"k, ГПа"}]

2=Table[{t-1,(3*G2[[t,2]]-2*2[[t,2]])/(6*G2[[t,2]]+22[[t,2]])},{t,1,11}];

21=ListPlot[2,PlotRangeAll,PlotStyle{Black,PointSize[0.02]},AxesOrigin {0,0}];

22=Plot[{a FLi+(1-a-V3)*FCo}/.a0.1*t*(1-V3),{t,0,10},PlotRangeAll,PlotStyle{Black,Thick},AxesOrigin {0,0}];

(*23=Plot[{1/(a/RLi +((1-a-V3)/RCo))}/.a0.1*t*(1-V3),{t,0,10},PlotRangeAll,PlotStyle{Green,Thick},AxesOrigin {0,0}];*)

Show[21,22,PlotRangeAll,AxesOrigin {0,0},AxesStyle->{Directive[Black,13],Directive[Black,13]},AxesLabel{ StandardForm["V1/(V1+V2)"],""}]

E2=Table[{t-1,2 2[[t,2]](1+2[[t,2]])},{t,1,11}];

E21=ListPlot[E2,PlotRangeAll,PlotStyle{Black,PointSize[0.02]},AxesOrigin {0,0}];

E22=Plot[{a EFLi+(1-a-V3)*EFCo}/.a0.1*t*(1-V3),{t,0,10},PlotRangeAll,PlotStyle{Black,Thick},AxesOrigin {0,0}];

(*E23=Plot[{1/(a/ERLi +((1-a-V3)/ERCo))}/.a0.1*t*(1-V3),{t,0,10},PlotRangeAll,PlotStyle{Green,Thick},AxesOrigin {0,0}];*)

Show[E21,E22,PlotRangeAll,AxesOrigin {0,0},AxesStyle->{Directive[Black,13],Directive[Black,13]},AxesLabel{ StandardForm["V1/(V1+V2)"],"E, ГПа"}]

S0=FullSimplify[Inverse[C0]]

Sum[S0[[i,i]],{i,1,3}]+(1/2)*Sum[S0[[i,i]],{i,4,6}]//FullSimplify

Sum[S0[[i,i]],{i,1,3}]+2*(S0[[2,3]]+S0[[1,3]]+S0[[1,2]])//FullSimplify

a=2*(c11Co+c12Co+2*c13Co)+c33Co

b=c11Co+c12Co+2*c33Co-4*c13Co

c=(c11Co+c12Co)*c33Co-2*c13Co2

b[_]:= /6*(9*GRCo+8*)/( GRCo+2*)

q0[_]:= -8*b*2- 9*c*

q1[_]:= -48*2-6**(a-2 b)+6*c

q2[_]:= 4*(18*+a)

FindRoot[q2[]*b[]2+q1[]*b[]+q0[]0,{,FCo}]

29

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]
Помощь Обратная связь Вопросы и предложения Пользовательское соглашение Политика конфиденциальности