Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
БАКАЛАВР_Р.docx
Скачиваний:
8
Добавлен:
19.03.2015
Размер:
695.96 Кб
Скачать

Додаток

Числовий розв'язок задачі виконано за допомогою створеної програми в Maple 13:

> nx:= 0;

ny:=0;

nz:=1;

C11:=250;

C12:=160;

C44:=118.5;

a0:=3.5243;

L:=-0.073;

M:=58.7;

K:=(C11+2*C12)/3;

ksi:=(C11-C12-2*C44)/C44;

X:=nx^2*ny^2+nx^2*nz^2+ny^2*nz^2;

Y:=nx^2*ny^2*nz^2;

D=C11+ksi*(C11+C12)*X+ksi^2*(C11+2*C12+C44)*Y;

A:=(-9*K^2*a0^3*L^2*(1+2*ksi*X+3*ksi^2*Y))/(4*250);

omega:=53.72;

B:=(9*K^2*a0^3*L^2*M*omega^2)/(256*4*C11^2);

nnx:=0.9035;

nny:=0.4287;

nnz:=0;

XX:=nnx^2*nny^2+nnx^2*nnz^2+nny^2*nnz^2;

YY:=nnx^2*nny^2*nnz^2;

DD:=C11+ksi*(C11+C12)*XX+ksi^2*(C11+2*C12+C44)*YY;

Q:=9*K^2*a0^3*L^2*(1+2*ksi*XX+3*ksi^2*YY)/(4*DD);

r:=a0/sqrt(2);

S4:=25.338;

S8:=12.80193;

epsilon:=-4.45;

a:=(-epsilon*a0^4)/S4;

b:=(-epsilon*a0^8)/(8*S8);

Phi:=-a*r^(-4)+b*r^(-8);

k:=0.0349;

V:=A+B*k+Q;

W:=V+Phi;

31