Teknillinen korkeakoulu
Akustiikka ja
äänenkäsittelytekniikka
S-89.112 Akustiikan ja
äänenkäsittelytekniikan työt, syksy 1996
Panu Maijala, 37847W
24.1.1997
Erikoistyön tarkoituksena oli selvittää TKK:n akustiikan laboratorion seisovan aallon putken käyttökelpoisuus mittauskäyttöön ja tehdä putkeen tiettyjä parannuksia, joiden myötä putken ominaisuudet saataisiin mahdollisimman hyvin vastaamaan ISO/DIS 10534 -standardiehdotusta. Olennainen osa tätä erikoistyötä muodostui itse putken lisäksi siihen liittyvän mittausjärjestelmän osien valinnasta, suunnittelusta, asennuksesta ja mittauksesta.
Erikoistyön laajuus aiheutti ongelmia tehtävän rajaamisessa. Tämä näkyi erityisesti työn hitaana etenemisenä, sillä uusia huomioitavia seikkoja tuli jatkuvasti esiin työn aikana. Standardiehdotuksessa olevien suositusten osalta putken kaikkia ominaisuuksia ei edes päästy tutkimaan erikoistyön puitteissa.
Erikoistyön toteutuksen vaiheet on kuvattu lyhyesti seuraavassa. Erikoistyön laajuuden takia osa tehdystä työstä ei suoraan liity putkimittauksiin, vaan esimerkiksi mittausjärjestelmän suunnitteluun ja rakentamiseen. Selvitys näistä on sijoitettu selostuksen liitteisiin. On kuitenkin syytä huomata, että kyseiset liitteet ovat vaatineet merkittävän osan erikoistyöhön liittyvästä työmäärästä.
Työssä lähdettiin liikkeelle tekemällä sarja mittauksia alkuperäisessä tilassa olevalla putkella. Mittausten tarkoituksena oli selvittää mittaussignaalin suodattamisen vaikutuksia järjestelmän tuottamiin tuloksiin.
Mittausjärjestelmä koostui putken lisäksi mikrofonivahvistimesta ja kahdesta erillisestä suodinyksiköstä, joita käytettiin heräte- ja paluusignaalin käsittelyyn. Tarkoituksena oli korvata nämä erillislaitteet varta vasten suunnitellulla suodatin-vahvistin-yhdistelmällä, joka jää pysyvästi käyttöön putken kanssa.
Ensimmäisestä mittaussarjasta laskettujen tulosten pohjalta oli tarkoitus suunnitella erillinen suodattimet ja mikrofonivahvistimen sisältävä yksikkö putken mittausjärjestelmää varten.
Ensimmäinen mittaussarja ei kuitenkaan tuottanut niin selkeitä tuloksia, että suodattimien valinta olisi voitu määrätä yksikäsitteisesti niihin pohjautuen. Suodattimien rakenne ja arvot valittiin lopulta työn ohjaajalta saatujen, toisaalla hyväksi koettujen arvojen perusteella. Tästä syystä ensimmäiset mittaustulokset ja suodattimille valitut arvot eivät ole suoraan yhteydessä toisiinsa.
Erikoistyön yhtenä tavoitteena oli saavuttaa putken ja mittausjärjestelmän osalta sellaiset ominaisuudet, jotka mahdollisimman suurelta osin täyttäisivät ISO/DIS 10534 -standardiehdotuksissa määritellyt ominaisuudet.
Mittausjärjestelmään tehtiin muun muassa seuraavat muutostyöt:
Putkeen vaihdettiin laadukkaammat kaiutinelementit ja mikrofonikapseli, ja niiden kiinnitys suunniteltiin ja toteutettiin uudelleen. Mikrofonin sijainnin mittausjärjestelyyn suunniteltiin muutoksia. Putken rakenteelliset puutteet, kuten tiiviysongelmat, hitsaussaumojen raot ja näytepalan kiinnitysmahdollisuudet tutkittiin, ja niille pyrittiin löytämään ja toteuttamaan paremmat ratkaisut.
Uusien kaiutinelementtien THD-arvot mitattiin IEC-standardin 268-5 mukaan.
Lisäksi mittausjärjestelmään suunniteltiin ja toteutettiin laajahko erillinen suodatin- ja mikrofonivahvistinosa. Selvitimme myös muita mittausjärjestelmään liittyviä laitteistohankintoja, kuten tarkkojen ilmanpaine- ja lämpömittarien hankintaa.
Tämän lisäksi kirjoitettiin mittaustulosten käsittelyä ja esittämistä varten tarvittava Matlab-koodi. Tämä sisälsi myös MLLSA-järjestelmässä tallennettujen TIM-näytetiedostojen Matlab-ympäristöön tuomiseen ja viemiseen tarvittavan koodin. Erikoistyön tulosten laskentaan käytetyt Matlab-ohjelmat ovat työselostuksen liitteissä.
Toinen mittaussarja tehtiin sen jälkeen, kun putken muutostyöt ja suodatin-vahvistin-yksikkö saatiin valmiiksi. Koepaloina käytettiin seuraavia akustiikkalevyjä: Duratex, Partec EL100 ja Plain. Koepalat tiivistettiin putken päähän silikonimassalla ja ne mitattiin sekä sini- että MLS-herätteillä. Mittausjärjestelmän yksityiskohtainen kokoonpano on eritelty luvussa 3.
Toinen mittaussarja jouduttiin kertaalleen uusimaan säröytyneen signaalin takia. Ongelma havaittiin vasta tulosten käsittelyssä, ja virhelähteeksi paikannettiin kaiutinelementit, jotka olivat tuottaneet säröä herätesignaaliin. Virheen aiheutti se, ettei kaiuttimeen syötettävää tehoa ollut tarkistettu ennen mittausten aloittamista. Mittaussarja uusittiin, ja tulokset paranivat huomattavasti säröytyneisiin verrattuna. Samalla kuitenkin todettiin, että uusittu järjestelmä on edelleen sangen altis ylimääräisten virhetermien syntymiselle.
Tulosten hyvyyden selvittämiseksi samat koepalat (Duratex, EL100 ja Plain) mitattiin kahdella VTT:n järjestelmällä. Otaniemessä (VTT Rakennustekniikka/Akustiikka) mittaukset tehtiin B&K 4002 -laitteistolla käyttäen sini- ja MLS-herätteitä. Tampereella (VTT Valmistustekniikka/Turvallisuustekniikka) mittaukset tehtiin automatisoidulla B&K 4206 -laitteistolla, joka tuotti valmiit tulokset. Tuloksia vertaillaan selostuksen kohdassa .
Kompleksinen heijastuskerroin R määritellään kahden väliaineen rajapintaan tulevan ja siitä heijastuvan (kompleksisen) paineen summaksi. Heijastussuhteella r tarkoitetaan kompleksisen heijastuskertoimen itseisarvon neliötä.
Absorptiosuhde määritellään pintaan saapuneen ja pinnan absorboiman äänitehon suhteena. Absorptiosuhde voidaan ratkaista heijastuskertoimesta seuraavasti: = 1 - R2.
Pinnan kohtisuora ominaisimpedanssi määritellään paineen ja hiukkasnopeuden kohtisuoran komponentin suhteena: Zn(f) = p/un. Normalisoitu impedanssi saadaan normalisoimalla tämä ilman karakteristiseen impedanssiin: z(f) = Z(f)/c. Tässä erikoistyöselostuksessa käytetään lyhyesti termejä impedanssi ja normalisoitu impedanssi.
Seisova aalto on yksiulotteiseen tasoaaltokenttään syntyvä stabiili aaltomuoto, jonka solmut ja huiput säilyvät paikallaan ajan funktiona. SAS eli seisovan aallon suhde voidaan määrittää seisovan aallon painemaksimin ja -minimin itseisarvojen suhteena. Seisovan aallon putki on puolestaan rakenne, johon yksiulotteinen seisova aalto synnytetään putken päätteenä olevan materiaalin tutkimista varten.
MLS on menetelmä, jolla voidaan mitata suoraan järjestelmän impulssivaste käyttäen näennäissatunnaista herätesignaalia.
Seisovan aallon putki perustuu äänikentän yksiulotteiseen malliin. Putkeen syötetään herätesignaali, joka muodostaa seisovista aalloista koostuvan tasoaaltokentän heijastuessaan putken päästä. Mitattavat kappaleet sijoitetaan putken päätteeksi, jolloin ne vaikuttavat muodostuvaan äänikenttään. Äänikenttä mitataan kahdesta tai useammasta putken pituusakselilla olevasta pisteestä. Näiden sijainti ja painesignaaleista saatu informaatio riittävät putken sisällä olevan äänikentän määrittämiseen yksikäsitteisesti.
Yksiulotteisen seisovan aallon äänikenttä koostuu tulevasta ja palaavasta aallosta:
,
()
jossa x on sijainti putken pituusakselilla. Putken päästä heijastunut palaava aalto vaimentuu heijastuskertoimen R verran tulevaan aaltoon nähden:
;
()
Äänikenttää mitataan kahdella painemikrofonilla, jotka sijoitetaan putkeen etäisyydelle d toisistaan. Koska yksiulotteisessa tasoaaltokentässä taajuusvastefunktiot tunnetaan analyyttisesti
, ()
voidaan mikrofonien painesignaalien välisestä taajuusvastefunktiosta H12(f) määrittää heijastus- ja absorptiosuhde sekä impedanssi.
Kompleksinen heijastuskerroin mikrofonin 1 kohdalla on
. ()
Heijastuskerrointa voidaan laskennallisesti siirtää putkessa vapaavalintaiseen kohtaan x, esimerkiksi tutkittavan materiaalin pintaan.
()
Absorptiosuhde saadaan seuraavasti:
.
()
Heijastuskertoimesta voidaan laskea myös impedanssi, jonka arvoksi saadaan
. ()
Impedanssi voidaan normalisoida kertoimella 1/c.
Putken pituus L määrää alimman mahdollisen moodin eli menetelmän alarajataajuuden
.
()
Yksiulotteisen äänikentän ylärajataajuus fU määräytyy vastaavasti putken halkaisijan, tai suorakulmaisen putken tapauksessa poikkipinnan pidemmän reunan W mukaan:
. ()
Tätä suuremmilla taajuuksilla putkessa syntyy pitkittäisten moodien lisäksi myös poikittaisia moodeja, jolloin äänikenttää ei voida enää käsitellä yksiulotteisen mallin avulla.
SAS-mittaustuloksista laskettiin heijastussuhteet noudattaen ISO/DIS 10534 -standardiehdotuksessa annettuja ohjeita.
Mittauksissa määritettiin ensimmäisen paineminimin ja -maksimin amplitudit ja paikat putkessa. Nämä tiedot riittävät seisovan aallon suhteen määrittämiseen, josta voidaan puolestaan laskea heijastus- ja absorptiosuhde.
Seisovan aallon suhde taajuudella f saadaan kaavasta
. ()
Tästä saadaan suoraan kompleksisen heijastuskertoimen itseisarvo
. ()
Kun tunnetaan lisäksi ensimmäisen minimin etäisyys putken päästä (tutkittavasta kappaleesta), vaihe voidaan laskea.
()
Tulokseksi saadaan kompleksisen heijastuskertoimen reaali- ja imaginaariosat:
. ()
Absorptiosuhde saadaan heijastuskertoimesta samoin kuin edellä.
Seuraavana esitetty laskenta suoritetaan Matlab-rutiinien ekaputki.m ja sasputki.m avulla. Mitattujen impulssivasteiden aikatasossa olevat näytejonot luetaan MLSSA-ohjelmassa tallennetuista TIM-tiedostoista Matlabiin timread.m-makron avulla.
Matlab-makrojen listaukset ovat työselostuksen liitteissä.
Matlabissa näytejonoista lasketaan energia-autospektrit G11(f) ja energiaristispektrit G12(f) käyttämällä PSD (Power Spectral Density)- ja CSD (Cross Spectral Density)-funktioita, jotka löytyvät Matlabin Signal Processing Toolkitistä. Funktioille annetaan argumenteiksi aikatason näytejono tai näytejonot, näytejonon pituus (16 383), näytteenottotaajuus (4000) sekä taajuuspehmennykseen käytettävä ikkunointifunktio (256 pisteen Hanning-ikkuna).
Lasketuista auto- ja ristispektreistä muodostetaan taajuusvastefunktiot H12(f). Tarkasteltaessa heijastuskertoimien käyriä havaittiin, että osassa mittauksia vaste saatiin käyttämällä taajuusvastefunktiota H'12(f). On huomattava, että kaavat poikkeavat toisistaan vain mallin tarkastelusuunnan tai -järjestyksen osalta. [Lahti: Akustinen mittaustekniikka, s. 97]
. ()
Seuraavaksi Matlab-makrossa rajataan saaduista taajuusvastevektoreista haluttua ylärajataajuutta korkeammat taajuudet pois. Tämä nopeuttaa tulosten jatkokäsittelyä ja rajaa tulostettavien käyrien taajuusalueen halutuksi.
Tämän jälkeen muodostetaan aaltolukuvektori, joka saadaan seuraavasti pistetaajuusvektorista f :
()
Kompleksiset heijastuskertoimet ensimmäisissä mittauspisteissä saadaan kaavasta
,
()
jossa d on mittauspisteiden välinen etäisyys (0,15 m kaikissa mittauspareissa).
Koska TKK:n putkella mittaukset tehtiin varmuuden vuoksi kolmesta pisteestä, heijastuskertoimet lasketaan sekä ensimmäisen ja toisen että toisen ja kolmannen mittauspisteen väliltä.
Heijastuskertoimet siirretään materiaalin pintaan käyttäen kaavaa
,
()
jossa x on lähemmän mittauspisteen ja materiaalin pinnan välinen etäisyys.
Heijastuskertoimista oktaavikaistoittain lasketut absorptiosuhteet on taulukoitu kohdassa . Mittaustuloksista lasketut heijastuskertoimet löytyvät myös käyrinä selostuksen liitteistä.
Impedanssin lausekkeeksi saadaan [Lahti: Akustinen mittaustekniikka, s.122]
, ()
jossa d on mittauspisteiden välinen etäisyys (0,15 m kaikissa mittauspareissa).
Oktaavikaistoittain lasketut impedanssit on taulukoitu kohdassa . Mittaustuloksista lasketut impedanssikäyrät ovat selostuksen liitteissä.
Jotta impulssivasteen luotettavuus voidaan varmistaa, mittauksissa tulee huomioida eri mittauspisteistä saatujen näytejonojen keskinäiset koherenssit. MLS-herätteellä tehty mittaus tuottaa jo sinällään keskiarvoistetun impulssivasteen, joten yksi mittaussarja riittää tuottamaan tarvittavan tarkkuuden. MLS-menetelmän käyttö puolustaa paikkaansa, sillä N-pituisen MLS-jonon autospektri on (N+1)-kertainen verrattuna yksikköimpulssijonoon, jolla on sama amplitudi [Lahti: Akustinen mittaustekniikka, s. 137]. Tästä seuraa voimakas parannus signaali-kohinasuhteeseen, joten myös mittausten koherenssin pitäisi olla verraten hyvä.
Koska putken vaimennuskerrointa ei tunnettu eikä sitä ryhdytty erikseen määrittämään, laskuissa käytetään aaltoluvulle k arvoa [omega]0 / c0 . Väliaineen vaimennuksen huomioiminen vaikeuttaisi laskuja myös siksi, että tällöin kaava 1 ei enää olisi vain etäisyyden x funktio, vaan äänen etenemisnopeudessa pitäisi huomioida myös väliaineen kompleksinen tiheys ja puristuvuus [Toivanen: Teknillinen akustiikka, s. 46]. Laskuissa oletetaan siis häviötön väliaine.
Tärkeimmät heijastuskertoimen, akustisen impedanssin ja intensiteetin mittaustavat putkimenetelmän ohella ovat korvausmenetelmä, ikkunointimenetelmä ja erilaiset yksi- ja kaksikanavaiset intensiteettimenetelmät.
Korvausmenetelmässä mittaus suoritetaan ensin ilman mitattavaa pintaa joko vapaassa kentässä tai äärettömäksi katsottavalla pinnalla. Tutkittava pinta tuodaan äänikenttään ja mittaus uusitaan kappaleen pinnalla. Menetelmä tuottaa yleensä vain pinnan reaalisen tehoheijastussuhteen.
Ikkunointimenetelmässä tuleva ja heijastunut ääniaalto mitataan yhtäaikaa, ja rajataan aikatasossa toisistaan erilleen. Menetelmä vaatii suurehkon tilan, jotta suorasta ja heijastuneesta signaalista saadaan kyllin pitkät näytteet ennen niiden sekoittumista muihin heijastuksiin. Ikkunointimenetelmällä lasketaan pinnan tuottama taajuusvastefunktio, josta voidaan määrittää heijastuskerroin. Yleensä taajuusresoluutio jää kuitenkin heikoksi mittaustilojen rajallisen koon vuoksi.
Intensiteettimenetelmät perustuvat äänikentän paineen mittaamiseen samanaikaisesti kahdessa pisteessä, jolloin painesignaaleista voidaan muodostaa hiukkasnopeuden approksimaatio. Näistä voidaan muodostaa suoraan äänen kompleksinen intensiteetti ja kaikki muut halutut suureet. Intensiteettimenetelmä tuottaa muihin menetelmiin verrattuna tarkimmat tulokset, mutta vaatii suuren tarkkuuden lisäksi kaksikanavaisen mittausjärjestelmän tai erikoismenetelmiä yksikanavaisissa mittauksissa.
Alla on eritelty joitakin standardiehdotuksessa mainittuja keskeisiä vaatimuksia, jotka pyrittiin huomioimaan putkeen tehdyissä muutostöissä.
ISO/DIS 10534, kohta 5: Putken halkaisijan tulee olla vakio 0,2 % tarkkuudella. Seinämissä ei saa olla reikiä eikä rakoja mittausalueella, ja niiden tulee olla jäykät.
ISO/DIS 10534, kohta 5.1: Poikkileikkaukseltaan suorakulmaisissa putkissa kulmaliitosten tulee olla riittävän jäykät, jotta sivuseinämien muoto ei pääse vääristymään. Sivuseinämien suositeltava paksuus on noin 10 % putken halkaisijasta.
ISO/DIS 10534, kohta 5.5: Siirreltävän mikrofonin sijainti on voitava mitata vähintään 0,5 mm tarkkuudella.
ISO/DIS 10534, kohta A.2.1.1: Jos näytekappaleen pinnan ja lähimmän mikrofonin välinen etäisyys on suurempi kuin kolme kertaa putken sisähalkaisija, se täytyy huomioida tulosten laskennassa.
ISO/DIS 10534, kohta 5.9: Kaiutinelementin kalvon tulee peittää pinta-alaltaan vähintään 2/3 putken poikkileikkauksesta.
Erikoistyön mittaukset tehtiin TKK:n akustiikan laboratorion suurella seisovan aallon putkella, jonka sisähalkaisija on noin 250 x 250 mm ja pituus 8820 mm. Putken alarajataajuus on siis fL = 20 Hz ja ylärajataajuus fU = 690 Hz. Mittaustulosten käsittelyssä rajoitutaan tälle taajuusalueelle.


HP / LP -suodatin
|
KEMO
Dual Variable Filter, type VBF/4
|
| PINK-suodatin
|
B&K
Weighting Network, type WB0078
|
| Päätevahvistin
|
B&K
Power Amplifier, type 2706
|
| Mikrofonivahvistin
|
B&K
Microphone Amplifier, type 2636
|
Seuraavassa on esitetty järjestelmään kohdistuneet muutostyöt suhteessa ISO-standardiehdotuksessa asetettuihin vaatimuksiin.
Putki
* TKK:n putken sisämitat täyttävät 0,2 % tarkkuusvaatimuksen. Ehto aiheutti tiivisteen uusimisen näytekappaleen vaihdossa käytettävään "luukkuun" sekä putken kappaleiden välitiivisteen uusimisen. Putken näytepäässä havaittiin 10 mm:n reikä, joka tukittiin.
* Putken tiiviys testattiin paineilmalla. Tätä varten valmistettiin erillinen painetta pitävä päätylevy, jossa oli manometri kiinni. Putki piti yhden ilmakehän painetta kohtuullisen ajan. (Osoitin ei liikahtanut minuutin aikana.)
* Näytteen asennusaukkoon tulee kiinnittää erillinen tiivistyspanta, sillä putken suunnittelussa kyseisen aukon tiiviyttä ei ole otettu lainkaan huomioon.
* Putken seinämät eivät täytä standardin suorakulmaputkelle asettamaa 10 prosentin paksuusehtoa. Seinämien jäykkyys on siis standardiin verrattuna kyseenalainen. Olsiko vaihtoehtona siis koko putken upottaminen hiekkalaatikkoon??




Mikrofoni
* Mikrofonin etäisyyden määritys 0,5 mm:n tarkkuudella on kyseenalaista. Arviolta 3,5 mm:n tarkkuus on toistettavissa nykyisellä varustuksella. Mikrofonin ripustuksessa on hystereesiä varsin paljon, sen saa pysymään eri paikoissa 7 mm:n matkalla siirtotankoon nähden.
* Mikrofonin ripustuksen tukevuutta ei määritetty. On mahdollista, että siirtotankoa pitkin tulee värähtelyä mikrofoniin saakka.
* Etäisyys näytteeseen mikrofonin etäisimmässä mittauspisteessä pitää olla alle kolme putken halkaisijaa tai kolme kertaa suorakulmaputken suurin poikittainen ulottuvuus. Muussa tapauksessa täytyy määrittää putken vaimennuskerroin ISO/DIS 10534-2 liitteen A.2.1.2 mukaan.
* Mikrofonin siirtelysauva vaihdettu pidempään, sillä vanha ei ylettynyt päätyyn saakka!
Kaiutin
* Putken kaiutinelementeiksi vaihdettiin kaksi kappaletta 8 ohmin bassoelementtejä vastakkaisiin vaiheisiin rinnankytkettynä.
(19)
Kuten kaavasta voidaan havaita, kaiuttimen kalvon pinta-ala ei ole 2/3 putken poikkileikkauksen pinta-alasta, vaan huomattavasti vähemmän. Toisaalta rinnankytketty elementtipari antaa enemmän tehoa kuin yksittäinen elementti.
Järjestelmään suunniteltiin kaistanpäästösuotimet sekä herätteelle, että mikrofonilta palaavalle mittaussignaalille. Kummankin suotimen -3 dB rajataajuudet ovat 15 ja 800 Hz. Suunnittelussa käytettiin hyväksi Aplac-piirisimulointiohjelman optimointiominaisuuksia.
Herätteelle haluttiin kaistanpäästösuodatuksen
lisäksi -3 dB/okt laskeva suodatus. Tällä oli tarkoitus
varmistaa, että signaali-kohinasuhde pysyy järkevänä
sillä taajuusalueella, missä sitä tarvitaan. Optimointitavoite
oli että 15 Hz:n ja 800 Hz:n kohdalla tulisi siirtofunktion olla
-3
dB:ssä ja välillä 20-600 Hz tulisi pysyä 0.5 dB:n
sisällä.

GoalData Gain
+TYPE FIRST LAST PTS MIN MAX
+LOG 1 15 8 -100 -3
+LOG 20 600 15 -0.5 0.5
+LOG 800 3.2k 10 -100 -3





Mittaussignaali haluttiin myös kaistarajoittaa optimoidun tuloksen saamiseksi. Tavoitteeksi asetettiin samat -3 dB:n rajataajuudet kuin herätesignaalin suotimelle, eli - 15 Hz ja 800 Hz.





Seuraavassa kuvataan uusittu mittausjärjestelmä ja eritellään muutokset vanhaan järjestelmään verrattuna pääpiirteittäin.
Seuraavassa kuvassa on esitetty putken mittoihin tapahtuneet muutokset.

Putken sisätilan efektiivinen pituus päätyseinästä kaiutinkartioiden keskipisteeseen on 8963 mm.

Kaistanpäästösuodatin
|
Black
Box
|
|
PINK-suodatin
|
Black
Box
|
| Päätevahvistin
|
Yamaha
MX-70
|
| Mikrofonivahvistin
|
Black
Box
|
| Kaistanpäästösuodatin
|
Black
Box
|
Referenssijärjestelminä käytettiin kahta erilaista VTT:llä käytössä olevaa mittausjärjestelmää. Näillä mitattiin samat näytepalat (Duratex, Plain ja EL100), kuin TKK:n uusitulla järjestelmällä.
Otaniemessä (VTT Rakennustekniikka/Akustiikka) mittaukset tehtiin tutkija Raimo Euraston opastuksella B&K 4002 -laitteistolla, jossa on liikuteltava mikrofoni. Putkella voidaan määrittää akustinen impedanssi taajuusalueella 90 - 6500 Hz. Käytimme ainoastaan 10 cm:n halkaisijaltaan olevaa putkea, joten taajuusalue rajoittui yläpäästä 1800 Hz:iin. Käytössä oli lisäksi B&K 1022 -siniaaltogeneraattori. Amplitudin poikkeamat luettiin B&K 2203 -äänitasomittarilla, jossa B&K 1613 -oktaavisuodin rajoitti ei-tutkittavan alueen pois. Herätteinä käytettiin sekä normaaleja sinisignaaleja että MLS-signaalia.
Tampereella (VTT Valmistustekniikka/Turvallisuustekniikka) tutkijat Jukka Tanttari ja Kari Saarinen suorittivat mittaukset automatisoidulla B&K 4206 -laitteistolla. PC-pohjaisen mittausjärjestelmän mittausohjelmistona on B&K:n BZ 5050. Ohjelma laskee ja näyttää graafisesti absorptio- ja heijastuskertoimet sekä normalisoidun impedanssin testinäytteestä. Samoin taajuusvastefunktio voidaan tulostaa järjestelmän kahden mikrofonin välille.Järjestelmä käyttää kohinaherätettä.
Alkuperäisen järjestelmän mittauksiin käytettiin valitettavasti eri näytemateriaaleja kuin uuden järjestelmän mittauksiin. Alkuperäisellä järjestelmällä mitatut näytekappaleet ja mittauksissa käytetyt suodatinyhdistelmät on esitetty alla olevassa taulukossa.
Tutkittava kappale
|
HPF
|
LPF
|
Pink
|
Mic
LPF
| |
| 1
|
vuorivillakiila1
|
--
|
--
|
--
|
--
|
| 2
|
30
mm akustolevy
|
--
|
--
|
--
|
--
|
| 3
|
15
Hz
|
1
kHz
|
--
|
--
| |
| 4
|
50
mm akustolevy
|
--
|
--
|
--
|
--
|
| 5
|
15
Hz
|
1
kHz
|
--
|
--
| |
| 6
|
70
mm akustolevy
|
--
|
--
|
--
|
--
|
| 7
|
15
Hz
|
1
kHz
|
--
|
--
| |
| 8
|
70
mm akustolevy
|
50
Hz
|
1
kHz
|
--
|
--
|
| 9
|
50
Hz
|
700
Hz
|
--
|
--
| |
| 10
|
70
mm akustolevy
|
30
Hz
|
1
kHz
|
--
|
--
|
|
11
|
30
Hz
|
700
Hz
|
--
|
--
| |
| 12
|
70
mm akustolevy
|
15
Hz
|
1
kHz
|
-3
dB/okt.
|
--
|
| 13
|
15
Hz
|
1
kHz
|
-3
dB/okt.
|
22,4
Hz
| |
| 14
|
15
Hz
|
1
kHz
|
--
|
22,4
Hz
|
Alkuperäisen järjestelmän mittaustulokset ovat selostuksen liitteissä, koska niiden painoarvo ei osoittautunut suureksi työn kannalta. Tuloskäyrästöihin on laskettu näytemateriaalien kompleksiset heijastuskertoimet ja impedanssit.
Saadut tulokset vaikuttavat verraten yhtenäisiltä, ja sietävät vertailua ainakin keskenään. Myös tulosten kokoluokka vaikuttaa yleisesti ottaen järkevältä. Virhetermien läsnäolosta kertoo kuitenkin mm. heijastuskertoimen itseisarvossa paikoitellen näkyvät piikit, jotka saavat ykköstä suurempia arvoja. Tämä ei tietenkään ole fysikaalisesti perusteltavissa, vaan aiheutuu mittauksessa tapahtuneista virheistä.
Kaikki järjestelmän signaalitiellä olevat laitteet aiheuttavat omalta osaltaan systemaattisia virheitä. Suurimmat epälineaarisuudet esiintyivät varmasti mikrofonissa (hiilimikrofoni!), päätevahvistimessa ja kaiuttimessa. Lisäksi käyrästöistä saattoi helposti havaita, että suodattimet aiheuttavat levottomuutta vasteeseen rajataajuutensa ympäristöön.
MLSSA-järjestelmä sisältää luonnollisesti omat virheterminsä, mutta ne ovat luonnoltaan sangen pieniä edellä mainittuihin verrattuna.
Putken käyttökelpoisen taajuusalueen rajallisuutta on käsitelty jo aiemmin tässä selostuksessa. Alle 20 Hz ja yli 690 Hz taajuuksilla menetelmän malli kaatuu, eikä tuota oikeita tuloksia. 1 kHz suodattimen vaikutus ei tullut kovin selkeästi esiin käyrästöistä, koska niiden taajuusalue pyrittiin rajaamaan putken käyttökelpoiselle alueelle.
Kaikki mitattavat kappaleet eivät olleet aivan tarkasti putken päädyn kokoisia ja muotoisia. Tiivis ja aukoton sijoittelu ovat tärkeitä mittauksen luotettavuuden takia, joten jotkut käyrästöjen epämääräisyydet johtuivat varmasti pienistä raoista levyjen kylkien ja putken seinämien välillä.
Mikrofonin tarkan sijainnin mittaaminen osoittautui myös vaikeaksi, ja laskennassa tyydyttiin 0,5 cm mittaustarkkuuteen. Tämä on eräs virhetekijä, jota voisi pienentää rakenteellisilla ominaisuuksilla (valmis mittajana) ja huolellisella mittaustyöllä.
Seuraavat akustiikkalevyt mitattiin sekä uusitulla järjestelmällä että molemmilla referenssijärjestelmillä. Tulokset eivät siis ole vertailukelpoisia alkuperäisen järjestelmän tuottamien tulosten kanssa, mutta kylläkin keskenään.
materiaali
|
paksuus
|
| Villa
EL100/Paroc
|
100
mm
|
| Duratex
|
16
mm
|
| Plain
|
15
mm
|
Näytepalojen kiinnityksen tiiviyteen kiinnitettiin TKK:n putkessa erityistä huomiota tiivistämällä näytepalojen ja putken seinämien saumat silikonimassalla mittausten ajaksi.
TKK:n putkella ja VTT:n B&K 4002 -järjestelmällä (Otaniemi) mittaukset tehtiin sekä sini- että MLS-herätteellä. VTT:n B&K 4206 -järjestelmä (Tampere) käyttää mittauksiin kohinaherätettä.
Mittauksista lasketut absorptiosuhteet ja normalisoidut impedanssit on koottu seuraavassa kohdassa oleviin taulukoihin, joiden avulla TKK:n järjestelmällä saatuja tuloksia voi selkeästi vertailla referenssijärjestelmien tuottamiin tuloksiin.
Selostuksen liitteissä on myös mittaustuloksista lasketut kompleksisen heijastuskertoimen ja impedanssin käyrät sekä SAS-mittausten yksityiskohtaiset tulostaulukot.
Seuraaviin taulukoihin on koottu mittaustuloksista lasketut absorptiosuhteet ja normalisoidut impedanssit. Referenssimittausten tulokset ovat samassa taulukossa vertailun helpottamiseksi. VTT:n 4206-järjestelmän tulokset on arvioitu käyrästöiltä, koska tuloksia ei ollut saatavana numeerisesti käsiteltävässä muodossa. Selostuksen liitteissä on mittaustuloksista lasketut heijastuskertoimen ja impedanssin käyrästöt, joiden avulla tulosten yksityiskohtiin voi tarkemmin tutustua.


Alla olevissa taulukoissa on lueteltu valmistajien akustolevyille ilmoittamia arvoja:
Plain Duratex (ilmarako takana)
f / Hz
Taulukko . Valmistajien ilmoittamia arvoja
abs.kerroin
f
/ Hz
Lmin
/ dB
125
0,33
125
0,16
250
0,31
250
0,11
500
0,53
500
0,10
1000
0,68
1000
0,07
Toisistaan poikkeavat mittausjärjestelmät ja tulosten laskentatavat
vaikeuttavat tulosten vertailua keskenään. Esimerkiksi B&K 4002
-laitteistolla ei yleensä tehdä mittauksia käyttäen
MLS-herätettä, joten putken soveltuvuudesta tarkoitukseen ei
sinänsä voitu antaa pitävää tietoa. B&K 4206
-järjestelmä puolestaan käyttää
kohinaherätettä, ja laskee oheisohjelmistoineen mittaustulokset
valmiiksi käyriksi, joten jo erilainen laskentatapa, ikkunointi,
taajuuspehmennys ym. saa järjestelmän tulokset erottumaan muista.
Lisäksi voidaan olettaa, että kyseisen järjestelmän
tarkkuus ja mittaustulosten paikkansapitävyys ovat näissä
puitteissa omaa luokkaansa.
Yleisesti voidaan todeta, että ainakin TKK:n putkella ja VTT:n 4002-putkella tehtyjen mittausten tuloksissa on selvästi havaittavia virhetermejä, jotka sinällään vaikeuttavat tulosten arviointia. Parhaiten tämä käy ilmi liitteissä olevista heijastuskerrointen kuvaajista, sillä niissä näkyy paljon yli yhden olevia arvoja, jotka eivät sinänsä ole fysikaalisesti perusteltavissa.
Alkuperäisissä mittaustuloksissa esiintyvien suurten virhemarginaalien pääasialliseksi syyksi paljastui putken kaiutinpäädystä heijastuva ääni, joka vääristää mittausten tuottamat impulssivasteet. Alla olevassa kuvassa on esitetty EL100-villan mittauksessa TKK:n ja VTT:n 4002-putkilla MLS-herätteellä saadut aikatason impulssivasteet. Kummastakin käy erittäin selvästi ilmi, ettei materiaalin impulssivaste ole todellinen, vaan sen loppuosa muodostuu putkessa edestakaisin heijastelevasta äänestä. Heijastukset vaimenevat VTT:n putkessa melko nopeasti, mutta TKK:n putkessa vielä kymmeneskin heijastus erottui selvästi käyrästä jopa kuvissa käytetyllä lineaarisella amplitudiasteikolla.
On selvää, ettei tämän kaltaisesta mittausdatasta saada oikeita tai edes järkeviä tuloksia ainakaan erikseen käsittelemättä. Päädyimme ikkunoimaan TKK:n mittausdatat aikatasossa niin, että jäljelle jäi vain suoran äänen ja materiaalista saapuvan heijastuksen pituinen näytejono, joka katkaistiin pehmennetysti ensimmäisen putken kaiutinpäästä saapuvan heijastuksen kohdalta. Todellisten mittausten kannalta tämä lähestymistapa ei tietenkään ole käyttökelpoinen, sillä 310 pisteen mittausdatasta laskettu taajuusresoluutio on erittäin paljon epätarkempi kuin alkuperäisellä 16 000 pisteen näytejonolla.


Ensimmäinen häiriöheijastus alkaa kuvassa näytteestä 305, joten sitä seuraavat näytteet rajattiin ikkunoimalla pois. Ikkunan pituudeksi valittiin 310 pistettä, sillä häiriöheijastuksen alku vaihteli näytteiden 303-309 välillä eri mittauksissa eri paksuisten materiaalien johdosta. Ikkuna muodostettiin suorakaideikkunasta, jonka viimeiset 20 pistettä korvattiin 40 pisteen Hanning-ikkunan laskevalla puoliskolla. Täten impulssivasteen alkupää saatiin säilymään muuttumattomana, mutta samalla torjuttiin ensimmäisten häiriöheijastusten kohdaltaan vaihteleva alkuosa.
Häiriöheijastusten poistaminen on siis ensisijainen tavoite, jotta putkea voitaisiin hyödyntää todellisissa mittauksissa. Putken tämänhetkinen rakenne aikaansaa erittäin voimakkaita pitkittäisiä heijastuksia, ja tyhjänä putken sisus muistuttaakin lähinnä kaiuntahuonetta. Ratkaisuna voidaan luontevasti esittää heijastusten vaimentamista asentamalla putken kaiutinpäähän sopivia absorbentteja. Koska heijastukset vaimenevat hitaasti jopa yllä kuvatun 100 mm:n vuorivillanäytteen kanssa, on kaiutinpäätyyn löydettävä sangen tehokas vaimennusratkaisu.
Eräs vaihtoehto on asentaa kaiutinpäätyyn koko putken poikkipinnan täyttävä villatukos, jonka paksuus ja tarkka sijainti täytyy luultavasti kokeellisesti määrittää. Vaimentimen paksuus vaikuttaa tietenkin myös kaiuttimesta tulevan suoran äänen vaimennukseen, joten liian paksu vaimennin voi heikentää mittauksen signaalikohinasuhdetta tai aiheuttaa signaalin säröytymistä kaiuttimessa liian suurilla voimakkuuksilla. Vaimentimen sijainnin on oltava lähellä kaiutinelementtejä, jotta putken efektiivinen pituus ei lyhene ja sen myötä mittausten alarajataajuus ei nouse.
Seuraavassa on esitetty ikkunoiduista näytejonoista lasketut käyrät:








Plain, normalisoitu impedanssi, TKK:n mittaukset, ikkunoitu 310
näytettä aikatasossa
Eräs järjestelmässä esiintyvä, hankalasti mitattava virhetermi on herätesignaalin kytkeytyminen mikrofoniin putken rakenteita pitkin. Kaiutinelementtien tukeva kiinnitys on samalla myös jäykkä, eikä elementtien kiinnityslevyä ole eristetty putkesta kovinkaan paljon. Signaali voi siis kytkeytyä mekaanisesti putken runkoon ja edetä putken seinämiä pitkin mikrofonikapselin tukirakenteisiin. Myöskään kaiutinelementtien taustaa ei ole eristetty millään tavoin, joten akustinen takaisinkytkeytyminenkin on teoriassa mahdollista.
Eräs mittauksia vaikeuttava virhetermi on signaalin säröytyminen, mikäli kaiutinelementtejä ohjataan liian suurella tasolla. Jouduimme uusimaan yhden mittaussarjan, koska putkeen syötetty signaali oli säröytynyt ja tulokset olivat erittäin hajanaisia.
Mikrofonin ja mikrofonivahvistimen dynamiikka asettavat omat rajansa järjestelmän suorituskyvylle. Liian voimakkaalla signaalilla myös mikrofonivahvistin voidaan saada aiheuttamaan säröä.
Putken vaimennuskertoimen mittaukseen ei ryhdytty, vaikka se kuvattiinkin ISO-standardin vedoksessa. Toisaalta vaimennuksen aiheuttama virhe voidaan ainakin toistaiseksi katsoa melko mitättömäksi järjestelmän muiden virhetermien rinnalla.
Koska MLSSA-ohjelmalla mitatut näytejonot ovat pitkiä, PSD- ja CSD-funktioiden laskeminen kestää nopeallakin koneella melko kauan. Tulosten käyriä voisi mahdollisesti vielä tasoittaa onnistuneemmalla Hanning-ikkunan pituuden valinnalla. VTT:n 4206-järjestelmän tuottamiin tuloksiin verrattuna Matlab-laskennan tuottamat tulokset sisältävät paljon häiriöitä, jotka olisi vielä mahdollista peittää. Jossakin tulee kuitenkin vastaan raja, jossa keskiarvoistus alkaa peittää yksityiskohtia. Toisaalta vertailu referenssimateriaaliin voisi puolustaa voimakkaamman keskiarvoistuksen käyttöä.
Tärkein mittausjärjestelmään tarvittava parannus on putken kaiutinpäädyn vaimentaminen niin tehokkaasti, ettei putkessa synny merkittäviä jälkiheijastuksia. Tämä parantaisi mittaustuloksia todella merkittävästi.
Muita parannusehdotuksia ovat ISO-standardiehdotuksessa mainitut kaiutinelementtien mekaaninen ja akustinen eristäminen, mitattavan kappaleen kiinnityksen järkeistäminen ja mikrofonin etäisyyden mittaustarkkuuden parantaminen. Kukin näistä parannuksista vaatisi kuitenkin huomattavia muutoksia putken nykyiseen rakenteeseen.
Erikoistyön alkuperäiseksi tavoitteeksi annettiin parannetun mittausjärjestelmän suunnittelu ja toteutus. Tämä osoittautui kuitenkin laajuudeltaan niin mittavaksi, ettei merkittävää osaa ISO-standardiehdotuksessa mainituista koemittauksista ja parannuksista voitu toteuttaa. Järjestelmää ei siis saatu kovin tyydyttävään toimintakuntoon erikoistyön puitteissa. Toisaalta työn myötä kertyi runsaasti havaintoja järjestelmän puutteista ja mahdollisista parannuksista. Lisäksi käytännönläheinen tehtävänasettelu toi mukanaan myös käytännön kokemusta mittauksista ja niiden ongelmista. Kaiken kaikkiaan työ oli tekijöilleen sangen mielenkiintoinen, joskin aivan liian laaja kyseisen kurssin suoritukseen nähden.
1. IEC-standardi 268-5. Sound system equipment, Part 5: Loudspeakers, CEI 1989. 69 s.
2. ISO/DIS 10534 -standardiehdotus. Acoustics - Determination of sound absorption coefficient and impedance or admittance by the impedance tube method, ISO 1994. 27 s.
3. ISO/DIS 10534 -standardiehdotuksen laajempi vedos. Acoustics - Determination of sound absorption coefficient and impedance or admittance by the impedance tube method, ISO 1994. 33 s.
4. Lahti T., Akustinen mittaustekniikka. Raportti 38, joulukuu 1995. TKK 1995. 152 s.
5. Lahti T., Analysis methods for acoustical systems based on the FFT and intensity techniques. VTT publications 67, VTT 1990. 288 s.
6. Kivelä S., Matlab-opas. 4. painos, Otatieto 1991. 60 s.
7. Toivanen J., Teknillinen akustiikka. 4. painos, Otatieto 1976. 518 s.
8. Aplac 7.0 Manuals. TKK OFFSET 1996.
1. Mittaustulokset ja niiden tarkastelu
* Alkuperäisellä järjestelmällä saadut tulokset
* Uusitulla järjestelmällä saadut tulokset
* Referenssijärjestelmillä saadut tulokset
2. Tulosten laskentaan käytetyt Matlab-makrot
* ekaputki.m -- ensimmäisten mittaustulosten laskenta
* sasputki.m -- uusitun järjestelmän ja referenssitulosten laskenta
* timread.m -- MLSSA:n .tim-tiedostojen luku Matlabiin
* ikkputki.m -- ikkunoitujen tulosten laskenta
* absimp.m, minmaxka.m -- absorptiosuhteiden ja impedanssien laskenta oktaavikaistoittain
3. MLSSA-mittausjärjestelmän asetukset (.SET-tiedosto)
4. Kaiutinelementtien ripustukseen käytetty tukilevy ja niiden kiinnitys
5. Äänilähteen parannusten määrittäminen
6. Kaiutinelementtien THD-mittaustulokset
7. Aplac-tiedostot
Herätesignaalin suunnitteluohje
Mikrofoniesivahvistimen suunnitteluohje
8. B&K 4002- ja 4206-järjestelmien kuvaukset
9. Mikrofonikapselin valmistajan ilmoittamat tiedot