blob: 5d4c7f1183ea7d5ee7ee5e16ce96a97486a48ca7 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
|
clc
clear
mprintf('Mechanical vibrations by G.K.Grover\n Example 1.7.1\n')
//given data
Po=25//amplitude of force in N
Xo=0.05//ampliyude of displacement in m
W=20*%pi
//calculations
//case 1
t0=0
t1=1
v1=integrate('sin(W*t)*cos(W*t-%pi/6)','t',t0,t1)
WD1=Po*Xo*W*v1
//case 2
t0=0
t1=1/40
v2=integrate('sin(W*t)*cos(W*t-%pi/6)','t',t0,t1)
WD2=Po*Xo*W*v2
//output
mprintf(' (i)work done during the first second is %f N-m\n (ii)work done during the first 1/40th of second is %f N-m',WD1,WD2)
|