Name' Mnesh l'anda AEP- Amigpmen-1

Rou Ale!r SE23MCAS O10

q1. SS2
6ii:[
Sa S22
LS31 S32 S33

Q. Jo datermine tu trta uien

Vecton aleng aplane dekined by a Ngxmal
We can use fomla fon tration
bilb.
VI t8n +e3

-|5 3 -17

Unit norumal vett o N

cal utak]
2 - 3
3 -l9

22.1 + I2:0-I6.

-32
-|6.| +3:0-14|

So, the tnauien vten Olong the plane 1: a32

Vikual SKeh! Ohe tomnonutA N whih has (omponenth n he

and z dineetion nd 0 en y din?. Jh means lane aiented to

diogonay in nz plane.
8 Add a legend for the traction vector
legend ('Plane', Traction Vector');
8 View the plot from a good angle
view (45, 30) :
t Hold off to end plotting
hold off;

Visualization of Plane and Traction Vector

Plane
Tracion Vector

10

N0

-10
-10 10

X 10 -10

% Q-1(B) Evaluate eigen values and eigen vectors

clc
clear all

% Define the stress tensor

sigma = (22, 12, -15; 12, -7, 3; -15, 3, -17];
% Calculate eigenvalues and eigenvectors
[eigenvecs, eigenvals] = eig (sigma)
$ Display eigenvalues and eigenvectors
disp(Eigenvalues : ');
disp (eigenvals);
disp ('Eigenvectors: ') ;
disp (eigenvecs) ;

2
 Skereh ot p1one

Shadel anea ’ Jaatien

plane
’ ’ Jnation veten

-10

Jrauien Plane (MATLAS Fig )

b. rtan Cay
Tay Syg Cyz
-\5 3

chanutisties ogn
Stp-1 4inding he co-ticlut at

-? 3 22 -l6

3 -14

= 5363

Stp-2 6 262 + 187 6-5363 =0

6,: 20.085 (Maxm)

625-94
24-936

Sip-3 Jo Eind the prinunal plane eorneS ponduny Jo tu 4 n principie Aret

Eigenvalues :
-24.9359 0
-7.1488
30.0848

Eigenvectors :
-0.3691 -0.1257 -0.9208
0.3882 0.8794 -0.2757
-0.8444 0.4592 0.2758

8 Q-8(a) Use MATLAB directly to calculate c Ax and D AB.

8 Define tensors A, B, and x
A= [1, 2, 1; 1, 4, 3; 2, 1, 2] ;
B [1, 2, 1; 4, 1, 2; 3, 2, 1];
x=(2; 4; 1];
8 Calculate Ax
C A * x;

8 CalculateD AB
D=A * B;

4 Display results
disp ('c - Ax:'):
disp (c) ;
disp ('D AB:'):
disp (D) ;
CAx:
11
21
10

D AB:
12 6 6
26 12 12
12 9 6

Q-8(B)
& Define the tensors A and B
A = (1 2 1; 1 4 3; 2 1 2]:
B (1 2 1; 4 1 2; 3 2 1]1
8 Initialize the result tensor D
D zeros (size (A) ) ;

8 Perform the calculations using index notation

for 1 1:size(A, 1)
for k = 1:size (B, 2)
for j 1:size (A, 2)
D(i, k)D(i, k) + A(i, j) * B(, k) ;

3
 -8.086 12 -I5
a2-90.0&5 I2 -I6
12 -23-085 3
I2 -9-30. 0R6
-I6 3 13. 086L
- 15 -19 -30, 0R5 -0
3

-8-086 O t122-6n3 0
12n -23. O&5 n) 43n, 0 0
-I5n, t 3n, -13.0hna rt-)

-o.421
-0.246
0246

2 liven, vCR?y)
-M9,

Stness 4ensan Cy

24ulti hrium in y dineuien (tye)-o

equiibrium in z- dineution dz +d

We car npres divacive7 clt ,4z as

d -M
dy dz
 - 0
Ty Tyz: V(R'-y)

Cye -Tey 0
Uhe sAA) Jenso matrin Can be Writ 0 a.
6 vCR2. y2)
3I
V(R2y?)
3T (Symmutrit manin)

(onsiktunt with the Anuss bld provited for a berdny problum

n beam.

We can @rite en matria don,

-MII

zero be UUse th

Bimilarly,
b. 2quisubrium in y dmuui on 0 ooj -0

C. equilei briun n z-dinetion =

-ye
Jhe stnets enhe, givn satiAFiak eqculuibium quatian în au thre
dlneuieng,. Equitibricun ss denes in a- din= Sn : - My
Memnt in 4- dian
dM
dy tV=0
Ahu cenbirumthat he provde Arem tild r conwsiunt wth the Llemensny
Anength s& matUals theery Son a beam tder shean & bundin 9.
 I

23
3. a. CijK éxij 6
e i =S I 5 i ) u an ven pamutation (i,2,)
Ohe penmt0tion bymbo l )-i5 cijhx) an odd pormaukaion (12,3)
0, 6theruie (9% thene hae nepeatel indite)
permtazian wiL he
Jhene te -nekclth , t Shoud be 1y-1 on 0, then
9! 312y(- 6.
rpratd then ts zeno, ko we haw 10 ondey any 0on
% ony indites
repeat ed idices,

D ohen i,i, ) =(12,3) =1

have tijK Euije (-1) (-)
i) when (ii,) - (,3,2) , we

i) when (bi)-(,,3) hoe

e3.1) , we hawe eijK Ekij ci).(-1) = 4
i) when(iyix)
whnix)= (3,1,2), we hawe eijK éRij 1. 11
v)
(3,2,1), we have eiK Erij =(-1).(-1):1
vi) whnLi,)
(Ans),ds
eiiu euij = 1+1t1+1t1+4 -6

) when (i, i, <)-(1,2,3)

ea13 a,az -a, a,
i) ohen('5)=(9,1,3) then eiin ajak =

thy tiyk aj au= E212 a,a,' = lxajagaa

) awhen (ii,)= 9 , 2)
tijaj ax '= eg| 4,a, -aga,
vi)when j , ) = 6,2,1) heneqcatian.
samminy cp au sth
(ajaje aj ai)
NoW

(Ana )
4. when ie I, je1 &y Sue4:4- 1
9. when is 2, je2 Sa2632 It
1le 1
3 when is3, Je 3 &a3 Sa=
6cumming cuptu rdivdeual evaleialiary,

b. Sj SinSaj
when is K,SiK=4

KAjs&j eo
Simitary don k-j , &xj 1, and fon

&j
No@ we (onaidn tu knonKan deta
ohen is4 )Sije4

Jhenfone the snpnmion surmpubis further

when ie jek
Sij &j= Sitkii: 141

Summeny u 1t0+0 1.
d. &ij Aij
A A13
A2| A22 A3
S3 A32 A3
0

e. Ain &ij
Ain &ij: Aix * &ij

ixK AiK &i AlK: 0 O

Summiy up Aik+0= AiR

5. ax (byc) alt:b) - alb-c)
as (bre)
c,4
aibjCk xCe3xet) )
aibjex (i x eixér )

we Know that,
Exij frem = Sin Sim-SjuSim
.! ax Cbt) aib; a (SiL Sjm - &jm Si) Sj Sim - SjmSa
zai bjlkSit Sim -
aibjuSim Sit
 ) aib;c ( Sie Sim - &;L Sim)

obi
: @c)b- ab).c (Pevd )
6. Thotnepie Jnnen!- TAeBMenie JunGn in desined Us a Jenhen nocamng
rettlen G6 co-6ndination
tormpanznth that ane cnthangd by anbiang)
T [T]
Proot Jhe eyuatian 86 transtomatian Boa dind orden tenso Alj n

Gn

Sine ee have vauuas Gon

e132 E32| t213 = -|

A on Gu Qm Gng
 bfew aeng
1. 96 pauts haksimes ,mán equals thun desenminant valu egual t0
Zeno.
a. 9% the vals 85 lm &n ane
yeie peretaeien en the vatcue 85 detenminant
3. 96 the values s i,m an ane anti - tyelie pernutatians hen he vatue
& detuni nant in aal t0 - .

has the
Hente fer al posible values I,mgn thedeserunirant
ame Vlis a imn

A'mn =timn
eijK n an iotropie Jenser es anden 3.

2
8. Aij Bij
21 2 32 |

24 642
4tMt2
21 21 6
D: AB 42
43 26 12 l2
I 2 32| J

b. C: Aijj
Au Avn A3
A21 A22 A2
LA31 A32 A3
T iven, Angte bun primed anes
or ad a and cunpeimd ameA

a Primd eomponent 86 vet8N

Jiven by v: 2e,tC)at*3

Jo derunine prime emonnts ot dhe veten V 2ei -stla

wL can
ak }un l4 tocunden cloe wihe otatign ahot ng rn,
Use a nottien matin.

roraea andiceotKwih
cos (40-0) tes(qot0)
by an angie e', which
ces 2r010) esh6e-o)
Obtaned by retauirgth
- bin o )71-emain cinhange
(oSO
both

e, nemain cnchangd el=e,

ansicloexwine by abeut dhe ng ais
 No@ o fnd 4he paimd tomponunth ot you can maipy by
Jhe RotQuon matrin a.

3/6

So, the primed te mponzrds r the vetON V= 21-l2teg atten Tly

b. r 2 64
420

Rotaien matri QfI

26 4
0

Prime compenznt 6
 Eigenvalues:
-24.9359 0
-7.1488
30.0848
Eigenvectors:
-0.3691 -0.1257 -0.9208
0.3882 0.8794 -0.2757
-0.8444 0.4592 0.2758

$0-8 (a) Use MATLAB directly to

calculate c= Ax and D= AB.
% Define tensors A, B,
A = [1, 2, 1; 1, 4, 3; and x
B= [1, 2, 1; 4, 1, 2; 2, 1, 2];
3, 2, 1]:
x= [2; 4; 1];
8 CalculateC= Ax
C-A * x:

% CalculateD= AB
D= A * B;

t Display results
disp ('c = Axt ');
disp (c) ;
disp ("D = AB:');
disp (D) ;
C= Ax:
11
21
10

D= AB:
12 6 6
26 12 12
12 6

Q-8(B)
& Define the tensors A and B
A = [1 2 1; 1 4 3; 2 1 2]:
B= [1 2 1; 4 1 2; 3 2 1]:

% Initialize the result tensor D

D= zeros (size (A) );

& Perfo rm the calculations using index notation

for i = 1:size(A, 1)
for k = l:size (B, 2)
for j = 1:size (A, 2)
D(i, k)= D(i, k) + A(i, j) * B(j, k):

3
 end
end
end

Display the resultD

disp (D using index notation: "):
disp (D):
Dusing index notation:
12 6 6
26 12 12
12 6

AB= AjBjK
Published with MATLAB® R2023b