Visualizza wikitesto di Calibrazione dei trasduttori subacquei con il metodo di reciprocità

{{Risorsa
| tipo = lezione
| materia1 = Studio delle basi idrofoniche
| avanzamento = 100%
}}

La calibrazione dei trasduttori con il metodo della reciprocità consente di determinare con buona precisione la [[Rumori del mare e dei preamplificatori idrofonici del sonar|sensibilità di un idrofono]] sconosciuto in base a misure elettroacustiche opportunamente eseguite con l'aiuto di due trasduttori ausiliari non necessariamente calibrati.

Condizioni basilari per l'esecuzione della misura sono:

* Il trasduttore da calibrare deve essere '''reversibile''', deve cioè poter operare sia in ricezione, sia in trasmissione.
* Le onde acustiche che si muovono tra i trasduttori devono essere onde dirette (in assenza di riflessioni) <ref>Se l'ambiente non è privo di riflessioni le misure devono essere eseguite con l'emissione/ricezione di segnali impulsivi da discriminare dalle riflessioni tramite adatti sistemi a soglia temporale misurando le ampiezze in Vpp. con oscilloscopio ( i livelli in Vpp. devono poi essere trasformati in Veff. )</ref>.
* Il suono deve propagarsi per divergenza sferica.
* La distanza <math>(d)</math> stabilita per il campo di misura deve essere precisa e costante.
* Le misure devono essere condotte in acqua con strumentazione elettronica adatta.

In questa lezione svilupperemo un esercizio relativo alla calibrazione di un trasduttore simulando numericamente tutte le variabili che concorrono al risultato finale con riferimenti continui al campo di misura.

Per tale compito ci serviremo delle caratteristiche elettriche di tre trasduttori, uno dei
quali simulerà l'elemento da calibrare.

Mostreremo la formula per il calcolo della sensibilità, la tipologia dei trasduttori, la procedura di misura sul campo e, per meglio chiarire la metodologia, una simulazione numerica delle operazioni fisiche.

==La formula di calcolo per la calibrazione per reciprocità ==

La calibrazìone di un trasduttore consiste nella determinazione della sua sensibilità  in ricezione espressa  come: 

<math> S = dB / V / \mu Pa</math>.

La formula in per il calcolo di <math>S</math> è :

:<math>\text{Sens. in dB/Volt/microPa}=\text{20 Log}\left[1.14\cdot10^{-9}\sqrt{\frac{e_3}{it}\frac{e_2}{e_1}d\frac{1530}{F} }\right]</math>

dove le variabili sono:
*<math>e1 \  ; \ e2 \ ; \ e3</math>,  tensioni in Veff. (ricavate dalle misure in mare secondo le procedure di seguito indicate in base alle figure 1 e 2)
*<math>it</math>,   corrente in  Ampere (ricavata da una misura in mare secondo procedura)
*<math>F</math>, frequenza in Hz. (posta in fase di emissione)
*<math>d</math>,  distanza del campo di misura in <math>m</math>.  
 
L'applicazione della formula, disponendo a tavolino delle sei variabili menzionate, è semplice; resta invece la perplessità che una procedura di rilevamento dati come quella indicata per le figure 1 e 2 ,  possa non essere compresa a fondo.

Per la ragione sopra esposta assumeremo trasduttori con caratteristiche note
onde simulare i calcoli ed i rilievi relativi alle  misure in acqua, tale procedura è illustrata di seguito.

==Caratteristiche dei trasduttori impiegati nel sistema di misura==

Per l'impostazione del metodo di reciprocità per la calibratura di un trasduttore reversibile, al fine di simulare per via numerica le misure in acqua, è utile operare con dati relativi a trasduttori reali con i quali, applicando le formule opportune, si possa avere un riscontro immediato sulla validità della procedura adottata. 

Per far ciò ci serviamo delle caratteristiche di due trasduttori sferici della ITC. <ref>Vedi bibliografia</ref> 
(International Transducer Corporation) i cui modelli sono sotto indicati:

Come trasduttore che simula quello da calibrare si assume il modello ITC 1001

Per la coppia dei trasduttori ausiliari si assume il modello ITC 1042 ( nominati nel nostro contesto: 1042/d  e 1042/s ) 

I trasduttori ausiliari non necessariamente devono essere dello stesso tipo.

Le caratteristiche del modello 1001, secondo dati ITC, <ref>Si tratta di un trasduttore progettato per l'emissione di potenza.</ref> alla frequenza di <math>10000 \ Hz</math> sono:
*Sensibilità in ricezione : <math> -192 \ dB / V / \mu Pa </math> '''Questa è la sensibilità che dovrà risultare dai calcoli relativi alla simulazione delle misura in mare'''
*Risposta in trasmissione : <math> 135 \ dB / \mu Pa / V / 1 \ m</math>
*Modulo dell'impedenza : <math> | Z | = 207 \ Ohm </math>


Le caratteristiche del modello 1042, secondo dati ITC, <ref>Si tratta di un trasduttore progettato per l'emissione / ricezione a bassi livelli.</ref> alla frequenza di <math> 10000 \ Hz</math> sono:
*Sensibilità in ricezione : <math> -198 \ dB / V / \mu Pa </math>
*Risposta in trasmissione : <math> 110 \ dB / \mu Pa / V / 1 \ m</math>

==Disposizione dei trasduttori sul campo di misura==

La disposizione dei trasduttori in acqua prevede due configurazioni diverse (figure 1 e  2) per  due diverse misure da eseguire in successione:

*la misura <ref>Se possibile la misura del livello di tensione deve essere fatta con un voltmetro selettivo. </ref> di cui alla figura 1 mostra la prima fase del rilievo dati in cui il trasduttore 1042/s eccitato dal generatore gen. emette onde acustiche che, dopo il percorso <math>(d)</math>, colpiscono la coppia 1042/d e 1001 ai capi dei quali sono monitorate e misurate le rispettive tensioni in ca. 

[[File:542dtc.jpg|thumb|600px|left|figura 1]] 
{{clear}}

*la misura di cui alla figura 2 mostra la seconda fase del rilievo dati in cui il trasduttore 1001, eccitato dal generatore gen. emette onde acustiche che, dopo il percorso <math>(d)</math>, colpiscono il  1042/d ai capi del quale è  monitorata e misurata la tensione in ca.

In questa fase è inoltre misurata, in ca, la corrente assorbita dal trasduttore 1001.

[[File:543dtc.jpg|thumb|600px|left|figura 2]] 
{{clear}}

==Come eseguire le misure==
Le variabili rilevate sul campo nelle due fasi di misura sono ora nominate ed elencate con riferimento al numero di figura: 

*figura 1. (il 1042/s trasmette alla coppia 1001 e 1042/d)
Si regola la frequenza del generatore gen. (su 1042/s) per <math> F \ = \ 10000 \ Hz </math> <ref> Il valore della frequenza è quello per il quale sono date le caratteristiche del trasduttore</ref>, si aggiusta il livello di gen. affinché i trasduttori, disposti in ricezione alla distanza <math>(d)</math>, possano evidenziare segnali ben misurabili che saranno rilevati rispettivamente con le sigle:

<math>e1</math> = tensione efficace ai capi del 1042d

<math>e2</math> = tensione efficace ai capi del 1001 

*figura 2. (il 1001 prende il posto del 1042/s e trasmette al 1042/d) 
Si regola la frequenza del generatore gen. (su 1001) per <math> 10000 Hz </math> <ref> Il valore della frequenza è quello per il quale sono date le caratteristiche del trasduttore</ref>, si aggiusta il livello di gen. affinché il trasduttore, disposto in ricezione alla distanza <math>(d)</math>, possa evidenziare un segnale ben misurabile che sarà rilevato rispettivamente con la sigla:

<math>e3 </math>= tensione efficace ai capi del 1042/d

Contemporaneamente sarà misurata la corrente in ca assorbita dal 1001 indicandola
con la sigla:

<math>it</math> = corrente alternata circolante in 1001 durante la sua fase di trasmissione.

==Simulazione della misura relativa alla sensibilità del trasduttore 1001==
 
L'obiettivo della simulazione,sviluppata in sei passi, consiste nel determinare la sensibilità del trasduttore 1001 preso come elemento da calibrare; ovviamente di tale trasduttore si conoscono a priori le caratteristiche elettriche compresa la sensibilità in ricezione che, con la simulazione, si dovrà trovare prossima al dato calcolato.

Iniziamo l'esame della procedura secondo la figura 1; in essa compaiono: a destra 1042/d e 1001, a sinistra 1042/s.

Il trasduttore 1001 simula l'elemento in fase di calibrazione, i due 1042 sono accessori alla misura in sei passi che vediamo nell'ordine seguente.

== Fase preparatoria sul campo==
*I^) definizione della  distanza <math>(d)</math> del campo, ad esempio <math>  d \ = \ 10 \ m </math>, collocandovi i trasduttori secondo figura 1. 
*II^) accensione del  generatore gen.  impostando il valore <math>F</math> della frequenza di lavoro <math>F = 10000 \ Hz</math>, regolandolo per avere <math>25 \ Veff.</math> ai capi di 1042/s, quest'ampiezza può garantire, come vedremo, tensioni facilmente rilevabili all'uscita di 1042/d e 1001.

== Azioni sul campo simulate numericamente - prima fase -==

===Passo a ===
Misurando il valore di <math>e1 </math> ai capi di 1042d, essendo nota la risposta in trasmissione del 1042s ( <math> 110 \ dB / \mu Pa / V / \ 1 \ m</math> ) alla frequenza di <math>10000 \ Hz </math>) con  <math> 25 \ V </math> applicati da gen. pari a  <math> +28\  dB / V  </math>  si avrebbe, ad un metro di distanza, una pressione acustica di:

<math> +110 \ dB / \mu Pa / V / 1 \ m +  28 \ dB \ / V = + 138 \ dB / \mu Pa / 1\ m  </math>.


Detta pressione raggiungerebbe i trasduttori 1001 e 1042/d  dopo aver percorso <math>10 \ m</math>, subendo un'attenuazione pari a:

<math> Att = 20 \cdot log_{10} {10} = 20 \ dB  </math>

la pressione incidente sarebbe pertanto:

<math>Pi = + \ 138 \ dB / \mu Pa / 1 \ m - 20 \ dB </math> = <math> + \ 118 \ dB / \mu Pa </math> .

Dato che la sensibilità in ricezione a  <math> 10000 \ Hz </math> di 1042/d  è: 

<math> -198 \ dB / V / \mu Pa </math>  il valore della  tensione  e1 in uscita  sarebbe:

<math> e1 = -198 \ dB / V / \mu Pa </math>  <math> + \ 118 \ dB / \mu Pa </math>    = <math> -80 \ dB / V </math> pari a  <math> 100 \ \mu V </math>.


===Passo b ===
 
Misurando il valore di <math>e2 </math> ai capi di 1001, in modo analogo al passo a), si misurerebbe il livello all'uscita di 1001 che, in base alla sua sensibilità in ricezione di <math> - 192 \  dB / V / \mu Pa </math>, genererebbe una tensione <math> e2 </math> così computata:

<math> e2 = - 192 \ dB / V / \mu Pa \ + \ 118 \ dB / \mu Pa </math>  =  <math> - 74 \ dB / V  </math> pari a <math> 200  \ \mu V </math>. 

La prima fase di simulazione ha condotto alla determinazione delle variabili <math>e1 \ ed \ e2</math> così come la misura fisica  sul campo avrebbe potuto produrre.

== Azioni sul campo simulate numericamente - seconda fase -==

===Passo c ===

La simulazione procede secondo la disposizione di  figura 2; in essa compaiono: a destra il 1042/d e a sinistra il 1001.

Il trasduttore 1001, che simula l'elemento in fase di calibrazione, nella  geometria di figura 2  emetterebbe il segnale acustico verso 1042/d secondo i valori di seguito computati.

Per la misura del valore di e3 ai capi di 1042d,  confermato il valore <math>F </math> della frequenza del generatore gen. <math> F = 10000 \ Hz </math>; si regolerebbe gen. per avere <math>10 \ Veff.</math> ai capi del 1001, quest'ampiezza garantirebbe  tensioni facilmente rilevabili all'uscita di 1042/d . 

Essendo note le caratteristiche elettriche di trasmissione del 1001  in: <math> + \ 134 \ dB / \mu Pa / V / 1 \ m </math> alla frequenza di <math>10000 \ Hz</math> con <math>10 \ V </math> applicati da gen. pari a <math> + \ 20 \ dB / V </math> si avrebbe una pressione acustica di: 

<math> + \ 134 \ dB / \mu Pa / V / 1 \ m +  20 \ dB / V = +  \ 154 \ dB / \mu Pa / 1 \ m. </math> 

Detta pressione raggiungerebbe 1042/d  dopo aver percorso <math> 10 \ m</math>, subendo un'attenuazione pari a:

<math> Att = 20 \cdot log_{10} {10}   </math> = <math> 20 \ dB</math> , la pressione incidente sarebbe pertanto:

<math>Pi = + \ 154 \ dB / \mu Pa /  1 \ m - 20 \ dB = + 134 \  dB / \mu Pa </math>

Dato che la sensibilità in ricezione a  <math>10000 \ Hz </math>  di 1042/d  è: <math>- \ 198 \ dB / V / \mu Pa </math> il valore della tensione  e3 in uscita  sarebbe: 

<math> e3 = \ -198 \ dB \ / V / \mu Pa +  \ 134 \  dB / \mu Pa = - 64 \ dB / V</math>  pari a <math>631 \  \mu V</math>

===Passo d ===

Si deve ora simulare la misura della corrente it che scorre nel 1001, corrente determinata dal rapporto tra la tensione applicata, <math> 10 \ V </math>, e il modulo dell'impedenza già indicato nelle sue le caratteristiche:  <math> | Z | = 207 \ Ohm </math>; quindi:

<math> it =  10 \ V / 207 \ Ohm </math>  =  <math> 0.048 \ A </math>

Con quest'ultimo calcolo abbiamo completato l'acquisizione simulata delle variabili necessarie al
computo della sensibilità del 1001. 

Richiamando la formula:
:<math>\text{Sens. in dB/Volt/microPa}=\text{20 Log}\left[1.14\cdot10^{-9}\sqrt{\frac{e_3}{it}\frac{e_2}{e_1}d\frac{1530}{F} }\right]</math>

 con:

 e1 = 100 microVolt = 100/1000000 V    

 e2 = 200 microVolt = 200/1000000 V  

 e3 = 631 microVolt = 631/1000000 V   

 it = 0.048 A   

 F = 100000 Hz

 d = 10 m

 si ottiene:

Sens. = <math> - 193 \ dB / V / \mu Pa </math> contro il valore di targa di <math> - \ 194 \ dB / V  / \mu Pa.</math>

==  Conclusioni ==

Il contenuto di questa  lezione ha consentito di apprendere le modalità per la calibrazione di un trasduttore reversibile; il risultato del calcolo, in base ai valori forniti dal costruttore, porta ad un valore della sensibilità dedotta pari a <math> - 193 \ dB / V / \mu Pa </math> contro quella indicata da ITC di <math> - \ 194 \ dB / V / \mu Pa </math> .

Quanto sopra evidenzia la correttezza della simulazione che ha permesso di giungere al risultato con un errore del 10%.

Se i rilevi di: <math>  e1 \ ; \ e2 \ ; \ e3 \ ; \ it, </math> fossero stati eseguiti con cura sul campo  forse avremmo potuto contare su di un errore di quell'ordine purché la distanza <math>(d)</math> fosse rimasta costante e non fossero intervenute perturbazioni nel mezzo; non si deve dimenticare, infatti, che misure in acqua con errori del <math> \ +/- 10\% \ (+/- 1 \ dB) </math> sono da ritenersi ottime. 

Chiudiamo la lezione facendo notare che la formula per il calcolo della sensibilità del trasduttore discende da una teoria molto complicata per la quale si rimanda alla letteratura indicata in bibliografia.

==note==
<references/>

==Bibliografia==
G. Pazienza, Fondamenti della localizzazione marina, La Spezia, Studio grafico Restani, 1970.
	
International Transducer Corporation, ''Catalog of underwater sound transducers'', Califonia, 1970

*{{Cita libro|autore= Del Turco|titolo=''Sonar- Principi - Tecnologie – Applicazioni'' |editore=Tip. Moderna La Spezia|anno=1992|cid=Del Turco}}