source: Sophya/trunk/Poubelle/archediab.old/archediab.sources/c/dilution.c@ 649

Last change on this file since 649 was 649, checked in by ansari, 26 years ago

archediab 28 vol
File size: 15.4 KB
RevLine 
[637]1#include "diabolo.h"
2#include "dilution.h"
3#include "controle.h"
4
5
6//-------------------- fonction exec de la fenetre : dilution -------------------------------
7
8#define bit0 0x01
9#define bit1 0x02
10#define bit2 0x04
11#define bit3 0x08
12#define bit4 0x10
13#define bit5 0x20
14#define bit6 0x40
15#define bit7 0x80
16
17
[643]18#define cnt_temperature(i) (20+(int)val_temperature(blk->ADC_dil[i]))
[637]19
20#define ecrit_switch(val,ii,bit) {if(val) ii=ii | bit ;else ii=ii& (0xffffffff^bit) ;}
21#define calib(j) litD(fenetre_dilution_calibration,j,0L)
[639]22
23
[643]24void ecrit_moteur_dilution(int valeur); // valeur sur 10 bit: poid fort=on/off
[639]25
[643]26void ecrit_moteur_dilution(int valeur)
27{
28int i;
29char mot_tc[10];
30
31 i=litD(fenetre_dilution,dil_switch_1,0);
32i=i<<8; i+=litD(fenetre_dilution,dil_switch_2,0);
33i=i<<8; i+=litD(fenetre_dilution,dil_switch_3,0);
34i=i<<8; // les 3 octets en poid fort, rien dans le dernier octet
35i=i& 0x003fffff;
36i = i| ( (valeur & 0x3ff)<<22);
37ecritD(fenetre_dilution,dil_switch_1,(i>>24)&0xff);
38ecritD(fenetre_dilution,dil_switch_2,(i>>16)&0xff);
39ecritD(fenetre_dilution,dil_switch_3,(i>>8 )&0xff);
40
41mot_tc[0]=tc_switch_dil;
42mot_tc[1]=0;
43mot_tc[2]=0;
44mot_tc[3]=0;
45mot_tc[4]=(i>>8 )&0xff;
46mot_tc[5]=(i>>16)&0xff;
47mot_tc[6]=(i>>24)&0xff;
48mot_tc[7]=0;
49emission_telecommande(tc_dir_transputer,mot_tc);
50}
51
52
[639]53// fonction ecrit un bit des switch dilution
54// c'est une telecommande normale 64 bit a supprimer
55// la telecommande passe a travers le transputer sans etre traitee par lui
56// elle arrive directement sur l'EPLD
[649]57/*
[637]58void ecrit_bit_dilution(int valeur,int bit);
59
60void ecrit_bit_dilution(int valeur,int bit)
61{
62int i;
63char mot_tc[10];
64
65 i=litD(fenetre_dilution,dil_switch_1,0);
66i=i<<8; i+=litD(fenetre_dilution,dil_switch_2,0);
67i=i<<8; i+=litD(fenetre_dilution,dil_switch_3,0);
68i=i<<8; // les 3 octets en poid fort, rien dans le dernier octet
69ecrit_switch(valeur,i,bit);
70ecritD(fenetre_dilution,dil_switch_1,(i>>24)&0xff);
71ecritD(fenetre_dilution,dil_switch_2,(i>>16)&0xff);
72ecritD(fenetre_dilution,dil_switch_3,(i>>8 )&0xff);
73
74mot_tc[0]=tc_switch_dil;
75mot_tc[1]=0;
76mot_tc[2]=0;
77mot_tc[3]=0;
78mot_tc[4]=(i>>8 )&0xff;
79mot_tc[5]=(i>>16)&0xff;
80mot_tc[6]=(i>>24)&0xff;
81mot_tc[7]=0;
82emission_telecommande(tc_dir_transputer,mot_tc);
83}
[649]84*/
[639]85// fonction ecrit un des dac dilution
86// c'est une telecommande normale 64 bit a supprimer
87// la telecommande passe a travers le transputer sans etre traitee par lui
88// elle arrive directement sur l'EPLD
[637]89
[639]90void ecrit_dac_dilution(void);
[637]91
[639]92void ecrit_dac_dilution(void)
93{
94char mot_tc[10];
95int k;
96unsigned long i,j;
97i=litD(fenetre_dilution,dil_dac_1,0); // les 2 bits adresse du dac
98i=i<<14; i+=litD(fenetre_dilution,dil_dac_2,0); // les 12 bits valeur du dac
99// ranger les bits a l'envers;
100j=0; for(k=0;k<16;k++) {j=(j<<1) | (i&1);i=i>>1;}
101mot_tc[0]=tc_dac_dil;
102mot_tc[1]=(j)&0xff;
103mot_tc[2]=(j>>8)&0xff;
104mot_tc[3]=0;
105mot_tc[4]=0;
106mot_tc[5]=0;
107mot_tc[6]=0;
108mot_tc[7]=0;
109emission_telecommande(tc_dir_transputer,mot_tc);
110}
[637]111
[639]112
113
114
[637]115void exec_dilution(int fen,int item,double valeur,...)
116{
117int i;
118
119if(item>1000) item-=1000; // pour appeler le case pour tous les cara d'un edit texte
120
121switch(item)
[649]122 {
123 case dil_switch_helium : emission_tc_reduite(tc2_dilution,(valeur==1),nbit_switch_helium); break;
124 case dil_EVB : emission_tc_reduite(tc2_dilution,(valeur==1),nbit_vanne_EVB); break;
125 case dil_EVO : emission_tc_reduite(tc2_dilution,(valeur==1),nbit_vanne_EVO); break;
126 case dil_EVF : emission_tc_reduite(tc2_dilution,(valeur==1),nbit_vanne_EVF); break;
127 case dil_EVV : emission_tc_reduite(tc2_dilution,(valeur==1),nbit_vanne_EVV); break;
128 case dil_sw_pp5 : emission_tc_reduite(tc2_dilution,(valeur==1),nbit_switch_pile_par_5); break;
129 case dil_sw_pp15 : emission_tc_reduite(tc2_dilution,(valeur==1),nbit_switch_pile_par_15); break;
130
131 /*
132
[637]133 case dil_switch_helium : ecrit_bit_dilution((int)valeur,switch_helium); break;
134 case dil_EVB : ecrit_bit_dilution((int)valeur,vanne_EVB); break;
135 case dil_EVO : ecrit_bit_dilution((int)valeur,vanne_EVO); break;
136 case dil_EVF : ecrit_bit_dilution((int)valeur,vanne_EVF); break;
137 case dil_EVV : ecrit_bit_dilution((int)valeur,vanne_EVV); break;
138 case dil_sw_pp5 : ecrit_bit_dilution((int)valeur,switch_pile_par_5); break;
139 case dil_sw_pp15 : ecrit_bit_dilution((int)valeur,switch_pile_par_15); break;
[649]140*/
[637]141// ---- commande directe d'un mot avec les 3 case de codage hexadecimales
142
143 // case dil_switch_1 :
144 // case dil_switch_2 :
145 // case dil_switch_3 : ecrit_bit_dilution(0,0); break;
146
147// ---- telecommande du raz periodique des fets
148 case dil_raz_modul :
[649]149 case dil_fet_raz :
150 i =0;
151 if(litD(fen,dil_fet_raz,0)==2) i+=1;
152 if(litD(fen,dil_fet_raz,0)==3) i+=7;
153 emission_tc_reduite(tc2_raz_fet,0,i);
154/*
[637]155 {
156 char mot_tc[10];
157 mot_tc[0]=7; // commande la premiere bebo directe -> il faut ecrire tc reduite
158 i=litD(fen,dil_raz_modul,0);
159 i=i<<3;
160 if(litD(fen,dil_fet_raz,0)==2) i+=1;
161 if(litD(fen,dil_fet_raz,0)==3) i+=7;
162 i=i<<4;
163 mot_tc[1]=i;
164 for(i=2;i<8;i++) mot_tc[i]=0;
165 emission_telecommande(tc_dir_transputer,mot_tc);
166 }
[649]167*/
[637]168 break;
169
[649]170 case dil_ch1 : emission_tc_reduite(tc2_dilution,(valeur==1),nbit_chauffage1); break;
171 case dil_ch2 : emission_tc_reduite(tc2_dilution,(valeur==1),nbit_chauffage2); break;
172 case dil_ch3 : emission_tc_reduite(tc2_dilution,(valeur==1),nbit_chauffage3); break;
173 case dil_ch4 : emission_tc_reduite(tc2_dilution,(valeur==1),nbit_chauffage4); break;
174 case dil_ch5 : emission_tc_reduite(tc2_dilution,(valeur==1),nbit_chauffage5); break;
175 case dil_ch6 : emission_tc_reduite(tc2_dilution,(valeur==1),nbit_chauffage6); break;
176 case dil_ch7 : emission_tc_reduite(tc2_dilution,(valeur==1),nbit_chauffage7); break;
[639]177
178// case dil_dac_1 :
179// case dil_dac_2 :
[643]180// case dil_envoi_dac : ecrit_dac_dilution(); break;
181 case dil_envoi_dac : ecrit_moteur_dilution((int)litD(fenetre_dilution,dil_dac_2,0)); // les 12 bits
182 break;
[637]183
184 default : break;
185 }
186}
187
188
189
190
191void traite_block_dilution(block_type_dilution* blk)
192{
193int i;
194
195// ouvre la fenetre de calibration
196if(!fenetre(fenetre_dilution_calibration))
197 nouveauD(fenetre_dilution_calibration,dilution_calibrations_id,"dilution_calibration",exec_cache);
198
199if( (blk->switch_dil&0x3f ) != 6 ) {
200// printf("**************** affiche erreur \n");
201 changecontrole(fenetre_dilution,dil_erreur,idem,idem,idem,idem,idem,"erreur");
202 }
203else {
204 static int q;
205 char ss[10]="____";
206 q=(q+1)&3;
207 ss[q]='/';
208 changecontrole(fenetre_dilution,dil_erreur,idem,idem,idem,idem,idem,ss);
209// printf("**************** pas d' erreur \n");
210 }
211
212
213// lit les switchs
214// -- les voyants de fin de course vanne principale (logique negative)
215ecritD(fenetre_dilution,dil_voyant_EVO,((blk->switch_dil&switch_EVO)?0:1));
216ecritD(fenetre_dilution,dil_voyant_EVF,((blk->switch_dil&switch_EVF)?0:1));
217
218// les retour de commande des electrovannes
219ecritD(fenetre_dilution,dil_EVB_retour,((blk->switch_dil&vanne_EVB)?1:0));
220ecritD(fenetre_dilution,dil_EVO_retour,((blk->switch_dil&vanne_EVO)?1:0));
221ecritD(fenetre_dilution,dil_EVF_retour,((blk->switch_dil&vanne_EVF)?1:0));
222ecritD(fenetre_dilution,dil_EVV_retour,((blk->switch_dil&vanne_EVV)?1:0));
223
[639]224
[645]225#define c(j,i) (1e-4*(double)parametr.nom_coef[parametr.bolo[j].numero_nom_coef].coef[i])
[639]226
227// les temperatures sur les cartes modifiées
228{
[645]229double I,V,R[8],llR;
[639]230int j,k;
231def_gains;
232k=0;
[645]233
[639]234for(j=0;(j<nb_max_bolo) && (k<4);j++)
235 {
236 if(parametr.bolo[j].bolo_code_util==bolo_thermo_simplifie)
237 {
238 I = 1e-3 * (double)dac_V(gg->reglage.bolo[j]) * 2441. / parametr.bolo[j].bolo_capa; // I en µA
239 V=0.001*bol_micro_volt(blk->temperature[k],(double)parametr.bolo[j].bolo_gain*gain_ampli(gg->reglage.bolo[j]));
240 if(I>0.0000001) R[k]=V/I; else R[k]=0; // R en ‡
[645]241 //------ calcul de l'etalonnage en temperature des cartes temperature simplifiées
242 // j=numero de bebo T[k] = R[k+4] en Kelvin
243 if ((R[k]-c(j,6))>1.)
[649]244 {if((log(R[k]-c(j,6))-c(j,0))>0.001) llR= log(log(R[k]-c(j,6))-c(j,0)) ; else llR=0;
245// printf("\nk=%d j=%d R=%g c2=%g c3=%g llR=%g",k,j,R[k],c(j,2),c(j,3),llR);
246 }
[645]247 else llR=0;
248 R[k+4] = exp( c(j,1) + c(j,2)* llR + c(j,3)* llR* llR + c(j,4)* llR* llR* llR + c(j,5)* llR* llR* llR* llR) ;
249 if(R[k+4]>9999) R[k+4]=9999;
250 k++;
[639]251 }
252 }
253
[649]254gg->temperature_cryo=R[4]; // Charge la temperature_cryo (0.1K) pour les calculs futurs avec les bolos
255/* La germanium doit toujours etre en premier par rapport aux autres thermos absolus */
[645]256
257if(R[0]<90.) ecritD(fenetre_dilution,dil_temp
258 , "ge=%6.2f‡ 4K=%6.2f‡ 1.6K=%6.2f‡ 10K=%6.2f‡\rge= %6.3fK 4K=%6.2fK 1.6K=%6.2fK 10K=%6.2fK "
259 ,R[0],R[1],R[2],R[3],R[4],R[5],R[6],R[7]);
260else ecritD(fenetre_dilution,dil_temp
261 ,"ge=%6.3fk‡ 4K=%6.2f‡ 1.6K=%6.2f‡ 10K=%6.2f‡\rge= %6.3fmK 4K=%6.2fK 1.6K=%6.2fK 10K=%6.2fK "
262 ,R[0]*0.001,R[1],R[2],R[3],R[4]*1000.,R[5],R[6],R[7]);
263
264
[639]265if(litD(fenetre_dilution,dil_trace_tempe,0))// trace les sondes de niveau helium
266 {
267 int temps_cntl=numero_block(blk)*nb_per_block*2;
268 double secondes,minutes;
269 secondes=gg->periode_echantillonage*(double)temps_cntl;
270 minutes=secondes/60.;
271 if(!fenetre(fenetre_temperature_dil))
272 {
273 nouveauD(fenetre_temperature_dil,0,"temperature_dil",0);
274 selectgra(fenetre_temperature_dil);
275 strcpy(graph->xtitre,"minutes");
276 }
[645]277 tracen(fenetre_temperature_dil,8,minutes,R);
[639]278 }
279}
280
[637]281// les pressions et debits metres des injections de la dilution
282ecritD(fenetre_dilution,dil_p_d_3He,"3He: %4.1fb -> %5.2fµm/s -> %4.1fb"
[643]283 ,40. * val_multiplex(blk->ADC_dil[ p_R3]) -1.6 // 200 bars pour 5V et 1.6 bar d'offset
284 ,2. * val_multiplex(blk->ADC_dil[ d_3He]) // 10 MICRO MOLES pour 5V
285 ,20. * val_multiplex(blk->ADC_dil[ p_C3]) // 100 bars pour 5V
[637]286 );
287ecritD(fenetre_dilution,dil_p_d_4He,"4He: %4.1fb -> %5.2fµm/s -> %4.1fb"
[643]288 ,40. * val_multiplex(blk->ADC_dil[ p_R4]) // 200 bars pour 5V
289 ,8. * val_multiplex(blk->ADC_dil[ d_4He]) // 40 MICRO MOLES pour 5V
290// ,4. * val_multiplex(blk->ADC_dil[ d_4He]) // 20 MICRO MOLES pour 5V
291 ,20. * val_multiplex(blk->ADC_dil[ p_C4]) // 100 bars pour 5V
[637]292 );
293
294ecritD(fenetre_dilution,dil_p_air,"Van=%4.1fb charb=%4.1fb mmb=%4.3fb haut=%4.3fb"
[643]295 ,20.*val_multiplex(blk->ADC_dil[ p_air]),20.*val_multiplex(blk->ADC_dil[ p_charb])
296 ,0.2*val_multiplex(blk->ADC_dil[ p_memb]),0.2*val_multiplex(blk->ADC_dil[ p_haut]) );
[637]297
298ecritD(fenetre_dilution,dil_piles," Trp=%4.1fV Dil=%4.1fV / %4.1fV \rBebo=%4.1fV / %4.1fV / %4.1fV Ch=%4.1fV"
[643]299 ,2.03*val_multiplex(blk->ADC_dil[ p_10T]),3.90*val_multiplex(blk->ADC_dil[ p_p18D]),3.90*val_multiplex(blk->ADC_dil[ p_m18D])
300 ,2.03*val_multiplex(blk->ADC_dil[ p_10B]),3.90*val_multiplex(blk->ADC_dil[ p_p18B]),3.90*val_multiplex(blk->ADC_dil[ p_m18B])
301 ,3.8*val_multiplex(blk->ADC_dil[ p_Ch]));
[637]302
303/* valeur seuil au mini (5.3 ou 15.3V) , *17 = +5.9 V au maxi */
304
305
[643]306ecritD(fenetre_dilution,dil_p_10T,(int)(17*(2.03*val_multiplex(blk->ADC_dil[ p_10T])-5.3)));
307ecritD(fenetre_dilution,dil_p_p18D,(int)(17*(3.90*val_multiplex(blk->ADC_dil[ p_p18D])-15.3)));
308ecritD(fenetre_dilution,dil_p_m18D,(int)(17*(-3.90*val_multiplex(blk->ADC_dil[ p_m18D])-15.3)));
309ecritD(fenetre_dilution,dil_p_10B,(int)(17*(2.03*val_multiplex(blk->ADC_dil[ p_10B])-5.3)));
310ecritD(fenetre_dilution,dil_p_p18B,(int)(17*(3.93*val_multiplex(blk->ADC_dil[ p_p18B])-15.3)));
311ecritD(fenetre_dilution,dil_p_m18B,(int)(17*(-3.93*val_multiplex(blk->ADC_dil[ p_m18B])-15.3)));
312ecritD(fenetre_dilution,dil_p_Ch,(int)(10*(3.8*val_multiplex(blk->ADC_dil[ p_Ch])-25))); /* 25 a 35 V */
[637]313
314// -- les voyants de switch de piles
315ecritD(fenetre_dilution,dil_voyant_pp5,((blk->switch_dil&switch_pile_par_5)?1:0));
316ecritD(fenetre_dilution,dil_voyant_pp15,((blk->switch_dil&switch_pile_par_15)?1:0));
317
318
[649]319ecritD(fenetre_dilution,dil_temperatures,"temp-Haut = %4.1f° %4.1f° Base = %4.1f° %4.1f° \r tube He: %4.1f° piles %4.1f° driver %4.1f°"
320 ,val_temperature(blk->ADC_dil[ t_h2]),val_temperature(blk->ADC_dil[ t_h4])
321 ,val_temperature(blk->ADC_dil[ t_b1]),val_temperature(blk->ADC_dil[ t_b2]),val_temperature(blk->ADC_dil[ t_b3])
322 ,val_temperature(blk->ADC_dil[ t_pile]),val_temperature(blk->ADC_dil[ t_a1])
[637]323 );
324
[649]325ecritD(fenetre_dilution,dil_t_h1,cnt_temperature(t_h2));
326ecritD(fenetre_dilution,dil_t_h2,cnt_temperature(t_h4));
327ecritD(fenetre_dilution,dil_t_h3,cnt_temperature(t_b1));
328ecritD(fenetre_dilution,dil_t_h4,cnt_temperature(t_b2));
329ecritD(fenetre_dilution,dil_t_a1,cnt_temperature(t_b3));
330ecritD(fenetre_dilution,dil_t_a2,cnt_temperature(t_pile));
331ecritD(fenetre_dilution,dil_t_a3,cnt_temperature(t_a1));
[637]332
333
334// les afficheurs analogiques de pression helium et pirani
335
[643]336ecritD(fenetre_dilution,dil_RP_He,(int)(2500*val_multiplex(blk->ADC_dil[ RP_He])-7500));
[637]337// je veux 10000 pour 0.8 bars soit alors que l'on a 5V pour 1 bar ou 1V = 0.2 bar
338// soit multiplier par 2500
339// puis j'enleve 0.6 bars soit 7500
[643]340ecritD(fenetre_dilution,dil_pirani,(int)(1000*val_multiplex(blk->ADC_dil[ pirani])));
[637]341
342
343
344
345
346
347// ecriture des valeurs brutes de la dilution
348if(litD(fenetre_dilution,dil_sortie_brute,0))
349 {
350 if(!fenetre(fenetre_dilution_lecture_brute)) nouveauT(fenetre_dilution_lecture_brute,0,"lecture dilution");
351
352 ecritT(fenetre_dilution_lecture_brute,fin_f,"dil_switch= %x \n",blk->switch_dil);
353
[643]354 for(i=0;i<16;i++) ecritT(fenetre_dilution_lecture_brute,fin_f,"%2d=%6.3f ",i<8?i+11:i+13,val_multiplex(blk->ADC_dil[ i]));
[637]355 ecritT(fenetre_dilution_lecture_brute,fin_f,"\n");
[643]356 for(i=16;i<32;i++) ecritT(fenetre_dilution_lecture_brute,fin_f,"%2d=%6.3f ",i<24?i+15:i+17,val_multiplex(blk->ADC_dil[ i]));
[637]357 ecritT(fenetre_dilution_lecture_brute,fin_f,"\n");
[643]358 for(i=32;i<48;i++) ecritT(fenetre_dilution_lecture_brute,fin_f,"%2d=%6.3f ",i<40?i+19:i+21,val_multiplex(blk->ADC_dil[ i]));
[637]359 ecritT(fenetre_dilution_lecture_brute,fin_f,"\n");
360 }
361// lecture des sondes de niveau helium
362
[639]363
[637]364 {
365 double y[8];
[639]366 static int niveau_pre;
[637]367 int niveau,fin;
[639]368// printf("dilution : numero de block = %d \n",numero_block(blk));
[637]369 changecontrole(fenetre_dilution,dil_niveau_helium,idem,idem,idem,idem,calrefcon(0,0,rouge,blanc,blanc),"");
[643]370 y[0]=val_multiplex(blk->ADC_dil[ j_he1])-calib(dil_j_he1);
371 y[1]=val_multiplex(blk->ADC_dil[ j_he2])-calib(dil_j_he2);
372 y[2]=val_multiplex(blk->ADC_dil[ j_he3])-calib(dil_j_he3);
373 y[3]=val_multiplex(blk->ADC_dil[ j_he4])-calib(dil_j_he4);
374 y[4]=val_multiplex(blk->ADC_dil[ j_he5])-calib(dil_j_he5);
375 y[5]=val_multiplex(blk->ADC_dil[ j_he6])-calib(dil_j_he6);
376 y[6]=val_multiplex(blk->ADC_dil[ j_he7])-calib(dil_j_he7);
377 y[7]=val_multiplex(blk->ADC_dil[ j_he8])-calib(dil_j_he8);
[639]378 if(blk->switch_dil & switch_helium)
[637]379 {
[639]380 changecontrole(fenetre_dilution,dil_niveau_helium,idem,idem,idem,idem,calrefcon(0,0,rouge,blanc,blanc),"");
381 niveau=0;fin=0;
382 for(i=0;i<8;i++)
383 {
384 if(y[i]>0.1) {
385 niveau++;
386 if(fin) niveau=9;
387 }
388 else fin=1;
389 }
390// if(niveau==9) changecontrole(fenetre_dilution,dil_niveau_helium,idem,idem,idem,idem
391// ,calrefcon(0,0,jaune,blanc,blanc),"");
392 if(niveau<niveau_pre) niveau_pre=0;
393 else {niveau_pre=niveau; ecritD(fenetre_dilution,dil_niveau_helium,niveau);}
394// printf("niveau=%d \n",niveau);
[637]395 }
[639]396 else changecontrole(fenetre_dilution,dil_niveau_helium,idem,idem,idem,idem,calrefcon(0,0,vert,blanc,blanc),"");
[637]397
398 if(litD(fenetre_dilution,dil_trace_helium,0))// trace les sondes de niveau helium
399 {
400 int temps_cntl=numero_block(blk)*nb_per_block*2;
401 double secondes,minutes;
402 secondes=gg->periode_echantillonage*(double)temps_cntl;
403 minutes=secondes/60.;
404 if(!fenetre(fenetre_niveau_helium))
405 {
406 nouveauD(fenetre_niveau_helium,0,"niveau_helium",0);
407 selectgra(fenetre_niveau_helium);
408 strcpy(graph->xtitre,"minutes");
409 }
410 tracen(fenetre_niveau_helium,8,minutes,y);
411 }
412 }
413
414}
415
416
Note: See TracBrowser for help on using the repository browser.