Identificando as forças que agem sobre o bloco: mz̈ = −Fat − Fel − E + P Sendo l0 o comprimento rekaxada da mola, a o comprimento da aresta do bloco cúbico, ρ sua densidade e ρ0 a densidade do fluido no qual o bloco está mergulhado, a expressão anterior se torna: mz̈ = −ρż − k(z − l0)− ρa3g +mg Reescrevendo a EDO na forma usual: z̈ + ρρa3ż + kmz = g − ρ0ρg + kl0 A frequência natural do sistema é: ω0 = √ √ √ √ √ k ρa3 ︸︷︷︸ =m = √ 40 8 ≈ 2.23s−1 E o fator de amortecimento γ vale: γ = ρ m = 2 8 = 0.25s−1 Temos que γ/2 < ω0, portanto o movimento é subamortecido e a solução geral da EDO homogênea correspondente é da forma: xh = e− γ 2 t(a cos (ωt) + b sin (ωt)) Com ω: ω = √ ω2 0 − γ2 4 = √ 2.232 − 0.252 4 ≈ 2.23s−1 Substituindo na solução homogênea: xh = e−0.125t(a cos (2.23t) + b sin 2.23t) Para encontrar a solução particular iremos supor que ela é da forma xp = C, onde C representa uma constante, pois a expressão do lado direito da EDo também é constante. Substituindo na EDO: z̈ ︸︷︷︸ =0 + ρ ρa3 ż ︸︷︷︸ =0 + k m C = g − ρ0 ρ g + kl0 k m C = g − ρ0 ρ g + kl0 =⇒ C = g ω2 0 (1− ρ0 ρ ) + l0 Substituindo pelos valores numéricos: zp = C = 9.81 2.232 (1− 1.25 8 ) + 0.5 = 2.15m A solução geral é então: z = zh + zp = e−0.125t(a cos (2.23t) + b sin 2.23t) + 2.15 Com as constantes a e b a serem determinadas a partir das condições inicias. De acordo com o enunciado o bloco é solto a partir do repouso e a 1cm para baixo da posição de equiĺıbrio. A posição de equiĺıbrio está a 2.15m do teto, e é a mesma distância encontrada na solução particular, portanto a primeira condição inicial é z(0) = 2.15 + 0.01 = 2.16m: z(0) = 2.16 = e−0.125×0(a cos (2.23× 0) + b sin 2.23× 0) + 2.15 =⇒ a = 0.01m E pela segunda condição inicial ż(0) = 0, pois o bloco parte do repouso. Derivando z(t) obtém-se: ż(t) = −γ 2 z(t) + e− γ 2 t(−ωa sin (ωt) + ωb cos (ωt)) ż(0) = −γ 2 z(0) + ωb =⇒ b− ż(0) + γ 2 z(0) ω b = 0.01× 0.125 2.23 = 0.01× 0.056 Substituindo as constantes encontradas na solução da EDO: z = zh + zp = e−0.125t(0.01 cos (2.23t) + 0.01× 0.056 sin 2.23t) + 2.15 Simplificando: z(t) = 2.15 + 0.01e−0.125t(cos (2.23t) + 0.056 sin (2.