O005 · C120 · Bab 4

Pemantulan Batu di Permukaan Air

Joceline Lega · Edisi Bahasa Indonesia

Tujuan Pembelajaran

Pada akhir bab ini, Anda akan mampu melakukan hal-hal berikut.

  • Memahami pemantulan batu di permukaan air sebagai pergantian antara fase terbang bebas dan fase tumbukan.
  • Menerapkan konsep-konsep mekanika klasik untuk mendeskripsikan dinamika batu selama setiap fase.
  • Menjelaskan alasan penggunaan aproksimasi untuk menyederhanakan persamaan gerak.
  • Merumuskan kondisi awal untuk setiap fase berdasarkan hasil fase sebelumnya.
  • Menganalisis dinamika setiap fase.
  • Mengevaluasi prediksi teoretis dengan membandingkannya terhadap pengamatan eksperimen.

Sekarang kita akan membahas masalah pemodelan pemantulan batu di permukaan air. Tujuannya ialah memahami bagaimana batu meluncur sambil memantul di permukaan air dan seberapa jauh batu itu bergerak. Dengan informasi ini, kita diharapkan dapat menyarankan cara melempar batu yang optimal untuk meningkatkan jumlah pantulan. Salah satu perbedaan utama dari masalah bandul yang telah dibahas sebelumnya ialah bahwa batu tidak dapat dianggap sebagai massa titik. Oleh karena itu, kita akan menggunakan dua perangkat persamaan model: hukum kekekalan momentum untuk mendeskripsikan gerak pusat massa (sebuah titik), dan persamaan gerak benda tegar untuk mendeskripsikan rotasi batu terhadap pusat massanya.

Data eksperimen

Kita mungkin pernah mencoba memantulkan batu di permukaan air, dan secara intuitif mengetahui bahwa agar batu dapat memantul berkali-kali, batu sebaiknya cukup pipih serta dilempar dengan sudut datang yang kecil dan dengan putaran. Kita juga tahu bahwa jika batu tidak dilempar dengan benar, batu mungkin memantul sekali, tetapi kemudian mulai terguling-guling di udara lalu tenggelam. Menariknya, tersedia data eksperimen terkendali untuk masalah ini. Sebuah artikel tahun 2004Catatan: C. Clanet, F. Hersen, and L. Bocquet, Secrets of successful stone-skipping, Nature 427, 29 (2004). karya C. Clanet, F. Hersen, dan L. Bocquet melaporkan pengukuran laboratorium yang memberikan kecepatan awal minimum batu, sudut antara vektor kecepatan batu dan permukaan air, serta waktu tumbukan, sebagai fungsi sudut antara batu dan permukaan air.

Asumsi dan persamaan model

Pembahasan berikut didasarkan pada artikel tahun 2003Catatan: L. Bocquet, The physics of stone skipping, Am. J. Phys. 71, 150-155 (2003). karya L. Bocquet yang berjudul The physics of stone skipping. Kita menggunakan notasi yang sama agar pembahasan lebih mudah diikuti.

Kita akan mengasumsikan bahwa gerak batu berganti-ganti antara terbang bebas, ketika batu berada di udara, dan bertumbukan dengan permukaan air. Setiap tahap ini dapat dideskripsikan secara terpisah dengan menggunakan persamaan mekanika klasik. Sebagai contoh, kondisi awal suatu fase terbang bebas adalah posisi dan kecepatan batu ketika batu lepas dari air pada akhir fase tumbukan sebelumnya, dan demikian pula sebaliknya. Kita akan mulai dengan mendeskripsikan proses tumbukan, kemudian membahas bagian gerak yang berupa terbang bebas.

Fase tumbukan

Sketsa proses tumbukan. Deskripsi panjang tersedia.
Gambar 4.1. Penampang batu pipih selama proses tumbukan (diadaptasi dari artikel L. Bocquet tahun 2003). Deskripsi Gambar
Perhatikan sebuah batu tipis, pipih, homogen, dan simetris (misalnya berbentuk persegi atau lingkaran) dengan massa total MM, dan asumsikan bahwa selama proses tumbukan batu tidak berputar pada sumbu normalnya ataupun terguling-guling, sehingga orientasinya tetap sejajar dengan orientasi semula.Catatan: Apakah asumsi ini masuk akal? Mengapa atau mengapa tidak? Oleh karena itu, masalah ini dapat disederhanakan dengan meninjau penampang batu yang melalui pusat massanya dan sejajar dengan bidang yang ditentukan oleh arah vertikal dan vektor kecepatan batu (lihat Gambar 4.1). Dengan demikian, semua titik pada penampang ini (dan sebenarnya pada seluruh batu) bergerak dengan vektor kecepatan yang sama.

Parameter masalah ini adalah ukuran batu aa (aa adalah panjang sisi batu persegi atau diameter batu berbentuk lingkaran), ketebalannya hah \ll a, vektor kecepatannya v(t)\vec v(t), dan sudut θ\theta (yang dianggap konstan selama setiap fase tumbukan) antara batu dan permukaan air. Kita akan menyebut β(t)\beta(t) sebagai sudut antara v(t)\vec v(t) dan arah horizontal, serta mendefinisikan dua basis ortonormal: (x,y,z)(\vec x, \vec y, \vec z) dan (n,y,t)(\vec n, \vec y, \vec t), sedemikian sehingga y\vec y tegak lurus terhadap bidang penampang, z\vec z mengarah ke atas, n\vec n normal terhadap permukaan batu, dan t\vec t menyinggung permukaan batu. Lihat Gambar 4.1 untuk ilustrasinya.

Karena kita telah membuat asumsi bahwa batu tidak berotasi selama proses tumbukan, geraknya dapat direduksi menjadi gerak pusat massanya dan sepenuhnya dideskripsikan oleh hukum Newton. Gaya-gaya yang bekerja pada batu ialah gaya gravitasi, Mgz- M g\, \vec z, dan gaya F=Fxx+Fzz\vec F = F_x\, \vec x + F_z\, \vec z yang dikerahkan air pada batu. Jika kita tuliskan v=vxx+vzz\vec v = v_x\, \vec x + v_z\, \vec z, kita memperoleh

Mdvxdt=Fx,Mdvzdt=Mg+Fz.(4.1)M\ \frac{d v_x}{d t} = F_x, \qquad M\ \frac{d v_z}{d t} = - M g + F_z. \qquad (4.1)

Tentu saja, sekarang kita memerlukan model untuk gaya F\vec F. Pertama, gaya ini dapat dianggap sebagai jumlah gaya angkat yang sejajar dengan n\vec n dan gaya gesek yang sejajar dengan t\vec t. Kedua, mari kita tinjau dinamika air selama proses tumbukan. Fluida bergerak menurut persamaan Navier-Stokes,

ρw(Vt+(V)V)=p+μ2V+f,\displaystyle \rho_w \left( \frac{\partial \vec V}{\partial t} + \left(\vec V \cdot \vec \nabla \right) \vec V \right)= - \nabla p + \mu \nabla^2 \vec V + \vec f,

di mana ρw\rho_w adalah massa jenis air, V\vec V adalah vektor kecepatan sebuah partikel fluida, pp adalah tekanan, μ\mu adalah viskositas dinamis air, dan f\vec f menyatakan gaya-gaya volume (seperti gravitasi) yang bekerja pada air. Pada antarmuka antara air dan batu, tekanan pp diberikan oleh

p||F||Sim,\displaystyle p \simeq \frac{||\vec F||}{S_{im}},

di mana SimS_{im} adalah luas permukaan bagian batu yang terendam (lihat Gambar 4.1). Sekarang kita akan melakukan analisis dimensi untuk membandingkan suku inersia dan suku difusi dalam persamaan ini. Kita mempunyai

[ρw(V)V][μ2V]=[ρw||V||Lμ]=[||V||Lνw]Re,\displaystyle \frac{\left[ \rho_w \left( \vec V \cdot \vec \nabla\right) \vec V \right]}{\left[\mu \vec \nabla^2 \vec V\right]} = \left[ \frac{\rho_w ||\vec V|| L}{\mu} \right] = \left[ \frac{||\vec V|| L}{\nu_w} \right] \equiv Re,

di mana LL adalah panjang karakteristik, νw=μ/ρw\nu_w = \mu / \rho_w adalah viskositas kinematik air, dan ReRe disebut bilangan Reynolds. Jika kita mengasumsikan bahwa ||V||||\vec V|| berorde beberapa meter per detik, bahwa LL berorde ukuran batu—misalnya beberapa sentimeter—dan karena untuk air νw=106\nu_w = 10^{-6} m2 s1^{-1}, kita memperoleh

Re=10102106=1051.Re = \frac{10 \cdot 10^{-2}}{10^{-6}} = 10^5 \gg 1.

Akibatnya, suku-suku viskos dapat diabaikan, dan kita dapat menyeimbangkan suku-suku inersia dengan suku-suku tekanan. Hal ini memberikan

[pL][ρw(V)V]=[ρw||V||2L],\left[ \frac{p}{L} \right] \simeq \left[ \rho_w \left( \vec V \cdot \vec \nabla\right) \vec V \right] = \left[ \frac{\rho_w ||\vec V||^2}{L} \right],

yang bersama dengan p||F||/Simp \simeq ||\vec F||/S_{im} menghasilkan

||F||ρwSim||V||2.||\vec F|| \propto \rho_w\, S_{im}\, ||\vec V||^2.

Karena ||V||||\vec V|| dan ||v||||\vec v|| sebanding, kita akan menganggap bahwa besar gaya F\vec F sebanding dengan massa jenis air ρw\rho_w, luas permukaan yang terendam SimS_{im}, dan kuadrat besar vektor kecepatan batu, ||v||2||\vec v||^2. Dengan demikian, kita dapat menuliskan

F=(12Cln+12Cft)ρwSim(vx2+vz2),\vec F = \left(\frac{1}{2} C_l\, \vec n + \frac{1}{2} C_f \, \vec t \right) \rho_w \, S_{im} \left(v_x^2 + v_z^2\right),

di mana koefisien tak berdimensi ClC_l dan CfC_f mendeskripsikan bagaimana komponen normal (angkat) dan tangensial (gesek) dari F\vec F berubah terhadap ρwSim||v||2\rho_w \, S_{im}\, ||\vec v||^2. Kita akan mengasumsikan bahwa koefisien-koefisien ini konstan. Sekarang kita dapat menyubstitusikan ekspresi ini ke dalam Persamaan (4.1). Karena n=sin(θ)x+cos(θ)z\vec n = - \sin(\theta) \, \vec x + \cos(\theta)\, \vec z dan t=cos(θ)xsin(θ)z\vec t = - \cos(\theta)\, \vec x - \sin(\theta) \vec z, kita memperoleh

Mdvxdt=12ρwSim(z)(vx2+vz2)(Clsin(θ)+Cfcos(θ)),Mdvzdt=Mg+12ρwSim(z)(vx2+vz2)(Clcos(θ)Cfsin(θ)).(4.2)\begin{array}{ll} \displaystyle M\ \frac{d v_x}{d t} = - \frac{1}{2} \rho_w \, S_{im}(z) \left(v_x^2 + v_z^2\right) \left(C_l \sin(\theta) + C_f \cos(\theta)\right),\\ \displaystyle M \ \frac{d v_z}{d t} = - M g + \frac{1}{2} \rho_w \, S_{im}(z) \left(v_x^2 + v_z^2\right) \left(C_l \cos(\theta) - C_f \sin(\theta)\right). \end{array}\qquad(4.2)

Persamaan-persamaan ini nonlinear terhadap vxv_x dan vzv_z. Selain itu, kita telah menyatakan secara eksplisit kebergantungan SimS_{im} pada koordinat vertikal zz dari bagian batu yang berada paling dalam di bawah air (lihat Gambar 4.1). Persamaan (4.2) berlaku selama batu belum terendam seluruhnya. Jika batu terendam seluruhnya, gaya angkat yang dikerahkan air pada batu akan diberikan oleh hukum Archimedes, dan batu akan tenggelam. Pada awal fase tumbukan, vz<0v_z < 0. Agar batu memantul di permukaan air, kita harus mempunyai dvz/dt>0d v_z / d t > 0, sehingga vzv_z bernilai positif pada akhir fase tumbukan. Kita melihat bahwa sudut θ\theta harus kecil agar hal ini terjadi, sebab Clcos(θ)Cfsin(θ)C_l \cos(\theta) - C_f \sin(\theta) harus bernilai positif jika kita ingin

12ρwSim(z)(vx2+vz2)(Clcos(θ)Cfsin(θ))\frac{1}{2} \rho_w \, S_{im}(z) \left(v_x^2 + v_z^2\right) \left(C_l \cos(\theta) - C_f \sin(\theta)\right)

mengimbangi Mg- M g.

Fase terbang bebas

Selama fase terbang bebas, gerak pusat massa batu dideskripsikan oleh persamaan Newton dengan gravitasi sebagai satu-satunya gaya:

Mdvxdt=0,Mdvzdt=Mg.(4.3)\displaystyle M\ \frac{d v_x}{d t} = 0, \qquad M\ \frac{d v_z}{d t}= - M g. \qquad (4.3)

Jika kita menyatakan koordinat batu dengan XX dan ZZ, kita mempunyai

X(t)=X(0)+vx(0)tZ(t)=vz(0)t12gt2,(4.4)X(t)=X(0) + v_x(0) t \qquad Z(t) = v_z(0) t -\frac{1}{2} g t^2, \qquad (4.4)

di mana t=0t=0 ditetapkan pada awal fase terbang bebas. Posisi awal (X(0),Z(0)=0)(X(0), Z(0)=0) dan kecepatan awal (vx(0),vz(0))(v_x(0),v_z(0)) batu diberikan oleh posisi dan kecepatan akhir pada fase tumbukan sebelumnya, kecuali, tentu saja, jika ini merupakan fase terbang bebas yang pertama.

Gerak rotasi batu sebagai benda tegar terhadap pusat massanya diberikan oleh persamaan Euler,

Iydωydtωnωp(InIp)=NyIndωndtωpωy(IpIy)=NnIpdωpdtωyωn(IyIn)=Np,(4.5)\begin{array}{ll}&I_{-y} \displaystyle \frac{d \omega_{-y}}{d t} - \omega_n \omega_p (I_n - I_p) = N_{-y}\\& \displaystyle I_n \frac{d \omega_n}{d t} - \omega_p \omega_{-y} (I_p - I_{-y}) = N_n \\ & \displaystyle I_p \frac{d \omega_p}{d t} - \omega_{-y} \omega_n (I_{-y} - I_n) = N_p,\end{array} \qquad (4.5)

di mana IαI_\alpha, ωα\omega_\alpha, dan NαN_\alpha masing-masing adalah momen inersia, kecepatan sudut, dan momen gaya terhadap sumbu utama benda yang sejajar dengan α\vec \alpha. Di sini, y-\vec y dan p\vec p adalah vektor satuan yang menyinggung permukaan batu dan sejajar dengan dua sumbu utamanya. Kerangka ortonormal (y,n,p)(-\vec y, \vec n, \vec p) melekat pada pusat massa batu dan berotasi bersamanya. Dalam penampang pada Gambar 4.1, arah p\vec p berimpit dengan arah tangensial t\vec t. Perangkat persamaan serupa dapat dituliskan untuk fase tumbukan jika kita ingin memperhitungkan kemungkinan batu berputar pada sumbu normalnya dan terguling-guling saat bertumbukan dengan air.

Karena batu berbentuk persegi atau lingkaran, berlaku Iy=IpJ1I_{-y} = I_p \equiv J_1, dan kita akan menggunakan notasi InJ0I_n \equiv J_0. Selain itu, satu-satunya gaya yang bekerja pada batu, Mgz-M g\, \vec z, bekerja di pusat massanya, sehingga semua momen gaya bernilai nol. Persamaan (4.5) kemudian dapat disederhanakan menjadi

J1dωydtωnωp(J0J1)=0J0dωndt=0J1dωpdtωyωn(J1J0)=0.(4.6)\begin{array}{ll}& J_1 \displaystyle \frac{d \omega_{-y}}{d t} - \omega_n \omega_p (J_0 - J_1) = 0 \\ & J_0\, \displaystyle \frac{d \omega_n}{d t} = 0 \\& J_1 \displaystyle \frac{d \omega_p}{d t} - \omega_{-y} \omega_n (J_1 - J_0) =0. \end{array} \qquad (4.6)

Analisis

Kita akan mulai dengan menganalisis persamaan untuk fase terbang bebas, lalu kembali membahas proses tumbukan.

Fase terbang bebas

Dari Persamaan (4.6), kita melihat bahwa ωn\omega_n bersifat kekal sehingga nilainya sama dengan kecepatan sudut putar awal batu, Ω0\Omega_0, terhadap sumbu yang melalui pusat massanya dan tegak lurus terhadap permukaannya. Dengan demikian, tersisa dua persamaan linear,

 J1dωydtΩ0ωp(J0J1)=0J1dωpdtωyΩ0(J1J0)=0,(4.7)\begin{array}{ll} \displaystyle J_1 \frac{d \omega_{-y}}{d t} - \Omega_0\, \omega_p (J_0 - J_1) &= 0 \\ \displaystyle J_1 \frac{d \omega_p}{d t} - \omega_{-y}\, \Omega_0 (J_1 - J_0) &= 0, \end{array} \qquad (4.7)

yang dalam bentuk matriks dapat ditulis sebagai ddt(ωyωp)=(0δδ0)(ωyωp)\displaystyle \frac{d}{dt} \begin{pmatrix}\omega_{-y} \\ \omega_p \end{pmatrix} = \begin{pmatrix} 0 & \delta \\ -\delta & 0 \end{pmatrix} \begin{pmatrix}\omega_{-y} \\ \omega_p \end{pmatrix} dengan δ=Ω0(J0J1)J1.\delta = \frac{\Omega_0 (J_0 -J_1)}{J_1}. Jika tidak ada putaran awal, yaitu jika Ω0=0\Omega_0 = 0, kita memperoleh

ωy=ωy(0)ωp=ωp(0),\omega_{-y} = \omega_{-y}(0) \qquad \omega_p = \omega_p(0),

yang, karena ωy=dθ/dt\omega_{-y}=d \theta/dt, menyiratkan

θ(t)=θ(0)+ωy(0)t.\theta(t) = \theta(0) + \omega_{-y}(0)\, t.

Dengan kata lain, jika selama fase tumbukan batu memperoleh sedikit saja rotasi terhadap sumbu y-\vec y, yaitu jika ωy(0)\omega_{-y}(0) tidak nol, batu akan mulai berjungkir selama fase terbang bebas. Pada tumbukan berikutnya, batu kemungkinan besar menghantam air dengan sudut datang θ\theta yang besar, lalu tenggelam.

Karena itu, Ω0\Omega_0 harus tidak nol. Dalam keadaan ini, ωy\omega_{-y} dan ωp\omega_p sama-sama berosilasi terhadap waktu dengan frekuensi |δ||\delta| dan diberikan oleh

ωy(t)=ωy(0)cos(δt)+ωp(0)sin(δt)ωp(t)=ωp(0)cos(δt)ωy(0)sin(δt).\substack{\omega_{-y}(t) = \omega_{-y}(0) \cos(\delta t) + \omega_p(0) \sin(\delta t) \\ \omega_p(t) = \omega_p(0) \cos(\delta t) - \omega_{-y}(0) \sin(\delta t). }

Akibatnya,

θ(t)=θ(0)+ωy(0)δsin(δt)+ωp(0)δ[1cos(δt)]\theta(t) = \theta(0) + \frac{\omega_{-y}(0)}{\delta} \sin(\delta t) + \frac{\omega_p(0)}{\delta} \left[1-\cos(\delta t)\right]

dan nilainya akan tetap dekat dengan θ(0)\theta(0) jika koefisien-koefisien dalam persamaan di atas kecil. Jadi, putaran terhadap arah n\vec n menstabilkan batu dan mencegahnya berjungkir. Gejala ini disebut efek giroskopik.

Sekarang kita beralih ke gerak pusat massa batu selama fase terbang bebas. Persamaan (4.4) menunjukkan bahwa durasi τf\tau_f setiap fase terbang bebas memenuhi Z(τf)=0Z(\tau_f)=0, yaitu

τf=2vz(0)g,vz(0)>0,\tau_f = \frac{2 v_z(0)}{g}, \qquad v_z(0) \gt 0,

dan panjang setiap lompatan, yakni jarak yang ditempuh batu sepanjang arah x\vec x, adalah

λ=vx(0)τf=vx(0)2vz(0)g.\lambda = v_x(0) \tau_f = v_x(0) \frac{2 v_z(0)}{g}.

Dari Persamaan (4.2), jelas bahwa ketika θ\theta kecil, vxv_x berkurang selama setiap fase tumbukan. Karena itu, nilai vx(0)v_x(0)—yakni nilai vxv_x pada awal setiap fase terbang bebas—akan berkurang dari satu fase terbang bebas ke fase berikutnya. Karena kita tidak mengharapkan vzv_z meningkat di antara dua fase terbang bebas, panjang setiap lompatan akan berkurang dari satu lompatan ke lompatan berikutnya, sebagaimana yang kita harapkan.

Fase tumbukan

Analisis kita terhadap Persamaan (4.2) pertama-tama harus menunjukkan dalam kondisi apa batu akan muncul kembali dan lepas dari air pada akhir fase tumbukan. Perhatikan persamaan untuk vzv_z,

Mdvzdt=Mg+12ρwSim(z)(vx2+vz2)(Clcos(θ)Cfsin(θ)).(4.10)M \displaystyle \frac{d v_z}{d t} = - M g + \frac{1}{2} \rho_w \, S_{im}(z) \left(v_x^2 + v_z^2\right) \left(C_l \cos(\theta) - C_f \sin(\theta)\right) . \quad (4.10)

Kita dapat memahami informasi yang terkandung dalam persamaan ini dengan menggunakan pendekatan vxvzv_x \gg v_z, yang cukup masuk akal, dan juga mengasumsikan bahwa vx2+vz2=v2v02v_x^2 + v_z^2 = v^2 \equiv v_0^2 konstan selama fase tumbukan. Pernyataan terakhir ini tidak sepenuhnya benar karena gesekan dengan air memperlambat batu. Namun, jika batu memantul sekitar 20 kali, kehilangan relatif (v2(0)v2(τc))/v2(0)(v^2(0) - v^2(\tau_c))/v^2(0) selama setiap tumbukan kecil, setidaknya untuk sebagian besar tumbukan. Di sini, τc\tau_c menyatakan durasi tumbukan. Dengan syarat-syarat ini, Persamaan (4.10) menjadi

Mdvzdt=Mg+12ρwSim(z)v02(Clcos(θ)Cfsin(θ))(z).(4.11)M \displaystyle \frac{d v_z}{d t} = - M g + \frac{1}{2} \rho_w S_{im}(z) v_0^2 \left(C_l \cos(\theta) - C_f \sin(\theta) \right) \equiv - {\mathcal F}(z). \quad (4.11)

Ruas kanan, (z)- {\mathcal F}(z), hanya merupakan fungsi dari zz (ingat bahwa sejak awal kita mengasumsikan θ\theta konstan selama setiap fase tumbukan). Karena vz=dz/dtv_z = d z / d t, kita dapat mengalikan kedua ruas persamaan ini dengan vzv_z dan mengintegralkannya terhadap tt untuk memperoleh

M2(dzdt)2=(z)dzdtdt=(z)dz𝒰(z)+E,(4.12)\displaystyle \frac{M}{2} \left(\frac{d z}{d t}\right)^2 = \int - {\mathcal F}(z) \frac{d z}{d t} \, dt = \int - {\mathcal F}(z) \, dz \equiv -\, {\mathcal U}(z) + E, \quad (4.12)

di mana 𝒰(z)=0z(s)ds,{\mathcal U}(z) = \int_0^z {\mathcal F}(s)\, ds, dan konstanta EE bersifat sembarang. Ketika z=0z=0, Sim(z)=0S_{im}(z)=0; selain itu, Sim(z)S_{im}(z) bertambah ketika zz menjadi semakin negatif. Oleh karena itu, potensial 𝒰(z){\mathcal U}(z) berbentuk seperti yang digambarkan pada Gambar 4.2 untuk nilai zz negatif, asalkan θ\theta cukup kecil.

Bentuk potensial U; tersedia deskripsi panjang.
Gambar 4.2. Bentuk potensial U(z) sebagai fungsi z untuk z < 0. Deskripsi Gambar
Dari bentuk 𝒰{\mathcal U}, kita dapat melihat bahwa lintasan batu sedemikian rupa sehingga zz memiliki titik balik pada z=zmz=z_m. Pada akhir fase tumbukan, batu akan muncul kembali dengan kecepatan vertikal vz(τc)=vz(0)v_z(\tau_c) = - v_z(0), berlawanan dengan kecepatan vertikalnya pada awal fase tumbukan. Hal ini hanya terjadi jika energi EE memungkinkan batu belum terbenam seluruhnya ketika z=zmz=z_m. Jika EE dipertahankan tetap, keadaan ini memberlakukan syarat pada v02v_0^2 karena semakin besar v02v_0^2, semakin kecil |zm||z_m|. Jadi, batu tidak akan tenggelam selama fase tumbukan apabila v02>vc2v_0^2 > v_c^2, dengan kecepatan kritis vcv_c yang bergantung pada parameter-parameter masalah. Karena vxvzv_x \gg v_z, syarat pada v0v_0 juga dapat ditulis kembali sebagai syarat pada kecepatan awal batu dalam arah x\vec x, yaitu vx(0)>vcv_x(0) > v_c. Sebagaimana ditunjukkan dalam artikel L. Bocquet tahun 2003 dan dalam latihan pada akhir bab ini, kita dapat menghitung Sim(z)S_{im}(z) untuk batu dengan bentuk tertentu, menyelesaikan Persamaan (4.11) secara eksplisit, dan memperoleh ekspresi untuk vcv_c.

Seperti disebutkan di atas, kita mengasumsikan bahwa dinamika sistem bersifat konservatif (yaitu v2=v02v^2 = v_0^2 dianggap konstan). Hal ini jelas tidak sepenuhnya benar karena batu tidak terus memantul tanpa henti. Energi hilang selama proses tumbukan akibat gesekan batu dengan permukaan air. Semua yang telah kita nyatakan sebelumnya tetap benar secara kualitatif—khususnya, batu akan muncul kembali dengan kecepatan vertikal yang hanya sedikit lebih kecil daripada |vz(0)||v_z(0)|—selama kehilangan energi relatif dalam satu fase tumbukan kecil. Kondisi ini tentu terpenuhi pada pantulan-pantulan awal, tetapi tidak lagi berlaku pada pantulan-pantulan terakhir. Karena kehilangan energi inilah penyebab utama batu melambat, sekarang kita akan mencoba memperkirakan besarannya. Perubahan energi selama tumbukan, yang dinyatakan dengan 𝒲\mathcal W, sama dengan usaha yang dilakukan oleh gaya gesek air pada batu dan bernilai negatif untuk disipasi. Secara kasar, nilainya dapat diperkirakan sebagai

𝒲[gaya sepanjang x][jarak yang ditempuh batu sepanjang x]Fxl,\begin{align} {\mathcal W} & \simeq \left[ \text{gaya sepanjang } \vec x \right] \cdot \left[ \text{jarak yang ditempuh batu sepanjang } \vec x \right] \\ & \simeq \vec F \cdot \vec x l, \end{align}

di mana l=vx(0)τcl = v_x(0)\ \tau_c. Dari Persamaan (4.11), kita melihat bahwa

[ρwMv02]=L1T2,\left[\frac{\rho_w}{M} v_0^2 \right] = L^{-1} T^{-2},

sehingga kita dapat memperkirakan

τcMaρwv02Maρw1vx(0),\tau_c \propto \displaystyle \sqrt{\frac{M}{a \rho_w v_0^2}} \simeq \sqrt{\frac{M}{a \rho_w}} \frac{1}{v_x(0)},

di mana aa adalah ukuran batu (panjang sisi atau diameter). Oleh karena itu, pada orde terendah, l=vx(0)τcl = v_x(0) \tau_c tidak bergantung pada vx(0)v_x(0). Selama tumbukan, ruas kanan Persamaan (4.10) bernilai kecil, sehingga kita dapat menulis

MgFz.M g \simeq \vec F \cdot \vec z.

Selain itu, karena

FxFz=Clsin(θ)+Cfcos(θ)Clcos(θ)Cfsin(θ)μ,\displaystyle \frac{\vec F \cdot \vec x}{\vec F \cdot \vec z} = - \frac{C_l \sin(\theta) + C_f \cos(\theta)}{C_l \cos(\theta) - C_f \sin(\theta)} \equiv - \mu,

(rasio μ\mu tidak boleh tertukar dengan viskositas dinamis air), kita memperoleh

FxMgμ,\vec F \cdot \vec x \simeq - M g\, \mu,

dan

𝒲Mgμl.{\mathcal W} \simeq - M g\, \mu \, l.

Jadi, karena selisih energi kinetik antara awal dan akhir fase tumbukan sama dengan 𝒲\mathcal W, kita memperoleh

12Mvx2(τc)12Mvx2(0)=𝒲=Mgμl,\frac{1}{2} M v_x^2(\tau_c) - \frac{1}{2} M v_x^2(0) = {\mathcal W} = - M g\, \mu\, l,

dan kecepatan awal vx(0)v_x(0) harus memenuhi

vx(0)>2gμlVc.v_x(0) > \sqrt{2\, g\, \mu\, l} \equiv V_c.

Dengan menggabungkan seluruh hasil tersebut, kecepatan dalam arah x\vec x pada awal fase tumbukan harus memenuhi

vx(0)>max(vc,Vc),v_x(0) > \max(v_c,V_c),

di mana vcv_c diperoleh dengan menetapkan bahwa batu tidak boleh terbenam seluruhnya selama tumbukan, sedangkan VcV_c memastikan batu bergerak cukup cepat untuk mengimbangi kehilangan energi akibat gaya gesek. Pembahasan kita tentang fase terbang bebas menunjukkan bahwa panjang setiap lompatan sebanding dengan kecepatan horizontal vxv_x. Kita tahu bahwa vxv_x bersifat kekal selama fase terbang bebas (lihat Persamaan (4.3)) dan bahwa vx2v_x^2 berkurang sebesar 2gμl2\, g\, \mu\, l selama setiap fase tumbukan. Karena |vz||v_z| sama pada awal dan akhir setiap fase terbang bebas, dan karena vzv_z hanya berubah sedikit selama fase tumbukan, panjang lompatan setelah NN fase tumbukan diberikan oleh

λN=2|vz(0)|gvx,022Ngμl=2|vz(0)|vx,0g12Ngμlvx,02=λ01NNc,λ0=2|vz(0)|vx,0g,Nc=vx,022gμl,\begin{align}\lambda_N &= \frac{2 |v_z(0)|}{g} \sqrt{v_{x,0}^2 - 2 \, N \, g\, \mu\, l} = \frac{2 |v_z(0)| v_{x,0}}{g} \sqrt{1 - \frac{2 \, N \, g\, \mu\, l}{v_{x,0}^2}}\\ &= \lambda_0 \sqrt{1 - \frac{N}{N_c}}, \qquad \lambda_0 = \frac{2 |v_z(0)| v_{x,0}}{g}, \ N_c = \frac{v_{x,0}^2}{2 \, g\, \mu\, l}, \end{align}

di mana vx,0v_{x,0} adalah kecepatan batu dalam arah x\vec x ketika mula-mula dilempar. Jadi, panjang setiap lompatan berkurang ketika NN membesar dan berskala seperti 1N/Nc\sqrt{1-N/N_c}.

Ringkasan

Model ini menjelaskan adanya kecepatan minimum untuk melempar batu; menjelaskan mengapa batu perlu diberi putaran; dan menggambarkan bagaimana panjang setiap lompatan berkurang ketika jumlah pantulan bertambah. Model ini juga menunjukkan bahwa sudut θ\theta antara batu dan permukaan air harus kecil. Eksperimen yang diuraikan dalam artikel C. Clanet et al.Catatan: Gambar 1.b dan 1.c dalam artikel C. Clanet et al. tahun 2004 menunjukkan daerah keberhasilan pemantulan batu pada bidang (θ\theta,||v||||v||) dan (θ\theta,β\beta); Gambar 1.d menunjukkan waktu tumbukan sebagai fungsi θ\theta untuk berbagai nilai β\beta (lihat Gambar 4.1 untuk definisi parameter-parameter tersebut). menunjukkan bahwa θ=20\theta = 20^\circ merupakan nilai optimal. Nilai θ\theta ini meminimumkan waktu tumbukan; ini juga merupakan sudut datang yang memungkinkan pemantulan batu berhasil dengan laju awal paling kecil. Untuk menjelaskan pengamatan ini, misalnya, kita perlu memperhitungkan ketergantungan ClC_l dan CfC_f pada θ\theta, lalu mencari nilai θ\theta yang meminimumkan kehilangan energi 𝒲\mathcal W.

Deskripsi Gambar

Gambar 4.1: Diagram teknis fisika yang memperlihatkan penampang batu pipih (digambarkan berwarna abu-abu) ketika menghantam permukaan air pada suatu sudut. Diagram ini mencakup sistem koordinat Kartesius (x,y,zx, y, z), vektor satuan normal (n\vec{n}) dan tangensial (t\vec{t}) yang terkait dengan batu, serta sudut kemiringan batu (θ\theta), vektor kecepatan batu (v\vec{v}), dan sudut datang (β\beta). [Kembali ke Gambar 4.1]

Gambar 4.2: Grafik yang menunjukkan kurva energi potensial sebagai fungsi kedalaman zz. Garis horizontal EE menyatakan energi mekanik total sistem. Perpotongan garis EE dengan kurva pada zmz_m mendefinisikan sebuah “titik balik”, yaitu kedalaman maksimum yang dapat dicapai benda sebelum gaya apung atau gaya hidrodinamik mendorongnya kembali ke atas. [Kembali ke Gambar 4.2]

Bahan Renungan

Soal 1

Perhatikan Persamaan (4.5).

  1. Gunakan persamaan-persamaan tersebut untuk menentukan dimensi ωα\omega_\alpha, dengan α\alpha masing-masing bernilai y,n-y,\,n, dan pp.
  2. Hubungkan dimensi IαI_\alpha dengan dimensi NαN_\alpha.
  3. Diketahui bahwa torsi berbentuk N=r×f\vec N = \vec r \times \vec f, dengan r\vec r sebagai vektor posisi dan f\vec f sebagai gaya. Tentukan dimensi NαN_\alpha.
  4. Gunakan hasil-hasil di atas untuk menentukan dimensi IαI_\alpha. Apakah hasil ini sesuai dengan definisi Iα||r||2ρdV,I_\alpha \equiv \int ||\vec r||^2 \rho\,dV, dengan ρ\rho sebagai massa jenis dan dVdV sebagai elemen volume?

Soal 2

Soal ini didasarkan pada artikel L. Bocquet tahun 2003.

Perhatikan Persamaan (4.11) dan anggap bahwa batu berbentuk persegi. Ambil t=0t=0 pada awal tumbukan, dengan kondisi awal z(0)=0z(0)=0 dan ż(0)=vz(0)<0\dot z(0)=v_z(0)\lt 0. Rumus luas terendam pada soal ini berlaku selama perendaman parsial, yaitu asin(θ)z0-a\sin(\theta)\leq z\leq 0.

  1. Buktikan bahwa Sim=a|z|/sin(θ)S_{im} = a |z| / \sin(\theta).
  2. Dengan menggunakan bentuk SimS_{im} tersebut, buktikan bahwa Persamaan (4.11) menjadi persamaan diferensial linear untuk zz.
  3. Carilah solusi umum persamaan diferensial biasa ini.
  4. Terapkan syarat batas dan buktikan bahwa z(t)z(t) diberikan oleh

z(t)=gω02+gω02cos(ω0t)+vz(0)ω0sin(ω0t),\displaystyle z(t) = \frac{-g}{\omega_0^2} + \frac{g}{\omega_0^2} \cos(\omega_0 t) + \frac{v_z(0)}{\omega_0} \sin(\omega_0 t),

dengan ω02=(Clcos(θ)Cfsin(θ))ρwvx(0)2a2Msin(θ)\displaystyle \text{dengan } \qquad \omega_0^2 = \frac{\left(C_l \cos(\theta)-C_f \sin(\theta)\right) \rho_w v_x(0)^2 a}{2 M \sin(\theta)}.


Soal 3

Soal ini didasarkan pada artikel L. Bocquet tahun 2003.

Dinamika sisi batu persegi yang terendam selama fase tumbukan dideskripsikan oleh (lihat Soal 2). Untuk bagian-bagian berikut, abaikan ketebalan hh, tetapkan CClcos(θ)Cfsin(θ)>0C \equiv C_l\cos(\theta)-C_f\sin(\theta)\gt 0, gunakan vz(0)=vx(0)tan(β)v_z(0)=-v_x(0)\tan(\beta), dan asumsikan ruas dalam tanda kurung siku pada hasil akhir positif. Dengan konvensi z<0z\lt 0 di bawah permukaan, batu mulai terendam seluruhnya ketika |z|asin(θ)|z|\geq a\sin(\theta).

z(t)=gω02+gω02cos(ω0t)+vz(0)ω0sin(ω0t),\displaystyle z(t) = \frac{-g}{\omega_0^2} + \frac{g}{\omega_0^2} \cos(\omega_0 t) + \frac{v_z(0)}{\omega_0} \sin(\omega_0 t),

dengan

ω02=(Clcos(θ)Cfsin(θ))ρwvx(0)2a2Msin(θ).\displaystyle \omega_0^2 = \frac{\left(C_l \cos(\theta)-C_f \sin(\theta)\right) \rho_w v_x(0)^2 a}{2 M \sin(\theta)}.

  1. Gunakan bentuk ini untuk membuktikan bahwa sisi batu mencapai kedalaman |zm||z_m|, dengan zm=gω02[1+1+ω02vz(0)2g2].z_m = - \frac{g}{\omega_0^2} \left[1+\sqrt{1+\frac{\omega_0^2 v_z(0)^2}{g^2}}\right].
  2. Syarat apa yang harus dipenuhi zz agar batu terendam seluruhnya?
  3. Buktikan bahwa batu tidak akan terendam seluruhnya selama fase tumbukan jika

    vx(0)2>4MgCρwa2[12tan2βMa3Cρwsin(θ)]1,C=Clcos(θ)Cfsin(θ).\displaystyle v_x(0)^2 > \frac{4 M g}{C \rho_w a^2} \left[1 - \frac{2 \tan^2 \beta \, M}{a^3 C \rho_w \sin(\theta)}\right]^{-1}, \quad C = C_l \cos(\theta)-C_f \sin(\theta).


Soal 4

Soal ini didasarkan pada artikel L. Bocquet tahun 2003.

Dengan menggunakan persamaan gerak Euler (4.5) bersama bentuk gaya angkat yang sesuai pada batu selama fase tumbukan, turunkan Persamaan (20) dalam artikel L. Bocquet tahun 2003,

d2θdt2+ν2(θθ(0))=θJ1,\displaystyle \frac{d^2 \theta}{d t^2} + \nu^2 \left(\theta - \theta(0)\right) = \frac{{\mathcal M}_\theta}{J_1},

dengan ν=((J0J1)/J1)Ω0\nu = ((J_0-J_1)/J_1)\,\Omega_0, sedangkan θNy{\mathcal M}_\theta \equiv N_{-y} adalah komponen torsi yang diberikan air pada batu dalam arah peningkatan sudut kemiringan (sumbu y-\vec y). Bentuk ini mempertahankan kuadrat frekuensi ν2\nu^2 pada Persamaan (20) artikel sumber.

Petunjuk, pemeriksaan, dan pembahasan

Bagian tambahan ini ditulis untuk edisi Bahasa Indonesia. Bukalah seperlunya setelah Anda berusaha menyelesaikan soal secara mandiri.

Dukungan untuk Soal 1

Petunjuk

Sudut tidak mempunyai dimensi fisik, sehingga setiap komponen kecepatan sudut berdimensi waktu invers. Setelah itu, samakan dimensi semua suku pada satu persamaan Euler dan gunakan [f]=MLT2[\vec f]=MLT^{-2}.

Periksa jawaban
Pemeriksaan akhir

Untuk α{y,n,p}\alpha\in\{-y,n,p\}, berlaku [ωα]=T1[\omega_\alpha]=T^{-1}. Persamaan Euler memberi [Nα]=[Iα]T2[N_\alpha]=[I_\alpha]T^{-2}. Karena N=r×f\vec N=\vec r\times\vec f, maka [Nα]=L(MLT2)=ML2T2[N_\alpha]=L(MLT^{-2})=ML^2T^{-2}, sehingga [Iα]=ML2[I_\alpha]=ML^2. Definisi Iα=||r||2ρdVI_\alpha=\int ||\vec r||^2\rho\,dV menghasilkan dimensi yang sama: L2(ML3)(L3)=ML2L^2(ML^{-3})(L^3)=ML^2.

Cek cepat

Suku Iαω̇αI_\alpha\dot\omega_\alpha, suku IωωI\omega\omega, dan NαN_\alpha semuanya harus berdimensi ML2T2ML^2T^{-2}.

Pembahasan atau rubrik
Langkah penyelesaian
  • Penjelasan

    Kecepatan sudut adalah laju perubahan sudut; radian merupakan besaran tak berdimensi.

    Rumus

    [ωy]=[ωn]=[ωp]=T1,[ω̇α]=T2[\omega_{-y}]=[\omega_n]=[\omega_p]=T^{-1},\qquad [\dot\omega_\alpha]=T^{-2}

  • Penjelasan

    Dalam setiap persamaan Euler, suku turunan dan suku kopling harus memiliki dimensi yang sama.

    Rumus

    [Iαω̇α]=[Iωω]=[Iα]T2=[Nα][I_\alpha\dot\omega_\alpha]=[I\omega\omega]=[I_\alpha]T^{-2}=[N_\alpha]

  • Penjelasan

    Torsi adalah hasil kali silang lengan posisi dan gaya. Hasil kali silang mengubah arah, tetapi tidak mengubah perkalian dimensinya.

    Rumus

    [Nα]=[r][f]=L(MLT2)=ML2T2[N_\alpha]=[\vec r][\vec f]=L(MLT^{-2})=ML^2T^{-2}

  • Penjelasan

    Gabungkan dua relasi sebelumnya untuk memperoleh dimensi momen inersia.

    Rumus

    [Iα]=[Nα]T2=ML2[I_\alpha]=[N_\alpha]T^2=ML^2

  • Penjelasan

    Periksa secara independen melalui definisi kontinu, dengan massa jenis berdimensi massa per volume.

    Rumus

    [||r||2ρdV]=L2(ML3)L3=ML2\left[\int ||\vec r||^2\rho\,dV\right]=L^2(ML^{-3})L^3=ML^2

Simpulan

Analisis persamaan gerak dan definisi integral memberikan hasil yang identik. Tanda pada nama sumbu y-y hanya menyatakan orientasi dan tidak memengaruhi dimensi.

Dukungan untuk Soal 2

Petunjuk

Selama perendaman parsial, panjang bagian sisi batu yang berada di bawah permukaan air adalah |z|/sinθ|z|/\sin\theta. Karena z<0z<0, gunakan |z|=z|z|=-z, tetapkan C=ClcosθCfsinθC=C_l\cos\theta-C_f\sin\theta, dan kumpulkan koefisien zz sebagai ω02\omega_0^2.

Periksa jawaban
Pemeriksaan akhir

Luas terendam ialah Sim=a|z|/sinθS_{im}=a|z|/\sin\theta. Dengan pendekatan v02vx(0)2v_0^2\simeq v_x(0)^2, C=ClcosθCfsinθ>0C=C_l\cos\theta-C_f\sin\theta>0, dan z<0z<0, Persamaan (4.11) menjadi z̈+ω02z=g\ddot z+\omega_0^2z=-g, dengan ω02=Cρwvx(0)2a/(2Msinθ)\omega_0^2=C\rho_wv_x(0)^2a/(2M\sin\theta). Solusi umumnya z=g/ω02+Acos(ω0t)+Bsin(ω0t)z=-g/\omega_0^2+A\cos(\omega_0t)+B\sin(\omega_0t). Kondisi z(0)=0z(0)=0 dan ż(0)=vz(0)\dot z(0)=v_z(0) memberi A=g/ω02A=g/\omega_0^2, B=vz(0)/ω0B=v_z(0)/\omega_0, sehingga diperoleh rumus yang diminta.

Cek cepat

Substitusi t=0t=0 harus memberi z(0)=0z(0)=0, sedangkan diferensiasi lalu substitusi t=0t=0 harus memberi ż(0)=vz(0)\dot z(0)=v_z(0).

Pembahasan atau rubrik
Langkah penyelesaian
  • Penjelasan

    Penampang batu membentuk sudut θ\theta terhadap air. Kedalaman vertikal |z||z| bersesuaian dengan panjang terendam |z|/sinθ|z|/\sin\theta sepanjang batu; kalikan dengan lebar persegi aa.

    Rumus

    Sim=a|z|sinθS_{im}=a\frac{|z|}{\sin\theta}

  • Penjelasan

    Selama fase masuk, z<0z<0, jadi |z|=z|z|=-z. Masukkan luas tersebut ke persamaan gaya vertikal dengan besar kecepatan dianggap tetap dan didominasi oleh komponen horizontal.

    Rumus

    Mz̈=MgCρwvx(0)2a2sinθz,C=ClcosθCfsinθM\ddot z=-Mg-\frac{C\rho_wv_x(0)^2a}{2\sin\theta}z,\qquad C=C_l\cos\theta-C_f\sin\theta

  • Penjelasan

    Untuk C>0C>0, definisikan frekuensi karakteristik real dan bagi persamaan dengan massa.

    Rumus

    ω02=Cρwvx(0)2a2Msinθ,z̈+ω02z=g\omega_0^2=\frac{C\rho_wv_x(0)^2a}{2M\sin\theta},\qquad \ddot z+\omega_0^2z=-g

  • Penjelasan

    Jumlahkan satu solusi partikular konstan dengan solusi homogen osilator harmonik.

    Rumus

    z(t)=gω02+Acos(ω0t)+Bsin(ω0t)z(t)=-\frac{g}{\omega_0^2}+A\cos(\omega_0t)+B\sin(\omega_0t)

  • Penjelasan

    Gunakan kedua kondisi awal pada saat sisi batu pertama kali menyentuh air.

    Rumus

    A=gω02,B=vz(0)ω0A=\frac{g}{\omega_0^2},\qquad B=\frac{v_z(0)}{\omega_0}

  • Penjelasan

    Masukkan konstanta untuk memperoleh lintasan selama asumsi perendaman parsial masih berlaku.

    Rumus

    z(t)=gω02+gω02cos(ω0t)+vz(0)ω0sin(ω0t)z(t)=-\frac{g}{\omega_0^2}+\frac{g}{\omega_0^2}\cos(\omega_0t)+\frac{v_z(0)}{\omega_0}\sin(\omega_0t)

Simpulan

Solusi ini hanya berlaku sebelum batu terendam seluruhnya, untuk sudut tetap, koefisien gaya tetap, dan perubahan besar kecepatan yang dapat diabaikan selama satu tumbukan.

Dukungan untuk Soal 3

Petunjuk

Tulis bagian berosilasi dari z(t)+g/ω02z(t)+g/\omega_0^2 sebagai satu sinusoid; amplitudonya langsung memberi nilai minimum. Setelah menetapkan |zm|<asinθ|z_m|<a\sin\theta, kuadratkan hanya sesudah memastikan kedua ruas yang relevan positif.

Periksa jawaban
Pemeriksaan akhir

Untuk vz(0)<0v_z(0)<0, kedalaman maksimum ialah zm=gω02[1+1+ω02vz(0)2/g2]z_m=-\frac{g}{\omega_0^2}\left[1+\sqrt{1+\omega_0^2v_z(0)^2/g^2}\right]. Batu persegi terendam seluruhnya ketika |z|asinθ|z|\ge a\sin\theta, sehingga pantulan parsial memerlukan |zm|<asinθ|z_m|<a\sin\theta. Dengan C=ClcosθCfsinθ>0C=C_l\cos\theta-C_f\sin\theta>0, vz(0)=vx(0)tanβv_z(0)=-v_x(0)\tan\beta, dan D=12Mtan2β/(a3Cρwsinθ)>0D=1-2M\tan^2\beta/(a^3C\rho_w\sin\theta)>0, syarat akhirnya adalah vx(0)2>4MgCρwa2D1v_x(0)^2>\frac{4Mg}{C\rho_wa^2}D^{-1}. Jika D0D\le0, tidak ada kecepatan horizontal berhingga yang memenuhi syarat tersebut dalam model dengan β\beta tetap.

Cek cepat

Faktor 4Mg/(Cρwa2)4Mg/(C\rho_wa^2) berdimensi kecepatan kuadrat, sedangkan DD tak berdimensi; batas untuk β=0\beta=0 mereduksi menjadi vx(0)2>4Mg/(Cρwa2)v_x(0)^2>4Mg/(C\rho_wa^2).

Pembahasan atau rubrik
Langkah penyelesaian
  • Penjelasan

    Bagian sinus dan kosinus mempunyai amplitudo resultan yang diperoleh dari jumlah kuadrat koefisiennya.

    Rumus

    R=(gω02)2+(vz(0)ω0)2=gω021+ω02vz(0)2g2R=\sqrt{\left(\frac{g}{\omega_0^2}\right)^2+\left(\frac{v_z(0)}{\omega_0}\right)^2}=\frac{g}{\omega_0^2}\sqrt{1+\frac{\omega_0^2v_z(0)^2}{g^2}}

  • Penjelasan

    Nilai minimum lintasan adalah pergeseran konstan dikurangi amplitudo tersebut.

    Rumus

    zm=gω02R=gω02[1+1+ω02vz(0)2g2]z_m=-\frac{g}{\omega_0^2}-R=-\frac{g}{\omega_0^2}\left[1+\sqrt{1+\frac{\omega_0^2v_z(0)^2}{g^2}}\right]

  • Penjelasan

    Jika sisi terendah berada pada z<0z<0, sisi tertinggi berada pada z+asinθz+a\sin\theta. Seluruh batu berada di bawah air tepat ketika sisi tertinggi tidak lagi positif.

    Rumus

    z+asinθ0|z|asinθz+a\sin\theta\le0\quad\Longleftrightarrow\quad |z|\ge a\sin\theta

  • Penjelasan

    Minta agar titik balik terjadi sebelum batas perendaman penuh. Setelah mengisolasi akar dan menguadratkan kedua ruas positif, syarat ini menyederhana menjadi batas pada ω02\omega_0^2.

    Rumus

    |zm|<asinθω02>vz(0)2a2sin2θ+2gasinθ|z_m|<a\sin\theta\quad\Longleftrightarrow\quad \omega_0^2>\frac{v_z(0)^2}{a^2\sin^2\theta}+\frac{2g}{a\sin\theta}

  • Penjelasan

    Gunakan hubungan geometris kecepatan datang dan ekspresi frekuensi dari model gaya angkat.

    Rumus

    vz(0)=vx(0)tanβ,ω02=Cρwvx(0)2a2Msinθ,C>0v_z(0)=-v_x(0)\tan\beta,\qquad \omega_0^2=\frac{C\rho_wv_x(0)^2a}{2M\sin\theta},\qquad C>0

  • Penjelasan

    Kumpulkan semua suku yang mengandung vx(0)2v_x(0)^2 pada satu ruas. Pembagian mempertahankan arah pertidaksamaan hanya bila faktor dalam kurung positif.

    Rumus

    vx(0)2>4MgCρwa2[12Mtan2βa3Cρwsinθ]1v_x(0)^2>\frac{4Mg}{C\rho_wa^2}\left[1-\frac{2M\tan^2\beta}{a^3C\rho_w\sin\theta}\right]^{-1}

Simpulan

Batas ini merupakan syarat dalam model ideal: batu masih harus terendam sebagian, C>0C>0, sudut dan koefisien gaya dianggap tetap, vx|vz|v_x\gg|v_z|, dan faktor penyebut harus positif.

Dukungan untuk Soal 4

Petunjuk

Pilih sumbu utama (1,2,3)=(y,n,p)(1,2,3)=(-y,n,p) dengan I1=I3=J1I_1=I_3=J_1, I2=J0I_2=J_0, dan ω1=θ̇\omega_1=\dot\theta. Ambil N1=θN_1={\mathcal M}_\theta, N2=N3=0N_2=N_3=0, lalu integralkan persamaan ketiga setelah persamaan kedua memberi ω2=Ω0\omega_2=\Omega_0.

Periksa jawaban
Pemeriksaan akhir

Dengan I1=I3=J1I_1=I_3=J_1, I2=J0I_2=J_0, ω1=θ̇\omega_1=\dot\theta, hanya torsi angkat N1=θN_1={\mathcal M}_\theta, serta N2=N3=0N_2=N_3=0, persamaan Euler kedua memberi ω2=Ω0\omega_2=\Omega_0. Jika ω3(0)=0\omega_3(0)=0, persamaan ketiga memberi ω3=((J1J0)/J1)Ω0(θθ0)\omega_3=((J_1-J_0)/J_1)\Omega_0(\theta-\theta_0). Substitusi ke persamaan pertama menghasilkan θ̈+[(J0J1)Ω0/J1]2(θθ0)=θ/J1\ddot\theta+[(J_0-J_1)\Omega_0/J_1]^2(\theta-\theta_0)={\mathcal M}_\theta/J_1, yaitu θ̈+ν2(θθ0)=θ/J1\ddot\theta+\nu^2(\theta-\theta_0)={\mathcal M}_\theta/J_1.

Cek cepat

ν\nu berdimensi T1T^{-1}, sehingga ν2(θθ0)\nu^2(\theta-\theta_0) dan θ/J1{\mathcal M}_\theta/J_1 sama-sama berdimensi T2T^{-2}. Koefisien tanpa Ω02\Omega_0^2 dan tanpa kuadrat rasio inersia tidak lolos pemeriksaan ini.

Pembahasan atau rubrik
Langkah penyelesaian
  • Penjelasan

    Gunakan sumbu 1 untuk perubahan kemiringan, sumbu 2 untuk normal batu, dan sumbu 3 untuk arah tangensial pada bidang batu. Simetri batu pipih memberi dua momen inersia yang sama.

    Rumus

    (1,2,3)=(y,n,p),I1=I3=J1,I2=J0,ω1=θ̇(1,2,3)=(-y,n,p),\qquad I_1=I_3=J_1,\quad I_2=J_0,\quad \omega_1=\dot\theta

  • Penjelasan

    Abaikan torsi gesek rotasional dan pertahankan hanya torsi akibat gaya angkat pada sumbu kemiringan. Jika gaya angkat bekerja pada pusat daerah terendam PP, torsinya adalah proyeksi hasil kali silang terhadap sumbu tersebut.

    Rumus

    N1=θ=(OP×Fl)e1,N2=N3=0N_1={\mathcal M}_\theta=(\overrightarrow{OP}\times\vec F_{\!l})\cdot\vec e_1,\qquad N_2=N_3=0

  • Penjelasan

    Persamaan Euler kedua, bersama I1=I3I_1=I_3, menyatakan bahwa putaran terhadap normal batu tetap.

    Rumus

    J0ω̇2=0ω2=Ω0J_0\dot\omega_2=0\quad\Longrightarrow\quad\omega_2=\Omega_0

  • Penjelasan

    Persamaan Euler ketiga menghubungkan kecepatan sudut tangensial dengan laju perubahan kemiringan.

    Rumus

    J1ω̇3=(J1J0)Ω0θ̇J_1\dot\omega_3=(J_1-J_0)\Omega_0\dot\theta

  • Penjelasan

    Integralkan dengan θ0=θ(0)\theta_0=\theta(0) dan asumsi ω3(0)=0\omega_3(0)=0, yakni tidak ada komponen sudut tangensial awal yang terpisah.

    Rumus

    ω3=J1J0J1Ω0(θθ0)\omega_3=\frac{J_1-J_0}{J_1}\Omega_0(\theta-\theta_0)

  • Penjelasan

    Masukkan hasil tersebut ke persamaan Euler pertama dan bagi dengan J1J_1. Dua faktor selisih inersia menghasilkan kuadrat positif.

    Rumus

    θ̈+(J0J1)2J12Ω02(θθ0)=θJ1\ddot\theta+\frac{(J_0-J_1)^2}{J_1^2}\Omega_0^2(\theta-\theta_0)=\frac{{\mathcal M}_\theta}{J_1}

  • Penjelasan

    Definisikan frekuensi giroskopik untuk memperoleh bentuk Persamaan (20) yang ringkas.

    Rumus

    ν=J0J1J1Ω0,θ̈+ν2(θθ0)=θJ1\nu=\frac{J_0-J_1}{J_1}\Omega_0,\qquad \ddot\theta+\nu^2(\theta-\theta_0)=\frac{{\mathcal M}_\theta}{J_1}

Simpulan

Bentuk ini mengasumsikan batu simetris, putaran normal Ω0\Omega_0 konstan, hanya torsi angkat pada sumbu kemiringan yang dipertahankan, dan ω3(0)=0\omega_3(0)=0. Jika komponen tangensial awal tidak nol, integrasi menghasilkan suku pemaksa konstan tambahan. Rumus yang tercetak di bab sumber tanpa Ω02\Omega_0^2 dan tanpa kuadrat rasio inersia bukan Persamaan (20) berdimensi konsisten dari artikel Bocquet.