diff options
Diffstat (limited to '3875/CH1/EX1.2/Ex1_2.sce')
| -rw-r--r-- | 3875/CH1/EX1.2/Ex1_2.sce | 12 |
1 files changed, 6 insertions, 6 deletions
diff --git a/3875/CH1/EX1.2/Ex1_2.sce b/3875/CH1/EX1.2/Ex1_2.sce index 44e020473..c8b805dbc 100644 --- a/3875/CH1/EX1.2/Ex1_2.sce +++ b/3875/CH1/EX1.2/Ex1_2.sce @@ -12,20 +12,20 @@ c=m*2*damp_omega omega_n=sqrt(K/m)
//for (c)
-damping_ratio=damp_omega/omega_n
-Q=1/(2*damping_ratio)
-omega_d=omega_n*sqrt(1-damping_ratio^2) //in radian/s
+eta=damp_omega/omega_n
+Q=1/(2*eta)
+omega_d=omega_n*sqrt(1-eta^2) //in radian/s
//for (d)
-fract_change=0.5*(damping_ratio^2) //fractional change in frequency
+fract_change=0.5*(eta^2) //fractional change in frequency
percent_change=fract_change*(10^2)
mprintf("\n(a)\n")
-mprintf("Resistive force constant c = %.0e newton/s/meter\n",c)
+mprintf("Resistive force constant c = %.0e newton-s/meter\n",c)
mprintf("(b)\n")
mprintf("Natural angular frequency omega_n = %d rad/s\n",omega_n)
mprintf("(c)\n")
-mprintf("Damping ratio damping_ratio = %.0e\n",damping_ratio)
+mprintf("Damping ratio eta = %.0e\n",eta)
mprintf("Q factor = %d\n",Q)
mprintf("(d)\n")
mprintf("percent change in frequency = %.0e\n",percent_change)
|