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

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

archediab 27
File size: 14.5 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
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}
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)
122 {
123 case dil_switch_helium : ecrit_bit_dilution((int)valeur,switch_helium); break;
124 case dil_EVB : ecrit_bit_dilution((int)valeur,vanne_EVB); break;
125 case dil_EVO : ecrit_bit_dilution((int)valeur,vanne_EVO); break;
126 case dil_EVF : ecrit_bit_dilution((int)valeur,vanne_EVF); break;
127 case dil_EVV : ecrit_bit_dilution((int)valeur,vanne_EVV); break;
128 case dil_sw_pp5 : ecrit_bit_dilution((int)valeur,switch_pile_par_5); break;
129 case dil_sw_pp15 : ecrit_bit_dilution((int)valeur,switch_pile_par_15); break;
130
131// ---- commande directe d'un mot avec les 3 case de codage hexadecimales
132
133 // case dil_switch_1 :
134 // case dil_switch_2 :
135 // case dil_switch_3 : ecrit_bit_dilution(0,0); break;
136
137// ---- telecommande du raz periodique des fets
138 case dil_raz_modul :
139 case dil_fet_raz :
140 {
141 char mot_tc[10];
142 mot_tc[0]=7; // commande la premiere bebo directe -> il faut ecrire tc reduite
143 i=litD(fen,dil_raz_modul,0);
144 i=i<<3;
145 if(litD(fen,dil_fet_raz,0)==2) i+=1;
146 if(litD(fen,dil_fet_raz,0)==3) i+=7;
147 i=i<<4;
148 mot_tc[1]=i;
149 for(i=2;i<8;i++) mot_tc[i]=0;
150 emission_telecommande(tc_dir_transputer,mot_tc);
151 }
152
153 break;
154
155 case dil_ch1 : ecrit_bit_dilution((int)valeur,chauffage1); break;
156 case dil_ch2 : ecrit_bit_dilution((int)valeur,chauffage2); break;
157 case dil_ch3 : ecrit_bit_dilution((int)valeur,chauffage3); break;
158 case dil_ch4 : ecrit_bit_dilution((int)valeur,chauffage4); break;
159 case dil_ch5 : ecrit_bit_dilution((int)valeur,chauffage5); break;
160 case dil_ch6 : ecrit_bit_dilution((int)valeur,chauffage6); break;
161 case dil_ch7 : ecrit_bit_dilution((int)valeur,chauffage7); break;
[639]162
163// case dil_dac_1 :
164// case dil_dac_2 :
[643]165// case dil_envoi_dac : ecrit_dac_dilution(); break;
166 case dil_envoi_dac : ecrit_moteur_dilution((int)litD(fenetre_dilution,dil_dac_2,0)); // les 12 bits
167 break;
[637]168
169 default : break;
170 }
171}
172
173
174
175
176void traite_block_dilution(block_type_dilution* blk)
177{
178int i;
179
180// ouvre la fenetre de calibration
181if(!fenetre(fenetre_dilution_calibration))
182 nouveauD(fenetre_dilution_calibration,dilution_calibrations_id,"dilution_calibration",exec_cache);
183
184if( (blk->switch_dil&0x3f ) != 6 ) {
185// printf("**************** affiche erreur \n");
186 changecontrole(fenetre_dilution,dil_erreur,idem,idem,idem,idem,idem,"erreur");
187 }
188else {
189 static int q;
190 char ss[10]="____";
191 q=(q+1)&3;
192 ss[q]='/';
193 changecontrole(fenetre_dilution,dil_erreur,idem,idem,idem,idem,idem,ss);
194// printf("**************** pas d' erreur \n");
195 }
196
197
198// lit les switchs
199// -- les voyants de fin de course vanne principale (logique negative)
200ecritD(fenetre_dilution,dil_voyant_EVO,((blk->switch_dil&switch_EVO)?0:1));
201ecritD(fenetre_dilution,dil_voyant_EVF,((blk->switch_dil&switch_EVF)?0:1));
202
203// les retour de commande des electrovannes
204ecritD(fenetre_dilution,dil_EVB_retour,((blk->switch_dil&vanne_EVB)?1:0));
205ecritD(fenetre_dilution,dil_EVO_retour,((blk->switch_dil&vanne_EVO)?1:0));
206ecritD(fenetre_dilution,dil_EVF_retour,((blk->switch_dil&vanne_EVF)?1:0));
207ecritD(fenetre_dilution,dil_EVV_retour,((blk->switch_dil&vanne_EVV)?1:0));
208
[639]209
[645]210#define c(j,i) (1e-4*(double)parametr.nom_coef[parametr.bolo[j].numero_nom_coef].coef[i])
[639]211
212// les temperatures sur les cartes modifiées
213{
[645]214double I,V,R[8],llR;
[639]215int j,k;
216def_gains;
217k=0;
[645]218
[639]219for(j=0;(j<nb_max_bolo) && (k<4);j++)
220 {
221 if(parametr.bolo[j].bolo_code_util==bolo_thermo_simplifie)
222 {
223 I = 1e-3 * (double)dac_V(gg->reglage.bolo[j]) * 2441. / parametr.bolo[j].bolo_capa; // I en µA
224 V=0.001*bol_micro_volt(blk->temperature[k],(double)parametr.bolo[j].bolo_gain*gain_ampli(gg->reglage.bolo[j]));
225 if(I>0.0000001) R[k]=V/I; else R[k]=0; // R en ‡
[645]226 //------ calcul de l'etalonnage en temperature des cartes temperature simplifiées
227 // j=numero de bebo T[k] = R[k+4] en Kelvin
228 if ((R[k]-c(j,6))>1.)
229 {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;}
230 else llR=0;
231 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) ;
232 if(R[k+4]>9999) R[k+4]=9999;
233 k++;
[639]234 }
235 }
236
[645]237
238if(R[0]<90.) ecritD(fenetre_dilution,dil_temp
239 , "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 "
240 ,R[0],R[1],R[2],R[3],R[4],R[5],R[6],R[7]);
241else ecritD(fenetre_dilution,dil_temp
242 ,"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 "
243 ,R[0]*0.001,R[1],R[2],R[3],R[4]*1000.,R[5],R[6],R[7]);
244
245
[639]246if(litD(fenetre_dilution,dil_trace_tempe,0))// trace les sondes de niveau helium
247 {
248 int temps_cntl=numero_block(blk)*nb_per_block*2;
249 double secondes,minutes;
250 secondes=gg->periode_echantillonage*(double)temps_cntl;
251 minutes=secondes/60.;
252 if(!fenetre(fenetre_temperature_dil))
253 {
254 nouveauD(fenetre_temperature_dil,0,"temperature_dil",0);
255 selectgra(fenetre_temperature_dil);
256 strcpy(graph->xtitre,"minutes");
257 }
[645]258 tracen(fenetre_temperature_dil,8,minutes,R);
[639]259 }
260}
261
[637]262// les pressions et debits metres des injections de la dilution
263ecritD(fenetre_dilution,dil_p_d_3He,"3He: %4.1fb -> %5.2fµm/s -> %4.1fb"
[643]264 ,40. * val_multiplex(blk->ADC_dil[ p_R3]) -1.6 // 200 bars pour 5V et 1.6 bar d'offset
265 ,2. * val_multiplex(blk->ADC_dil[ d_3He]) // 10 MICRO MOLES pour 5V
266 ,20. * val_multiplex(blk->ADC_dil[ p_C3]) // 100 bars pour 5V
[637]267 );
268ecritD(fenetre_dilution,dil_p_d_4He,"4He: %4.1fb -> %5.2fµm/s -> %4.1fb"
[643]269 ,40. * val_multiplex(blk->ADC_dil[ p_R4]) // 200 bars pour 5V
270 ,8. * val_multiplex(blk->ADC_dil[ d_4He]) // 40 MICRO MOLES pour 5V
271// ,4. * val_multiplex(blk->ADC_dil[ d_4He]) // 20 MICRO MOLES pour 5V
272 ,20. * val_multiplex(blk->ADC_dil[ p_C4]) // 100 bars pour 5V
[637]273 );
274
275ecritD(fenetre_dilution,dil_p_air,"Van=%4.1fb charb=%4.1fb mmb=%4.3fb haut=%4.3fb"
[643]276 ,20.*val_multiplex(blk->ADC_dil[ p_air]),20.*val_multiplex(blk->ADC_dil[ p_charb])
277 ,0.2*val_multiplex(blk->ADC_dil[ p_memb]),0.2*val_multiplex(blk->ADC_dil[ p_haut]) );
[637]278
279ecritD(fenetre_dilution,dil_piles," Trp=%4.1fV Dil=%4.1fV / %4.1fV \rBebo=%4.1fV / %4.1fV / %4.1fV Ch=%4.1fV"
[643]280 ,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])
281 ,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])
282 ,3.8*val_multiplex(blk->ADC_dil[ p_Ch]));
[637]283
284/* valeur seuil au mini (5.3 ou 15.3V) , *17 = +5.9 V au maxi */
285
286
[643]287ecritD(fenetre_dilution,dil_p_10T,(int)(17*(2.03*val_multiplex(blk->ADC_dil[ p_10T])-5.3)));
288ecritD(fenetre_dilution,dil_p_p18D,(int)(17*(3.90*val_multiplex(blk->ADC_dil[ p_p18D])-15.3)));
289ecritD(fenetre_dilution,dil_p_m18D,(int)(17*(-3.90*val_multiplex(blk->ADC_dil[ p_m18D])-15.3)));
290ecritD(fenetre_dilution,dil_p_10B,(int)(17*(2.03*val_multiplex(blk->ADC_dil[ p_10B])-5.3)));
291ecritD(fenetre_dilution,dil_p_p18B,(int)(17*(3.93*val_multiplex(blk->ADC_dil[ p_p18B])-15.3)));
292ecritD(fenetre_dilution,dil_p_m18B,(int)(17*(-3.93*val_multiplex(blk->ADC_dil[ p_m18B])-15.3)));
293ecritD(fenetre_dilution,dil_p_Ch,(int)(10*(3.8*val_multiplex(blk->ADC_dil[ p_Ch])-25))); /* 25 a 35 V */
[637]294
295// -- les voyants de switch de piles
296ecritD(fenetre_dilution,dil_voyant_pp5,((blk->switch_dil&switch_pile_par_5)?1:0));
297ecritD(fenetre_dilution,dil_voyant_pp15,((blk->switch_dil&switch_pile_par_15)?1:0));
298
299
300ecritD(fenetre_dilution,dil_temperatures,"temp-Haut = %4.1f° %4.1f° %4.1f° %4.1f° \rtemp-Bas = %4.1f° %4.1f° %4.1f° %4.1f° "
[643]301 ,val_temperature(blk->ADC_dil[ t_h1]),val_temperature(blk->ADC_dil[ t_h2]),val_temperature(blk->ADC_dil[ t_h3]),val_temperature(blk->ADC_dil[ t_h4])
302 ,val_temperature(blk->ADC_dil[ t_b1]),val_temperature(blk->ADC_dil[ t_b2]),val_temperature(blk->ADC_dil[ t_b3]),val_temperature(blk->ADC_dil[ t_b4])
[637]303 );
304
305ecritD(fenetre_dilution,dil_t_h1,cnt_temperature(t_h1));
306ecritD(fenetre_dilution,dil_t_h2,cnt_temperature(t_h2));
307ecritD(fenetre_dilution,dil_t_h3,cnt_temperature(t_h3));
308ecritD(fenetre_dilution,dil_t_h4,cnt_temperature(t_h4));
309ecritD(fenetre_dilution,dil_t_a1,cnt_temperature(t_a1));
310ecritD(fenetre_dilution,dil_t_a2,cnt_temperature(t_a2));
311ecritD(fenetre_dilution,dil_t_a3,cnt_temperature(t_a3));
312ecritD(fenetre_dilution,dil_t_a4,cnt_temperature(t_a4));
313ecritD(fenetre_dilution,dil_t_b1,cnt_temperature(t_b1));
314ecritD(fenetre_dilution,dil_t_b2,cnt_temperature(t_b2));
315ecritD(fenetre_dilution,dil_t_b3,cnt_temperature(t_b3));
316ecritD(fenetre_dilution,dil_t_b4,cnt_temperature(t_b4));
317
318
319// les afficheurs analogiques de pression helium et pirani
320
[643]321ecritD(fenetre_dilution,dil_RP_He,(int)(2500*val_multiplex(blk->ADC_dil[ RP_He])-7500));
[637]322// je veux 10000 pour 0.8 bars soit alors que l'on a 5V pour 1 bar ou 1V = 0.2 bar
323// soit multiplier par 2500
324// puis j'enleve 0.6 bars soit 7500
[643]325ecritD(fenetre_dilution,dil_pirani,(int)(1000*val_multiplex(blk->ADC_dil[ pirani])));
[637]326
327
328
329
330
331
332// ecriture des valeurs brutes de la dilution
333if(litD(fenetre_dilution,dil_sortie_brute,0))
334 {
335 if(!fenetre(fenetre_dilution_lecture_brute)) nouveauT(fenetre_dilution_lecture_brute,0,"lecture dilution");
336
337 ecritT(fenetre_dilution_lecture_brute,fin_f,"dil_switch= %x \n",blk->switch_dil);
338
[643]339 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]340 ecritT(fenetre_dilution_lecture_brute,fin_f,"\n");
[643]341 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]342 ecritT(fenetre_dilution_lecture_brute,fin_f,"\n");
[643]343 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]344 ecritT(fenetre_dilution_lecture_brute,fin_f,"\n");
345 }
346// lecture des sondes de niveau helium
347
[639]348
[637]349 {
350 double y[8];
[639]351 static int niveau_pre;
[637]352 int niveau,fin;
[639]353// printf("dilution : numero de block = %d \n",numero_block(blk));
[637]354 changecontrole(fenetre_dilution,dil_niveau_helium,idem,idem,idem,idem,calrefcon(0,0,rouge,blanc,blanc),"");
[643]355 y[0]=val_multiplex(blk->ADC_dil[ j_he1])-calib(dil_j_he1);
356 y[1]=val_multiplex(blk->ADC_dil[ j_he2])-calib(dil_j_he2);
357 y[2]=val_multiplex(blk->ADC_dil[ j_he3])-calib(dil_j_he3);
358 y[3]=val_multiplex(blk->ADC_dil[ j_he4])-calib(dil_j_he4);
359 y[4]=val_multiplex(blk->ADC_dil[ j_he5])-calib(dil_j_he5);
360 y[5]=val_multiplex(blk->ADC_dil[ j_he6])-calib(dil_j_he6);
361 y[6]=val_multiplex(blk->ADC_dil[ j_he7])-calib(dil_j_he7);
362 y[7]=val_multiplex(blk->ADC_dil[ j_he8])-calib(dil_j_he8);
[639]363 if(blk->switch_dil & switch_helium)
[637]364 {
[639]365 changecontrole(fenetre_dilution,dil_niveau_helium,idem,idem,idem,idem,calrefcon(0,0,rouge,blanc,blanc),"");
366 niveau=0;fin=0;
367 for(i=0;i<8;i++)
368 {
369 if(y[i]>0.1) {
370 niveau++;
371 if(fin) niveau=9;
372 }
373 else fin=1;
374 }
375// if(niveau==9) changecontrole(fenetre_dilution,dil_niveau_helium,idem,idem,idem,idem
376// ,calrefcon(0,0,jaune,blanc,blanc),"");
377 if(niveau<niveau_pre) niveau_pre=0;
378 else {niveau_pre=niveau; ecritD(fenetre_dilution,dil_niveau_helium,niveau);}
379// printf("niveau=%d \n",niveau);
[637]380 }
[639]381 else changecontrole(fenetre_dilution,dil_niveau_helium,idem,idem,idem,idem,calrefcon(0,0,vert,blanc,blanc),"");
[637]382
383 if(litD(fenetre_dilution,dil_trace_helium,0))// trace les sondes de niveau helium
384 {
385 int temps_cntl=numero_block(blk)*nb_per_block*2;
386 double secondes,minutes;
387 secondes=gg->periode_echantillonage*(double)temps_cntl;
388 minutes=secondes/60.;
389 if(!fenetre(fenetre_niveau_helium))
390 {
391 nouveauD(fenetre_niveau_helium,0,"niveau_helium",0);
392 selectgra(fenetre_niveau_helium);
393 strcpy(graph->xtitre,"minutes");
394 }
395 tracen(fenetre_niveau_helium,8,minutes,y);
396 }
397 }
398
399}
400
401
Note: See TracBrowser for help on using the repository browser.