Skip to content

Commit

Permalink
Browse files Browse the repository at this point in the history
calculs sur le CS
dans ce dossier programmes pour calculer le CPC, et les angles et
paramètres géométriques du miroir de Fresnel
  • Loading branch information
ooriane committed Jun 27, 2015
1 parent 0b45ef4 commit 3937aa4
Show file tree
Hide file tree
Showing 2 changed files with 68 additions and 0 deletions.
35 changes: 35 additions & 0 deletions géométrie CS matlab/CPC.m
@@ -0,0 +1,35 @@
% test equation param CPC ferme these veynandt
thetaC=45*pi/180; % angle max d'acceptance (rad)
r=48/2; % rayon du tube
thetaMin=pi/2+thetaC;
thetaMax=3*pi/2-thetaC;
theta=[0:0.1:thetaMax];
Mtheta=theta;
% Mtheta2=theta;
u=length(theta);
for i=1:u
if theta(i)>thetaMin
Mtheta(i)= (thetaMin+theta(i)-cos(theta(i)-thetaC))/(1+sin(theta(i)-thetaC));
end
end
figure
x=r*(sin(theta)-Mtheta.*cos(theta));
y=r*(-cos(theta)-Mtheta.*sin(theta));
% x2=r*(sin(theta)-Mtheta2.*cos(theta));
% y2=r*(-cos(theta)-Mtheta2.*sin(theta));
plot(x,y)
hold on
plot(-x,y)
theta_tube=[0:0.1:2*pi];
xtube=r*cos(theta_tube);
ytube=r*sin(theta_tube);
plot(xtube,ytube);
% plot(x2,y2)
% for k=1:(u-1)
% dx(k)=x(k+1)-x(k);
% dy(k)=y(k+1)-y(k);
% df(k)=dy(k)./dx(k);
% end
% hold on
% plot(x(1:u-1),df)

33 changes: 33 additions & 0 deletions géométrie CS matlab/paramCS.m
@@ -0,0 +1,33 @@
% parametrage miroir
% angle d'incidence solaire (thtrans et thi)
% orientation N-S
% jour 269 (alternatiba paris 26/09) à midi
S=soleil_ah(269,15,45); % azimut S(1) et hauteur S(2) solaire
thtrans=atan(sin(S(1))/tan(S(2)))
thi=asin(cos(S(1))*cos(S(2)));
% 2 cadre
n=10; % nb facettes par cadre
lm=10; % largeur facette miroir
dx=14.3; % intervalle centres deux facettes consécutives
% centre facette m1 = 14.3
H=150; % hauteur absorbeur (cm)
i=-n:n; % indice facette
xi=dx*i; % abscisse centre facette i
thini=atan(xi/H); % angle entree CPC facette i
alphai=(thini+thtrans)/2; % angle d'inclinaison facette i
figure
hold on
plot(i,thini*180/pi)
plot(i,alphai*180/pi)
% largeur faisceau réflechi intercepté par l'horizontale facette i
li=(cos(alphai)+xi/H.*sin(alphai))*lm;
% en tenant compte en plus du diametre apparant du soleil
gamma=16/60*pi/180; %16 min d'angle
li2=li+2*H*tan(gamma)./(cos(thi).*(cos(thini).^2));
figure
hold on
plot(i,li)
plot(i,li2)
% ajouter les effets d'ombrage!!


0 comments on commit 3937aa4

Please sign in to comment.