Bab 11 Persamaan Diferensial Parsial - unsri.ac.id differensial parsial.pdf · Kita akan menurunkan persamaan itu sebagai model dari sistem fisik dan ... parsial linear jika persamaan
Post on 01-Feb-2018
309 Views
Preview:
Transcript
Bab 11
Persamaan Diferensial Parsial
Persamaan diferensial parsial dijumpai dalam kaitan dengan berbagai masalah fisik dan geometris bila
fungsi yang terlibat tergantung pada dua atau lebih peubah bebas. Tidak berlebihan jika dikatakan bahwa
hanya sistem fisik yang paling sederhana yang dapat dimodelkan dengan persamaan diferensial biasa
mekanika fluida dan mekanika padat, transfer panas, teori elektromagnetik dan berbagai bidang fisika
lainnya penuh dengan masalah-masalah yang harus dimodelkan dengan persamaan differensial parsial.
Yang sesungguhnya, kisaran penerapan persamaan diferensial parsial sangatlah besar, dibandingkan dengan
kisaran penerapan persamaan diferensial biasa. Peubah-peubah bebas dapat berupa waktu dan satu atau
lebih koordinat di dalam ruang. Bab ini akan ditujukan untuk beberapa persmaan diferensial parsial paling
penting yang dijumpai di dalam penerapan rekayasa. Kita akan menurunkan persamaan itu sebagai model
dari sistem fisik dan mengupas cara-cara untuk memecahkan masalah nilai awal dan masalah nilai batas,
dengan kata lain metode untuk memperoleh solusi bagi persamaan yang berkaitan dengan masalah fisik
yang dihadapi.
Di dalam pasal 11.1, kita akan mendefinisikan pengertian solusi persamaan diferensial Parsial.
Pasal 11.2 – 11.4 akan ditujukan untuk persamaan gelombang berdimensi-satu, yang mengatur gerak seutas
dawai yang bervibrasi. Persamaan panas akan dibahas dalam pasal 11.5 dan 11.6; persamaan gelombang
berdimensi-dua (membran bervibrasi) dalam pasal 11.7 – 11.10, dan persamaan Laplace dalam pasal 11.11
dan 11.12.
Di dalam pasal 11.13 dan 11.14 kita akan melihat bahwa persamaan diferensial parsial dapat juga
dipecahkan melalui transformasi Lapace (lihat Bab 5) atau transformasi Fourier (lihat pasal 10.10-10.12).
Metode numerik untuk persamaan diferensial parsial akan disajikan dalam pasal 20.4-20.7.
Prasyarat untuk Bab ini: persamaan diferensial biasa (Bab 2) dan deret Fourier (Bab 10).
Pasa-pasal yang dapat dilewati untuk kuliah yang lebih singkat : 11.6, 11.9, 11.10.
11.1 Konsep-konsep Dasar
Persamaan yang mengandung satu atau lebih turunan parsial suatu fungsi (yang
diketahui) dengan dua atau lebih peubah bebas dinamakan persamaan diferensial
parsial. Ordo turunan tertinggi dinamakan ordo persamaan tersebut. Seperti pada
persamaan diferensial parsial biasa, kita katakan bahwa suatu persamaan diferensial
parsial linear jika persamaan itu berderajat satu dalam peubah biasanya dan turunan
parsialnya. Jika setiap suku persamaan demikian ini mengandung peubah takbebasnya
atau salah satu dari turunannya, persamaan itu dikatakan homogen; bila tidak, persamaan
itu dikatakan tak homogen.
Teladan 1. Persamaan diferensial parsial linear ordo kedua yang
penting
2 22
2 2
22
2
2 2
2 2
2 2 2
2 2 2
persamaan gelombang dimensi-satu
persamaan panas berdimensi-satu
( , ) persamaan poisson berdimensi-dua
u uc
u x
u uc
t x
u uf x y
x y
u u u
x y z
0 persamaan Laplace berdimensi-tiga
Dalam hal ini c adalah konstanta, t adalah waktu, x y z adalah adalah koordinat Kartesius.
Persamaan 4 dengan ( f ≠ 0) adalah tak homogen, sedangkan persamaan-persamaan
lainnya homogen.
Yang dimaksud dengan solusi suatu persamaan diferensial pada suatu daerah R di
dalam ruang peubah (-peubah) bebasnya ialah fungsi yang memiliki turunan parsial yang
muncul di dalam persamaan itu, yang didefenisikan pada suatu domain mengandung R
dan yang memenuhi persamaan itu dimana-mana di dalam R. Ada kalanya orang hanya
menyaratkan bahwa fungsi tersebut kontinu pada batas daerah R, mempunyai turunan-
turunan tersebut di dalam interior daerah R, dan memenuhi persamaan itu di dalam
interior daerah R.)
Secara umum, keseluruhan solusi suatu persamaan diferensial adalah sangat besar.
Misalnya, fungsi-fungsi
(6) 2 2 2 2, cos . ln( ),xu x y u e y u x y
Yang berbeda sama sekali satu sama lain, semuanya merupakan solusi bagi (3), sebagai
pembaca dapat menverifikasi sendiri. Kita akan melihat nanti bahwa solusi tunggal suatu
persamaan diferensial parsial yang berasal dari suatu masalah fisik tertentu akan
diperoleh dengan memanfaatkan informasi tambahan dari dari situasi fisik tersebut.
Misalnya, sering kali nilai solusi yang diinginkan pada batas domainnya diketahui
(“syarat atau kondisi batas”); dalam kasus lain, bila t menyatakan waktu, nilai solusi
pada t = 0 adakalanya diberikan (”syarat awal”).
Kita tahu bahwa jika suatu persamaan diferensial biasa bersifat linear dan
homogen, maka dari solusi yang diketahui dapat diperoleh solusi-solusi lain melalui
superposisi. Pada kasus persamaan diferensial parsial linear, keadaannya sangat serupa.
Dan memang, teorema berikut ini berlaku.
Teorema Dasar 1 (Prinsip Superposisi)
Jika u1 dan u2 adalah solusi bagi suatu persamaan diferensial parsial homogen linear
pada sutu daerah, maka
1 1 2 2 ,u c u c u
dengan c1 dan c2 sembarang konstanta, juga merupakan solusi bagi persamaan itu dalam
daerah tersebut.
Bukti teorema penting ini mudah dan sangat mirip dengan bukti untuk Teorema 1 Pasal
2.1 dan disediakan untuk pembaca.
Di dalam pasal berikut, kita akan memulai pembahasan dengan persamaan
penting pertama yang dicantumkan pada Teladan 1, yaitu persamaan gelombang
berdimensi satu. Istilah ”berdimensi satu” menunjukkan bahwa persamaan itu
mengandung hanya satu peubah ruang, yaitu x. Persamaan ini mengatur gerak seutas
dawai elastik, misalnya dawai biola.
Soal-soal Latihan untuk Pasal 11.1
1. Buktikan Teorema Dasar 1 untuk persamaan diferensial ordo-kedua dalam
dua dan tiga peubah bebas.
2. Verifikasi bahwa fungsi-fungsi (6) merupakan solusi bagi ( 3).
Verifikasi bahwa fungsi-fungsi berikut merupakan solusi bagi persamaan
gelombang (1) untuk nilai konstan c di dalam (1) yang sesuai.
3. 2 2u x t 4.
2 29u x t 5. sin sinu t x6. cos sinu t t 7. cos sinu ct x 8. sin sinu ct x Verifikasi bahwa fungsi-fungsi berikut merupakan solusi bagi persamaan panas (2) untuk nilai konstanta c yang sesuai.
9. 9cos 10. sin 2 11. cost t tu e x u e x u e x
912. sin 3tu e x 13. 16 cos 2tu e x 14.
2 2
sinc tu e x
2 2 2 2
2
15. 16. 3 17. cos
18. sin 19. / 20. sin cosh
xu x y u x y u e y
u e y u arct y x u x y
21.Tunjukkan bahwa u(x, y) = v(x + c) + w(x - ct) merupakan solusi bagi persamaan gelombang (1); dalam hal ini v dan w terdiferensialkan dua kali. Verfikasi bahwa u(x, y) =
2 2lna x y b memenuhi persamaan laplce (3) dan tentukan a dan b sedemikian rupa
sehingga u memenuhi syarat batas u = 0 pada lingkaran 2 2 4.x y
22. Tunjukkan bahwa 2 2 21/u x y z merupakan suatu solusi bagi persamaan
Laplace (5).23. Tentukan potensial elektrostatik [soluis bagi (5)] antara dua bola konsentrik
S1: 2 2 2 2 2 2
21 dan : 4x y z S x y z bila S1 bertegangan 110 volt sedangkan S2
dihubungkan dengan tanah [0 volt]. Petunjuk. Ingat soal 23.
25. Tunjukkan bahwa 2 2/u y y x y memenuhi persamaan Laplace (3) dan bernilai
0 pada sumbu-x maupun pada lingkaran 2 2 1.x y
Hubungan dengan persamaan diferensial biasa
Jika suatu persamaan diferensial parsial mengandung turunan terhadap hanya salah satu
peubahnya saja, maka kita dapat memecahkannya seperti suatu persamaan diferensial
biasa, dengan memperlakukan peubah (-peubah) bebas lainnya sebagai prameter.
Pecahkan persamaan berikut jika u merupakn fungsi dua peubah x dan y.
Pecahkan persamaan berikut jika u merupakan fungsi dua peubah x dan y.
26. 0 27. 0 28. 0
29. 30. 0 31. 2 0
y x xx
xx z yy xx x
u u u u
u u u u u u u
Dengan memisalkan ux = p, pecahkan32. 0, 0 36. 0, 0 37. 0, 0
38. 0, 0
x y xx yy xx xy
zz yy
u u u u u u
u u
39. Perlihatkan jika kurva datar z = konstanta dari suatu permuklaan z = (x, y) berupa
garis-garis lurus yang sejajar sumbu-x, makla z merupakan solusi bagi persamaan
diferensial zx = 0. Berikan beberapa cotoh sebagai ilustrasi.
40. Tunjukkanbahw solusi , bagi 0x yz x y yz xz merepresentasikan
permukaan benda-benda putar. Berikan contoh ilustrasinya. Petunjuk. Ambil cos , sinx r y r dan kemudian tunjukkan bahwa persamaan itu menjadi 0z c .
11.2Pemodelan: Dawai Bervibrasi, Persamaan Gelombang
Berdimensi-Satu
Berbagai persamaan diferensial parsial penting yang pertama, Marilah kita turunkan
persamaan yang mengatur vibrasi pada seutas dawai elastis, yang diregangkan sampai
panjang L dan diikatkan pada kedua ujungnya. Misalkan kemudian dawai itu ditarik dan
diganggu dan kemudian dilepaskan pada t = 0 agar bergetar. Masalahnya adalah
menentukan vibrasi dawai tersebut, dengan kata lain menentukan penyimpangannya u(x,
t) pada sembarang titik x dan waktu t > 0; lihat gambarar 265.PQ
aT
1
T2
bxX+ΔxL
T2
T1
uPQ
ab
0
Gambar 265 Dawai yang bervibrasi
Bila menurunkan suatu persamaan diferensial yang bersumber pada suatu masalah
fiska tertentu, biasanya kita harus menyederhanakan asumsi-asumsi untuk menjamin agar
persamaan yang dihasilkan tidak menjadi terlalu rumit. Kenyataan penting telah kita
peroleh ketika kita mempelajari persamaan diferensial biasa, dan hal yang sama pun
berlaku untuk persamaan diferensial parsial.
Pada kasus kita sekarang ini, kita memberlakukan asumsi-asumsi berikut :
1. Massa dawai persamaan satuan panjang adaah konstan (”dawai homogen”)
Dawainya elastis sempurna dan tidak memberikan perlawanan terhadap
pelengkungan.
2. Tegangan yang disebabkan oleh peregangan dawai itu sebelum pengikatan kedua
ujungnya lebih besar diandingkan dengan gaya gravitasi, sehingga yang terakhir
ini dapat diabaikan.
3. Dawai itu mengalami gerak tranversal (melintang) kecil pada suatu bidang vertikal;
artinya, setip partikel dawai itu bergerak secara vertikal dengan defleksi
(penyimpangan) dan kemiringan di setiap titik dawai tetap kecil nilai mutlaknya.
Asumsi-asumsi itu sedemikian rupa sehingga sehingga kita dapat berharap bahwa solusi u(x, t) bagi persamaan diferensial yang diperoleh dapat menerangkan dengan cukup baik vibrasi kecil dawai ”nonideal” yang bermassa kecil dan homogen yang mengalami
tegangan besar.Untuk memperoleh persamaan diferensialnya, kita simak gaya-gaya yang bekerja pada suatu bagian kecil dawai tersebut (gambar 265). Karena dawai itu tidak memberi perlawanan terhadap pelengkungan (does not offer resistance to bending), maka tegangan bersifat tangensial terhadap kurva dawai itu pada setiap titik. Misalkan T1 dan T2 adalah tegangan di kedua titik ujng P dan Q dari bagian kecil tersebut. Karena tidak ada gerak dalam arah horisontal, berarti komponen horisontal tegangan ini pasti konstan. Dengan menggunakan notasi seperti pada gambar 265, kita memperoleh
(1) 1 2cos cosT T T konstanta
Pada arah vartikal terdapat dua gaya, yaitu komponen vertikal –T1 sin α dari T1 dan T2 sin
β dari T2; tanda minus disini berarti bahwa komponen yang bersangkutan di P arah ke
bawah. Menurut hukum kedua Newton, resultan kedua gaya tersebut sama dengan massa
bagian itu, px, kali percepatannya, 2u/ t2, dihitung pada suatu titik daerah x dan x + Δx;
dalam hal ini p adalah massa dawai yang tidak terdefleksi persatuan panjang, sedangkan
Δx adalah panjang bagian dawai yang tidak terdefleksi persatuan panjang, sedangkan Δx
adalah panjang bagian dawai yang tidak terdefleksi.
(2)
2
2 1 2sin sin
uT T x
t
Dengan menggunakan (1) kita memperoleh
22 1
22 1
sin sintan tan
cos cos
T T x u
T T T t
Sekarang tan α dan tan β adalah kemiringan kurva dawai itu di titk x dan x + Δx :
tan dan tanx x x
u u
x x
Di sini kita harus menuliskan turunan parsial sebab u juga tergantung pada t.Pembagian (2) dengan Δx menghasilkan
2
2
1
x x x
u u u
x x x T t
Jika Δx mendakati nol, kita memperoleh persamaan diferensial parsial linear
2 22
2 2
u uc
t x
2 Tc
.
Inilah yang dinamakan Persamaan gelombang berdimensi-satu, yang mengatur
masalah kita. Kita lihat bahwa persamaan homogen dan berorde dua. Notasi c2 (alih-alih c
untuk konstanta fisik T/p diambil unutuk menunjukkan bahwa konstanta positif. Solusi
persamaan ini akan diperoleh dalam pasal berikut.
11.3 Metode Pemisahan Peubah (Metode Hasil Kali)
Di dalam pasal sebelumnya kta telah menunjukkan bahwa vibrasi seutas dawai, misalnya
dawai bola, mengikuti persamaan gelombang berdimensi-satu
(1)
2 22
2 2
u uc
t x
,
dengan u(x ,t) adalah defleksi dawai tersebut. Untuk mengetahui bagaimana dawai itu bergerak, kita harus memecahkan persamaaan ini; lebih tepatnya, kita tentukan solusi u bagi (1) yang juga memenuhi syarat yang dikenakan oleh sistem fisik bersangkutan. Karena dawai itu diikat pada kedua ujung x = 0 dan x = L, kita kenai dua syarat batas
(2) (0, ) 0, ( , ) 0 untuk semua u t u L t t
Bentuk gerak dawai itu akan tergantung pada defleksi awalnya (defleksi pada t =0) dan
pada kecepatan awalnya (kecepatan pada t = 0). Jika f(x) adalah defleksi awalnya dan
g(x) adalah kecepatan awalnya, maka kita akan memperoleh kedua syarat awal
( , 0) ( )u x f x
(3)
dan
0
( ).t
ug x
t
(4)
Masalah kita sekarang adalah menemukan solusi bagi (1) yang memenuhi syarat (2)-(4).
Langkah-langkah yang akan kita tempuh adalah sebagai berikut.
Langkah Pertama. Dengan menerapkan apa yang dinamakan metode pemisahan
peubah atau metode hasil kali, kita akan memperoleh dua persamaan diferensial
biasa.
Langkah Kedua. Kita akan menetukan solusi bagi kedua persamaan itu yang
memenuhi syarat batas yang ditentukan.
Langkah Ketiga. Solusi-solusi itu akan digabungkan sedemikian rupa sehingga
hasilnya merupakan solusi persamaan gelombang (1) memenuhi syarat awal yang
diberikan.
Rinciannya adalah sebagai berikut.
Langkah Pertama. Dua Persamaan Diferensial Biasa
Metode hasil kali ini menghsailkan solusi bagi persamaan gelombang (1) yang berbentuk
( , ) ( ) ( )u x t F x G t
(5)
yang merupakan hasil kali dua fungsi, masing-masing tergantung pada salah satu peubah
x atau t. Akan kita lihat nanti bahwa metode ini mempunyai banyak penerapan di dalam
matematika rekayasa. Dengan mendiferensialkan (5), kita memperoleh
2 2''
2 2 dan
u uF G F G
t x
Dengan tanda dua titik di atas melambangkan turunan terhadap t, sedangkan tanda aksen
melambangkan turunan terhadap x. Dengan menyisipkan ini ke dalam persamaan
diferensial (1), kita memperoleh
2 ''F G c F G
Dengan membagi dengan c2FG kita akan memperoleh
''
2
G F
c G F
.
Ruas kiri hanya mengandung fungsi yang tergantung hanya pada t sedangkan ruas kanan mengandung fungsi yang tergantung hanya pada x. Ini berarti ruas kanan maupun ruas kiri sama dengan suatu konstanta, maka pengubahan t akan mengubah nilai ruas ini namun jelas tidak mengubah nilai ruas kanan, sebab ruas kanan tidak tergantung pada t begitu pula, jika ruas kanan tidak sama dengan konstanta, maka pengubahan x akan mengubah nilai ruas ini namun jelas tidak mengubah nilai ruas kiri. Jadi,
''
2
G Fk
c G F
.
Ini menghasilkan dua persamaan diferensial biasa linier, yaitu :
(6) '' 0F kF
dan
(7) 2 0G c kG
Di sini, k masih merupakan sembarang bilangan.
Langkah Kedua. Memenuhi Syarat Batas
Kita sekarang akan menentukan solusi Fungsi dan G bagi (6) dan (7) sedemikian rupa
sehingga u = FG memenuhi syarat batas (2); dengan kata lain
(0, ) (0) ( ) 0, ( , ) ( ) ( ) 0u t F G t u L t F L G t
untuk semua t.
Jelaslah, jika G ≡ 0, maka u ≡ 0, sama sekali tidak menarik. Jadi, kita meginginkan G ≡
0, sehingga haruslah
(8) (a) (0) 0 (b) ( ) 0F F L
Untuk k = 0, solusi umum bagi (6) adalah F = ax + b, sehingga dari (8) kita
memperoleh a = b = 0. Dengan demikian F ≡ 0, juga tidak menarik sebab akibatnya G ≡
0. Untuk k = m2 positif, solusi umum bagi (5) adalah
,x xF Ae Be
dan dari (8) kita memperoleh F ≡ 0, seperti sebelumnya. Sekarang, kita tinggal
kemungkinan mengambil k negatif k = -p2 . persamaan (6) sekarang mempunyai bentuk
2" 0.F p F
dengan solusi umum adalah
( ) cos sin .F x A px B px
dari sini dan (8) kita memperoleh
(0) 0 sehingga ( ) sin 0F A F L B pL
Kita harus mengambil B 0 sebab jika tidak demikian F ≡ 0. ini berarti sin pL = 0.
Akibatnya ,
(9) atau ( bilangan bulat).
npL n p n
L
Dengan mengambil B = 1, kita memperoleh tak hingga banyaknya solusi F(x) = Fn(x),
dengan
(10)( ) sin n = 1,2,...n
nF x x
L
Solusi-solusi itu memenuhi (8). Untuk bilangan bulat negatif n kita memperoleh
solusi yang pada hakekatnya sama, kecuali tandanya, sebab sin (-α) =
- sin α.
Bilangan k sekarang terbatas pada nilai-nilai k = -p2 = -(np/L)2, yang berasal dari
(9). Untuk nilai-nilai k tersebut, persamaan (7) mempunyai bentuk
2 0 dengan n n
cnG G
L
Suatu solusi umumnya adalah
*( ) cos sin n n n n nG t B t B t .
Jadi, fungsi-fungsi un(x, t) = Fn(x)Gn(t), yakni
(11)
*( , ) ( cos sin )sin ( =1,2,...)n n n n n
nu x t B t B t x n
L
,
Merupakan solusi-solusi bagi (1) yang memenuhi syarat batas (2). Fungsi-fungsi itu
dinamakan fungsi-eigen, atau fungsi ciri, dan nilai λn = cnp/L dinamakan nilai-eigen,
atau nilai ciri, dari dawai yang bervibrasi tersebut. Himpunan (λ1, λ2, ...)
dinamalan spektrum.
Kita lihat bahwa setiap un merepresentasikan suatu gerak harmonik
yang memiliki frekuensi λn /2 λn = cn/2L siklus per satuan waktu. Gerak ini
dinamakan modus normal ke-n dawai tersebut. Modus normal pertama dikenal sebagai
modus dasar (fundemental mode) (n=1), sedangkan yang lain dikenal sebagai overtone,
pada musik, modus-modus tersebut menghasikan oktaf, oktaf plus seperlima, dan
seterusnya. Karena di dalam (11)
2 1sin 0 pada , ,...,
n L L nx x L
L n n n
,
Modus normal ke-n memiliki n-1 apa yang dinamakan simpul (nodes), yaitu pada dawaii yang tidak bergerak (Gambar 266).0L0L0L0L
N = 1N = 2N = 3N = 4
Gambar 266. Modus normal dari dawai yang bervibrasi
Gambar 267 memperlihatkan modus normal kedua untuk berbagai nilai t. Pada tiap saat
dawai itu berbentuk gelombang sinus. Bila bagian kiri dawai itu bergerak dibawah, maka
bagian yang bergerak ke atas, dan sebaliknya. Untuk modus yang lain, keadaannya
serupa. L
x
Gambar 267. Modus normal kedua untuk berbagai nilai t
Langkah Ketiga. Solusi Masalah Keseluruhannya
Jelaslah, satu solusi un(x, t) pada umunya tidak memenuhi syarat awal (3) dan (4).
Sekarang, karena persamaan (1) linier dan hmogen, maka menurut Teorema Dasai 1
Pada pasal 11.1 jumlah terhingga banyaknya solusi un juga merupakan solusi bagi (1).
Untuk memperoleh solusi yang memenuhi (3) dan (4), kita simak
(12)
*
1
( , ) ( cos sin )sin n n n n n
n
nu x t B t B t x
L
Dengan λn = cnp/L seperti sebelumnya. Dari sini dan syarat awal (3), kita memperoleh
(13) 1
( , 0) sin ( ). n
n
nu x x f x
L
Jadi, agar (12) memenuhi (3), koefisien-koefisien Bn harus diambil sedemikian rupa
sehingga u(x, 0) menjadi uraian separuh-kisaran (half-range expansion) bagi f(x); yaitu
deret sinus Fourier bagi f(x). Jadi [lihat (4) pasal 10.5],
(14) 0
2( )sin n=1,2,...
L
n
nB f x dx
L L
Begitu pula, dengan mendiferensialkan (12) terhadap t dengan menggunakan syarat awal
kedua, yaitu (4), kita memperoleh
*
01 0
( sin cos )sinn n n n ntn t
u n xB t B t
t L
*
1
sin ( ).n nn
n xB g x
L
Oleh karena itu, agar (12) memenuhi (4), koefisien-koefisien Bn* harus dipilih sedemikian
rupa sehingga, untuk t=0, turunan u/ t menjadi deret sinus Fourier bagi g(x); jadi menurut (4) pasal 10.5.
*
0
2( )sin
L
n n
nB g x dx
L L
Atau, karena λn = cnp/L,
(15)
*
0
2( )sin
L
n n
nB g x dx
cn L
n = 1, 2, ...
Akibatnya, u(x, t) yang diberkan leh (12) dengan koefisien-koefisien (14) dan (15)
merupakan salah satu solusi bagi (1) yang memenuhi syarat (2)-(4), asalkan deret (12)
tersebut konvergen dan juga asalkan deret yang diperoleh meleui pen diferensialan (12)
suku-demensi-suku terhadap x dan i konvergen dan mempunyai jumlah bertut-turut, yang
bersifat kontinu.
Jadi, solusi (12) masih berupa ekspresi formal belaka, kita masih harus
membuktikannya. Untuk kemudahan, kita hanya akan menyimak kasusu bila kecepatan
awalnya g(x) sama dengan nol. Ini berimplikasi Bn* sama dengan nol dan (12) tereduksi
menjadi
(16) 1
( , ) cos sin , = . n n nn
n cnu x t B t
L L
Sekarang kita akan mejumlahkan deret ini, artinya menuliskan hasilnya dalam bentuk tertutup atau terhingga. Untuk maksud ini, kita gunakan rumus [lihat (11dalam apendiks 3.1]
1cos sin x = sin ( ) sin ( ) .
2
cn n n nt x ct x ct
L L L L
Akibatnya, kita dapat menuliskan (16) dalam bentuk
1 1
1 1( , ) sin ( ) sin ( ) .
2 2n nn n
n nu x t B x ct B x ct
L L
Kedua deret itu yang diperoleh melalui subsitusi, berturut-turut, x - ct dan x + ct peubah x
dalm deret sinus Fourier (13) untuk f(x). Oleh karena itu,
(17) * *1
( , )2
u x t f x ct f x ct
Dengan f* adalah perluasan periodik ganjil fungsi f dengan periode 2L (Gambar 268).
Karena defleksi awal f(x) kontinu pada selang 0 x L dan bernilai nol di kedua titik
ujngnya, maka dari (17) kita dapat menyimpulkan bahwa u(x, t) merupakan suatu fungsi
kontinu untuk semua niai x dan t. Dengan mendiferensialkan (17) kita lihat bahwa u(x, t)
merupakan salah satu solusi bagi (1), asalkan f (x) dideferensialkan dua kali pada selang 0
< x < L, dan memiliki turunan kedua satu-arah (one sided second derivatives) di x = 0,
yang bernilai nol. Di bawah kondisi tersebut, terbuktilah bahwa u(x, t) merupakan salah
satu bagi (1) yang memenuhi 1-(4)
0Lx
Gambar 268. perluasan periodik ganjil bagi f(x)
Jika f’(x) dan f”(x) hanya kontinue sepotong-sepotong (Lihat Pasal 1.5), atau jika
turunan-turunan satu arah itu tidak nol, maka untuk setiap t akan ada terhingga
banyaknya nilai x yang pada nilai-nilai itu turunan kedua fungsi u yang muncul di dalam
(1) tidak ada. Kecuali pada titik-titik tersebut, persamaan gelombang masih terpenuhi,
dan kita dapat memandang u(x, t) sebagai solusi bagi masalah kita dalam pengertian yang
lebih luas. Misalnya, dalam kasus defleksi awal yang berbentuk segi tiga (teladan 1 di
bawah ini) menghasilkan solusi jenis ini.
Kiranya perlu dikemukakan suatu tafsiran fisik yang sangat menarik bagi (17).
Grafik f*(x + ct) diperoleh dari grafik f*(x) dengan cara menggeserkan yang terakhir ini
ct satuan ke kanan (gambar 269). Ini berarti bahwa f*(x – ct ) (c > 0) merepresentasikan
suatu gelombang yang bergerak ke kanan sejalan dengan naiknya t. Begitu pula, f* (x +
ct) merepresentasikan suatu gelombang yang bergerak ke kiri, dan u(x, t) merupakan
superposisi kedua gelombang tersebut.
ctx
f *(x)f *(x - ct)
Gambar 269.Tafsiran bagi (17)
Teladan 1. Dawai bervibrasi jika defleksi awalnya berupa segi tiga
Tentukan solusi persamaan gelombang (1) jika defleksi awalnya berbentuk segi tiga
2
2( )
02
2
( )K
xL
KL x
L
Ljika x
Ljika x L
f x
dan kecepatan awalnya nol. (gambar 270 memperlihatkan f(x) = u(x, 0) di puncaknya.)
jawab. Karena g(x) ≡ 0, kita memperoleh Bn* = 0 di dalam (12), dan dari Teladan 1
dalam pasal 10.5 kita melihat bahwa Bn ditentukan oleh (5) Pasal 10.5. Jadi, (12)
mempunyai bentuk
2 2 2
8 1 1 3 3( , ) sin cos sin cos ..... .
1 3
k c cu x t x t x t
L L L L
Untuk memplotkan grafik solusi itu, kita dapat memanfaatkan u(x, 0) = f(x) dan tafsiran kedua fungsi di dalam representasi (17). Ini menghasilkan grafik dalam gambar 270.
Sangat menarik bahwa solusi (17) juga dapat diperoleh dengan cepat melalui transformasi yang sesuai terhadap persamaan gelombang itu, mengikuti gagasan cerdik d’Alembert yang kita akan bahas dalam pasal berikutnya. Soal-soal Latihan untuk Pasal 11.3
1. Bagaimanakah pengaruh pelipatduaan tegangan terhadap titik nada (pitch) hada dasar
(fundamental tone) dawai tersebut ?
2. Tentukan defleksi , u x t dawai yang bervibrasi (panjang L = p, kedua ujungnya
terikat , dan c2 = T/p = 1) jika kecepatan awalnya nol dan defleksi awalnya adalah:
4 4
2 2
2. 001 sin 3. sin 2 4. sin sin 3
1 15. 0.01 6. 0.01 7. k
2 2
x k x k x x
x x x x x
8. 9. 10.
Gambar 270. Solusi ,u x t
dalam Teladan 1 untuk bebagai nilai t (bagan kanan gambar) yang diperoleh melalui superposisi gelombang yang bergerak ke kanan (garis putus-putus) dengan gelombang yang bergerak ke kiri (bagian kiri gambar).11. Berapakah nisba atau rasioa antara amplitudo modus dasar dengan amplitudo
overtune kedua didalam Soal 5 ? Nisbah 2 21 1 2/ ...a a a
? Petunjuk. Gunakan (7) dalam Pasal 10.8 dengan tanda sama dengan.
12. Bagaimanakah ketergantungan frekuensi modus dasar dawai yang bervariasi pada
panjang dawai, tegangan, dan massa dawai persamaan satuan panjang ?
Tentukan defleksi u(k, t) dawai yang bervibrasi (Panjang L = p, kedua ujungnya terikat,
c2 = 1) jika defleksi awalnya f(x) dan kecepatan awalnya g(x) adalah :
13. f = 0, g(x) = 0.1 sin 2x
14. f = 0, g(x) = 0.01x jika 1 1
0 , 0.01 jika .2 2
x g x x x
15. f(x) = 0,1 sin x, g(x) = -0,2 sin xDengan memisahkan peubah-peubahnya tentukan solusi u(x, y) bagi persamaan-
persamaan berikut.
16. 17. 0 18.
19. 2 =0 20. 21.
x y x y x y
x y x y x y
u yu u u u u
u u x y ayu bxu xu yu
2 2
22. 3 0 23. 24. 0
xy xy xx yyx u y u u u u u
25. Tunjukkan bahwa suatu masalah (1) – (4) dengan syarat batas lebih rumit, katakanlah
0, 0, , , u t u L t b t dapat direduksi menjadi suatu masalah pencarian fungsi baru
v yang memenuhi syrat 0, , 0,v t v L t 1, 0 ,v x f x 1, 0tv x g x
, namun suatu persamaan gelombang yang tak homogen. Petunjuk. Ambil u = v + w dan tentukan w yang cocok.
11.4 Solusi D’Alembert bagi Persamaan Gelombang
Menarik untuk dicatat bahwa solusi (17) Pasal 11.3 bagi persamaan gelombang
(1)
2 22 2
2 2
u u Tc c
t x
dapat diperoleh dengan cepat mentransformasikan (1) secara cerdik, yaitu dengan
memperkenalkan dua peubah bebas baru1
(2) , - .v x ct z x ct
Dengan demikian u menjadi fungsi dari v dan z, dan turunan-turunan di dalam (1) dapat
diucapkan dalam turunan terhadap v dan z melalui penerapan kaida rantai dalam Pasal
8.7. Dengan melambangkan turunan parsial dengan subskrip, kita lihat dari (2) bahwa vx =
1 dan zx = 1 . Untuk kemudahan marilah lambangkan u(x, t), sebagai suatu fungsi dari v
dan z, dengan huruf yang sama u . dengan demikian
x v z z x v zu v v u z u u
Dengan menerapkan kaidah rantai terhadap ruas kanan dan memanfaatkan kenyataan
bahwa vx = 1 zx = 1 , kita memperoleh
2x x v z x v z x vv vz zzv zu u u v u u z u u u
Kita transformasikan turunan satu lagi di dalam (1) dengan cara yang sama , kita peroleh
2 2 .tt vv vz zzu c u u u
Dengan menyisipkan kedua hasil itu kedalam (1) kita memperoleh (lihat catatan kaki 1
dalam apendiks 3.1)
(3)
2
0.vz
uu
z v
jelaslah, maksud pendekatan ini adalah agar (3) yang dihasilkannya dapat segera dipecahkan dengan mengintegralkan dua kali berturut-turut. Dengan mengintegralkan terhadap z , kita memperoleh1perlu dikemukakan bahwa teori umum persaman diferensial parsial memberikan suatu cara bersistem untuk menemukan transformasi ini, yang dapat menyderhanakan persamaannya. Lihat acuan [C 14] dalam Apendiks 1
uh v
v
dengan h(v) adalah suatu fungsi dari v . pengintegralan terhadap v menghasilkan
u h v dv z
dengan y(z) adalah suatu fungsi dari z. Karena integral itu merupakan fungsi dari v,
gunakanlah (v), maka solusi u berbentuk u = (v) + y(z). Karena (2), maka
(4) u (x, t) = f(x + ct) + y(x - ct).
Ini dikenal sebagai solusi d’Alembert2 bagi persamaan gelombang (1).
Fungsi dan y dapat ditentukan dari kondisi awal. Mari kita ilustrasikan ini
untuk kasus kecepatan awal nol dan defleksi awal, u(x, 0) = f(x). Dengan
mendifrensialkan (4) kita memperoleh
uc
t
f’ 'x ct c x ct
Dengan tanda aksen melambang turunan masing-masing terhadap keseluruhan argumen x
+ ct dan x – ct. Dari (4) dan (5) dan syarat awal, kita memperoleh
( , 0) ( ) ( ) ( )
( , 0) '( ) '( ) 0t
u x x x f x
u x c x c x
Dari persamaan terakhir, kita peroleh y’ = ’. Oleh karena itu, y = + k, dan dari sini
dan persamaan pertama, kita peroleh 2 + k = f = (f - k)/2.Dengan fungsi –fungsi
dan y tersebut, solusi (4) menjadi
u(x, t) =
1
2 f [(x + ct) + f(x – ct) ],sesuai dengan (17) dalam Pasal 11.3. Pembaca dapat mennjukkan bahwa karena syarat
batas (2) dala pasal tersebut, fungsi f haruslah ganjil dan berperiode 2L, jika kecepatan
awal g(x) tidak identik dengan nol, alih-alih (6) kita akan memperoleh
(7)
1 1( , ) [ ( ) ( )] ( ) .
2 2
x ct
x ct
u x t f x ct f x ct g s dsc
Hasil kita ini menunjukkan bahwa syarat awal dan syarat batas menentukan solsinya
secara tunggal.
Solusi persamaan gelombang mealui transformasi Laplace dan transformasi linier
Solusi persamaan gelombang melalui transformasi Laplace dan transformasi Fourier akan
dibahas dalam Pasal 11.13 dan 11.14.
Di dalam pasal berikut, kita akan menengok pada persamaan panas, yang merupakan
persamaan diferensial sangat penting lainnya.
Soal-soal Latihan untuk Pasal 11.4
Dengan menggunakan (6), buatlah gambar (sejenis Gambar 2.70 dalam Pasal 11.3)
defleksi u(x, t) dari waktu bervibrasi (panjang L = 1, kedua ujungnya terikat) jika
kecepatan awalnya nol dan defleksi awalnya f(x), k sebuah bilangan kecil, misalkan k =
0.01.
3
2 4 2
1. sin 2. 1 3.
4. 5. 1 cos2 6. sin
f x k x f x kx x f x k x x
f x k x x f x k x f x k x
Dengan menggunakan transformasi yang diberikan, pecahkan persamaan diferensial berikut.
7. ,
8. ,
9. 2 0 ,
10. , z 2
11. 4 0 , 3
xy xx
xy
xy
xy xx
x
xx yy
yy yy
xx
u u v y z x y
yu xu u v y z xy
u u u v y z x y
u u v y x y x
u u v x y z x y
Jenis dan bentuk-bentuk normal persamaan diferensial parsial linear. Suatu persamaan yang berbentuk
(8) 2 , , , xx xy yy x yAu Bu Cu F x y u u
dikatakan eliptik jika2 0,AC B hiperbolik jika
2 0.AC B Di sini A, B, C, mungkin merupakan fungsi dari x dan y, dan jenis (8) mungkin berbeda di bagian-bagian bidang –xy yang berbeda.
13. Tunjukkan bahwa
Persamaan Laplace 0 adalah eliptik,xx yyu y
persamaan panas 2 adalah parbolik,t xxu c u
persamaan gelombang 2 adalah hiprbolik,tt xxu c u
persamaan Tricomi 0 xx yyyu u
adalah berjenis campuran (eliptik di separuh-bidang bagian atas dan hiperbolik di separuh bagian bawah).14. Jika persaman (8) hiperbolik, maka persaman itu dapat ditarnsformasikan ke dalam
bentuk nomal uvz = , , , , v zF v z u u u
dengan mengambil , , , ;v x y z x y
dalam hal ini konstanta dan konstanta merupakan solusi y y x bagi
persamaan 2' 2 ' 0Ay By C (lihat acuan [C12] ). Tunjukkan bahwa dalam kasus
persamaan gelombang (1)., .x ct x ct
15. jika (8) parabolik, subsitusi , , ,v x z x y dengan didefinisikan seperti pada
soal (14), mengubahnya ke bentuk normal , , , , .vv v zu F v z u u u Verifikasi fakta ini
untuk persamaan 2 0xx xyu u
. 16. (Persamaan Airy) Tunjukkan bahwa melalui pemisahan peubah kita dapat
memperoleh persamaan Airy " 0G yG dari persamaan Tricomi.(untuk solusinya, lihat hlm. 446 Acuan [1] yang dicantumkan dalam Apendiks 1, lihat juga Soal Ulangan 30 Bab 4.)
uyx
x=L
Gambar 271. Balok yang tidak mengalami perubahan bentuk dalam Soal 17
Vibrasi sebuah balok. Dapat diperlihatkan bahwa vibrasi vertikal bebas dari suatu balok
seragam (Gambar 271) mengikuti persamaan orde-empat
(9)
2 42
2 40
u uc
t x
(Acuan [C14]).
Dengan dalam hal ini 2 /c EI pA (E = modulus elastisitas Young, I = momen inersia
irisan melintang terhadap sumbu –y di dalam gambar, p = kerapatan atau densitas, A = luas irisan melintang ).
17. Dengan mensubsitusikan u F x G y ke dalam (9) dan memisahkan peubah-
peubahnya, tunjukkan bahwa
(4) 2 4
2 2
/ / konstanta.
cos sin cos sin ,
cos sin .
F F G c G
F t A x bB x C x D x
G t a c t b c t
18. Tentukan solusi n n nu F x G t bagi (9) jika kecepatan awalnya nol dan memenuhi
syarat-syarat batas (lihat Gambar 272)
0, 0, , 0u t u L t (kedua ujungnya diganjal untuk semua t),
0, 0, , 0xx xxu t u L t (di kedua ujungnya, momennya nol, sehingga
lengkungannya juga nol).
ux=L
x
Gambar 272 . Balok dalam Soal 18
19. Tentukan solusi bagi (9) yang memenuhi syarat-syarat dalam Soal 18 dan syarat awal , 0 .u x f x x L x
20. Bandingkan hasil-hasil dari Soal 19 dengan Soal 5 Pasal 11.3.
Apakah perbedaan pokok antara frekuensi modus normal dawai yang bervibrasi
dengan frekuensi modus normal balok yang bervibrasi?
21. Apakah syarat-syarat batasnya jika balok itu ditahan/dijepit di kedua ujungnya ? (lihat
Gambar 273.)
Gambar 273. Balok dalam Soal 21
Gambar 274. Balok dalam Soal 21
22. Tunjukkan bahwa F x dalam Soal 17 memenuhu syarat-syarat dalam Soal 21 jika
L adalah akar persamaan
(10) cosh cos 1L L
23. Tentukan solusi hampiran bagi (10).
24. Jika balok itu dijepit di ujung kirinya dan dibiarkan bebas di ujung kanannya
(Gambar 274), maka syarat-syarat batasnya adalah
0, 0, 0, 0, , 0, , 0.x xx xxxu t u t u L t u L t
Tunjukkan F(x) dalam Soal 17 memenuhi syarat-syarat itu jika L adalah akar persamaan
(11) cosh cos 1L L .
25. Tentukan soluisi hampiran bagi (11).
11.5 Aliran Panas
Aliran panas di dalam suatu benda homogen mengikuti persamaan panas (lihat pasal 9.7)
2 2uc u
t
2 kc
Dengan u(x, y, z, t) adalah suhu di dalam benda tersebut, K adalah konduktivias termal, s
adalah panas jenis dan r adalah kerapatan benda, v2 u adalah Laplacian dari u, dan relatif
terhadap koordinat Kartesius x, y, z
22 2 2
2 2 2.
u u uu
x y z
Sebagai salah satu penerapan penting, marilah kita simak suhu pada suatu batang
atau kawat tipis panjang, yang irisan melintangnya konstan dan terbuat dari bahan yang
homogen, dan terletak pada sumbu-x. Maka u tergantung hanya pada x dan waktu t dan
persamaan panasnyan menjadi apa yang dinamakan persamaan panas berdemensi-satu.
22
2
u uc
t x
xL
Gambar 275. Batang yang dibahas
Kita akan memecahkan (1) untuk beberapa jenis syarat batas dan syarat awal yang
penting. Prosedur penyeleseiannya akan mirip dengan pada kasus persamaan gelombang.
Perilaku solusinya akan sama sekali berbeda dibandingkan dengan perilaku persamaan
gelombang, sebab (1) melibatkan uu. (Klasifikasi dalam soal 13 pada pasal sebelumnya
bukan cuma untuk formalitas belaka, namun mempunyai konsekuensi yang dalam relatif
terhadap perilaku umum solusinya.)
Marilah kita mulai dengan kasus kedua ujung batangnya (x = 0 dan x = L )
dipertahankan pada suhu nol. Maka syarat-syarat batasnya adalah
(0, ) 0, ( , ) 0u t u L t untuk semua t.
Perhatikan bahwa ini mempunyai bentuk yang sama dengan (2) Pasal 11.3. Jika f(x)
adalah suhu awal batang tersebut, maka syarat awalnya adalah
( , ) ( ) [ ( ) diketahui].u x t f x f x
Kita akan menentukan solusi u(x, t) bagi (1) yang memenuhi (2) dan (3).
Langka pertama. Dengan menerapkan metode pemisahan peubah, mula-mula kita
tentukan solusi bagi (1) yang memenuhi syarat batas (2). Kita mulai dengan
( , ) ( ) ( )u x t F x G t .Pensubstitusian ini ke dalam (1) menghasilkan persamaan FG = c2
( , ) ( ) [ ( ) dikrtahui]u x t f x f x F”G, dengan tanda petik di atas melambangkan turunan terhadap t, sedangkan aksen menandakan turunan terhadap x. Untuk memisahkan peubah, kita bagi persamaan ini dengan c2F”
"
2.
G F
c G F
g
Ruas kiri hanya tergantung pada t, sedangkan ruas kanan hanya tergantung pada x.
Seperti pada pasal 11.3, kita simpulkan bahwa kedua ruas itu pasti sama dengan suatu
konstanta, katakanlah k. Pembaca dapat memperlihatkan bahwa untuk k 0, satu-satunya
solusi u = -p2 yang negatif, kita memperoleh dari (5)
"2
2. .
G Fp
c G F
g
Kita lihat bahwa ini menghasilkan dua buah persamaan diferensial biasa
(6)" 2 0F p F
(7) 2 2 0G c p G
Langkah kedua. Kita perhatikan (6). Solusi umumnya adalah
(8) ( ) cos sin .F x A px B px Syarat batas atas (2) berakibat bahwa
(0, ) (0) ( ) 0 dan ( , ) ( ) ( ) 0U t f G t u L t F L G t Karena G ≡ 0 berimplikasi u ≡ 0, maka kita menginginkan F(0) = 0 dan F(L) = 0.
Menurut (8), F(0) = A. jadi, A = 0, sehingga
F(L) = b sin pL
Haruslah b 0, sebab jika tidak F ≡ 0. oleh karena itu syarat F(L) = 0 membawa pada
sin 0, sehingga , 1, 2,...n
pL p nL
Dengan mengambil B = 1, kita memperoleh solusi bagi (6) yang memenuhi (2):
( ) sin , = 1,2,...n
n xF x n
L
(Seperti pada Pasal 11.3, kita tidak perlu menyimak nilai-nilai n yang bulat negatif)
Sekarang kita simak persamaan diferensial (7). Untuk nilai-nilai p = np/L yang baru
diperoleh, persamaan ini mempunyai bentuk
2 0 dengan .n n
cnG G
L
g
Solusi umumnya adalah 2
( ) , = 1,2,...,n tn nG t B e n
dengan Bn adalah konstanta. Jadi, fungsi –fungsi
2
( , ) ( ) ( ) sin = 1,2,...,n tn n n n
n xu x t F x G t B e n
L
Merupakan solusi bagi persamaan panas (1) yang memenuhi (2).
Seperti halnya pada pasal sebelumnya, untuk memperoleh solusi bagi masalah
perubahannya, kita akan menggunakan deret tak hingga.
Langkah ketiga. Untuk memperoleh solusi yang juga memenuhi (3), kita perhatikan
(10)
2
n=1 1
( , ) ( , ) sin , .n tn n n n
n
n x cnu x t u x t B e
L L
Disini dan (3), maka
1
( ,0) sin ( ).n nn
n xu x B f x
L
Agar (10) memenuhi (3), koefisien-koefisien Bn harus dipilih sedemikian rupah hingga
u(x, 0) merupakan uraian separuh-kisaran bagi f(x), yaitu deret sinus bagi f(x); dengan
kata lain [lihat (4) Pasal 10.5)],
(11) 0
2 ( )sin 1, 2,..., L
L
n
n xB f x dx n
L
3ihat catatan kaki 7 di dalam Pasal 10.2.
Solusi bagi masalah kita ini dapat dibuktikan dengan mengasumsikan bahwa f(x) itu
sepotong-sepotong pada selang 0 x L (lihat pasal 5.1), dan memiliki satu arah3 di
semua ttik interior selang tersebut; artinya, di bawah asumsi-asumsi tersebut, deret (10)
dengan koefisien-koefisien (11) merupakan solusi bagi masalah fisik yang kita hadapi.
Buktinya, yang membutuhkan pengetahuan tentang konvergenan seragam, akan diberikan
pada kesempatan nanti (soal 21, 22 di akhi pasal 14.8).
Karena adanya faktor eksponensial, semua suku di dalam (10) mendekati nol jika
t mendekati tak hingga. Laju penurunan suhu semakin cepat sejalan dengan naiknya n.
Teladan 1. Suhu awal sinusiodal
Misalkan suhu u(x, t) di dalam sebatang tembaga yang telah diisolasi yang panjang 80
cm suhu awalnya adalah 100 sin (px/80)o C dan ujung-ujungnya dipertahankan pada suhu
0oC berapa lama sampai suhu maksimum di dalam batang tembaga itu turun menjadi
50oC? Mula-mula kira-kiralah, dan kemudian hitunglah. Data fisik untuk tembaga:
kerapatan 8.92 gm/cm3, panas jenis 0,092 kal/gm oC konduktifitas termal 0.95 kal/cm det
oC.
Jawab. Syarat awal menghasilkan
1
( , 0) sin ( ) 100sin . 80 80n
n
n x xu x B f x
Dengan demikian, dengan melihat begitu saja atau dari (10) kita memperoleh B1 = 100,
B2 = B3 = ... = 0. Di dalam (10) kita membutuhkan λ12 = c2p2/ L2, dengan c2 = K/sp
= 0.95/(0.092)(8.92) = 1.158 [cm2/det]. Jadi, kita memperoleh l12 =
(1.158)(9.870)/6400 = 0.001785 [det-1]. Solusi bagi (10) adalah 0.0017853
( , ) 100sin .80
txu x t e
Selanjutnya, 100 -0.001785t = 50 jika t = (ln 0.5)/(-0.001785) = 388 [detik] 6.5 [menit].
Teladan 2. Kecepatan turunnya suhu
Pecahkan soal pada Teladan jika suhu awalnya adalah 100 sin (3px/80) oC dan data
lainnya sama
Jawab. Di dalam Pasal (10), alih-alih n = 1 sekarang n =3, dan l32 = 32l1
2 = 9(0.001785) =
0.01607, sehingga solusinya adalah
0.016073( , ) 100sin .
80tx
u x t e
Jadi, suhu maksimunnya turun menjadi 50 oC dalam waktu t = (ln 0.5)/(0.01607) 43
[detik], jauh lebih cepat (9 kali lebih cepat dibandingkan pada teladan 1).
Kalau kita mengambil n yang lebih besar, penurunan suhu akan lebih cepat lagi,
dan di dalam dewret yang suku-sukunya demkian, setiap suku mempunyai laju
penurunannya sendiri, dan suku dengan n besar praktis bernilai nol, dalam waktu singkat.
Teladan kita berikut adalah semacam ini, dan kurva dalam Gambar 276 untuk t = 0.5
sangat mirip dengan kurva sinus; artinya, praktis ini merupakan grafik suku pertama
solusi tersebut.
0=lxπu
5.0=lxπu
2=lxπu
1.0=lxπu
Gambar 276. Solusi pada Teladan 3
Untuk L = p, c = 1 dan bagian nilai t
Aliran Panas Berdemensi-Dua Keadaan-Stabil
Persamaan bedemensi-dua (lihat di awal Pasal ini)
Jika aliran panas itu stabil atau stasioner (tidak tergantung waktu), maka /dy t = 0 sehingga persamaan panas ini berubah menjdi persamaan Laplace4
2u =
2 2
2 20
u u
x y
Jadi suatu masalah panas rediri atas persamaan ini yang harus disimak pada suatu daerah
R pada bidang-xy dan suatu syarat batas pada kurva yang menjadi batas daerah R tesebut.
Masalah ini dinamakan masalah nilai batas (boundary value problem). Orang
menyebutnya
masalah Dirichlet jika u ditentukan pada C.
masalah Neumann jika turunan normalnya un = /u n ditentukan pada .masalah campuran jika u ditentukan pada sebagian dari C sedangkan un pada
sebagian lain dari C.
Sekarang kita simak suatu masalah Dirichlet pada suatu persegi panjang R(Gambar 77),
dengan mengasumsikan bahwa suhu u(x, y) sama dengan suatu fungsi f(x) tertentu di sisi
atas dan sama dengan 0 di ketiga sisi lainnya.
4Persamaan sangat penting muncul dalam pasal 9.7 dan akan dibahas lebih lanjut dalam pasal-pasal 11.9, 11.11, 11.12, 12.5 dan dalam Bab 17
Gambar 277. Persegi Panjang R dan nilai batasnya
Kita pecahkan masalah ini melalui pemisahan peubah. Subtitusi
( , ) ( ( )u x y F xG y
ke dalam (13) dan pembagian dengan FG menghasikan
2 2
2 2
1 1. . .d F d G
kF dx G dy
Dari sini menghasilkan syarat batas pada sisi kiri dan sisi kanan persegi panjang R, kita
memperoleh
2
20
d FkF
dx
0 0,F 0.F a
Ini menghasilkan k = (np/a)2 dan solusi bukan-nol
(14) F(x) = Fn(x) = sin ,
nx
a
n = 1, 2, . . .
Dengan demikian, persamaan untuk G menjadi4Persamaan sangat penting ini muncul dalam pasal 9.7 akan dibahas lebih lanjut dalam Pasal-Pasal 11.9, 11.11, 11.12, 12.5, dan dalam bab 17.
22
20
d G nG
dy a
Solusinya adalah
/ /( ) ( ) n y a n y an n nG y G y A e b e .
Sekarang syrat batas u = 0 di sisi bawah persegi panjang R berimplikasi bahwa Gn(0);
sehingga Gn(0) = An + Bn = 0 atau Bn = -An . ini menghasilkan
/ /( ) ( ) 2 sinh .n y a n y an n n n
n yG y A e b e A
a
Dari sini dan (14), dengan menuliskan 2an
*, kita memperoleh ”fungsi –eigen” bagi
masalah kita
(15) ( , ) ( ) ( )n n n nu x y f x G y A sin sin .
n x n y
a a
Fungsi-fungsi ini memenuhi syarat batas u = 0 pada sisi kiri, sisi kanan, dan sisi bawah
persegi panjang R. Untuk memperoleh solusi yang juga memenuhi syarat batas
(16) ( , ) ( )u x b f x Pada sisi atas persegi panjang R, kita perhatikan deret tak hingga
( , y)u x 1n
( , )nu x y
dari ini, (16) dan (15) dengan y = b, kita memperoleh
u(x, b) = f(x) = 1n
sin sinhn
n x n bA
a a
kita dapat menuliskan ini dalam bentuk
1
, sinh sinhnn
n b n xu x b A
a a
Ini menunjukkan bahwa eksperesi di dalam tanda kurung harus merupakan koefisien
Fourier bn bagi f(x); dengan kata lain, menurut (4) dalam pasal 10.5.
0
2sinh ( )sin
a
n n
n b n xb A f x dx
a a a
dari ini dan (15) kita memperoleh solusi bagi masalah kita, yaitu
(17)
1
, sin sinhnn
n x n yu x y A
a a
Dengan
(18) 0
2sin
sinh /
a
n
n xA f x dx
a n b a a
Solusi ini, yang diperoleh tanpa memperhatikan konvergen dan jumlah deret bentuk u, uxx dan uyy, dapat dibuktikan bila f dan f’ kontinu dan f” kontinu sepotong-sepotong pada
selang 0 x a. Buktinya agak rumit dan didasarkan pada konvrgenan seragam (uniform konvegen); bukti ini dapat anda temukan di dalam Acuan [C2] yang di cantumkan di dalam Apendiks 1.
Tegangan Elektrostatis
Persamaan Laplace (13) juga mengatur tegangan elektrostatis muatan listrik di
gambarang daerah yang bebas dari muatan tersebut. Jadi, masalah panas keadaan-stabil
dapat juga ditafsirkan sebagai masalah tegangan elektrostatis, sehingga (17), (18)
merupakan tegangan di dalam persegi panjang R jika sisi atas persegi panjang R
mempunyai tegangan f(x)sedangkan ketiga sisi lainnya dihubungkan dengan tanah.
Sesungguhnya, di dalam keaadan-stabil, persamaan gelombang berdemensi-dua
(akan dibahas dalam pasal 11.7, 11.8) tereduksi menjadi (13) dan (17); dengan demikian
(18) merupakan persegi (displacement) membran elastis persegi panjang (lembaran karet,
permukaan genderang) yang terikat (fixed) sepanjang batasnya, dengan ketiga sisinya
terletak pada bidang-xy dan sisi keempat mengalami pergeseran f(x). ”Mode karet” bagi
medan tegangan elektrostatis ini dapat digunakan (dan telah digunakan) untuk
menemukan lintasan elektron: kita tinggal membiarkan sebutir kelereng kecil
menyelusuri turunan dan memotret lintasannya. Begitu pula, ini akan menghsilkan kurva
yang merupakan lintasan aliran panas di dalam model panas.
Ini merupakan petunjuk yang mengesankan tentang matematika sebagai kekuatu
penyatu: sistem fisik yang sama sekali berbeda ternyata dapat diterangkan matemtis
yang sama-dan analogi demikian ini sering di gunakan di dalam pemodelan.
Di dalam pasal berikutnya, kita akan menunjukkan bahwa pembahasan kita
tentang persamaan panas dapat diperluas ke batang tak hingga. Ini merupakan model
yang baik bagi batang atau kawat yang sangat panjang (misalnya, kawat yang panjangnya
100 m). Di dalam model tersebut, itegral Fourier (Pasal 10.9) akan menggantikan deret
Fourier.
Soal-soal latihan untuk Pasal 11.5
1. Bandingkan Gambar 270 dengan Gambar 276, dan jelaskan perbedaan perilaku solusi
kedua persaman itu .
2. Bagaimanakah ketergantungan laju penurunan suhu (9) untuk n tetap tertentu pada
panas jenisnya, kerapatannya, dan konduktifitas termalnya?
3. Gambarkan grafik u1, u2, u3, [lihat (9) dengan 1, 1, nB c L ] sebagai fungsi dari x untuk nilai-nilai t = 0, 1, 2, 3. Bandingkan perilaku fungsi-fungsi tersebut.
Tentukan suhu x(x, t) di dalam sebuah batang perak (panjang 10 cm, luas irisan melintang
1 cm2, kecepatan 10.6 gm/cm3, kondukvifitas termal 1.04 kal/cm detoC, panas jenis 0.056
kal/gm oC) yang terisolasi sempurna, ujung-ujungnya dipertahankan pada suhu 0 oC dan
yang suhu awalnya (dalam oC) adalah f(x), jika
2
4. 0.2 5. sin 0.1
jika 0 5 jika 0 56. 7. =
0 jika 5< 10 10 jika 5 10
8. 100 9.
f x sin x f x x
x x x xf x f x
x x x
f x x x f x
= 10
10. = jika 0 2.5, =2.5 jika 5 7.5, = 10
jika 7.5 10
x
f x x x f x x f x x
x
11. Misalkan sebuah batang memenuhi asumsi-asumsi yang dinyatakan di dalam buku ini dan kedua ujungnya dipertahankan pada suhu konstan yang berbeda
1 20, dan , .u t U u L t U Tentukan suhu U1(x) di dalam batang setelah waktu yang
lama (secara teoritis, untuk t ).12. Di dalam Soal 11, misalkan suhu awalnya adalah u(x, 0) = f(x). Tunjukkan bahwa
suhunya pada sembarang waktu 1 110 adalah , ( ) , , t u x t u x u x t dengan u1
sepertisebelumnya dan
2/ ,
111
sin cn L tn
n
n xu B e
L
Dengan
1
0
2 1
2sin
2 2 sin [ 1 ].
L
n
n
n xB f x u x dx
L L
n xf x dx U U
L L n
13. (Ujung0ujung terisolasi, syarat batas diabatik) Tentukan suhu u(x, t) di dalam sebuah batang yang panjang L yang terisolasi sempurna, juga kedua ujungnya pada x = 0 dan x = L, dsg mengasumsikan bahwa u(x, 0) = f(x). infomasi fisik: fluks panas yang melalui permukaan kedua ujungnya sebanding dengan nilai /u x di kedua titik tersebut. Tunjukklan bahwa situasi ini dapat diucapkan secara matematis dalam bentuk syatrat-syarat berikut:
0, 0, , 0 0, .x x xu t u L t u t f x
Tunjukkan bahwa metode pemisahanpeuabahmenghasilkansolusi
2/
01
, cos c n L tn
n
n xu x t A A e
L
dengan, menurut (2) Pada Pasal 10.5.
0
0 0
1 2, cos , 1, 2,...
L L
n
n xA f x dx A f x dx n
L L L
14. Dalam Soal 13, u . Apakah ini sesuai dengan intiusi fisik anda?
Tentukan suhu di dalam bidang pada Soal 13 jika L = p, c = 1, dan
2
15. ( ) konstanta 16. ( )
17. ( ) cos 2 18. ( ) 1 1 /
19. ( ) cos 20. ( ) sin
21. ( ) jik
f x k f x x
f x x f x x
f x x f x x
f x x
a 0 / 2, ( ) 0 jika / 2 x f x x
22. Perhatikan batang dalam Soalp-Soal 4-10. Jika ujung-ujungnya dipertahankan pada
suhu 100 oC untuk waktu yang lama, maka pada suatu saat, katakanlah pada t = 0,
suhu di x = L tiba-tiba perubahan menjadi 0 oC dan tetap pada suhu ini, tengah batang
pada t = 1, 2, 3, 10, 50, detik?
23. (Pemancaran di ujung batang) Perhatikan sebuah batang yang terisolasi secara
memanjang, yang panjangnya p dan sedemikian rupa sehingga c = 1 di dalam (1),
yang ujung kirinya dipertahankan pada suhu 0 oC dan yang ujung kanannya dibiarkan
bebas memancarkan panas ke udara yang suhunya dijaga 0 oC. Informasi fisik. ”syarat
batas radiasi” ini adalah
0, [ , ],zu t k u t u dengan u0 = 0 adalah suhu udara di sekitar dan k adalah konstanta, misalkanlah k = 1 untuk mudahnya. Tunjukkan bahwa suatu solusi yang memenuhi syarat batas tersebut
adalah 2
( , ) sin p tu x t px e , dengan p adalah solusi bagi tan p p . Tunjukkan secara grafis bahwa persamaan ini memiliki tak hingga banyaknya solusi positif
2 2 3, , ,..., dengan p p p
1 1 dan 0.
2 2n np n p n
24. (Persamaan panas tak homogen) Perhatikan masaah yang terdiri atas persamaan
2 xt xxu c u Ne
dan syarat-syarat (2), (3). Di sini suku di atas kanan dapat merepresentasikan hilangnya panas akibat memancarnya panas keluar. Tunjukkan bahwa masalah dapat diubah menjadi suatu masalah yang dengan persamaan homogen sebagai berikut. Ambil
, , ( )u x t v x t w x dan tentukan ( )w x sehingga v memenuhi persamaan homogen
itu dan syarat 0, , 0, , 0 ( ).v t v L t v x f x w x
25. (Pemancaran) Jika batang yang dinyatakan di dalam buku dibiarkan memancarkan
panas ke medium di sekitarnya yang dijaga pada suhu konstan nol, maka persamaan
menjadi
2 .t xxv c v v Tunjukkan bahwa persamaan ini dapat direduksi ke bentuk (1) dengan mengambil
, , ( )v x t u x t w t.
Masalah Berdemensi-Dua
26. (Persamaan Laplace) Tentukan tegangan di dalam bujur sangkar
0 3, 0 3x y jika sisi atas dijaga pada tegangan
1sin
3x
sedangkan sisi yang lain dipertahankan pada tegangan nol.
27. Tentukan tegangan di dalam persegi panjang 0 20,0 40x y yang sisi atasnya dipertahankan pada tegangan 110 volt dan yang sisi-sisi lainnya dihubungkan dengan tanah.28. (Aliran Panas pada plat) Permukaan suatu plat tipis berbentuk bujur sangkar (Gambar 278, dengan a = 12) diisolasi sempurna. Sisi atas dipertahankan pada suhu 100 oC sedangkan sisi lainnya dipertahankan pada 0 oC . Tentukan suhu keadaan-stabil
, u x y pada plat tersebut.
xyaa
Gambar 278. Plat bujur sangkar]
29. Tentukan solusi u bagi persamaan panas berdemensi-dua
2t xx yyu c u u
pada plat tipis dalam Gambar 278 dengan a sedemkian rupa sehingga u = 0 di kedua sisi vertikalnya, jika diasumsikan bahwa kedua permukaan dan kedua sisi horizontal plat itu diisolasi sempurna.
29. Tentukan rumus yang serupa dengan (17) dan (18) untuk distribusi suhu di dalam
persegi panjang R yang digambarkan di dalam buku bila sisi bawah R diperhatikan
pada suhu f(x) sedangkan sisi lainnya dipertahankan pada 0.
30. Tentukan rumus yang serupa dengan (17) dan (18) untuk distribusi suhu di dalam
persegi panjang R yang digambarkan dalam buku bila sisi bawah R dipertahankan
pada suhu f(x) sedangkan sisi lainnya dipertahankan pada 0.
11.6 Aliran Panas pada Sebuah Batang Tak Hingga
Sekarang kita akan mempelajari solusi bagai persamaan panas
(1)
22
2
u uc
t x
Pada suatu batang yang memanjang tak hingga pada keda ujungnya (dan terisolasi
mebujur, seperti sebelumnya). Dalam kasus demikian, tidak ada syarat batas, yang ada
hanyalah syrat awal
(2) ( , 0) ( ) [ ]u x f x x dengan f(x) adalah suhu awal batang tersebut.
Untuk memecahkan masalah ini, kita mulai pada pasal sebelumnya; kita mensubsitusikan
u(x, t) = F(x)G(t) ke dalam (1). Ini menghasilkan persamaan diferensial biasa
(3) " 2 0 [lihat (6) Pasal 11.5]F p F
Dan
(4) 2 2 0 [lihat (7) Pasal 11.5]G c p G
Solusinya masing-masimg
(5) 2 2
cos + sin dan ( ) c p tF x A spx B px G t e
dalam hal ini A dan B adalah sembarang komstanta. Jadi, salah satu solusi bagi (1) adalah
2 2
, ; cos sin .c p tu x t p FG A px B px e
[Seperti pada soal sebelumnya, kita haurs mengambil konstanta pemisah k yang negatif, k
= -p2, sebab nilai k positif menghasilkan fngsi eksponensial yang naik di dalam (5),
sesuatu yang tidak ada makna fisikanya.]
Sembarang deret untuk fungsi (5), yang diperoleh melalui cara bisa dengan
mengambil p sebagai kelipatan sebuah bilangan tertentu, akan menghasilakan sesuatu
fungsi periodik, di dalam x bila t = 0. Akan tetapi, karena f(x) di dalam (2) tida
diasumsikan periodik, maka amat wajarlah untuk menggunakan integral Fourier sebagai
ganti deret Fourier.
Karena A dan B di dalam (5) adalah sembarang, kita dapat memandang kedua
kuantitas itu sebagai fungsi dari p dan menuliskan A = A(p) dan B = B(p). Karena
persamaan panas di atas linier dan homogen, maka fungsi
(6)
2 2
0 0
( , ) ( , ; ) [ ( )cos ( )sin ] c p tu x t u x t p dp A p px B p px e dp
merupakan solusi bagi (1), asalkan integral ini ada dan dapat didiferensialkan dua kali terhadap x dan sekali terhadap t.
Dari represintasi (6) dan syarat awal (2), kita memperoleh
(7) 0
( , ) [ ( )cos ( )sin ] ( )u x t A p px B p px dp f x
.
Dengan menggunakan (4) dan (5) dalam pasal 10.9, kita memperoleh
(8) 0 0
1 1( ) ( )cos , ( ) ( )sin . A p f v pvdv B p f v pv dv
Menurut (1 ) dalam pasal 10.11, integral Fourier ini dapat dituliskan sebagai
0
1, 0 cosu x f v px pv dv dp
,
sehingga (6) pasal ini menjadi
2 2
0
1, cos .c p tu x t f v px pv e dv dp
Dengan mengasumsikan bahwa kita dapat membalik urutan integrasinya, kita
memperoleh
(9) 2 2
0
1, cos .c p tu x t f v e px pv dp dv
kita hitung integral sebelah dalam dengan menggunakan rumus
(10)
2 2
0
cos 2 2
s be bs ds e
Yang akan diturunkan dalam pasal 15.4 (soal 13). Untuk maksud ini, kita dikenalkan di
dalam (10) suatu peubah integrasi baru p dengan mengambil s = t dan menetapkan
2
x vb
c t
Dengan demikian 2 dan ds = c t ,bs x v dp sehingga (10) menjadi
2 22 2 / 4
0
cos .2
x v c tc p te px pv dp ec t
dengan menyisipkan ini ke dalam (9) kita memperoleh reprentasi
2
2
1, exp .
42
x vu x t f v dv
c tc t
Akhirnya kita perkenalkan peubah integrasi z = (v - x)/(2c t ), sehingga
21
, 2 .zu x t f x cz t e dz
Jika f(x) terbatas (bounded) untuk semua nilai x dan terintegralkan pada setiap yang
terhingga, maka dapat diperlihatkan (lihat Acuan [C18] ) bahwa fungsi (11) dan (12).
Jadi, fungsi ini merupakan solusi yang kita cari.
Teladan 1. suhu pada batang yang panjangnya tak hingga
Tentukan suhu pada batang tak hinggga jika suhu awalnya adalah (Gambar 279)
0 jika 10 konstanta
0 jika 1
( ) xU
x
f x
Jawab. Dari (11) kita memperoleh1 2
02
1
( )( , ) exp .
42
U x vu x t dv
c tc t
Gambar 279. suhu awal pada Teladan 1
Gambar 28.Solusi U(x,t)pada teladan 1 untuk U0 100 oC, C2 = 1 cm2/det, dan beberapa nilai tJika kita memperkenalkan peubah integrasi z di atas, maka integrasi terhadap v dari -1
sampai 1 sama dengan integrasi terhadap z dari (-1-x)/2c t sampai (1-x)/2c t ,dan(13) (1 )2
0
(1 )2
( , )x c t
x c t
Uu x t
Telah di singgung bahwa integral ini bukan sebuah fungsi elementer, namun dengan
mudah dapat diucapkan sebagai fungsi galat, yang nilai-nilainya telah ditabelkan. (Tabel
A5 dalam Apendiks 4 mencantumkan beberapa niali; tabel yang lebih besar diberikan
dalam acuan []1 pada Apendiks 1. lihat juga Soal-Soal 11-20 di bawah.) Gambar 280
memperlihatkan u(x, t) untuk U0 = 100 oC, c2 = 1 cm2/det, beberapa nilai t
Persamaan-persamaan yang dibahas sejauh sampai ini (kecuali persamaan
Laplace pada separuh kedua Pasal 11.5) bersifat ”berdemensi-satu”; artinya, hany
melibatkan satu peubah ruang x (dan waktu t). Pada Pasal berikut, kita akan memulai
membahas persamaan ”bedemensi-dua” (persamaan gelombang berdemensi-dua), yang
melibatkan dua peubah ruang x dan y dan yang mengatur vibrasi membran elastis.
Soal-soal untuk Pasal 11.6
1. Gambrkan grafik suhu pada Teladan 1 (dengan U0 = 100 oC dan c2 = 1 cm2/detik)
pada titk-titik x = 0.5 1 dan 1.5 sebagai fungsi dari t. Apakah hasil sesuai dengan
intuisi fisik anda?
2. Jika f(x) = 1 untuk x > 0 dan f (x) = 0 untuk x < 0?, tunjukkan bahwa (12) tereduksi
menjadi
2
/ 2
1( , ) ( 0) z
x c tu x t e dz t
3. Perlihatkan bahwa untuk x = 0, solusi pada Soal 2 tidak tergantung pada t. Apakah
hasil ini sudah dapat diperkirakan sebelumnya bedasarkan alasan fisik?
4. Jika batangnya semi tak-hingga, memanjang dari nol sampai , ujungnya pada 0x
dipertahankan pada suhu 0 dan suhu awalnya adalah ( ),f x tunjukkan bahwa
(14)
2 2
/
1( , ) ( ) ( ) ,
xu x t f x w e dw f x w e dw
4. Turunkan (14) dari (11) dengan mengasumsikan bahwa f(v) di dalam (11) adalah
ganjil .
5. Jika f(x) = 1 dalam Soal 4, tunjukkan bahwa
/
0
2( , ) ( 0)
xu x t e dw t
6. Bagaimanakah bentuk (14) akan berubah jika f(x) = 1 untuk ( 0)a x b a dan f(x) = 0 untuk selainnya?7. Tunjukkan bahwa hasil pada Soal 6, dapat diperoleh dari (11) dan (12) dengan
menggunakan f(x) = 1 bila x > 0, dan f(x) = -1 bila x < 0. Apakh alasan yang
mendasarinya?
8. Tunjukkan bahwa pada Soal 6, waktu yang dibutuhkan agar dua titik mencapai suhu
yang sama sebanding dengan kuadrat jarak keduanya pada batas x = 0.
9. (Sebaran normal) Dengan mengasumsikan 2w z sebagai peubah integrasi yang baru, tunjukkan bahwa (12) menjadi
2 / 21( , ) ( 2 ) .
2u x t f x cw t e dw
Pembaca yang sudah mengenal teori peluagdapat melihat bahwa ini melibatkan
kepekatan (density) 2 / 2 2e
sebaran atau distribusi normal (lihat Pasal 23.7).Fungsi galat. Fug\ngsi galat (error function) didefenisikan oleh integral
2
0
2erf
xx e dw
Fungsi ini sangat penting bagi matematika rekayasa. Agar terbiasa dengan fungsi ini,
pembaca dipersilahkan memecahkan soal-soal berikut. (Berapa rumus diperikan pada
apendiks 3.1; lihat (35)-(37). Lihat juga teladan 1 dalam Pasal 18.6.)
10. Gambar integral fungsi erf x (yang dinamakan kurva bentuk-bel).
11. Tunjukkan bahwa erf x ganjil.
12. Tunjukkan bahwa erf x ganjil
13. Tunjukkan bahwa (13) dapat dituliskan dalam bentuk
0 1 1( , ) erf erf ( 0).
2 2 2
U x xu x t t
c t c t
14. Tunjukkan bahwa dalam Soal 6,
15. Dapat diperlihatkan bahwa ref ( ) = 1. Dengan menggunakan ini, tunjukkan bahwa
dalam Soal 2,
16. Tunjukkan bahwa
2 2
(erf erf ), = erf . 2
b b
a ae dw b a e dw b
17. Tunjukkan bahwa erf
22 / 2
0.
2
xx e d
18. Tentukan deret Maclaurin bagi erf x pengintegralan suku-demi-suku terhadap deret
Maclaurin bagi integral erf x tersebut.
19. Tunjukkan bahwa 2
erf xy e x memenuhi persamaan diferensial ' 2 2 / .y xy
20. Hitunglah erf x untuk x = 0(0.1)0.5 dari Soal 18 dan bandingkan dengan nilai-3D
0.000, 0.112, 0.223, 0.239, 0.428, 0.520.
11.7 Pemodelan: Membran Bervibrasi. Persaman Gelombang
bwedemensi-dua
Sebagai masalah penting lain yang masih berkaitan dengan vibrasi, marilah kita simak
gerak suatu membran yang teregang, misalnya pada sebuah gendang. Pembaca akan
mengetahui bahwa pembahasan kita ini sangat mirip dengan pada kasus dawai bervibrasi
pada pasal 11.2.
Kita memberlakukan asumsi-asumsi beriku:
1. Massa memberan per satuan luas adalah konstan (”membran homogen”).
Memberannya feleksibel sempurna dan begitu tipisnya sehingga tidak melawan
terhadap pelengkungan.
2. Membrang itu diregangkan dan kemudian di-fixed-kan sepanjang batasnya
pada bidang x-y. Tegangan per satuan panjang T yang disebabkan oleh
peregangan membran itu sama di semua tritik dan dalam semua arah, dan
tidak berubah selama gerak tersebut.
3. Deflexi u(x, y, t) membran itu selama gerak tersebut relatif kecil dibandingkan
dengan besarnya membran, dan semua sudut inkilinasi adalah kecil.
Meskipun asumsi-asumsi itu tidak mungkin terpenuhi di dalam praktek, vibrasi tranversal
kecil dari suatu membran tipis akan memenuhi asumsi-asumsi itu relatif cukup teliti.
Untuk menurunkan persamaan diferensial yang mengatur gerak membran
demikian ini, kita perhatukan gaya-gaya yang bekerja pada suatu bagian kecil membran
itu, seperti pada Gambar 281. karena defleksi membran itu dan sudut inklinasinya kecil,
sisi-sisi bagian kecil tersebut kira-kira sama dengan Dx dan Dy. Tegangan T adalah gaya
per satuan panjang. Oleh karena itu, gaya-gaya yang bekerja pada sisi-sisi bagian kecil
tersebut kira-kira sama dengan T Dx dan T Dy. Karena membrannya feleksibel sempurna,
gaya-gaya itu berarah tangen (menyinggung) membran tersebut.
Pertama-pertama marilah kita simak komponene horisontal gaya-gaya tersebut.
Komponen ini diperoleh melalui penggandaan gaya-gaya itu dengan kosinus sudut
inklinasinya. karena sudutnya kecil, kosinus dekat pada 1, sehingga komponen horisontal
gaya ini pada sisi yang berhadapan kira-kira sama. Akibatnya, gerak partikel memberan
tersebut dalam arah horisontal dapat diabaikan. Dari sini kita bisa menyimpulkan bahwa
gerak membran tersebut dapat dianggap tranversal; artinya, setiap partkel bergerak
vertikel.
Gambar 281. Membran bevibrasi
Gambar 281. Membran bervibrasi
Komponen vertikel gaya tersebut sepanjang sisi-sisi yang sejajar dengan bidang-yu
adalah5 (gambar 281)
sin dan sin ; T y T y di sini tanda minus muncul sebab gaya di sisi kiri terarah ke bawah. Karena sudutnya
kecil, kita dapat mengganti sinusnya dengan tangennya. Dengan demikian, resultan kedua
komponen vertikal itu adalah
(1) 1 2
(sin sin ) (tan tan )
= [ ( , ) ( )] x x
T y T y
T y u x x y u x y
Dalam hal ini subskrip x melambangkan turunan parsial, sedangkan parsial, sedangkan y1 dan y2 mempunyai nilai antara y dan y + Dy. Begitu pula, resultan komponen vertikal dari gaya-gaya yang bekerja pada kedua sisi lainnya adalah
1 2(2) [ ( , ) ( , )]
y yT x u x y y u x y
Dengan x1 dan x2 terletak antara x dan x + Dx.Menurt hukum kedua newton (lihat pasal 2.6), jumlah gaya-gaya yang diberikan
oleh (1) dan (2) sama dengan massa p DA dari bagian kecil itu kali percepatan 2 u/ t2;
disini p adalah massa membran yang yang tak terdefleksi per satuan luas, sedangkan DA
= Dx Dy adalah luas bagian kecil itu bila tidak terdefleksi. Jadi
2
1 22, , x x
ux y T y u x x y u x y
t
1 2 y , .y yT x u x y u x y ]
dalam hal ini turunan dari ruas kiri dihitung pada titik tertentu (??) pada bagian kecil tersebut. Pembagian dengan p Dx Dy menghasilkan
21 2 22
2
, , y yx x u x y u x yu x x u x yu T
t x y
Jika Dx dan Dy mendekati nol, maka ini memperole persamaan diferensial
(3)
2 2 22 2
2 2 2 .
u u u Tc c
t x y
Persamaan ini dinamakan persamaan gelombang berdimensi-dua. Ekspresi di dalam
tanda kurung adalah Laplacian D2 u bagi u (lihat pasal 8.8), dan (3) dapat dituliskan
sebagai
(3’)
22 2
2.
uc u
t
Solusi bagi persamaan diferensial akan dibahas dalam pasal berikut.
11.8 Membran Persegi Panjang
Untuk memecahkan masalah membran yang bevibrasi, kita harus menentukan suatu
solusi u(x, y, t) bagi persamaan gelombang berdemensi-dua
(1)
2 2 22
2 2 2
u u uc
t x y
Syarat memenuhi syarat batas (2) u = 0 pada batas membran tersebut untuk semua t 0jika kedua syrat awal
(3) ( , , 0) ( , ) [ ( , ) adalah pergeseran awal]u x y f x y f x ydan
0
, t
ug x y
t
[g(x, y) adalah kecepatan awal]
Syarat-syrat ini sangat mirip dengan syrat pada kasus dawai yang bervibrasi. Sebagai
kasus penting pertama, marilah kita membahas memberan persegi panjang R yang
ditunjukkan dalam Gambar 282.
Rxyab
Gambar 282. Membran persgi panjang
Langkah Pertama. Tiga Persamaan Diferensial Biasa
Dengan menerapkan metode pemisahan peubah, mula-mula kita tentukan solusi bagi
yang memenuhi syarat batas (2). Kita mulai dari
(5) ( , , ) ( , ) ( )u x y t F x y G t .
Dengan menerapkan ini ke dalam persamaan gelombang (1) kita memperoleh
2
1.xx yy
GF F
c G F
gg
Dalam hal ini subskrip melambangkan turunan persial, sedangkan dua titik di atas huruf
melambangkan turunan terhdap t. Untuk memisahkan peubah kita bagi ruas dengan
c2FG:
2
1.xx yy
GF G F Y
c G F
karena fungsi –fungsi ini di ruas kiri tergantung hanya pada t sedangkan fungsi –fungsi di
ruas kanan tidak tergantung pada t, berarti masing-masing ruas itu sama dengan suatu
konstanta. Dengan sedikit upaya dapay ditunjukkan bahwa hanya nilai-nilai negatif
sajalah yang akan menghasilkan solusi yang memenuhi (2) tanpa harus identik dengan
nol; ini sama dengan prosedur Pasal 11.3. dengan melambangkan konstanta negatif ini
dengan –v2, kita memperoleh
22
1.xx yy
GF F v
c G F
Dari sini kita memperoleh dua persamaan, diferensial biasa
(6) 2 0 dengan G G cv
Dan persamaan diferensial parsial
(7) 2 0.xx yyF F v F
Persamaan (7) dinamakn persamaan Helmholtz berdemensi-dua. Selanjutnya kita simak
(7) dan kita terapkan metode pemisahan peubah sekali lagi, dengan kata lain, kita
tentukan solusi bagi (7) yang berbentuk
(8) ( , ) ( ) ( )F x y H x Q y
yang bernilai nol pada batas memberan tersebut. Kita substitusikan (8) ke dalam (7):
2 22
2 2
1.
d H d QQ H v HQ
H dx dy
Fungsi di sebelah kiri hanya tergantung pada x, sedangkan fungsi di sebelah kanan hanya
tergantung pada y. Karena itu, ekspresi dari kedua ruas harus sama dengan sebuah
konstanta. Konstanta ini harus negatif, misalnya, -k2, karena hanya nilai negatif yang
akan memberikan solusi yang memenuhi (2) dan tidak identik dengan nol. Jadi,
2 22 2
2 2
1 1.
d H d Qv HQ k
H dx Q dy
Ini menghasilkan dua persamaan diferensial biasa linear untuk H dan Q:
(9)
22
20
d Hk H
dx
dan
(10)
22
20
d Hk H
dx
dengan 2 2 2.p v k
Langkah kedua. Memenuhi Syarat BatasSolusi umum bagi (9) dan (10) adalah
( ) cos sin dan ( ) cos sinH x A kx B kx Q y C py D py
Dengan A, B, C, dan D adalah konstanta. Dari (5) dan (2) haruslah fungsi F = HQ sama
dengan nol pada batas tersebut, yaitu x = 0, x = a, y = 0 dan y = b; lihat gambar 282. Ini
menghasilkan syrat
(0) 0, ( ) 0, (0) 0, ( ) 0H H a Q Q b Oleh karena itu, H(0) = A = 0, sehingga
( ) sin 0.H a B ka Kita harus mengambil B 0, sebab jika tidak H 0 dan F 0 dengan demikian sin ka =
0 sehingga ka = mp, atau
mk
a
(m bilangan bulat)
Dengan cara yang persis sama, kita akan memperoleh C = 0 dan membatasi nilai-nilai p =
np/b dengan n bilangan bulat. Adi, kita memperoleh solusi
( ) sin dan ( ) sin 1,2, ,
1,2, ,
m n
m x m xH x Q y m
a bn
.
(Seperti halnya pada kasus dawai yang bervibrasi, kita perlu memprhatikan m, n = -1,
-2, ... sebab solusi yang dihasilkan pada hakekatnya sama dengan solusi untuk m dan n
positif, kecuali perbedaan tanda.) Ini berarti bahwa fungsi
(11)
( , ) ( ) ( ) sin sin 1,2,3, ,
1,2,3, ,
mn m n
m x m xF x y H x Q y m
a bn
Merupakan solusi persamaan (7) yang bernilai nol pada batas membran persegi panjang
tersebut.
Karena p = v2 – k2 di dalam(10) dan l = cv di dalam (6), maka
2 2 .c k p
Dengan demikian untuk k = mp/b kita memperoleh nilai
(12)
2 2
2 2 mn
m nc
a b
m = 1, 2, ... n = 1, 2, ...
Di dalam persamaan diferensial (6) solusi umum bagi (6) adalah
( ) cos sin .mn mn mn mn mnG t B t B t
Ini berakibat bahwa fungsi –fungsi umn (x, y, t) = Fmn(x, y)Gmn(t) yaitu
(13)( , , ) ( cos sin )sin sinmn mn mn mn mn
m x m yu x y t B t B t
a b
dengan lmn dihitung menurut (12) merupakan solusi bagi persamaan gelombang (1) yang
bernilai nol pada batas memberan persegi panjang dalam gambar 282. Fungsi –fungsi
tersebut dinamakan Fungsi-eigen (eigenfunction). Atau fungsi-ciri (characteristik
funnction), sedangkan bilangan-bilangan lmn dinamakan nilai eigen (eigenvalue) atau
nilai-ciri (characteristik value) dari membran yang bervibrasi. Frekuensi dari umn adalah
lmn/2p.
Menarik untuk dicatat bahwa, tergantung antara a dan b, berapa fungsi Fmn
mungkin saja berasal dari nilai-eigen yang sama. Ditinjau dari segi fisik, ini berarti
bahwa mungkin saja ada vibrasi yang memiliki frekuensi yang sama namun dengan garis
simpul (nodal lines) (kurva titik pada memberan yang tidak bergerak) .
Marilah kita ilustrasikan ini dengan sebuah teladan.
Teladan 1.Nilai-eigen dan fungsi-eigen bagi menbran bujur sangkar
Perhatikan membran bujur sangkar dengan a = b = 1. Dari (12) kita memperoleh nilai-
eigen .
(14)2 2 .mn c m n
Dengan demikian
mn nm namun untuk m n fungsi-fungsi
sin sin dan sin sinmn nmF m x n y F n x m y
jelas berbeda. Misalnya, l12 = l21 = cp 5 masing-masing berpadanan dengan fungsi
12 21sin sin 2 dan sin 2 sin .F x y F x y Dengan demikian, solusinya
12 12 12 12( cos 5 sin 5 )u B c t B c t F
dan
21 21 21 21( cos 5 sin 5 )u B c t B c t F
masing-masing mempunyai garis simpul y = 1/2 dan x = 1/2(lihat gambar 283). Dengan
mengambil B12 = 1 dan B12 = B21 = 0, kita memperoleh
(15) 12 21 12 21 21cos 5u u c t F B F
yang ternyata mempresentasikan vibrasi lain yang berpadanan dengan nilai-eigen cp 5 . Garis simpul fungsi ini merupakan solusi bagi persamaan
12 21 21 21sin sin 2 sin 2 sin 0.F B F x y B x y
Atau, karena sin 2a = 2 sin a cos a,
(16) 21sin sin (cos cos ) 0.x y y B x
Solusi ini bergantung pada nilai B21 (Gambar 824)Dari (14) kita melihat bahwa bahkan mungkin saja lebih dari dua fungsi
bepedanan dengan nilai lmn yang sama. Misalnya, keempat fungsi F18, F81, F47, dan F74
berpadanan dengan nilai l18 = l81 = l47 = l74 = py 65 e, sbab1 2 2 21 8 4 7 65.
Ini terjadi karena 65 dapat dinyatakan sebagai jumlah dua kuadrat bilangan asli dalam
beberapa cara. Menurut sebuah teorema oleh Gauss, ini berlaku untuk setiap jumlah dua
kuadrat bilangan asli yang diantara faktor-faktor primanya terdapat sedikitnya dua faktor
yang berbeda yang berbentuk 4n + 1, n bilangan bulat positif. Dalam kasus kita, 65 =
(5)(13) = (4 + 1)(12 +1).
U11
U12
U21
U22
U13
U31
21 10B
21 1B
21 05B
21 0B
21 0.5B
Gambar 283. Garis simpul bagi solusi- Gambar 284. Garis-garis simpul
solusi 11, 12 21 22 13 31, , , , U U U U U U pada solusi (15) untuk beberapa nilai 21B
kasus membran berbentuk bujur sangkar.
Langkah Ketiga. Solusi Masalah Keseluruhannya
Untuk memperoleh solusi yang juga memenuhihi syarat batas awal (3) dan (4), kita akan
menempuh jalan seperti pada pasal 11.3. marilah kita simak deret ganda dua6
(17) 1 1
( , , ) ( , , )mnm n
u x y t u x y t
1 1
( cos sin )sin sin .mn mn mn mnm n
m x mB t B t
a b
sin
n y
b
6Kita tidak dapat menyinggung masalah kekonvergenan dan ketunggalan.
Dari sini dan (3) kita memperoleh
1 1
( , , 0) sin sin ( , ).mnm n
m x m yu x y B f x y
a b
)
deret ini dinamakan deret Fourier ganda. Jika f(x, y) dapat diuraikan menjadi deret
demikian ini7, maka koefisien-koefisien Fourier Bmn bagi f(x, y) di dalam (18) dapat
ditentukan sebagai berikut. Jika kita ambil
(19) 1
( ) sinm mnm
m yK y B
b
Maka (18) dapat kita tulis dalam bentuk
1
( , ) ( )sin .mm
m xf x y K y
a
untuk nilai y tertentu, ini merupakan deret sinus Fourier bagi f(x, y), dipandang sebagai
sutu fungsi dari x, dan dari (4) dalam pasal 10.5 kita dapat melihat bahwa koefisien-
koefisien uraian ini adalah
(20) 0
2( ) ( , )sin .
a
m
m xK y f x y dx
a a
Lebih lanjut, (19) merupakan deret sinus Fourier bagi Km(y), dan dari (4) pasal 10.5 kita
memperoleh koefisien-koefisien tersebut, yaitu
0
2( )sin .
a
mn m
m yB K y dy
b b
Dari sini dan (20) kita memperleh rumus Euler umum
(21)
0 0
4( , )sin sin 1,2, ,
1,2, ,
b a
mn
m x m yB f x y dx dy m
ab a bn
bagi koefisien-koefisien Fourier untuk f(x, y) di dalam deret Fourier ganda (18) Bmn di
dalam (17) ternyata sekarang ditentukan dalam f(x, y). Untuk menentukan B*mn, kita
diferensialkan (17) suku-demi-suku terhadap t; dengan menggunakan (4) kita
memperoleh
0 1
sin sin = , .
t
mn mnm n i
u m x m yB g x y
t a b
jika g(x, y) dapat diuraikan menjadi deret Fourier ganda, maka dengan menempuh
langkah-langkah sebelumnya kita memperoleh
0 0
4(22) , sin sin 1, 2, ,
b a
mnmn
n x n yB g x y dxdy m
ab a b
1, 2, ,n
Kesimpulannya adalah agar (17) memenuhi syarat-ayarat awalnya, maka koefisien-
koefisien Bmn dan mnB harus ditentukan menurut (21) dan (22).
Teladan 2. vibrasi sebuah memberan persegi panjang
Gerakan vibrasi sebuiah memberan persegi panjang dengan sisi-sisi a = 4 ft dan b = ft
jika tegangan adalah 12.5 lb/ft, kerapatannya adalah 2.5 slug/ft2 (seperti pada karet
ringan), kecepatan awalnya adalah 0 dan pergeseran awalnya adalah 7Syarat cukup:f, f/ x, f/ y, 2f/ x y kontinu pada daerah persegi panjang R yang dipekajari.
(23) 2 2( , ) 00.1(4 )(2 )ft.f x y x x y y
Jawab c2 = Tlp = 12.5 = 5(ft2/sec2) juga B*mn = 0 dari (22). Selanjutnya dari (21) dan
(23)
2 42 2 2
0 2
4 22 2
0 0
40.14 )(2 )sin (2 )sin
4 2 4
1 = (4 )sin (2 )
20 4 2
mn
m x m yB x x y y y y dx dy
m y
m x m yx x dx y y dy
Dua kali pengintegralan bagian-demi-bagian terhadap integral yang pertama
menghasilkan
3 3 3 3
128 2561 1 ganjil
mm
m m
dan terhadap integral yang kedua menghasilkan
3 3 3 3
16 321 1 ganjil
nn
m m
untuk m dan n genap kita memperoleh nol. Setelah digabungkan dengan faktor 1/20, kita
memperoleh Bmn genap dan
3 3 3 3
256 3 0.246 050 dan ganjil
20mnB m nm m
Dari sini dan (17) kita memperoleh solusi
2 23 3
.
1 5 ( , , ) 0.426050 cos 4 sin sin
4 4 2
5 5 1 5 37 3(24) 0.426050 cos sin sin cos sin sin
4 4 2 27 4 4 2
1 5 13 3 1 cos sin sin cos
27 4 4 2 729
m n vdd
m x m xu x y t m n
m n
x y x y
x y
5 45 3 3sin sin ...
4 4 2
x y
Perhatikan betapa cepat turunnya koefisien-koefisien deret ini.Untuk menyelidiki membran bervibrasi yang berbentuk lain, kita perlu menggunakan
suatu sistem koordinasi sehingga batas memberan itu dapat dinyatakan dalam arah rumus
yang sederhana; ini diperlukan karena adanya syarat batas (u = 0 sepanjang batas itu). Di
dalam praktek, kasus yang paling penting adalah membran lingkaran (permukaan
tambur, gendang). Untuk ini kita perlu menggunakan koordinat kutub dan mengucapkan
Laplacian dalam persamaan gelombang ke dalam koordinat kutub. Ini kita akan lakukan
pasal berikut.
Soal-Soal Latihan Untuk Pasal 11.8
1. Bagaiman frekuensi solusi (13) akan berubah jadi tegangan membrannya
ditingkatkan?
2. Tentukan dan gambarkan grafik garis simpul bagi solusi (130 untuk 1,2,3,4m dan 1,2,3,4 jika 1.n a b
3. Kerjakan tugas yang sama seperti pada Soal dua namun jika 2 dan 1.a b 4. Tentukan nilai-eigen lain bagi membran bujur sangkar dengan sisi 1 sedemikian rupa sehingga setiap nilai-eigen bagi membran bujur sangkar dengan sisi 1 sedemikian rupa sehingga setiap nilai-eigen menghasilkan emapat fungsi-eigen yang bebeda.
5. Tentukan nilai-eigen bagi membran perasegi panjang dengan sisi 2, , 1a b sedemikian rupa sehingga setiap nilai-eigen itu menghasilkan dua atau lebih fungsi-eigen yang berbeda6. Tunjukkan bahwa diantara semua semua membran persegi panjang dengan luas yang
sama 11 dan yang sama, A ab c u [lihat(13)] dari membran bujur sangkar mempunyai frekuensi terendah.7. Carilah hasil seperti pada Soal 6 bigi frekuensi solusi (13) dengan m dan n tetap tertentu.
8. Dapatkah andah memvisualisasikan ( , )f x y pada teladan 2? Dapatkah anda melihat
bahwa suku pertama pada deret (24) telah menghampiri ( , )f x y dengan sangat baiknya? Apa kira-kira alasannya?
9. Dengan mengintegralkan bagian-demi-bagian, verifikasi perhitungan mnB pada teladan 2.
10. mnB pada Teladan 2 merupakan hasil kali dua integral. Carilah ke fungsi-fungsi mana integral itu bepedanan pada verifikasi nilai-nilai integral itu dengan metode loncatan (Pasal 10.6) .
11. Dalam kasus apakah mnB di dalam (21) merupakan hasil kali dua integral (seperti pada Teladan 2)?
Deret Fourier ganda. nyatakanlah ( , )f x y sebagai suatu deret Fourier ganda yang
berbentuk (18) dengan 0 , 0 .x y 12. ( , ) 1 13. ( , )
14. ( , ) 15. ( , )
f x y f x y x
f x y y f x y xy
16.
jika 0 / 21
0 jika /2<
( , )x
x
f x y
17.
jika /2< , 1( , )
selainnya.0
x yf x y
Nyatakan fungsi-fungsi ( , )F x y berikut (0 , 0 )x a y b sebagai suatu deret Fourier ganda yang bentuk (18).
18. 1 19. ( , )
20. 21. ( )( )
f f x y x
f x y f xy a x b y
2 2 2 222. ( 1)( 1) 23. ( )( )f x y f xy a x b y
Tentukan defleksi ( , , )u x y t pada membran bujur sangkar dengan 1a b dan 1c
jika kecepatan awalnya nol dan defleksi awalnya adalah ( , )f x y :
2 2 2 2
24. ( , ) sin sin 25. ( , ) sin sin 2
26. ( , ) sin3 sin 4 27. ( , ) 0.01 (1 )(1 )
28. ( , ) sin sin 29. ( , ) (1 )(1 )
f x y k x y f x y k x y
f x y k x y f x y xy x y
f x y k x y f x y kx x y
30. (Vibrasi paksa terhadap membran) Tunjukkan bahwa vibrasi paksa terhadap suatu
membran mengikuti persamaan
2 2 /ttu c u P
dengan ( , , )P x y t adalah gaya persamaan satuan luas yang bekerja dengan arah normal terhadap bidang-xy. 11.9 Laplacian Dalam Koordinat Kutub
Dalam kaitan dengan masalah nilai batas untuk persaman diferensial parsial, menjadi
suatu perinsip umum untuk menggunakan suatu sistem koordinat yang dapat
mengucapakan batas daerah yang dipelajari menjadi rumus sederhana. Pada pasal berikut,
kita akan membahas membran lingkaran. Dalam kaitan dengan ini, koordinat kutub r dan
q yang didefenisikan sebagai
cos , sinx r y r cepat untuk digunakan, sebab batas memberannya dapat dinyatakan dalam persamaan
sederhana r = konstanta
Bila menggunakan r dan q, kita harus mengubah Laplacian
2 22
2 2
u uu
x y
di dalam persamaan gelombang ke dalam koordinat baru tersebut.
Transformasi bentuk-bentuk diferensial dari sistem koordinat yang satu ke sistem
koordinat yang lain sering digunakan di dalam penerapan. Oleh karena itu pembaca harus
menyimak pembahasan berikut ini dengan penuh perhatian.
Seperti pada pasal 11.4, kita akan menggunakan kaidah rantai. Untuk kemudahan,
kita akan melambangkan turunan parsial dengan subskrip dan u(x, y, t) sebagai fungsi
dari r, q, t dengan huruf yang sama u. Dengan menerapkan kaidah rantai (4) pasal 8.7,
kita memperoleh
(1) z r x xu u r u
Dengan mendiferensialkan sekali lagi terhadap x kita akan memperoleh
( )
=( ) ( ) .xx r x x x
r x x r xx x x xx
u u r u
u r u r u u
Sekarang, dengan menerapkan kaidah rantai sekali lagi, kita memperoleh
( ) dan ( )r x rr x r x x r x xu u r u u u r u .
Untuk menentukan turunan parsial rx dan qx, kita harus mendiferensialkan
2 2 y dan = arc tan ,
xr x y
Dan hasilnya adalah
x 2 2 22 2
1 .
1 /x
x x y yr
r x ry xx y
Dengan mendiferensialkan kedua rumus itu sekali lagi, kita memperoleh
2 2
2 3 3 3 4
1 2 2, .x
xx xx x
r xr x y xyr y r
r r r r r r
Sekarang kita subsitusikan semua ekspresi itu ke dalam (1). Jika turunan-turunan parsial
pertama dan kedua itu kontinu, kita memperoleh r ru u , dan setelah penyederanaan,
(2)
2 2 2
2 3 4 3 42 2 .xx rr r r
x xy y y xyu u u u u u
r r r r r
Dengan cara yang sama akan kita peroleh
(3)
2 2 2
2 3 4 3 42 2 .yy rr r r
y xy x x xyu u u u u u
r r r r r
Dengan menjumlahkan (2) dan (3) kita memperoleh Laplacian dari u dalam koordinat kutub
(4)
2 22
2 2 2
1 1.
u u uu
x r r r
Soal-soal latihan untuk pasal 11.91.Berikut rincian perhitungan yang menghasilkan (2) dan (3).
2.Tunjukkan bahwa (4) dapat dituliskan sebagai
22
2 2
1 1.
u u uu r
r r r r
3. Jika u tdak tergantung pada q, maka (4) tersederhanakan menjadi 2 /rr ru u u r . Turunkan hasil ini langsung dari laplacain dalam koordinat Kartesius
dengan mengasumsikam bahwa u tidak tergantung pada q.4.Transformasikan (4) kemali ke dalam koordinat Kartesius.
5. Tunjukkan bahwa satu-satunya solusi bagi 2 0u yang tergantung hanya pada
2 2 adalah ln .r x y u a r b
6. Jika x dan y adalah koordinat Kartesius, perhatikan bahwa cos sinx x y
dan sin cosy x y adalah koordinat Kartesius, dan vrifikasi melalui perhitungan
bahwa 2 .x x x y y yu u u
7. Ucapkan 2u dalam koordinat-koordinat , , jika , x ax b y cy d x y
adalah koordinat Kartesius.
8. Tunjukkan bahwa cos , sin , 0,1...,n nn nu r n u r n n adalah solusi bagi
2 0u dengan 2u dinyatakan oleh (4).
9. Tentukan solusi bagi persamaan laplcian di dalam cakram yang berjari-jari-jari 1
dengan nilai batas 2( ) 10cos .f
10. Dengan mengasumsikan bahwa pendiferensial suku-demi-suku dibolehkan, tunjukkan bahwa solusi persamaan bagi Laplace di dalm cakram R < 1 yang memenuhi
syarat batas ( , ) ( )( diketahui) adalahu R f f
01
1( , ) cos sin
n n
n nn
ru r a a n a
R R
dalam hal ini am bn adalah koefisien Fourier bagi f (lihat Pasal 10.3).
Tegangan elektrostatis. Masalah keadaan stabil. Tegangan elektrostatis u memenuhi
persamaan Laplace 2 0u di dalam setiap daerah yang bebas muatan. Selain itu,
persamaan panas 2 2
tu c u (lihat Pasal 1.5) tereduksi menjadi persamaan Laplace jika suhu tidak tergantung pada waktu t (”kasus keadaan-stabil”). Tentukan tegangan elektrostatis (sama artinya dengan: distribusi suhu-stabil) di dalm cakram R < 1 nilai-nilai batasnya adalah sebagai berikut
1 jika 0 110 jika / 2 011. ( ) 12. ( )
1 jika 0< < 0 jika /2< <3 / 2
jika / 2 / 213. ( ) 14. ( )
0 jika /2< <3 / 2 - jika /2<
u u
u u
2
jika / 2 / 2
<3 / 2
15. ( ) 15. ( )
u u
17. Tentukan tegangan elektrostatis di dalam separuh cakram 1, 0 ,r yang sama dengan 110q ( ) pada separuh lingkaran r = 1 yang merupakan batasnya dan sama dengan 0 pada ruas garis 1 1.x
18. Tentukan suhu keadaan-stabil u pada olat tipis separuh lingkaran , 0 < < ,r a jika separuh lingkaran r = a dijaga pada suhu tetap uo sedangkan ruas garis a x a dijaga pada suhu tetap u = 0. (Gunakan pemisahan peubah.) 19. Tentukan rumus bagi tegangan u pada sunbu-x di dalam Soal 15. gunakan empat suku
pertama deret ini untuk menghitung di 0.75, 5, 0.25,u x 0, 0.25, 0.5, 0.75 (dua desimal).20. Tentukan rumus bagi tegangan u pada sumbu-y dalam Soal 15.
21. Misalkan , ,r adalah koordnat bola, yang didefinisikan oloeh cos sin , sin sin , cos . x r y r z r
Jika u(x, y, z) hanya merup[akan fungsi dari 2 2 2 ,r x y z tunjukkan bahwa
2 2.rr ru u u
r
22. (Laplacian di dalam koordinat bola) Tunjukkan bahwa Laplacian di dalam
koordinat bola , ,r yang didefinisikan di dalam Soal 21 adalah 2
2 2 2 2
2 1 cot 1.
sinrr ru u u u u ur r r r
Tunjukkan bahwa ini juga dapat dituliskan sebagai
2 22 2
1 1 1sin .
sin sinru r u u ur r
23. Jika permukaan bola 2 2 2 2r x y z R dipertahankan pada suhu nol dan suhu di
dalam bola itu adalah f(r), tunjukkan bahwa suhu u(r, t) di dalam bola adalah solusi 2 2( ),t rr ru c u u
r
yang memenuhi syarat ( , ) 0, ( , 0) ( ).u R t u r f r 24. Tunjukkan bahwa dengan mengambil v ru , rumus-rumus dalam Soal 23 mengambil
bentuk 2 , ( , ) 0, ( , 0) ( ).t rrv c v v R t v r rf r sertakan syarat (0, ) 0v t (yang berlaku
sebab u harus terbatas pada r = 0), dan pecahkan masalah ini melalui pemisahan peubah.25. (Laplacian dalam koordinat tanbung) Tentukan bahwa Laplacian dalam koordinat
tabung (cylindrical coordinates) yang didefinisikan oleh cosx r , sin , y r z z adalah
22
1 1.rr r zzu u u u u
r r
11.10. Membran lingkaran. Persaman Bessel
Sekarang kita membahas vibrasi membran berbentuk bundar yang berjari-jari R (gambar
285). Dengan menggunakan koordinat kutub, yang didefenisikan oleh x = r cos q = y r
sin q, kita melihat dari (4) pasal 11.9 bahwa persamaan gelombang (3’) dalam pasal 11.7
mengambil bentuk
2 2 22
2 2 2 2
1 1.
u u u uc
x x r r r
yxR
Gambar 285. membran lingkaran
Kita akan mengupas solusi u(r, t) bagi persamaan ini yang setangkup secara radial
(radially symmetric)8, artinya tidak tergantung pada q. Dengan demikian, persamaan
gelombang itu tereduksi menjadi
(1)
2 22
2 2
1.
u u uc
t t r r
Karena membran itu diletakkan (fikxed) sepanjang batasnya r = r, kita mempunyai syrat
batas
(2) , 0 untuk semua 0u R t t kondisi yang tidak
tegantung pada q akan diperoleh jka syarat awal tidak tergantung pada q, dengan kata lain syarat awal itu berbentuk
(3) , 0 defleksi awal fu r f r r
dan
(4)
0
( ) kecepatan awal t
ug r g r
t
Langkah Pertama. Dua Persaman Diferensial Biasa
Dengan menggunakan metode pemisahan peubah, mula-mula kita menentukan solusi
bagi (1) yang memenuhi syarat batas (2). Kita mulai dari
( , ) ( ) ( ).u r t W r G t
8Untuk solusi yang tergantung pada q, lihat soal-soal latihan.
dengan mendefinisikan dan menyisipkan (5) ke dalam (1) dan membagi persamaan yang
dihasikannya dengan c2WG kita memperoleh.
2
1 1" '
GW W
c G W r
gg
dalam hal ini dua-titik di atas huruf menandakan turunan terhadap t, sedangkan tanda
aksen menandakan turunan terhadap t. Eksperesi di kedua ruas itu pasti sama dengan
suatu konstanta, dan konstanta ini harus negatif, katakanlah –k2, agar diperoleh solusi
yang memenuhi syarat batas dan dan tidak idetentik dengan nol. Jadi, Ini menghasilkan
dua persamaan diferensial linear biasa
(6) 2 0 dengan = G G ck
gg
Dan
(7)
21" ' 0W W k W
r
Langkah Kedua. Pemenuhan Syarat batas
Pertam-tama kita simak (7). Dengan memperkenalkan peubah bebas baru s = kr, kita
memperoleh 1/r = k/s, sehingga menurut kaidah rantai, turunan-turunan itu menjadi.
22
2' dan W"= .
dW dW ds dW d WW k k
dr ds dr ds ds
Dengan mensubsitusikan ini ke dalam (7) dan menghapus faktor sekutu k2, kita
memperoleh
2
2
1W 0
d W dW
ds s ds
Ini tidak lain adalah persamaan Bessel (1), Pasal 4.5, dengan v = 0. solusi umum adalah
(lihat pasal 4.6)
1 0 2 0W c j s c y s dengan J0 dan Y0 masing-masing adalah fungsi Bessel jenis pertama dan kedua yang
berordo nol. Karena defleksi membran selalu terhingga, sedangkan Y0 menjadi tak hingga
jika s mendekati 0, berarti kita tidak mungkin menggunakan Y0 sehingga harus
mengambil C2 = 0. Jelaslah, C 0, sebab jika tidak demikian W 0. kalau kita ambil c1 =
1, maka kita peroleh
(8) 0 0 .W r J s J kr
Pada batas r = R haruslah u(R, t) = W(R)G(t) = 0. karena G 0 berimplikasi
u 0, maka kita syaratkan
0 0.W R J kR Fungsi Bessel J0 mempunyai tak hingga banyaknya ke-nol-an atau nilai nol (zeros) yang
nyata. Marilah kita lambangkan nilai nol 1, J0(s) Dengan 1 2, ,...s (lihat Gmbar 286). Nilai numeriknya adalah (teliti sampai 4 desimal)
1 2 3 4 52.4048, 5.5201, 8.6537, 11.7915, 14.9309. Kita lihat bahwa nilai-nilai nol itu tidak berjarak sama, dan begitu pula dengan nilai-nilai
nol lainnya 6 7, ... (Tabel yang lebih lengkap diberikan di dalam acuan [1] pada Apindeks 1.) Sekarang, persamaan (8) berimplikasi
(9) atau , 1, 2, .m
m mkR k k mR
Gambar 286. Fungsi Bessel J0(s)
Dengan demikian fungsi-fungsi
(10) 0 0( ) ( ) 1, 2, ,m
m mW r J k r J mR
merupakan solusi bagi (7) yang brnilai nol pada r = R.
Solusi umum padanannya bagi (6) dengan m mck adalah ( ) cos sinm m m m mG t a t b t .
Jadi, fungsi-fungsi
(10) 0( , ) ( ) ( ) ( cos sin ) ( )m m m m m m mu r t W r G t a t b t J k r
Gambar 287. Modus normal membran lingkaran pada kasus vibrasi yang tidak tergantung pada sudutnya
dengan 1, 2, ,m adalah solusi bagi persamaan gelombang (1) yang memenuhi syarat batas (2). Fungsi-fungsi itu tidak lain adalah fungsi-aigen bagi masalah kita dengan
nilai-eigennya m
Vibrasi membran yang menghasilkan mu dinamakan modus normal ke-m; vibrasi
ini memilki frekuensi / 2m siklus per satuan waktu. Karena nilai nol dari J0 tidak berjarak sama (not regularly spaced) pada sumbu tersebut (berbeda dengan nilai nol fungsi sinus yang muncul dalam kasus vibrasi dawai), ini berakibat suar drum/gendang sama sekali bebeda dengan suara biola. Bentuk modus normal dapat diperoleh dengan mudah dari Gambar 286 dan ditunjukkan dalam Gambar 287. Untuk m = 1, semua titik membran bergerak ke atas (atau ke bawah) pada saat yang sama. Untuk m = 2, keadaannya sebagai berikut. Fungsi
22 0W r J r
R
bernilai nol jika 2 1 1 2/ atau = / .r R r R Dengan demikian, lingkaran 1 2= /r R merupakan garis simpul, sehingga ketika pada suatu saat bagian tengah membran itu
bergerak ke atas, maka bagian luar 1 2( / )r R bergerak ke bawah, dan begitu
sebaliknya. Solusi ( , )mu r t mempunyai m – 1 garis simpul, yang berupa lingkaran-
lingkaran konsentris (Gambar 287).
Langkah Ketiga. Solusi Masalah Keseluruhan
Untuk memperoleh solusi yang juga memenuhi syarat-syarat awal (3) dan (4), kita dapat
menempuh cara seperti pada kasus dawai yang bervibrasi, yaitu menyimak deret9
(12) 0
1 1
( , ) ( ) ( ) ( cos sin ) .mm m m m m
m m
u r t W r G t a t b t J rR
Dengan mengambil t = 0 dan menggunakan (3), kita memperoleh
(13) 0
1
,0 mm
m
u r a j r f rR
9Kita tidak akan menyinggung maalah konvergenan dan ketunggalan
Jadi, memenuhi (3), am haruslah merupekan koefisien bm deret Fourier-Bessel
mempresentasikan f(r) dalamJ0 (amr/R); dengan kata lain [lihat (9) pasal 4.9],
02 21 0
2 m =1, 2, ...
Rm
mm
a rf r j r drR j R
Keterdiferensialan f(r) pada selang 0 r R sudah cukup untuk menjamin keberadaan uraian (13), Lihat acuan [A 14]. Koefisien bm di dalam (12) dapat ditentukan dari (4)
dengan cara yang serupa. Untuk memperoleh ini numerik bagi dan ,ma b kita dapat menerapkan salah satu metode integrasi hampiran, dengan menggunakan tabel-tabel J0 -
dan J1. Adakalanya kita tidak perlu menggunakan integrasi numerik, sebagai diilustrasikan di bawah ini.
Teladan 1. vibrasi sebuah membran lingkaran
Vibrasi permukaan genderang yang berjari-jari 1 ft dan berkerapatan 2 slug/ft2 jika
tegangannya adalah 8 lb/ft, kecepatan awalnya 0, dan pergeseran awalnya adalah
2( ) 1 [ft].f r r
Jawab. c2 T/p = 8/2 = 4 [ft2/det2]. Karena kecepatan awalnya 0, maka bm = 0. Dari (14)
dan dalam pasal 4.9, karena R = 1,
122
02 21 10
421 m
m mm m
ja r r j r dr
j j
Menggunakan tabel-tabe dalam Acuan (1) hlm. 390 – 394, kita menghitung
m am J1(am) J2(am) am
1 2.4048 0,51915 0.43170 1.108
2 5.5201 -0.34026 -0.12316 -0.140
3 8.6537 0.27145 0.06262 0.045
4 11.7915 -0.23246 -0.03940 -0.021
5 14.9309 0.20655 0.02762 0.012
Jadi,
f(r) = 1.108J0(2.4048r) – 0.140J0(5.5201r) + 0.045J0(8.6537r) - ...
Kiat lihat bahwa koefisien-koefisien itu turunan relatif lambat. Jumlah koefisien-
koefisien yang diuraikan dalam tabel di atas adalah 1.004. jumlah semua koefisien-
koefisien adalah 1.(Mengapa?)
Karena lm = ckm = cam/R = 2am, maka dari (12) kita memperoleh solusi (dengan r
diukur dalam kaki (feet) dan t dalam detik)
u (r, t) = 1.108J0(2.4048r)cos 4.8096t - 0.140J0(5.5201r)cos 11.0402t…
+ 0.045J0 (8.6537r) cos 1.73074t - ...
Di dalam pasal berikutnya, kita akan memperluas pembahasan tentang
persamaan Laplace, yang telah kita singgung sebelumnya (pasal 9.7 11.5), yang
mungkin merupakan persamaan diferensial parsial paling penting di dalam fisika dan
penerapan rekayasanya.
Soal-Soal Latihan Untuk Pasal 11.10
1. Apakah alasan penggunaan koordinat kutub yang disebutkan dalam pasal ini?
2. Jika tegangan membran ditingkatkan, bagaimanakah frekuensi masing-masing modus
normal (11) berubah?
3. Jika tegangan dan kerapatannya sama, genderang yang lebih kecil memiliki modus
dasar yang lebih tinggi dibandingkan dengan genderang yang lebih besar.
Bagaimanakah ini kita dapat simpulkan dari rumus-rumus kita di atas?
4. Tentukan nilai numerik bagi jari-jari kedua garis simpul u2 dan u3 [lihat (11)] bila R =
1.
5. Buatlah Gambar seperti 287, untuk 4 5 6, , .u u u
6. Pertanyaan yang sama pada Soal 4 untuk 4 5 6, , .u u u
7. Tunjukkan bahwa agar (12) memenuhi (4) maka,
(15) 0
0
2( ) / , 1,2 ,
R
m mm m
b rg r J r R dr mc RJ
8. Mengapa pada Teladan 1 harus berlaku 1 2 1?a a
9. Verifikasi perhotungan ma pada Teladan 1.
Tentukan defleksi ( , )u r t bagi membran yang berjari-jari-jari 1R jika 1c , kecepatan
awalnya nol dan defleksi awalnya adalah ( )f r :2 4
0 210. 0.1 ( ) 11. (1 ) 12. (1 )f J r f k r f k r Petunjuk. Ingat Soal-Soal 25-32, Pasal 4.9.
13. Mungkinkah untuk nilai c dan R tertentu, dua atau lebih fungsi mu [lihat (11) yang memiliki garis simpul yang berbeda berpedanan dengan nilai-eigen yang sama? Vibrasi suatu Membran Lingkaran yang Bergantung pada jari-jari r
dan
14. Tnjukkan bahwa substitusi ( , ) ( )u F r G t ke dalam persamaan gelombang
(16)
22
1 1tt rru c u u
r r
menghasilkan
(17) 2 0, dengan ,G G ck
(18)
22
1 10.rr rF F F k F
r r
15. Tunjukkan bahwa subtitusi ( ) ( )F W r Q ke dalam (18) menghasilkan
(19) 2" 0Q n Q
(20) 2 2 2 2 2' ( ) 0.r W rW k r n W
16. Tunjukkan bahwa Q(q)memounyai perioda 2 , sehingga oleh karena itu 0, 1,...n di dalam (19) dan (20). Tunjukkan bahwa ini menghasilkan solusi
sin , ( ), 0, 1,...n n nQ n W J kr n 17. Tunjukkan bahwa syarat batas
(21) ( , , ) 0u r t 18. Tunujukkan bahwa solusi bagi (16) yang memenuhi (21) adalah
( cos sin ) ( )cosmn mn mn mn mn n mnu A ck B ck t J k r n (22)
( cos sin ) ( )sin .mn mn mnmn mn n mnu A ck B ck t J k r n
19. Tunjukkan bahwa 0mou dan mou identik dengan (11) dalam Pasal ini.
20. Tunjukkan bahwa mou mempunyai 1m n garis simpul.Gambar garis simpul solusi-solusi berikut ini.
3 421. , 1, 2, 3 22. , 1, 2, 3 23. , 1, 2, 3n n mnu n u n u m n
24. Tunjukkan bahwa syarat awal ( , , 0) 0 tu r menghasilkan 0, 0mn mnB B dalam (22)
25. Tunjukkan bahwa 11u mempresentasikan modus dasar (fundamental mode) suatu
membran seperuh lingkaran dan tentukan frekuensinya jika 2 1c dan R = 1.
11u
12u
13u
Gambar 288. Garis simpul sebagian solusi (22)
11.11 Persaaman Laplace, Potensial
Salah satu persaman diferensial parsial paling penting di dalam fisika adalah persaman
Laplace
(1)2 0u
Di sini 2u adalah Laplacian fungsi u. dalam koordinat kartesius x, y, z di dalam ruang,
(2)
2 2 2 22
2 2 2 2.
u u u uu
x x y z
Teori tentang solusi persaman Laplace dinamakan teori potensial. Solusi bagi (1) yang
memiliki turunan parsial orde-kedua yang kontinu dinamakan fungsi harmonik.
Kasus berdimensi dua, bila u bergantung hanya pada dua peubah saja, paling
cocok dipecahkan melalui metode analisis kompleks dan akan dibahas dalam pasal 12.5
dan Bab 17.
Untuk mengilustrasikan pentingnya persamaan Laplace bagi matematika rekayasa
marilah kita singgung beberapa penerapan dasar.
Persamaan Laplace dijumpai dalam kaitan dengan gaya gravitasi. Pada teladan
pasal 8.8, kita telah melihat bahwa jika sebuah partikel A yang bermassa M berada pada
titik (X, Y, Z) dan sebuah partikel lain B yang bermassa m berada pada titik (x, y, z), maka
A dan B saling tari-menarik sebesar gradien fungsi skalar
u(x, y, z) =
c
r , c = GMm = const,
r = 2 2 2( ) ( ) ( )x X y Y z Z (> 0).
Fungsi x, y, z ini dinamakan potensial medan gravitasi, dan momentum persaman
Laplace.
Perluasan ke potensial dan gaya yang disebabkan oleh distribussi massa yang
kontinu tidaklah sulit. Jika suatu massa yang berkerapatan p(X, Y, Z) didistribusikan di
seluruh daerah T di dalam ruang, maka potensialnya u pada suatu titik (x, y, z) yang tidak
diduduki oleh massa itu didefenisikan sebagai
(3) , , >0
T
u x y z k dx dy dz kr
Dalam hal ini r diberikan oleh rumus sebelum ini. Karena 1/r (r > 0) merupakan solusi
bagi (1) , artinya 2(1/r) = 0, dan tidak bergantung pada x, y, z maka kita memperoleh
2 2 1 =0;
T
u k dx dy dzr
Ini berarti potensial gravitasi yang didefenisikan oleh (3) memenuhi persamaan Laplace
pada sembarang titik yang ditempati oleh massa tersebut.
Di dalam elektrostatik, gaya tarik-menarik atau tolak-menolak listrik antara
partikel-partikel yang bermuatan mengikuti hukum Coulomb (lihat Pasal 8.8), yang
mempunyai bentuk matematis yang sama dengan hukum gravitasi Newton. Ini berarti
medan yang diciptakan oleh suatu distribusi muatan listrik dapat diterangkan secara
matematis oleh sebuah fungsi potensial yang memenuhi persamaan Laplace pada
sembarang titik yang tidak diduduki oleh muatan.
Di dalam bab 17, kita akan melihat bahwa persamaan Laplace juga dijumpai
dalam teori aliran fluida yang tak termanpatkan (theory of incompressible fluid flow).
2 2tu c u
Lebih jauh, persaman pokok di dalam konduksi panas adalah persaaman panas
(lihat pasal 9.7 dan pasal 11.5). jika suhunya u tidak tergantung pada waktu t (“keadaan-
stabil”), maka persaman ini tereduksi menjadi persaman Laplace; lihat pasal 11.5.
Pada sebagian besar penerapan yang membawa pada persaman Laplace, kita
diharuskan memecahkan sesuatu masalah nilai batas, artiny menentukan solusi bagi (1)
yang memenuhi syarat batas yang diberikan pada permukaan S yang menjadi batas
daerah T yang merupakan daerah defenisi persaman tersebut. Ini dinamakan:
(I) Masalah nilai batas pertama atau masalah Dirichlet jika u ditentukan
pada S
(II) Masalah nilai batas kedua atau masalah Neumann jika turunan normal /nu u n ditentukan pada S.
(III) Masakah nilai batas ketiga atau campuran jika u ditentukan pada suatu
bagian dari S dan un ditentukan pada bagian sisanya dari S.
Selanjutnya kita perlu menggunakan koordinat di dalam ruang yang membuat S dapat
diucapkan dalam rumus yang sederhana. Ini mengharuskan kita melakukan transformasi
terhadap Laplace (2) ke dalam sistem koordinat lain. Tentu saja,
xzyrq(r, q, z)xzyrq(r, q, f)
Gambar 289. Kordinat tabung Gambar 290. Koordinat bola
transformasi demikian ini sangat serupa dengan transformasi pada kasus Laplacian suatu
fungsi dua peubah (lihat Pasal 11.9
Dari (4) Pasal 11.9 dapat ditunjukkan bahwa Laplacian sebuah fungsi u dalam
koordinat tabung (clinderical koordinates)10 (lihat gambar 289)
(4) r =2 2x y . q = arc tan
y
x , z = zadalah
(5)
2 2 22
2 2 2 2
1 1.
u u u uu
r r r r z
Koordinat penting lain adalah koordinat bola (spherical coordinates) r, q, , yang
berhubungan dengan koordinat kartesius sebagaiberikut (lihat gambar 290):
(6) cos sin , y sin sin , z cos .x r r r
Laplacian suatu fungsi u dalam koordinat bola adalah
10Perhatikan bahwa tidak ditentukan sepenuhnya oleh rasio y/x, namun harus diperhitungkan juga tanda x dan tanda y.
(7)
2 2 22
2 2 2 2 2 2 2
2 1 cot 1.
sin
u u u u uu
r r r r r r
Ini juga dituliskan dalam bentuk
(7’)
22 2
2 2 2
1 1 1sin .
sin sin
u uu r
r r r
Rumus ini diturunkan dengan cara yang sama seperti pada Pasal 11 rinciannya disediakan sebagai latihan bagi pembaca.
Di dalam Pasal berikut, kita akan menunjukkan bahwa pemisahan peubah pada koordinat bola akan menghasilkan persamaan Legendre, yang telah kita bahas pada Pasal 43.
Soal-soal Latihan untuk Pasal 11.11
1. Tentukan Laplacian dalam koordinat persegi panjang , y , ,x ax by z cz jika , ,x y z adalah koordinat-koordinat Kartesius.
2. Dengan mulai dari (2), verefikasi melalui perhitungan bahwa relative terhadap koordinat Kartesius yang didefinisikan oleh (7) dan (8) halaman 562.
2 .x x y y z zu u u u
Tunjukkan bahwa fungsi-fungsi ( , )u f x y berikut memenuhi persamaan Laplace dan buatlah ploy beberapa garis ekipotensial u = konstanta
2 2 2 2
3 2 2 2 2 2 2 2 2
3. 4. 5. y( )
6. - 3 7. / ( ) 8. ( )( )
x y xy x y
x xy x x y x y x y
9. Verifikasi bahwa /u c r memenuhi persaman Laplace dalam koordinat bola.
10. Tunjukkan bahwa satu-satunya solusi 2 0u yang tergantung hanya pada peubah
2 2 2r x y z adalah / ;u c r k dalam hal ini c dan k adalah konstanta.11. Tentukan c dan k dalam soal 10 sedemikian rupa sehingga u merepresentasikan
tegangan elektrostatis antara dua bola konsentris yang bejari-jari 1 2 cmr dan
2 10 cmr yang masing-masing dipertahankan pada tegangan 1 110 voltU dan
2 0 voltU .12. Tunjukkan bahwa satu-satunya solusi persamaan Laplace berdemensi-dua yang
tergantung hanya pada 2 2 2r x y z adalah ln .u c r k
13. Tentukan tegangan elektrostatis antara dua tabung koaksial yang berjari-jari 1 2 cmr
dan 2 10 cmr yang masing-masing dipertahankan pada tegangan 1 110 voltU dan
2 0 voltU . Gambarkan grafiknya dan bandingkan solusi Soal 11 dengan solusi Soal 13.14. Ucapkan koordinat bola yang didefinisikan oleh (6) dalam koordinat Kartesius.
15. Verifikasi (5) dengan mentransformasikan 2 0u kembali ke koordinat Kartesius.
Tentukan distribusi suhu keadaan-stabil (tidak bergantung waktu):
16. Antara dua plat paralel 0 1 dan x x x x yang masing-masing dijaga pada suhu
1 0 dan u u .
17. Antara dua tabung lingkaran koaksial yang berjari-jari-jari 0 1 dan r r yang masing-
masing dijaga pada suhu 1 0 dan u u .
18. Antara dua balok konsentris yang berjari-jari-jari 0 1 dan r r yang masing-masing
dijaga pada suhu 1 0 dan u u .
19. (persamaan Helmholts 11) Tunjukkan bahwa substitusi ( , , ) i tu U x y z e ( 1)i ke dalam persamaan gelombang berdemensi-tiga
2 2ttu c u menghasilkan
apa yang dinamakan persamaan Helmholts
2 2 0 /U k U k c 20. Misalkan , ,r adalah koordinat-kordinat bola. Jika (r, q, ) memenuhi
2 0u , tunjukkan bahwa v(r, q, ) = r-1u(r-1, q, ) memenuhi 2 0u .
11.12. persaman Laplace dalam Koordinat Bola. Persaman
Legendre
Marilah kita simak suatu masalah nilai batas tipikal yang melibatkan persamaan Laplace
dalam koordinat bola. Misalnya sebuah bola S yang berjari-jari r dipertahankan pada
suatu distribusi tegangan listrik tertentu
(1) 2 " '2 1 0r G rG n n G
dengan r, q, f adalah koordinat yang didefenisikan dalam Pasal sebelumnya dengan titik
asal di pusat S, dan f(f) adalah suatu fungsi yang diketahui. Kita ini mengetahui tegangan
u di semua titik di dalam ruang, yang diasumsikan bebas muatan-muatan lain. Karena
tegangan di S tidak tergantung q, maka begitu pula dengan tegangan di dalam ruang.
Oleh karena itu, 2u/ 2 = 0 dan dari (7’) dalam Pasal 11.11 maka kita melihat bahwa
persamaan Laplace tereduksi menjadi
(2)
2 1sin 0.
sin
u ur
r r
11HERMANN VON HELMHOLTZ (1821-1894),seoran fisikawan jerman, dikenal krena karya-karya pentingnya dalam bidan termodinamika, hidrodinamika dan akustika.
Lebih lanjut, di tak hingga tegangannya akan sama dengan nol; artinya, haruslah
(3) lim , 0.
ru r
Kita akan memecahkan masalah nilai batas yang terdiri atas persamaan (2) syarat
batas (1) dan syarat (3) di tak hingga dengan menggunakan metode pemisahan peubah.
Dengan mensubstitusikan solusi yang berbentuk
, u r G r H
ke dalam (2) dan membagi persamaan yang dihasilkannya dengan GH , memperoleh
21 1sin .
sin
d dG d dHr
G dr dr H d d
Melalui penalaran seperti biasa, kedua ruas persamaan harus sama dengan sebaran
Konstanta, katakanlah k, sehingga
(4)
21 1 sin 0 dan .sin G
d dh d dGkH r k
d d dr dr
Persamaan terakhir ini dapat dituliskan sebagai
2 " '2 0.r G rG kG
Ini tidak lain adalah persamaan Euler-Cauchy. Dari pasal 2.7 kita tahu bahwa persaman
ini mempunyai solusi yang berbentuk .G r Agar solusi itu mempunyai bentuk yang
sangat sederhana, kita ubah notasi kita dan kita ganti k dengan ( 1)n n .Dengan demikian
(5) 2 " '2 1 0.r G rG n n G
dengan n masih berupa sembarang bilangan. Dengan mensubstitusikan G r ke dalam (5), kita memperoleh
1 2 1 0.n n r
Nilai nol dari ekspresi yang di dalam kurung siku adalah 1.n jadi, kita memperoleh solusi
(6) n+1
1 dan r = .
rn
n nG r r G
Dengan menyisipkan 1k n n ke dalam (4) dan mengambil
cos w .
kita memperole 2 2sin 1 w dan
sin .
d d dw d
d dw d dw
Dengan demikian, (4) mengambil bentuk
2(7) 1 1 0.
d dhw n n H
dw dw
atau
2
22
(7 ') 1-w 2 1 0.
d H dHw n n H
dw dw
Ini adalah persamaan Legendre (lihat Pasal 4.3).
Solusi Persamaan Legendre
Untuk bilangan bulat12 n = 0, 1, ..., polinom Legendre
cos , n = 0, 1, ...,n nH P W p merupakan solusi persamaan Legendre (7). Jadi, kita memperoleh dua barisan solusi u GH di bawah bagi persamaan Laplace (2):
(8*) n 1
, cos , , cosn nn n n nn
Bu r A r p n r p
r
dengan n = 0, 1, ..., dan An dan Bn adalah konstanta.Untuk mendapatkan solusi bagi (2) yang sah di titik-titik di dalam bola tersebut
dan memenuhi (1), kita perhatikan deret13
(8)
0
, cos .nn n
n
u r A r P
Agar (8) memenuhi (1), haruslah
(9)
0
, cos ;nn n
n
u R A R P f
artinya (9) haruslah merupakan deret Fourier umum bagi f( ) dalam polinum Legendre
Dari (6) dan (7) Pasal 4.7 dan (2) Pasal 4.9, dapat diperlihatkan bahwa
1
1
2 1
2n n
n nnA R w P w dw
Rf
:
Dalam hal ini
wf:
menyatakan f ( ) sebagai fungsi dari w = cos karena dw= - d , dan batas-batas integrasi 1 dan -1 berpadanan dengan = p dan = 0 maka kita juga memperoleh
0
2 1cos sin , 0, 1, .
2n
n nnA P d n
Rf
Jadi, deret (8) dengan koefisien (10) merupakan solusi bagi masalah kita untuk titik
dalam bola.
Untuk memperoleh solusi eksterior terhadap bola (titik diluar bola), kita tidak
dapat menggunakan fungsi un(r, ) sebab dalam hal ini fungsi-fungsi itu tidak memenuhi
(3), namun kita dapat menggunakan fungsi-fungsi un (r, ), yang memenuhi (3), dan
dapat menempuh langkah-langkah seperti sebelumnya. ini menghasilkan solusi
(11) 1
0
, cos nnn
n
Bu r P r R
r
Dengan koefisien-koefisien
(12)
1
0
2 1cos sin .
2nn
n nB R P df
Teladan 1.kapasitor Bola Tentukan tegangan di dalam dan di luar seuah kapasitor berbentuk bola yang terdiri atas
dua separuh-bola yang berjari-jari-jari 1 ft yan dipisahkan oleh suatu celah kecil untuk
xyz110 volt
Gambar 291. kapasitor bola pada Teladan 1
insilusi,jika separuh bola bagian atas dipertahankan pada 110 volt sedangkan bagian
bawahnya dihubungkan ke tanah.(Gambar 291 ).
Solusi , syarat batasnya adalah (ingat gambar 290 di Pasal 11.11)
110 jika 0 / 2
0 jika /2< .
f
Karena R = 1, maka dari (10) kita memperoleh
/ 2
0
2 1110 cos sin
2n n
nA P dp
Jika kita ambil w = Cos , maka Pn (cos f) sin d = -Pn (w) dw dan kemudian kita
integralkan dari 1 sampai 0.Tanda minus dapat kita buang dengan mengintegralkan dari 0
sampai 1.Dengan demikian ,dari (11) dalam Pasal 4.3.
0
2 2 !55 2 1 1 .
2 ! ! 2 !
mM
n nm
n mA n
m n m n m
Dengan M =n/2 untuk n genap dan M =(n - 1)/2 unuk n ganjil. Integral diatas sama
dengan 1/(n - 2m + 1.Jadi
(13)
0
55 2 1 2 2 !1 .
2 ! ! 2 !
mM
n nm
n n mA
m n m n m
Dengan mengambil n = 0, kita memperoleh A0 = 55 (karena 0! = 1).
Untuk n = 1, 2, 3,….kita memperoleh
1
2
3
165 2! 165.
2 0!1!2! 2275 4! 2!
04 0!2!3! 1!1!1!
385 6! 4! 385
8 0!3!4! 1!2!2! 8
A
A
A
dan seterusnya. Jadi tegangan (8) di dalam bola itu adalah (karena Po = 1)
31 3
165 385, 55 cos ,
2 8u r rP r P coc
dalam hal ini P1, P3, ... diberikan oleh (11’) paal 4.3. Karena R = 1, kita peroleh dari (10)
dan (12) dalam pasal ini Bn = An, sehinggga (11) menghasilkan tegangan di luar bola
1 32 4
55 165 385, cos
2 8u r P P coc
r r r
Jumlah parsial kedua deret itu sekarang dapat digunakan untuk menghitung nilai
hampiran bagi tegangan tersebut. Selain itu, menarik untuk dilihat di bahwa jauh dari
bola tersebut, tegangannya kira-kira sama dengan satu muatan titik, yaitu 55/r (lihat
Teladan 4 pasal 8.8).
Di dalam dua pasal terakhir bab ini, kita akan menunjukkan bahwa masalah
melibatkan persamaan diferensial dapat juga dipecahkan melalui metode operasional
(transformasi Laplace dan trasformasi Fourier).
Soal-Soal Latihan untuk pasal 11.12
1. Verifikasi melalui substitusi bahwa nu (r, ), 0, 1, 2,n di dalam (8) adalah
solusi (2).
2. Tentukan pada permukaan mana fungsi-fungsi 1 2 3, , u u u bernilai nol.
3. Gambarkan grafik fungsi-fungsi Pn(cos ) untuk 0, 1, 2,n [ihat (11’) pasal 4.3].4. Gambarkan fungsi-fungsi P3(cos ) dan P4(cos ).
5. Hitunglah nilai-nila 0 1 3, , A A A pada Teladan 1 dari (13).6. Verifikasi nilai-nilai dalam Soal 5 dengan mengintegralkan rumus (11’) padanannya pada Pasal 4.3.Misalkan r, q, adalah koordinat boa yang digunakan di dalam buku ini. Tentukan
tegangan di dalam interior bola R = 0, jika di dalam interior bola ini tidak ada muatan dan
tegangan pada permukaannya adalah f( ), jika dalam hal ini
2 3
7. ( ) 1 8. ( ) cos 9. ( ) cos2
10. ( ) cos 11. ( ) cos 12. (
f f f
f f f
3 2
) cos3
13. ( ) 10cos 3cos 5cos 1f
14. Tunjukkan bahwa pada Soal 7, tegangan di eksterior bola sama dengan tegangan
suatu muatan titik di titik asal.
15. Gambarkan irisan permukaan-permukaan bertegangan sama (equipotential surfaces)
pada Soal 8 dengan bidang-xz.
16. Tentukan tegangan di eksterior pada Soal-Soal 7-13.
17. Verifikasi melalui substitusi bahwa jawaban Soal 10 memenuhi (2).
18. (Persamaan saluran teransmisi) Perhatikan suatu kabel panjang atau kawat telepon (Gambar 292) yang terisolasi tidak sempurna sehingga terjadi kebocoran
sepanjang kabel atau kawat tersebut. Sumber S dari arus ( , )i x t terletak pada x 0, sedangakan ujung penerimanya T terlekak pada x = l. Arus mengalir dari S ke T, melalui muatan, dan kembali ke tana. Misalkan konstanta-konstanta R ,L, C , dan G masing-masing melambangkan resistansi, induktansi, kapasitansi ke tanah, dan konduktansi ke
tanah kabel itu per satuan panjang. Tunjukkan bahwa
u iRi L
x t
(persamaan
saluran transmisi)
dengan ( , )u x t adalah tegangan di dalam kabel. Petunjuk. Terapkan hukum tegangan Kirchhoff pada suatu bagian kecil kabel itu antara dan x x x (selisih antara tegangan di dan x x x = turunnya resistansi + turunnya induktansi .)19. Tunjukkan untuk kabel dalam Soal 18,
i uGu c
x t
(Persamaan saluran transmisi kedua)
0=xlx =
TSBeban
Petunjuk. Gunakan hukum aru Kirchhoff (selisi arus di dan x x x = kehilangan akibat kebocoran ke tanah + kehilangan kapasitas).
Gambar 292. Saluran transmisi
20.Tunjukkan bahwa menghilangkan i dari persamaan saluran transmisi menghasilkan
( ) .xx tt tu LCu RC u RGu 20. (Persamaan telegraf) Untuk suatu kabel dasar laut, G dapat diabaikan dan
frekuensinya rendah. Tunjkkan bahwa ini menghasilkan apa yang dinamakan
persamaan kabel dasar laut (submarine cable equation) atau persaman telegraf
(telegraph equation)
, .xx v xx tu RCu i RCi
22. Tentukan tegangan pada satu kawat dasar laut yang ujung-ujungnya ( 0,x )x l dihubungkan ke tanah dan distribusi tegangan awalnya U0 = konstanta. 23. Masalah apakah pada konduksi panas yang analog dengan soal 22?
24. (Persamaan saluran frekuensi-tinggi ) Tunjukkan bahwa dalam kasus arus bolak-
balik berfrekuensi tinggo, persamaan dalam Soal 22 dan persamaan analoginya bagi
arus i dapat dihampiri dengan apa yang dinamakan persamaan saluran frekuensi-
tinggi (high-frekuncy line equation)
, .xx v xx ttu LCu i LCi 25. Pecahkan persamaan saluran frekuensi-tinggi yang pertama dengan mengasumsikan
bahwa tegangan awalnya adalah 0 sin( / ), ( ,0) 0tU x l u x dan 0u di kedua ujung 0 dan 1x x untuk semua t.
11.13 Penerapan Transformasi Laplace pada Persamaan
Differensial Parsial
Pembaca yang telah mengenal transformasi Laplace (Bab 5) mungkin bertanya-tanya
apakah metode itu dapat diterapkan untuk memecahkan persamaan differensial Parsial.
Jawabnya adalah ya, terutama jika salah sat dari peubah bebasnya mempunyai nilai-nilai
pada sumbu positif.Gagasan pokoknya akan dijelaskan dalam dua teladan tipikal dibawah
ini. langkah-langkah penyelesaiannya serupa pada Bab 5; untuk persamaan dengan dua
peubah bebas,langkah-langkah itu adalah sebagai berikut
1. Ambillah transform Laplace terhadap salah satu dari kedua peubah bebas itu
biasanya t.ini menghasilkan suatu persamaan diferensial biasa bagi transform
fungsi uang belum diketahui itu.ini disebabkan turunan fungsi terhadap peubah
yang satu lagi menyusup ke dalam persamaan yang telah ditransformasi .yang
terakhir ini juga memasukkan syarat batas dan syarat awal yang diberikan.
2. Pecahkan persamaan diferensial biasa itu untuk memperoleh tranformasi fungsi
yang tidak diketahui.
3. Terapkan transformasi kebalikan pada transformasi itu untuk memperoleh solusi
masalah yang dihadapi.
Jika koefisien-koefisien persamaan yang diberikan tidak tergantung pada
transformasi akan menyederhanakan masalahnya.
Teladan 1.Persamaan ordo-pertama
Pecahkan masalah
(1) 0, , 0, 0,
w wu x t w t t
x t
Kita lambangkan w sebab kita menggunakan u untuk melambangkan fungsi tangga satuan.(Pasal 5.3).
Jawab. Kita ambil transform Laplace bagi (1) terhadap t.menurut (1) dalam Pasal 5.
(2) L ,0 0,
wx sl w w x
x
Dalam hal ini w(x, o) = 0. Di dalam suku yang pertama kita mengasumsikan bahwa kita
boleh saling menukarkan perintegralan pendiferensialan dan pendiferensialan:
(3) L
0 0
, , .st stw we dt e w x t dt w x t
x x x x
Dengan melambangkan ( , 0) 0w x .maka kita memperoleh dari (2)
0.W
xsWx
Ini dapat dipandang sebagai suatu perasamaan diferensial biasa dengan x sebagai peubah
bebasnya ,sebab turunan terhadap s tidak muncul di dalam persamaan itu. solusi
umumnya adalah (Pasal 1.7).
2 / 2, .sxW x s c s e
Karena L (t) = 1/s2,syarat w (0, t) = t menghasilkan W (0, s) =1/s2.ini berarti
2
2
/ 22
0, 1/ .
Dengan demikian
1 , = .
sx
W s c s s
W x s es
Sekarang L -1{1/s2} = t, sehingga teorema pergeseran kedua (Pasal 5.3) dengan 2 / 2a x
(4)
sin jika 0 t 2
22
0 selainnya
1, ( )
2 2
t
xw x t t u t x
Karena kita mengikuti rumus-rumus saja, kita masih harus memverifikasi bahwa (4)
memenuhi (1). Penulis serahkan kepada pembaca.
Teladan 2. Dawai semi takhingga
Tentukan pergeseran (displacement) w(x, t) sebuah dawai elastis yang kenai syarat-syarat
berikut.
(i) Dawai ini mula-mula diam dan terletak pada sumbu –x dari x = 0 sampai
(”dawai semi-takhingga”).
(ii) Pada waktu t > 0, ujung kiri dawai itu digerakkan seperti pada Gambar
293
jika 0 2sin
0 selainnya
0, tt
w t f t
(iii) Lebih lanjut
lim , 0 untuk 0.x
w x t t
Tentu saja tidak ada dawai takhingga, namun model kini dapat diterapkan pada dawai
atau tali yang panjang sekali (yang beratnya dapat diabaikan) dengan ujung kanannya
diikatkan jauh pada sumbu –x (sekali lagi kita gunakan w sebab u kita perlukan untuk
melambangkan fungsi tangga satuan.
Jawab. Kita harus memecahkan persamaan gelombang (Pasal 11.2)
(5)
2
2
w
t
22 .
wc
t
Untuk x dan t positif, dengan ”syarat batas”
(6) 0, lim , = 0 ( 0).
xw t f t w x t t
dengan f seperti diberikan di atas ,dan syarat-syarat awal
1-1
f( t )p
2pt
Gambar 293. Gerak ujung kiri dawai pada
teladan 2 sebagai fungsi dari waktu t
(7) , 0 0w x
(8) 0
0t
w
t
Kita ambil transform Laplace terhadap t.menurut (2) Pasal 5.2
L
22
2
ws
t
L 2
0
, 0t
ww sw x c
t
L
2
.w
t
Berdasarkan (7) dan (8), dua suku hilang.diruas kanan kita mangasumsikan bahwa kita
boleh saling menukarkan perintegralan dan perdiferensial :
L
2 2 2 2
2 2 2 20 0
, st stW We dt e x t dt
x x x x
L
, .W x t
Dengan menuliskan, , W x s L
, ,w x t kita memperoleh
22 2
2
Ws W c
x
Atau
2 2
2 20.
W sW
x c
Karena ini mengandung hanya turunan terhadap x, maka persamaan ini dapat dipandang
sebagai sebuah persamaan diferensial biasa bagi , ,W x s dengan memandang terakhir
ini sebagai fungsi dari x. Solusi umumnya adalah
(9) / /, ( ) ( ) .sx c sx cW x s A s e B s e
Dari (6) kita memperoleh, dengan menulisnya F s L
,f t
0, W s L
0, w t L
f t F s
dan, dengan mengasumsikan bahwa uraian pengintegralan terhadap t dan pengambilan limt untuk 0x dapat saling dipertukarkan.
0 0
lim , lim , lim , 0x x
s st
x x xW x s e w x t dt e w x t dt
Ini berimplikasi 0A s di dalam (9) sebab 0,c sehingga untuk setiap s positif tertentu
fungsi /xs ce naik dengan naiknya x. Perhatikan bahwa kita dapat mengasumsikan 0s
sebab transforman Laplace ada untuk semua s yang lebih besar daripada suatu tertentu (Pasal 5). Oleh karena itu kita memperoleh
0, .W s B s F s Sehingga
/, .sx cW x s F s e
Dari teorema pergeseran kedua (Pasal 5.3) dengan /a x c kita memperoleh transformasi kebalikan (gambar. 294)
(10) , .
x xw x t f t u t
c c
(t = 0)(t = 2p)(t = 4p)(t = 6p)
2pcxxxx
yang berarti
Gambar 294. Gelombang berjalan pada Teladan 2
, sin jika < < 2 atau 2 x x x
w x t t t ct x t cc c c
dan nol solusinya. Ini merupakan gelombang sinus tunggal yang bergerak ke kanan
dengan kecepatan c. Perhatikan bahwa sembarang titik x tetap diam sampai t = x/c, waktu
yang dibutuhkan untuk mencapai x tersebut jika kita berangkat pada t = 0 (awal gerak
ujung kiri) dan bergerak dengan kecpatan c. Hasilnya sesuai dengan intuisi fisik kita.
Karena kita hanya mengikuti rumus-rumus saja, kita masih harus memverifikasi bahwa
(10) memenuhi syarat -syarat yang ditetapkan. Pembaca penulis persilahkan untuk
melakukannya sendiri .
Di dalam Pasal berikut, kita akan menunjukkan bahwa disamping transformasi
Laplace, metode oprasional lainnya, yaitu transformasi Fourier (lihat Pasal 10.10-
10.12), dapat juga digunakan untuk memecahkan persamaan diferensial parsial.
Soal-soal Latihan untuk Pasal 11.13
1. Buatlah Gambar seperti pada Gambar 294 jika 1c dan f adalah ”segi tiga” seperti pada Teladan 1 pasal 11.3 dengan / 2 1.k L 2. Bagaimanakah kecepatan gelombang pada Teladan 2 tergantung pada tegangan dan massa dawai tersebut?3. Verifikasi solusi pada Teladan 2. Gelombang berjalan apakah yang akan kita peroleh pada Teladan 2 jka kita menganakan suatu gerak sinusoidal (yang tidak henti-hentinya) di
ujung kiri mulai pada t = 0?Pecahkan melalui transformasi Laplace:
4. 2 2 , ( , 0) 1, (0, ) 1
u ux x u x u t
x x
5. , ( , 0) 0 jika 0, (0, ) 0 jika 0.
u ux xt u x x u t x
x t
6. Pecahkan Soal 5 dengan cara lain.
Tentukan suhu ( , )w x t di dalam sebuah batang semi-takhingga, yang di isolasi memanjang mulai dari 0x sepanjang sumbu-x sampai , jika diasumsikan bahwa suhu
awalnya adalah 0, ( , )w x t untuk setiap nilai t 0 tertentu, dan ( , ) ( ).w x t f x Tempulah langkah-langkah berikut ini.7. Buatlah modelnya dan tunjukkan bahwa transformasi Laplace menghasilkan
22
2( , )
wsW x s c W
x
L{w}
dan
/( , ) ( ) s x cW x s F s e F L{f}
8. Dengan menerapkan teorema konvolusi pada soal 7, tunjukkan bahwa
2 23/ 2 / 4
0
( , ) ( ) .2
tx cx
w x t f t e dc
9. Misalkan 0(0, ) ( ) ( )w t f t u t (Pasal 5.3). Lambangkan w, W dan F padanannya dengan w0, W0 dan F0. Tunjukkan bahwa dengan demikian di dalam Soal 8
2 23/ 2 / 40
0
( , ) 1 erf2 2
tx cx x
w x t e dc c t
dalam hal ini fungsi galat erf didefinisikan seperti pada Soal latihan untuk Pasal 11.6.
10 Rumus Duhamel14) Tunjukkan bahwa di dalam Soal 9
/0
1( , ) s x cw x t e
s
dan teorema konvolusi menghasilkan rumus Duhamel
0
0
( , ) ( ) .t w
w x t f t d
11.14 Penerapan Transformasi Fourier pada Persaman
Diferensial Parsial
14JEAN MARIE CONSTANT DUHAMEL (1797-1872), matematika Prancis.
Sekarang kita akan membahas sebuah cara lain untuk untuk memecahkan persamaan
diferensial parsial, yaitu penerapan salah satu transformasi Fourier. Ini berlaku untuk
masalah-masalah yang syrat-awal dan syarat batasnya diberikan pada sumbu-positif-
sehingga transformasi cosinus atau sinus Fourier (Pasal 10.10) dapat digunakan atau pada
seluruh sumbu-dalam hal demikian kita dapat menggunakan transformasi Fourier dalam
Pasal 10.11. Kita bahas cara-cara ini melalui penerapan tipikal.
Teladan 1. masalah panas pada sumbu-x
Tentukan suhu u(x, t) pada sebuah batang homogen yang irisan melintangnya konstan,
yang diisolasi memanjang dari x = sampai , untuk 0,x t dengan mengasumsikan bahwa suhu awalnya adalah
(1) , 0 u x f x x
Dan untuk semua 0,t solusi dan turunannya terhadap x memenuhi syarat
(2) , 0. , 0 untukxu x t u x t x
Sebagai kasus khusus, tentukan u(x, t) bila
(3) 0 f x U konstanta
jika 1 dan 0 jika >1.x f x x
Jawab. Kita harus memecahkan persamaan panas
(4) 2
t xxu c u
Yang memenuhi syarat (1), (2). Strategi kita adalah mengambil transformasi Fourier
terhadap x dan kemudian memecahkan persamaan diferensial biasa dalam t yang
dihasilkannya.rinciannya sebagai berikut.
Misalkan u
F u melambangkan transform Fourier bagi u, dipandang sebagai
suatu fungsi dari x. Dari (10) dalam Pasal 10.11, kita melihat bahwa (4) menghasilkan
F 2tu c F 2 2 .u c w u
Di ruas kiri, dengan mengasumsikan bahwa kita boleh saling menukarkan urutan
perdiferensialan dan perintegralan, kita memperoleh
F
1 1
2 2
i x i xt t
uu u e dx u e dx
t t
Sehingga
2 2 .
uc w u
t
Karena ini melibatkan hanya turunan terhadp t dan tidak melibatkan turunan terhadap w,
berarti persamaan ini merupakan persamaan diferensial biasa orde-pertama dengan t
sebagai peubah bebasnya dan w prameter. Dengan memishkan peubah (Pasal 1.2) kita
memperoleh solusi umum
2 2
, .c tu w t C w e
dengan dalam hal ini ”konstanta” C w tergantung pada prameter w . Syrat awal (1)
menghasilkan , 0 ( ) ( )u w C w f w
F ,f
sehingga kita memperoleh
2 2
, .c tu w t f w e
Rumus pembalikan (7) Pasal 10.11 sekarang menghasilkan solusi
(5) 2 21
, .2
c t i xu x t f w e e dw
Ke dalam ini kita dapat menyisipkan transform Fourier
1.
2i xf w f v e dv
Dengan mengasumsikan bahwa kita boleh mengubah urutan perintegralannya, maka kita
memperoleh
2 2
0
1, .
2i wx wvc tu x t f v e e dw dv
Menurt rumus Euler (3) Pasal 10.11 integral yang di dalam kurung siku sama dengan
2 2 2 2
cos sin .c t c te wy wv ie wx wv
Ini menunjukkan bahwa bagian khayalnya merupakan suatu fungsi ganjil dari w, sehingg
integral15 bagian ini sama dengan 0, sedangkan bagian nyatanya adalah genap, sehingga
integralnya sama dengan dua kali integral dari nol sampai ;
2 2
0
1, cos .
2c tu x t f v e wx wv dw dv
Ini sesuai dengan (9) Pasal 11.6, dan menghasilkan rumus-rumus (11) dan (13) dalam
Pasal 11.6
Teladan 2. Masalah pada Teladan 1 dipecahkan melalui metode
konvolusi
Pecahkan masalah panas pada teladan 1 dengan metode konvolusi.
IjawabI. Awalnya seperti sebelumnya dan menghasilkan (5), yaitu
(5) 2 21
, 2
c t i xu x t f w e e dw
Sekarang kita sampai pada gagasan metode konvolusi. Kita mengenali bahwa ini
bebentuk (13) Pasal 10.11, yaitu
(6 , i xu x t f g x f w g w e dw
)
dengan dalan hal ini
(7) 21
.2
c tg w e
Karena, menurut definisi konvolusi[(1) Pasal 10.11],
(8) .f g x f p g x p dp
maka sebagai langkah beriktnya kita harus menentukan transformasi Fourier kebalikan g
bagi g
. Untuk ini kita dapat menggunakan rumus 9 dalam Tabel III Pasal 10.12 (yang telah dibuktikan dalam Teladan 2 Pasal 110.11),15Sesungguhnya, bagian utama integral ini; lihat Pasal 15.3.
F
2 / 41( )
2ar u ae e
a
dengan a yang sesuai. Dengan menggunakan (7) dengan c2t = 1/4a atau a = 1/4c2t, kita
memperoleh
F 2 2 2 2/ 4 2 24 2 2 .x c t c te c t e c t g w
Jadi, g
mempunyai kebalikan 2 2/ 4
2
1? .
2
x c tec t
Mulalui penggantian x – p dan pensubstitusian ini ke dalam (8) akhirnya kita memperoleh
(9)
2
2
1, exp .
42
x pu x t f g x f p dp
c tc t
Hasil ini sesuai dengan (11) dalam pasal 11.6. perhatikan bahwa kita melambangkan f g x
tanpa menunjukkan prameter t, karena kita tidak melakukan pengintegralan terhadap peubah itu.
Teladan 3. Transform sinus Fourier diterapkan pada persaman panas
Tentukan solusi pada Teladan 1 jika batang itu memanjang dari 0 sampa , suhu awalnya
, 0u x t f x r
Dan syarat batasnya di ujung kiri adalah
0, 0 0u t t Jawab. Alih-alih transformasi Fourier, sekarang kita akan menerapkan transformasi sinus Fourier (Pasal 10.10), karena x berkisar antara 0 sampai . Dengan menemouh jalan
seperti pada Teladan 1, dari persamaan panas (9b) Pasal 10.10, karena (0) (0, 0) 0,f u kita memperoleh
F 2s
ts
uu c
t
F 2 2
xxu c w F 2 2 , .ss
u c w u w t
Solusi persamaan diferensial bias orde-pertama ini adalah
2 2
, .c tsu w t w e
Dari syarat awal , 0u x= f x
, kita memperoleh , 0 .s su w f w C w
Dengan
demikian
2 2
, .c ts su w t f w e
Dengan mengambil transform sinus kebalikan dan mensubstitusikan
0
2sin sf w f p wp dp
kita memperoleh solusi yang kita cari
(10) 2 2
0 0
2, sin sin c tu x t f p wpe wx dp dx
Teladan 4. Persamaan gelombang pada selang tak hingga. Solusi
d'Alembert
Pecahkan persamaan gelombang
, 0tt xxu cu x t Dengan syarat
, 0 defleksi awal diketahui
, 0 0 kecepatant
u x f x
u x
awal nol
0, 0 jika x untuk semua t. xu u dengan hal ini f dan g diasumsikan memiliki transfom Fourier.
Jawab. Kita transformasikan terhadap x, dengan menuliskan u F u
. Dengan
menggunakan (10) Pasal 10.11 kita memperoleh dari persamaan gelombang
F
2ttttu u c
F 2 2 .ttu c w u
2 2
sehingga
0.ttu c w u
Karena turunan terhadap w tidak muncul, berarti ini merupakan suatu persamaan
diferensial biasa ordo-kedua dengan koefisien c2w2 ”konstanta” (tidak tergantung pada t).
Solusi umumnya adalah
, cos sinu w t A w cwt B w cwt
Karena F , 0 , 0u x u w
, maka kita peroleh dari syarat awal, untuk
, 0u w A w f w
, 0 0tu w AcwB w
Jadi
, 0u w t f w
Ingat bahwa kita ingin mencari u. Ini perlu akal. Untuk ini kita ucapkan kosinus itu dalam
fungsi eksponensial dan kemudian kita terapkan rumus pergeseran
(11) F i ax a e
F .f x
yang merupakan akibat langsung dari definisi transformasi Fourier dengan mengambil
, , , sehinggax a p x p a dx dp
F 1 1
2 2
i p ai xf x a f p e dx f p e dp
11
.2
a i pe f p e dp
1
Jadi, ( ) cos ( )2
ic t ic tf w cwt f w e e
mempunyai transformasi Fourier kebalikan
(12) 1
, .2
u x t f x ct f x ct
Ini tidak lain adalah adalah solusi d'Alembert (6), Pasal 11.4
Ini akhir Bab 11. di dalam Bab ini kita telah memusatkan perhatian pada
persamaan-persamaan diferensial parsia paling penting di dalam fisika dan rekayasanya
sekaligus, ini juga merupakan akhir bagian C tentang analisa Fourier dan persamaan
diferensial parsial. Di dalam bagian D (Bab 12-Bab 17), kita akan menjelajahi suatu
cabang yang sama sekali berbeda sifatnya, analisis kompleks, yang juga sangat penting
bagi rekayasawan, sebagaian akan diperlihatkan oleh teladan-teladan dan soal-soal
latihan yang akan diberikan.
Soal-soal Latihan Ulangan untuk Bab 11.4
1. Apakah prinsip superposisi itu? Pada jeis persamaan bagaimanakah prinsip ini
berlaku?
2. Hukum fisik apakah yang menghasilkan persamaan gelombang, dan mengapakah
persamaan ini melibatkan turunan persiak kedua terhadap t? Turunan terhadap t
apakah yang apakah yang terdapat di dalam persamaan panas?
3. Verifikasi bahwa 2 2 34 dan 3u x t u x xt merupakan solusi bagi persamaan
gelombang bedemensi-satu dengan c yang sesuai.4. Apa yang dimaksud dengan solusi d'Alembert bagi persamaan gelombang?5. Apakah fungsi-eigen bagi dawai yang bervibrasi itu? Bagi membran yang bervibrasi?6. Di dalam pemisaan persamaan gelombang ber-demensisatu, kita hanyan memperole fungsi trigonometrik, sedangkan di dalam pemisahan persamaan gelombang kita juga memperoleh satu persamaan eksponensial. aaaaaaApak alasannya?
7. Bagaimanakah fungsi eksponensial di dalam Soal 6 akan berubah jika c2 di dalam
persamaan itu ditingkatkan? Apakah jawaban anda dapat dipahami dari sudut
pandang fisik?
8. Bagaimana mungkin beberapa syarat awal dikenakan pada kasus persamaan
gelombang? Di dalam kasus persamaan panas?
9. Verifikasi bahwa 4 2 2 46 dan sin sinhu x x y y u x y merupakan solusi bagi
persamaan Laplace.
10.Verifiaksi bahwa 2 22 2 2 2 2 2/ dan 2 / sin sinhu x y x y u x x y x y
merupakan solusi bagi persamaan Laplace.
11. Verifikasi bahwa 2 2 2u x y z merupakan solusi bagi persamaan Laplace.
12. Verifikasi 2 2 2u x y z memenuhi persamaan Possion.
13. Apa yang dimaksud dengan persamaan eliptik? Persamaan elptik manakah yang paling penting?
14. Apakah persamaan hiperbolik itu? Persamaan hiperbolik? Berikan teladan.
15. Berapa persamaan diferensial yang sederhana dapat dipecahkan dengan cara yang
biasa digunakan untuk persamaan diferensial biasa. Jelaskan. Berikan teladan.
16. Apakah fungsi galat itu? Dalam kaitan dengan apa fungsi ini muncul dalam bab ini?
17. Mengapa fungsi Bessel dan polinom Lagendre muncul di dalam bab ini?
18. Mengapa deret Fourier memainkan peranan pokok di dalam bab ini, walaupun fungsi
yang diberikan yang secara fisik menarik tidak bersifat periodik?
19. Atas pertimbangan apakah integral Fourier muncul di dalam bab ini?
20. Di dalam bab 5, persamaan pembantunya merupakan merupakan persamaan aljabar.
Mengapa kita memperoleh memperoleh persamaan diferensial biasa ketika
memecahkan persamaan diferensial parsial melalui metode transform (metode
oprasonal)?
21. Pecahkan 9 0.
22. Pecahkan 1 0.
23. Pecahkan 3 4 8.
xx
xy x
yy y
u u
u u x y
u u u
24. Carilah semua solusi ( , ) ( ) ( )u x y F x G y bagi persamaan Laplace dan dua peubah.
25. Tentukan semua solusi bagi x yyu xu melalui pemisah peubah.
Tentukan gerak dawai yang bervibrasi yang panjangnya 2,dengan / 1,c T p yang
mulai dengan kecepatan awal 0 dan defleksi
3
126. ( ) sin sin 2 27. ( ) 0.1sin3
2
28. ( ) / 2 / 2 29. ( ) sin
30. ( ) jika 0 / 3, ( ) ( ) / 2 jika / 3
f x x x f x x
f x x f x x
f x x x f x x x
Osilasi paksa terhadap suatu dawai elastis. Persamaan adalah
(1)
2 .xx
Pu c u
Simak dan pecahkan persamaan ini sebagai berikut. 31. Tunjukkan bahwa vibrasai paksa terhadap suatu dawai elastis mengikuti persamaan
(1), dengan dalam hal ini ( , )P x t adalah gaya eksternal persamaan satuan panjang yang bekerja tegak lurus pada dawai tersebut.32. Misalkan gaya eksternal itu sinusiodal, katakanlah sin .P A t Tunjukkan bahwa
1
/ sin ( )sinnn
n xP A t k t
L
dengan dalam hal ini ( ) (2 / )(1 cos )sin ;nk t A n n t akibatnya 0nk (untuk n
genap), dan ( ) (4 / )(1 cos )sin ;nk t A n n t (unutk n ganjil). Lebih lanjut, tunjukkan bahwa substitusi
••2 2
1
( , ) sin ke dalam menghasilkan 0n tt tt n nn
n xu x t G u c u G G
L
dengan
dalam hal ini .cn L
33. Tunjukkan bahwa dengan mensubstitusiakn ekspresi bagi dan /u P p di dalam Soal 32 ke dalam (1) kita memperoleh
2 2(1 cos )sin , .n n n n
A cnG G n t
n L
Tunjukkan bahwa jika 2 ,n solusinya adalah
2 2
2 (1 cos )( ) cos sin sin .
( )nn n n n nn n
n
A nG t B B G t t
n
34. Tunjukkan bahwa dalam kasus resonansi ( ),n
( ) cos sin (1 cos ) cos .n n n n
AG t B t B n t t
n
Dengan menggunakan transformasi yang ditunjukkan, pecahkan persamaan-persamaan
diferensial parsia berikut.
36. 2 0 ( , )
37. 2 0 ( , )
38. 2 5 2 0 ( 2 , 2 )
39. ( , )
40. 2 ( , 2 )
xx xy yy
xx xy yy
xx xy yy
xy yy y
xx xy yy
u u u v x z x y
u u u v x z x y
u u u v x y z y x
xu yu u v x z xy
u u u v x y z x y
Tentukan distribusi suhu di dalam batang tembaga yang diisolasi memanjang 2(c
2/ 1.158cm / det),K yang panjangnya 50 cm dan irisan melintannya constan, yang
ujung-ujungnya pada x = 0 dan x = 50 dipertahankan pada suhu 0 oC dan yang suhu awalnya adalah
3
41. ( ) sin( / 50) 42. ( ) 100sin( / 250)
43. ( ) sin ( /10) 44. ( ) (50 )
45. ( ) jika 0 25, ( ) 50 jika 25 50
f x x f x x
f x x f x x x
f x x x f x x x
Kita ingat dari Pasal 11.5 bahwa untuk syarat batas adiabatik, suhu ( , )u x t di dalam sebuah batang yang diisolasi memanjang dengan panjang L adalah
2
01
( , ) cos exp .nn
n x c xu x t A A t
L L
Jika dan 1,L c tentukan dari ini sousi yang memenuhi syarat awal 2
10 jika 0<x< /2 4x jika o<x< /2
0 jika /2<x< 4( ) jika /2<x<
46. ( ) 30 47. ( ) 95cos 2
47. ( ) 49. ( ) x
f x x f x x
f x f x
50. Dengan menggunakan tabel A2 di dalam Aapendiks 4 dan interpolasi linear dan koefisien-koefisien yang duberikan pada Teladan 1 Pasal 11.10, tunjukkan bahwa
(0.5) 0.7498.f beraapa nilai pastinya?
Tunjukkan bahwa membran berikut, yang luasnya 1 dan c2 = 1, mempunyai frekuensi
modus dasar (fundamental mode) sebagai berikut (dalam 4-desimal). Bandingkanlah.
51. Lingkaran: 1 / 2 0.6784
52. Bujursangkar: 1 2 0.7071
53. Seperempat lingkaran: 12 12/ 4 0.7244 ( 5.135652 nilai nol ( )zero positif yang pertama dari J2)
54. Setengah lingkaran: 3.832 / 8 0.7644
55. Persegi panjang (sisi-sisi 1 : 2): 5/8 0.7906
Ringkasan Bab 11
Persamaan Diferensial Parsial
Bila persamaan diferensial parsial biasa (Bab 1-Bab 5) merupakan model bagi masalah
rekayasa sederhana yang melibatkan suatu peubah bebas, masalah yang melibatkan dua
atau lebih peubah bebas (peubah ruang, atau waktu t dan satu atau lebih peubah ruangan)
menghasilkan persamaan diferensial parsial. Ini menunjukkkan bahwa pentingnya
persamaan-persamaan itu bagi rekayasawan dan fisikawan tidak dapat diabaikan.
Di dalam bab ini kita terutama mengupas persamaan-persamaan diferensial
parsial paling penting di dalam fisika dan rekayasa, yaitu:
(1) 2
u xxu c u Persamaan gelombang berdimensi-satu (Pasal 11.2-Pasal 11.4)
(2) 2 u xx yyu c u u
Persamaan gelombang berdemens-dua (Pasal 11.7-Pasal 11.10)
(3) 2
t xxu c u Persamaan panas bedemensi-satu (Pasal 11.5, Pasal 11.6)
(4) 0xx yyu u u
Persamaan Laplace berdemensi-dua (pasl 11.5, Pasal 11.9)
(5)2 0xx yy zzu u u u
Persamaan Laplace berdemensi-tiga (Pasal 11.11, Pasal 11.12)
(1) dan (2) adalah hiperbolik, (3) adalah parabolik, (4) dan (5) adalah elptik. (Lihat Soal-
soal Latihan untuk Pasal 11.4.)
Di dalam praktek, biasnya kita ingin memperoleh solusi persaman demikian ini
pada suatu daerah tertentu yang memenuhi syarat-syarat tambahan, misalnya syarat awal
(syarat pada waktu t = 0) atau syarat batas (nilai tertentu bagi solusi u atau turunannya
pada permukaan batas S, atau kurva batas C dari darah di atas) atau keduanya. Pada (1)
dan (2) kita menetapkan dua syarat awal (pergeseran awal dan kecepatan awal). Pada (3),
kita menetapkan distribusi suhu awal. Pada (4) atau (5) kita menetapkan suatu syarat
batas dan menamakan masalah tersebut
Masalah Dirichlet jika u ditetapkan pada S,
Masalah Neumann /nu u n ditetapkan pada S,Masalah campuran jika u ditetapkan pada sebagian lainnya dari S sedangkan un pada
sebagian lainnya dari S.
Lihat Pasal 11.5.
Suatu metode umum untuk memecahkan masalah demikian ini adalah metode
pemisahan peuba atau metode hasil kali. Di dalam metode ini kita mengasumsikan
solusinya berbentuk hasilkali dua fungsi, yang masing-masing tergantung pada satu
peubah saja. jadi, persamaan (1) dipecahkan dengan mengambil (Pasal 11.3)
, ;u x t F x G tbegitu pula untuk (3) dalam Pasal 11.5 substitusi ke dalam persaman yang diberikan
menghasilkan persaman diferensial biasa bagi F dan G, dan dari situ kita memperoleh tak
hingga banyaknya solusi F = dan G = Gn sedemikian rupa sehingga fungsi
, ;n nu x t F x G tmerupakan solusi bagi persamaan diferensial parsial di atas yang memenuhi syarat batas
yang diberikan. Itu adalah fungsi-eigen bagi masalah tersebut, sedangkan nilaiaigen-nya
menentukan frekuensi vibrasinya (atau kecepatan turunnya suhu dalam kasusu persamaan
panas, dan lain sebagainya). Agar juga memenuhi syarat awal, kita harus menyimak deret
takhinggan bagi un, yang koefisien-koefisiennya ternyata adalah koefisien Fourier fungsi
f dan g yang merepresentasikan syarat awal yang diberikan (Pasal-Pasal 11.3, 11.5). oleh
karena itu, deret Fourier (dan integral Fourier) sangat penting di sini (Pasal-Pasal 11.3,
11.5, 11.6, 11.8)
Masalah keadaan-stabil atau stasioner adalah masalah yang solusinya tidak
tergantung pada waktu t. Untuk masalah demikian ini, persamaan panas ut = c22u
menjadi persamaan Laplace (Pasal 11.5).
Sebelum memecahkan suatu masalah nilai awal atau masalah nilai batas, kita
sering mentransformasi persamaan itu ke dalam koordinat-koordinat sehingga batas
daerahnya dapat diucapkan dalam rumus yang sederhana. Dalam koordinat kutub yang
didefinisikan oleh x = r cos q, y = r sin q, Laplacian menjadi (pasal 11.9)
(6)
22
1 1rr r mru u u u
r r
Sedangkan dalam koordinat bola, lihat pasal 11.11. kalau kemudian kita terapkan metode
pemisahan peubah, maka akan kita peroleh persamaan Bessel dari (2) dan (6) (membran
lingkaran bervinrasi, pasal 11.10) dan persaman Legendre dari (5) yang ditransformasi
ke dalam koordinat bola (pasal 11.12).
Metode operasional (trnformasi Laplace, transformasi fourier)sangat membantu
untuk memecahkan persaman diferensial parsial dalam daerah tak hingga (pasal-pasal
11.12, 11.13), khususnya jika koefisien-koefisien persamaan itu adalah konstanta atau
tergantung hanya pada satu peubah.
110 jika 0 / 2
0 jika /2< .
f
sin jika 0 t 2
22
0 selainnya
1, ( )
2 2
t
xw x t t u t x
Publish By Awang, Himpunan mahasiswa Fisika (Himafi) Universitas Hasanuddin Makassar.
top related