Changeset 3632 in Sophya for trunk/Cosmo/SimLSS/geneutils.cc

Timestamp:

May 25, 2009, 1:57:27 PM (16 years ago)

Author:

cmv

Message:

add lobe N-S pour dipole avec regroupement analogique, cmv 25/05/2009

File:

• trunk/Cosmo/SimLSS/geneutils.cc (modified) (1 diff)

trunk/Cosmo/SimLSS/geneutils.cc

@@ -361 +361 @@
 
 //-------------------------------------------------------------------
+/*=========================================================
+Antenne de longueur totale d (lue au milieu cad 2 brins de d/2)
+see Jackson p401  (antenne selon Oz, centre z=0 en O)
+k = w/C = 2 Pi Nu / C = 2 Pi / Lambda
+- densite de courant: J = I sin(kd/2 - k|z|)*delta(x)*delta(y)
+                        (selon Oz) et pour |z|<d/2
+                        I current peak value if kd>=Pi
+- current at gap (center): I0 = I sin(kd/2)
+#
+Z=sqrt(Mu/Epsilon) is the impedance of the vacum (376.7 Ohms)
+  cad pour le vecteur de pointing: |S| = |E0|^2/Z
+- Time average radiated power by solid angle unit:
+ t is theta the polar angle, t=0 along antenna, Pi/2 perpendicular
+ dP/dOmega = Z*I^2/(8Pi^2) * (cos(kd/2*cos(t)) - cos(kd/2))^2 / sin(t)^2
+           if we define  L = d/lambda ,  kd/2= Pi * L
+ dP/dOmega = Z*I^2/(8Pi^2) * (cos(Pi*L*cos(t)) - cos(Pi*L))^2 / sin(t)^2
+#
+- pour t -> 0
+          dP/dOmega = Z*I^2/(8Pi^2) * ( 1/4 * (Pi*L*sin(Pi*L)*t)^2 )
+- pour t -> Pi
+          dP/dOmega = Z*I^2/(8Pi^2) * ( 1/4 * (Pi*L*sin(Pi*L)*(t-Pi))^2 )
+#
+- Polarisation in plane containing antenna and the direction of propagation
+#
+- pour kd << 1  (L->0)  -> approximation du dipole:
+          dP/dOmega = Z*I^2/(8Pi^2) * ( (Pi*L)^4/4 * sin(t)^2 )
+- pour kd=Pi  (L=1/2)  -> half-wave length:
+          dP/dOmega = Z*I^2/(8Pi^2) * cos(Pi/2*cos(t))^2/sin(t)^2
+- pour kd=2*Pi  (L=1)  -> full-wave length:
+          dP/dOmega = Z*I^2/(8Pi^2) * 4*cos(Pi/2*cos(t))^4/sin(t)^2
+#
+- Le programme retourne Val tel que: dP/dOmega = Z*I^2/(8Pi^2) * Val
+    L<0 (arg[0]<0) alors Val est normalise tel que le maximum = 1
+#
+La carte antenne fournie par Peterson:
+         d=8cm, distance entre 2 centres consecutifs 10cm
+#
+# ---- AntCentFed et AntDipole
+# Input:
+#   t = angle entre la direction de la radiation et le fil de l'antenne (radian)
+#   L = d/Lambda  (d=longueur totale de l'antenne)
+# Output:
+#   (8Pi^2) dP/dOmega / (Z*I^2)
+========================================================= */
+
+double AntCentFed(double L, double trad)
+{
+ double pil = M_PI*L;
+ double st = sin(trad), ct = cos(trad);
+ // L'antenne standing-wave
+ double v = 0.;
+ if(fabs(st)<1e-3) {
+   if(ct>0.) {  // theta ~= 0
+     v = pil*sin(pil)*trad/2.;  v *= v;
+   } else {  // theta ~= 180
+     v = -pil*sin(pil)*(trad-M_PI)/2.;  v *= v;
+   }
+ } else {
+   v = (cos(pil*ct) - cos(pil))/st;  v *= v;
+ }
+
+ return v;
+}
+
+double AntDipole(double L, double trad)
+{
+ double pil = M_PI*L;
+ double st = sin(trad);
+ // L'approximation du dipole
+ double vd = pil*pil*st/2.;  vd *= vd;
+ return vd;
+}
+
+//-------------------------------------------------------------------
 
 static unsigned short IntegrateFunc_GlOrder = 0;