Кучланиш ва дифференциаллар

 

---

27 
3.3-masalani yechish dasturi 
 
Программа 2. 
(P 3.3.mws) 
> restart;# Masala3.3 
> with(LinearAlgebra): 
Puasson koeffitsienti, Yung koeffitsienti (MPa) 
> nu:=0.4: E:=1e5: 
Berilgan kuchlanishlar sigma11,tau12,tau13,sigma22,tau23,sigma33 
> S:=Matrix([[11.,12,14],[13,1],[-33]], 
> shape=symmetric,scan=triangular [upper]); 












33
1
14
1
13
12
14
12
11
:
s
Bosh kuchlanishlar 
> S0:=Eigenvalues(S); 
S0 :=










690689272
25.9266482
323587076
2.27082811
823047852
37.1974763
-
> E1:=Vector(3,[1,1,1]);# Yordamchi vektor 
E1 :=










1
1
1
Bosh urinma kuchlanishlar 
> Tau:=CrossProduct(E1,S0)/2: Tau[2]:=-Tau[2]: Tau; 










450000008
19.7341522
249999994
31.5620623
800000006
11.8279100
> 2/3*Norm(Tau,2);# Oktaedrik kuchlanish 
26.03843313 
> s:=sort([S0[1],S0[2],S0[3]],‘>‘): 
> with(ListTools): 
> for i to 3 do 
> eps[i]:=(s[1]-(s[2]+s[3])*nu)/E: 

---

28 
> s:=Rotate(s,1);# Doiraviy o‘rin almashtirish 
> od: 
> evalf(seq(eps[i],i=1..3),4); 
> Hajmning nisbiy o‘zgarishi =add(eps[i],i=1..3); 
Hisob natijalari: 
0.0003990, 0.00006779, −0.0004848 
Hajmning nisbiy o‘zgarishi = −0.0000179984000