17.04.2013 Views

In cazul rezolvǎrii prin modelare numericǎ apar câteva deosebiri ...

In cazul rezolvǎrii prin modelare numericǎ apar câteva deosebiri ...

In cazul rezolvǎrii prin modelare numericǎ apar câteva deosebiri ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

.................................. LUCRAREA Nr.3 Data ................<br />

Anul II ETc seria A G…. g… N…<br />

MODELAREA UNUI CÂMP ELECTROSTATIC PRODUS DE SARCINI DISTRIBUITE.<br />

REZOLVARE ANALITICĂ ŞI NUMERICĂ (QFIELD)<br />

1.INTRODUCERE<br />

Orice problemă reală necesită pentru rezolvare o <strong>modelare</strong>, <strong>prin</strong> care se reţin aspectele esenţiale<br />

ale problemei , echivalându-se unele aspecte complicate <strong>prin</strong> altele mai simple şi renuntându-se la aspecte<br />

cu influenţă minoră în stabilirea soluţiei. Cu cât un model este mai simplu, cu atât şi rezolvarea problemei<br />

este mai simpla şi mai rapidă (chiar <strong>prin</strong> prisma timpului de calcul pe ordinator). Pentru a avea o situaţie<br />

mai apropiată de realitate modelul nu poate fi excesiv de simplu. Evident că şi în acest caz, ca şi în<br />

majoritatea problemelor de tip ingineresc, trebuie găsit un echilibru între cele două tendinţe.<br />

<strong>In</strong> <strong>cazul</strong> <strong>rezolvǎrii</strong> <strong>prin</strong> <strong>modelare</strong> <strong>numericǎ</strong> <strong>apar</strong> <strong>câteva</strong> <strong>deosebiri</strong> esenţiale faţǎ de situaţia realǎ:<br />

- materia şi câmpul sunt caracterizate din punct de vedere macroscopic de noţiunea de continuu, iar<br />

modelele numerice sunt caracterizate de noţiunea de discret<br />

- asupra realitǎţii se poate face un calcul integro-diferenţial, pe când asupra modelelor numerice se<br />

acţioneazǎ cu metode numerice, în <strong>prin</strong>cipal iterative, de rezolvare.<br />

<strong>In</strong> cadrul laboratorului vom utiliza şi modele fizice, care la o anumitǎ scarǎ , cu o anumitǎ precizie<br />

şi cu anumite simplificǎri reproduc realitatea pe care vrem să o studiem. Asupra acestor modele se pot<br />

face masurǎri directe şi <strong>prin</strong> relaţii de analogie şi corespondenţǎ, din aceste masurǎri se pot obţine<br />

marimile ce ne intereseazǎ.<br />

2. PROBLEMA PROPUSA PENTRU REZOLVARE PRIN MODELARE<br />

Pentru explicare vom analiza o problemǎ, simplă, pe care o vom<br />

rezolva utilizând mai multe modele simplificate şi mai multe metode de<br />

rezolvare.<br />

Problema se enunţǎ dupa cum urmeazǎ:<br />

Sǎ se calculeze expresia intensitǎţii câmpului electrostatic,<br />

E = f(x), câmp produs de un fir conductor de lungime l = 3 m întrun<br />

punct oarecare, situat la distanţa x de fir, în planul mediator al<br />

firului. Firul este încǎcat cu o sarcinǎ electricǎ distribuitǎ uniform<br />

în lungul firului, cu densitatea lineară ρl, exprimată în nC/m.<br />

Sarcina totalǎ cu care este încǎrcat firul este Q = 6·G·N nC = 6·G·N<br />

·10 -9 C, în care N este numǎrul de ordine în semigrupǎ, iar G este<br />

Figura 1.<br />

numǎrul grupei. Geometria problemei este prezentată în figura 1.<br />

<strong>In</strong> cele ce urmeazǎ ne vom ocupa de obţinerea unor soluţii<br />

aproximative, pe baza unor modele, cu diferite grade de aproximare. Prin aceste modele cu diferite grade<br />

de precizie în redarea situaţiei reale dorim să arătăm felul în care <strong>prin</strong> modele şi metode aproximative<br />

putem să ne apropiem de soluţia exactă a unei probleme.<br />

2.1. Sarcina concentrată într-un punct.<br />

Un model foarte simplificat (Fig.2) este cel în care considerǎm toatǎ<br />

sarcina (de fapt distribuită uniform pe lungimea firului), ca fiind<br />

concentratǎ la mijlocul lungimii firului. <strong>In</strong> acest caz soluţia<br />

aproximativǎ va fi:<br />

E<br />

1P<br />

1<br />

x ⋅<br />

Q<br />

=<br />

2<br />

4π ⋅ε<br />

⋅<br />

Semnificaţia notaţiilor este următoarea:<br />

0<br />

E1 P este intensitatea câmpului electric produsă de sarcina electrică<br />

concentrată Q, în punctul P situat la distanţa x pe mediatoarea<br />

segmentului încărcat, x u Figura 2.<br />

este versorul pe direcţia distanţei x , orientat de la fir către punct<br />

2.2. Sarcina concentrată în 3 puncte.<br />

Un model ceva mai precis se obţine dacǎ se consideră sarcina concentrată în 3 pǎrţi egale, simetric plasate<br />

în punctele A, B şi C. Aceasta înseamnǎ cǎ am presupus firul iniţial impǎrţit în trei segmente egale, de<br />

ux<br />

(1)


lungime l/3, şi sarcina electricǎ Q/3 uniform distribuitǎ pe fiecare din acestea. Fiecare sarcină Q/3 o<br />

consideram concentratǎ în centrul fiecǎruia din segmentele respective (în punctele A,B şi C).<br />

Pentru acest caz se vor însuma vectorial cele trei<br />

componente, fiecare componentǎ fiind datǎ de câte o sacinǎ<br />

Q/3. Rezultatul compunerii vectorilor este cu<strong>prin</strong>s în planul<br />

ortogonal pe mijlocul firului (plan mediator), datoritǎ<br />

condiţiilor de simetrie a configuraţiei. Cele două componente<br />

determinate de sarcinile plasate în A şi B au aceeaşi valoare<br />

Q / 3<br />

E A = = E şi formează acelaşi unghi cu direcţia x,<br />

2 C<br />

4π<br />

ε 0 d<br />

astfel că suma lor vectorială va fi orientata doar după direcţia x.<br />

Figura 3.<br />

Distanţele AP şi CP sunt egale: d = [x 2 +(l / 3) 2 ] 1/2<br />

Q / 3<br />

EB = este orientat de asemenea după direcţia lui x, astfel că şi suma (rezultanta) vectorială<br />

2<br />

4π<br />

ε x<br />

0<br />

E = E + E + E<br />

3 P A B C este şi ea orientată după direcţia lui x.<br />

⎛ Q / 3 x Q / 3 ⎞<br />

⎛ 2x<br />

1 ⎞<br />

E ⋅u<br />

= Q ⋅ ⋅ + ⋅<br />

⎝ d d x ⎠<br />

⎝ d x ⎠<br />

9<br />

9<br />

9<br />

3P<br />

= ⎜2<br />

⋅9<br />

⋅10<br />

⋅ + 9 ⋅10<br />

⎟ x 3 10 ⎜ ⎟ u<br />

2<br />

2<br />

3 2 x (2)<br />

Se poate uşor înţelege cǎ <strong>prin</strong> continuarea operaţiunii de împǎrţire în segmente şi mai mici (spre<br />

E 9 ca sumă a acţiunii a nouă sarcini<br />

exemplu ficare treime se imparte din nou în trei, obţinându-se P<br />

punctiforme Q/9) se va obţine un rezultat ceva mai precis, mai apropiat de o distribuţie liniară uniformă.<br />

Se poate apoi continua tot asa şi se va ajunge la un model care se apropie tot mai mult de situaţia<br />

realǎ (distribuţie continuǎ a sarcinii), <strong>prin</strong> adunarea a 27 de contribuţii (mai precis), apoi a 81de mici<br />

segmente (şi mai precis), apoi 243 (încă şi mai precis) ş.a.m.d.<br />

Cu cât numǎrul subdiviziunilor este mai mare, cu atât se complicǎ şi calculul. Utilizând metode<br />

numerice de calcul se pot folosi şi modele mai complicate, cu un numǎr foarte mare de subdomenii.<br />

Diferenţa între modelul discontinuu (discret) şi situaţia realǎ (continuu) poate fi sugeratǎ şi <strong>prin</strong><br />

compararea unei imagini reale cu un mozaic care <strong>prin</strong> foarte multe bucǎţele mici (discrete) încearcǎ sǎ<br />

redea o imagine cât mai apropiată de cea realǎ (continuǎ).<br />

2.3. Sarcina distribuită. Calcul analitic, conductor filiform de lungime l .<br />

Se poate face un calcul analitic, folosind metode integro<br />

diferenţiale de calcul şi, în acest caz, se poate lua în calcul faptul<br />

că sarcina este distribuită uniform, cu densitatea liniară ρl = Q/l<br />

,pe un fir conductor ideal extrem de subţire, filiform, adică fără<br />

dimensiune radială, dar de lungime finită l. Calculăm câmpul în<br />

planul mediator (perpendicular pe mijlocul segmentului)<br />

Se poate observa, în figura 4, că pentru orice element<br />

infinitezimal dy plasat în partea superioară există un element<br />

simetric plasat în partea inferioară. Fiecare din aceste elemente îşi<br />

aduce contribuţia (dE1 şi, respectiv, dE2) la câmpul total. Cele<br />

două componente infinitezimale însumate vor da o rezultantă<br />

(dEiP), în planul ortogonal, orientată perpendicular pe fir.<br />

Q<br />

⋅ dy<br />

l ⋅sin<br />

dEiP = 2<br />

l<br />

ρ β<br />

⋅ cosα<br />

⋅u<br />

x = ⋅dy⋅<br />

u<br />

2<br />

2 x (3)<br />

4π<br />

⋅ε<br />

0 ⋅ r<br />

2π<br />

⋅ε<br />

0 ⋅ r<br />

Figura 4.<br />

Prin integrarea acestei rezultante, dE, pentru a obţine<br />

contribuţia întregului fir unghiul β baleiază valorile de la β = 90° la β = θ. θ este unghiul sub care se vede<br />

punctul P de la unul din capetele firului, iar r este distanţa de la elementul infinitezimal dy la punctul P.<br />

2


θ<br />

2<br />

x<br />

x<br />

ρl<br />

sin β ⎛ x ⎞<br />

r = ; y=<br />

x ⋅ctgβ<br />

; dy = − dβ<br />

; EiP<br />

= ∫ sin β ⋅ ⋅⎜<br />

− dβ<br />

⎟<br />

2<br />

2<br />

2<br />

sin β<br />

sin β<br />

2πε<br />

x ⎝ β<br />

π / 2 0<br />

sin ⎠<br />

Q ⋅ cosθ<br />

ρl<br />

⋅ cosθ<br />

EiP =<br />

⋅u<br />

x = ⋅u<br />

x (4)<br />

2π<br />

⋅l<br />

⋅ε<br />

0 ⋅ x 2π<br />

⋅ε<br />

0 ⋅ x<br />

Dacă se poate considera că firul este foarte lung (teoretic infinit), mult mai lung decât distanţa la<br />

care dorim calcularea câmpului (l>>x), putem face un calcul mai simplu, utilizând teorema lui Gauss.<br />

Particularizarea rezultatului precedent obţinut pentru conductor de lungime finită se poate particulariza<br />

pentru <strong>cazul</strong> unui conductor foarte lung (teoretic infinit) punând θ = 0.<br />

Q<br />

ρl<br />

E∞P =<br />

⋅u<br />

x = ⋅u<br />

x (5)<br />

2π ⋅l<br />

⋅ε<br />

0 ⋅ x 2π<br />

⋅ε<br />

0 ⋅ x<br />

De asemenea, se poate face un calcul analitic al componentelor şi rezultantei câmpului, produs de<br />

un segment încărcat cu sarcina electrică uniform distribuită, într-un punct oarecare plasat coplanar cu<br />

segmentul de fir încărcat.<br />

2.4. Sarcina uniform distribuită pe un conductor filiform, presupus infinit lung; calcul utilizând<br />

legea fluxului electric (teorema lui Gauss).<br />

Legea fluxului electric, aplicată pe suprafaţa închisă<br />

cilindrică Σ, alcătuită din cele două baze Sb1 şi Sb2, precum şi din<br />

suprafaţa laterală Slat se scrie :<br />

D ⋅ dS<br />

= Q<br />

(6)<br />

3<br />

∫ Σ<br />

∫ D ⋅ dS<br />

= ∫D⋅dS+<br />

Σ<br />

= 0 + 0 +<br />

Σ<br />

D ⋅ dS<br />

+<br />

Sb 2<br />

Slat<br />

∫ D ⋅ dS ⋅ 0 ∫<br />

Slat<br />

Sb1 ∫<br />

= ε E(<br />

r)<br />

dS = ε E(<br />

r)<br />

⋅ 2π<br />

⋅ r ⋅h<br />

Slat<br />

∫<br />

1 0<br />

D ⋅dS<br />

=<br />

În integrala pe suprafaţa laterală, Slat , E(r) s-a putut scoate<br />

în afara semnului de integrală pentru că, din condiţii de simetrie,<br />

Elat (ca modul) are o acceaşi valoare, E(r), în orice punct al<br />

suprafeţei laterale cilindrice ( r fiind acelaşi).<br />

Atât pe baza superioară, b1, cât şi pe cea inferioară, b2, D = ε 0 E şi<br />

dS sunt vectori perpendiculari, deci integralele respective sunt<br />

nule (cos90 ° = 0). Pe suprafaţa laterală a lui Σ, D = ε 0 E şi dS sunt vectori omoparaleli (cos 0 ° Figura 5.<br />

=1), iar pe<br />

Slat modulul lui D (deci şi al lui E) este constant, fiecare punct de pe suprafaţa laterală fiind în condiţii<br />

identice faţa de orice alt punct situat pe aceeaşi suprafaţă, la aceeaşi distanţă r faţă de axa de simetrie a<br />

figurii (axă care coincide de fapt cu conductorul filiform).<br />

QΣ este sarcina repartizată pe fir şi este acea parte a<br />

sarcinii care se găseşte în interiorul suprafeţei închise Σ. Este<br />

vorba de o porţiune de lungime h, încărcata cu o sarcină<br />

Figura 6.<br />

uniform distribuită ρl . Deci:<br />

Q l h ⋅ = Σ ρ şi egalând cu (3) ε 0 E( r)<br />

⋅2π ⋅r<br />

⋅h<br />

= l h ⋅ ρ ceea ce<br />

ρ l<br />

conduce la E ( r)<br />

= = E (8)<br />

G<br />

2π<br />

ε 0 r<br />

2.5. Sarcina uniform distribuită pe suprafaţa unui<br />

conductor cilindric foarte lung de rază „a”<br />

Rezolvarea seamănă cu cele prezentate la punctul 2.4<br />

(figura 5). Metoda aplicată şi aici este tot utilizarea legii<br />

fluxului electric (sau teorema lui Gauss pentru <strong>cazul</strong> vidului).<br />

Pentru a putea să facem comparaţie între rezultatele obţinute<br />

în cele două situaţii vom considera o echivalenţă privind<br />

încărcarea cu sarcină electrică . O aceeaşi lungime, h, din<br />

cele două conductoare trebuie să conţină o aceeaşi sarcină<br />

electrică, fie că ea este uniform distribuită pe fir cu densitatea<br />

(7)


lineară ρl sau este distribuită uniform, cu densitatea ρS , pe suprafaţa laterală a unui cilindru conductor de<br />

rază „a”, foarte lung (presupus infinit lung).<br />

Se poate stabili astfel o corespondenţă între ρl şi ρS.<br />

ρl<br />

ρl<br />

⋅ h = ρ S ⋅2π<br />

a ⋅ h ; ρ S =<br />

(9)<br />

2π<br />

a<br />

Aplicând legea fluxului electric: D⋅ dS<br />

= Q ∫Σ<br />

Σ<br />

r r<br />

r<br />

D<br />

r<br />

dS<br />

=<br />

r r r r r r<br />

D ⋅ dS<br />

+ D ⋅ dS<br />

+ D ⋅ dS<br />

= 0 + 0 + D ⋅ dS ⋅1<br />

= D dS = D ⋅ 2π<br />

⋅ r ⋅ h<br />

∫ ⋅ ∫ ∫ ∫ ∫ ∫<br />

Σ<br />

S1<br />

S 2<br />

Slat<br />

Slat Slat<br />

ρ S ⋅ 2πah<br />

a<br />

iar QΣ = 2 π ⋅ a ⋅ h ⋅ ρ S<br />

D( r)<br />

= = ρ S ⋅<br />

(9)<br />

2π<br />

r h r<br />

Dacă în această expresie se înlocuieşte densitatea de sarcină distribuită pe suprafaţă cu valoarea<br />

a ρ l a ρ l<br />

echivalentă (9) se obţine: E( r)<br />

= ρ S ⋅ = ⋅ = , expresie identică cu (8).<br />

ε 0r<br />

2π<br />

a ε 0r<br />

2π<br />

ε 0 r<br />

2.6. (SUPLIMENTAR) Sarcina electrică distribuită pe un conductor (metalic) cilindric de lungime<br />

finită “l” şi rază “a” Această punere a problemei nu poate fi rezolvată “exact” cu metode integrodiferenţiale,<br />

de calcul analitic.Rezolvarea se poate face utilizând MEF, programul Qfield. Pentru cei<br />

interesaţi un model al acestei probleme se găseşte în fişierele L306.<br />

3. MODELAREA UNOR CONDENSATOARE<br />

3.1.Considerente teoretice<br />

Condensatorul electric se defineşte ca fiind un ansamblu de două corpuri conductoare (armături metalice)<br />

încărcate cu sarcini electrice egale în modul, dar de semn opus. Între armături se stabileşte un câmp<br />

electrostatic complet (nu există decât linii de câmp ce pornesc de pe armătura pozitivă şi care, toate, se<br />

termină pe armătura negativă, adică nu există linii de câmp care să unească una din armături cu un alt<br />

corp sau cu „infinitul”).<br />

Un caz particular este cel în care<br />

condensatorul este alcătuit din armături<br />

plane şi paralele. În <strong>cazul</strong> când armăturile<br />

plane au forma identică (indiferent de<br />

formǎ) şi sunt plasate “faţă în faţă”,<br />

condensatorul se numeste condensator plan.<br />

Pentru acest condensator, cu anumite<br />

aproximări care vor fi puse în evidenţă în<br />

cadrul lucrării, este valabilă formula<br />

capacitaţii condensatorului plan (cunoscută<br />

Figura 7.<br />

încă din liceu);<br />

ε 0 ⋅ ε r ⋅ S<br />

C =<br />

(10)<br />

d<br />

formulă în care ε0 este permitivitatea absolută a vidului, εr este permitivitatea relativă a materialului<br />

dielectric care umple complet spaţiul dintre armături, S este suprafaţa (identică) a armăturilor plasate „faţă<br />

în faţă”, iar d este distanţa între cele două armături paralele. Condensatorul cilindric ideal este alcătuit din<br />

doi cilindri concentrici (figura 8), între care se stabileşte o diferenţă de potenţial V1 – V2 . Pentru ca între<br />

armătura centrală şi cea exterioară să se stabilească un câmp electric complet trebuie ca armătura<br />

exterioară să fie la potenţial 0. În caz contrar vor exista şi linii de câmp care vor pleca către ∞ (dacă V2 ><br />

0) sau se vor închide la armǎtura exterioară venind de la ∞ (dacă V2 < 0) şi, în ambele, cazuri acestea vor<br />

<strong>apar</strong>ţine mediului exterior condensatorului, existând câmp şi deci şi energie în afara spaţiului din<br />

interiorul condensatorului.<br />

În <strong>cazul</strong> general, indiferent de geometria condensatorului (dimensiuni diferite şi formă diferită<br />

pentru cele douǎ armǎturi, poziţia lor relativǎ oarecare) şi de proprietǎţile mediului dintre armǎturi<br />

capacitatea electică a sa se calculează cu formula (11) :<br />

4


C<br />

Q<br />

+ = cu V1 > V2 (11)<br />

V1<br />

−V2<br />

Q+ este sarcina electrică a armăturii pozitive, adică aici Q1 , unde V1 > V2<br />

V1 şi V2 sunt potenţialele armăturii pozitive şi,<br />

respectiv, negative a condensatorului. Reamintim că<br />

V1 nu trebuie să fie neapărat pozitivă şi nici V2 nu<br />

trebuie să fie, neapărat, negativă, ci doar V1 > V2 ,<br />

doar că, aşa cum am arătat anterior, între o armătură<br />

cu potenţial nenul si un alt corp metalic (sau ∞) se<br />

stabileşte un câmp exterior. Acestui câmp exterior îi<br />

corespunde o capacitate parazită (în general de zeci<br />

sau de sute de ori mai mică decât capacitatea<br />

<strong>prin</strong>cipală) , care poate fi neglijată în multe cazuri,<br />

dar nu intotdeauna.<br />

Formula generală de calcul a capacităţii (11)<br />

este evident valabilă şi în <strong>cazul</strong> particular al unui<br />

condensator cilindric. Pentru <strong>cazul</strong> particular, ideal,<br />

prezentat în figura 8, se poate determina o formulă de<br />

Figura 8<br />

Figura 9<br />

calcul a capacităţii ce ţine seama de geometria<br />

condensatorului şi de proprietăţile de material,<br />

presupunând cǎ întregul câmp este concentrat între armǎturi (se neglijeazǎ efectul de capǎt ).<br />

Pentru efectuarea calculelor ne referim la figura 8. Condensatorul cilindric, presupus ideal, are razele<br />

armăturilor R1 şi R2 , înăţimea H şi mediul dielectric dintre armături este aerul (ε = ε0 sau εr = 1). Ţinând<br />

seama de simetria condensatorului legea fluxului electric se aplică pe o suprafaţă cilindrică, Σ ,de o rază<br />

curentă, oarecare, dar R1 < r < R2 şi înălţime h.<br />

Aplicând legea fluxului electric : ∫ D ⋅ dS<br />

= QΣ<br />

Σ<br />

r r<br />

r r r r r r r r<br />

D ⋅ dS<br />

= D ⋅ dS<br />

+ D ⋅ dS<br />

+ D ⋅ dS<br />

=<br />

∫ ∫<br />

Σ<br />

S<br />

1<br />

∫<br />

Slat Slat<br />

S<br />

Slat<br />

= 0 + 0 + ∫ D ⋅ dS ⋅1<br />

= D ∫dS<br />

= D ⋅ 2π<br />

⋅ r ⋅ h<br />

2<br />

∫<br />

Q1<br />

h<br />

iar QΣ = ⋅ 2π<br />

⋅ R1<br />

⋅ h = Q+<br />

⋅<br />

(14)<br />

2πR1<br />

H<br />

H<br />

În relaţia 12 integrala pe suprafaţa închisă Σ a fost descompusă în trei integrale. <strong>In</strong>tegralele pe cele două<br />

baze ale suprafeţei cilindrice Σ sunt nule, pentru că D şi dS sunt ortogonali şi deci produsul lor scalar este<br />

0 (cos π/2 = 0).Pe suprafaţa laterală a cilindrului vectorii D şi dS sunt omoparaleli (aceeaşi direcţie şi<br />

acelaşi sens) şi, din motive de simetrie, în orice punct <strong>apar</strong>ţinând suprafeţei laterale a cilindrului modulul<br />

vectorului D este acelaşi, motiv pentru care a putut fi scos în afara semnului de integrare. În relaţia (14)<br />

sarcina din interiorul suprafeţei Σ se găseşte doar pe porţiunea din suprafaţa laterală a armăturii interioare<br />

“interceptată” de suprafaţa Σ. Din egalarea relaţiilor (13) şi (14) rezultă:<br />

Q+<br />

1<br />

Q+<br />

1<br />

D = ⋅ şi pentru că D = ε 0 ⋅ E , rezultă E = ⋅<br />

(15)<br />

2π<br />

⋅ H r<br />

π ⋅ε<br />

⋅ H r<br />

2 0<br />

Se poate acum calcula tensiunea electrică între cele două ărmături, în funcţie de sarcină:<br />

U<br />

1,<br />

2<br />

= V<br />

1<br />

−V<br />

2<br />

1<br />

1<br />

1<br />

0<br />

(12)<br />

(13)<br />

R2<br />

R2<br />

R<br />

r<br />

2<br />

r<br />

Q+<br />

1 Q+<br />

R<br />

= ∫E⋅dr= ∫E⋅dr = ∫ ⋅ dr = ⋅ ln<br />

2π<br />

⋅ε<br />

⋅ H r 2π<br />

⋅ε<br />

⋅ H R<br />

R<br />

R<br />

R<br />

şi pentru capacitatea electrică a condensatorului cu armăturile de formă cilindrică, coaxiale, rezultă<br />

2π<br />

⋅ε<br />

0 ⋅ H<br />

C =<br />

R<br />

(17)<br />

2<br />

ln<br />

R<br />

1<br />

5<br />

0<br />

2<br />

1<br />

(16)


4.PARTEA PRACTICA.<br />

4.1. Calculul unor valori ale câmpului cu ajutorul unor expresii determinate analitic<br />

x = 1 m x = 3 m Se completează datele din tabel utilizând formulele adecvate.<br />

E1P V/m<br />

Semnificaţiile şi expresiile intensităţii câmpului în punctul P au fost<br />

E3P<br />

EiP<br />

E∞P<br />

EGP<br />

V/m<br />

V/m<br />

V/m<br />

V/m<br />

precizate la paragrafele 2.1. ... 2.4.<br />

Aceste calcule se fac în cadrul laboratorului, măcar parţial.<br />

Fiecare student trebuie să aibă calculate cel puţin 4 valori, pe care le<br />

trece în tabel. Restul calculelor pot fi efectuate acasă, verificarea lor<br />

având loc odată cu lucrarea L7 (de încheiere a activităţii de laborator).<br />

4.2. Modelarea unui condensator plan ideal cu armături dreptunghiulare.<br />

Problema modelată este L301.<br />

a) Se fixează valoarea tensiunii între armături la valoarea de 10V (se intervine în Data). Valoarea<br />

potenţialului armăturii pozitive se ia V1 = +G volţi, G fiind numarul grupei din care face parte studentul.<br />

EA V/m VA V<br />

EB V/m VB V<br />

b) Se notează datele geometrice ale condensatorului<br />

modelat; S – aria armăturilor şi d – distanţa între<br />

armături<br />

EC V/m VC V<br />

c) Se preiau din programul Qfield valorile câmpului<br />

electric şi ale potenţialului electric în puncte situate pe linia mijlocie ce uneşte cele două armături. Punctul<br />

A este situat la jumătatea distanţei între armături, punctul B este situat mai aproape de armătura pozitivă,<br />

la ¼ din distanţa dintre armături şi apoi punctul C este situat la ¾ din distanţa dintre armături, distanţă<br />

măsurată faţă de armătura pozitivă. Se va prezenta poziţia punctelor pe un desen calitativ. Se scriu<br />

valorile EA, EB şi EC , respectiv, VA, VB B şi VC<br />

în referat.<br />

d) Se calculeaza teoretic valoarea condensatorului modelat<br />

SUPLIMENTAR<br />

e) Se modifică distanţa dintre armături, apropiindu-le cu câte 50mm faţă de centru fiecare (pentru a păstra<br />

simetria) şi se constată creşterea valorii câmpului electric.<br />

SE LUCREAZA IN L305. NU MODIFICATI L301!<br />

f) Se modifică condiţiile de frontieră, lasând câmpul electric să iasă din spaţiul dintre armături (L307). Se<br />

poate observa ca în interiorul armăturilor câmpul nu mai este riguros uniform şi că în exterior campul este<br />

cu atât mai important cu cât armăturile sunt mai îndepărtate<br />

4.3. Modelarea unui codensator coaxial cilindric cu armăturile de raze R1 şi R2.<br />

Modelarea acestei probleme se poate face , în 2D, în două moduri.<br />

Se poate face o secţiune în lungul<br />

axei de simetrie a problemei (axa<br />

cilindrilor coaxiali). Aşa este<br />

prezentat în figura 10. Se remarcă<br />

direcţia radială a liniilor de câmp ,<br />

precum şi faptul că modulul<br />

vectorului intensitate a câmpului<br />

electric este mai mare în apropierea<br />

armăturii interioare şi este mai mic<br />

pe măsură cene îndepărtăm de<br />

aceasta.<br />

De asemenea este foarte important<br />

de remarcat că intensitatea câmpului<br />

electrostatic este zero în interiorul<br />

corpurilor metalice (armăturile<br />

condensatorului). Din acest motiv grosimea armăturii exterioare nu contează, ea neinfluenţând asupra<br />

felului cum este distribuit câmpul electric, localizat între armături.<br />

a). Să se reprezinte, calitativ, variaţia câmpului electric E, în lungul razei, <strong>prin</strong> centrul<br />

condensatorului. Linia de reprezentare se alege între punctele de coordonate : z = 500 mm, r = 0 şi z =<br />

500 mm şi r = 500mm. Se notează valorile R1 şi R2 .<br />

6


Se reprezintă intensitatea câmpului şi în interiorul armăturilor metalice (E = 0).<br />

b) Se reprezintă, calitativ, variaţia potenţialului, în lungul aceleiaşi linii.<br />

c) Să se noteze valoarea intensităţii în punctul situat la mijlocul distanţei dinte armături şi valoarea<br />

potenţialului în acelaşi punct.<br />

d) Se calculează teoretic valoarea capacotăţii condensatorului modelat<br />

Modelarea se poate face şi <strong>prin</strong>tr-o secţiune perpendiculară pe axa de simetrie a condensatorului, secţiune<br />

efectuată în centrul (zona centrală) condensatorului<br />

SUPLIMENTAR.<br />

e) Se vizualizează elementele corepunzătoare problemei L303 (<strong>modelare</strong> a secţiunii perpendiculare pe axa<br />

de simetrie a condensatorului). Se pot urmari şi aici cele de la 4.3.a,b şi c.<br />

f) Se poate urmări în problema L304 felul cum a fost modelat un condensator sferic. Se poate înţelege că<br />

din considerente de simetrie a putut fi făcută reprezentarea doar a unui sfert din secţiunea <strong>prin</strong> sferă. Se<br />

pot urmări vectorii de cămp şi reprezentarea (Plot) în funcţie de rază a diferitelor mărimi.<br />

<strong>In</strong> toate cazurile modelate se poate calcula capacitatea electrică a condensatorului modelat.<br />

TEME OBLIGATORII<br />

6. Condensatorul plan reprezentat în problema L301 (punctul 4.2) are dielectric arerul şi este alimentat cu<br />

tensiunea de 25V. Distanţa dintre armături este 40 mm, iar armăturile au dimensiunile 0,4·π m x 72N cm<br />

6.1. Să se calculeze capacitatea condensatorului C0 exprimată în pF<br />

6.2. intensitatea câmpului electric E0 în punctul central al condensatorului<br />

6.3. densitatea de sarcină ρ0 de pe armătura pozitivă<br />

Dacă se introduce între armături un bloc dielectric papalelipipedic, cu permitivitatea relativă 4, care ocupă<br />

complet spaţiul dintre armături, păstrând aceeaşi diferenţă de potenţial între armături, să se calculeze:<br />

6.4. capacitatea condensatorului C4,<br />

6.5. intensitatea câmpului E4 în punctul central al condensatorului<br />

6.6. densitatea de sarcină ρ4 de pe armătura pozitivă<br />

Se deconectează de la sursă condensatorul cu dielectricul între armaturi şi, după deconectare, se scoate<br />

parţial dielectricul dintre armaturi. Considerând că acum dielectricul ocupă jumătate din volum, iar<br />

cealaltă jumătate este ocupată de aer, precum şi considerând că întregul câmp electrostatic este concentrat<br />

în volumul dintre armături, să se deseneze schematic condensatorul neomogen şi să se calculeze:<br />

6.7. capacitatea condensatorului C1/2<br />

6.8. tensiunea la bornele condensatorului U1/2<br />

6.9. densitatea de energie in volumul ocupat de aer<br />

Răspunsurile la întrebările 6.1. .....6.9. se predau cu ocazia efectuării următoarei lucrări<br />

PROBLEME SUPLIMENTARE<br />

7. Să se calculeze analitic capacitatea teoretică a condensatoarelor cilindric şi sferic, modelate, şi să se<br />

compare cu rezultate obţinute folosind Qfield.<br />

8. Să se calculeze densitatea de sarcină pe armăturile condensatorului cilindric ideal, modelat<br />

9. Să se calculeze energia înmagazinată în fiecare din condensatoarele modelate, daca tensiunea aplicată<br />

este de 10V.<br />

7

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!