blob: 51a34d218f534f559275426d59b46d35ae251684 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
|
clc;
// Resolving forces 125 N and 90 N
f1=125;//N
f2=90;//N
f3=200;//N
theta=45;// degree, angle made by f1 with x axis
theta=theta*%pi/180;// Conversion of angle into radian
alpha=30;// degree, angle made by force f2 with y axis
alpha=alpha*%pi/180;// Conversion of angle into radian
f1x=f1*cos(theta);//N, X component of 125 N.
f1y=-f1*sin(theta);//N, Y component of 125 N
f2x=-f2*sin(alpha);//N, X component of 90N
f2y=-f2*cos(alpha);//N, Y component of 90N
// Equivalent force at A
Rx=f1x+f2x;//N
Ry=f1y+f2y-f3;//N
R=[Rx,Ry];//N
R=norm(R);//N, magnitude of resultant
theta=atan(Ry/Rx);//radian, angle of resultant with X axis
theta=theta*180/%pi;
printf("Equivalent force at A is R=%.2fN with angle %.2f Degrees \n ",R,theta);
//equivalent couple at A
//Clockwise moments are negative
Meq=-550*f2*sin(35*%pi/180)-800*f3*sin(65*%pi/180)-1200*f1*sin(65*%pi/180);//N.mm , sum of all moments
Meq=Meq/1000;// N.m conversion into N.m
printf("Equivalent couple at A is Meq= %.2f N.m \n",Meq);
|
