root / UCalcSimulT.pas
Historique | Voir | Annoter | Télécharger (37,434 ko)
| 1 | 3 | avalancogn | unit UCalcSimulT ;
|
|---|---|---|---|
| 2 | |||
| 3 | interface
|
||
| 4 | |||
| 5 | uses
|
||
| 6 | UVariables ; |
||
| 7 | |||
| 8 | procedure CalcSimulT (Simulation : integer ;
|
||
| 9 | Variable : integer ; Variation : double ; |
||
| 10 | CycleInitVar, CycleFinVar, StadeVar : integer ; |
||
| 11 | PResSimulT : PTabResSimulT) ; |
||
| 12 | // Simulation = N? de la simulation
|
||
| 13 | // Variable = N? de la variable ? moduler
|
||
| 14 | // Variation = Coefficient de variation
|
||
| 15 | // CycleInitVar = Port?e de d?but de variation
|
||
| 16 | // CycleFinVar = Port?e de fin de variation
|
||
| 17 | // StadeVar = Stade sur lequel s'applique la variation
|
||
| 18 | // PResSimulT = Tableau de r?sultats
|
||
| 19 | |||
| 20 | implementation
|
||
| 21 | |||
| 22 | uses
|
||
| 23 | Math, UFindRec, UCalcul ; |
||
| 24 | |||
| 25 | procedure CalcSimulT (Simulation : integer ;
|
||
| 26 | Variable : integer ; Variation : double ; |
||
| 27 | CycleInitVar, CycleFinVar, StadeVar : integer ; |
||
| 28 | PResSimulT : PTabResSimulT) ; |
||
| 29 | var
|
||
| 30 | i, j : integer ; |
||
| 31 | ok : boolean ; |
||
| 32 | Stade, StadeInit, StadeFin, TypeStade, Cycle : integer ; |
||
| 33 | Jour, Duree, Age, Ecart : integer ; |
||
| 34 | RuleSeqAli : array[1..MAX_RULE] of RecRuleSeqAliT ; |
||
| 35 | RuleRation : array[1..MAX_RULE] of RecRuleRationT ; |
||
| 36 | RuleLoge : array[1..MAX_RULE] of RecRuleLogeT ; |
||
| 37 | Unite : integer ; |
||
| 38 | NumRuleLoge, NumRuleSeqAli, NumRuleRation : integer ; |
||
| 39 | RuleSeqAliInit, RuleRationInit, RuleLogeInit : integer ; |
||
| 40 | RecCC1, RecCC2 : CompositionChimique ; |
||
| 41 | TabAAtotal1, TabAAtotal2, TabCUDAA1, TabCUDAA2 : array[0..12] of double ; |
||
| 42 | PctAli1, PctAli2, Quantite, Ingere, IngereSec, IngSec1, IngSec2 : double ; |
||
| 43 | TabAAdig : array[0..12] of double ; |
||
| 44 | NRAA : array[0..14] of double ; |
||
| 45 | EqLysd : double ; |
||
| 46 | Temperature, Activite : integer ; |
||
| 47 | PVTot, PVVTot, PVNet, PVVNet, P2, P2Var : double ; |
||
| 48 | p, l, PDNet, PDTot, LD, R1, R2, RA : double ; |
||
| 49 | NesTotaux, Sevres, PdsNais, PdsSev, PdsPt, GMQPt, GainLact, GainProfil, GainSimul : double ; |
||
| 50 | CorrPP, CorrCU, CorrJ, CorrNR, TeneurLysPot : double ; // Facteur de correction (poids port?e, contenus ut?rins, jour, r?tention azot?e)
|
||
| 51 | CoeffEntGest, CoeffNR, CoeffEntLact, RatioLP : double ; |
||
| 52 | EMIng, EMEnt, EMEntInit, EMThe, EMAct, EMProd, EMProt, EMLip, EMRes : double ; |
||
| 53 | ENProd, ENProdCum, ENProdTot, ENProt, ENLip, ENRes, ENResCum : double ; |
||
| 54 | NRPotAA, NRPot, NR, NRProd, NRProdCum, NRProdTot, NRRes : double ; |
||
| 55 | PProd, PProdCum, PProdTot : double ; |
||
| 56 | NExcF, NExcU : double ; // Azote excr?t? f?c?s, urine
|
||
| 57 | Chaleur : double ; |
||
| 58 | MATIng, MATdIng : double ; |
||
| 59 | begin
|
||
| 60 | // Initialisation du tableau de r?sultat
|
||
| 61 | for i := 1 to NB_COL_TRUIE do |
||
| 62 | for j := 1 to MAX_LIG_TRUIE do |
||
| 63 | PResSimulT.TabResult[i, j] := 0 ;
|
||
| 64 | // Initialisation de la simulation
|
||
| 65 | PResSimulT.NbJSim := 0 ;
|
||
| 66 | if Simulation = -1 |
||
| 67 | then // Profil |
||
| 68 | begin
|
||
| 69 | // D?but et fin de simulation
|
||
| 70 | StadeInit := 0 ;
|
||
| 71 | StadeFin := 23 ;
|
||
| 72 | // Age initial
|
||
| 73 | Age := PProfilT.Truies[1].AgeSail ;
|
||
| 74 | end
|
||
| 75 | else // Simulation |
||
| 76 | begin
|
||
| 77 | // Initialisation des param?tres de simulation
|
||
| 78 | PSimulT := ListSimulT[FindIdxSimulT (FindNomSimulT (Simulation))] ; |
||
| 79 | // Initialisation du profil animal
|
||
| 80 | PProfilT := ListProfilT[FindIdxProfilT (FindNomProfilT (PSimulT.Profil))] ; |
||
| 81 | // D?but et fin de simulation
|
||
| 82 | StadeInit := PSimulT.StadeInit ; |
||
| 83 | StadeFin := PSimulT.StadeFin ; |
||
| 84 | // Age initial
|
||
| 85 | if StadeInit mod 3 = 0 |
||
| 86 | then // Gestation : age saillie |
||
| 87 | Age := PProfilT.Truies[StadeInit div 3 + 1].AgeSail |
||
| 88 | else // Lactation : age mise-bas |
||
| 89 | Age := PProfilT.Truies[StadeInit div 3 + 1].AgeSail + DureeGest ; |
||
| 90 | // Sensibilit? (Age initial)
|
||
| 91 | if (Variable = 0) |
||
| 92 | then
|
||
| 93 | Age := Round (Age * Variation) ; |
||
| 94 | end ;
|
||
| 95 | // Coefficients de calibrage
|
||
| 96 | CoeffEntGest := PProfilT.CoefEntretienGest ; |
||
| 97 | CoeffEntLact := PProfilT.CoefEntretienLact ; |
||
| 98 | RatioLP := PProfilT.RatioLipProt ; |
||
| 99 | // Boucle stades
|
||
| 100 | for Stade := StadeInit to StadeFin - 1 do |
||
| 101 | begin
|
||
| 102 | Cycle := Stade div 3 + 1 ; |
||
| 103 | TypeStade := Stade mod 3 ; |
||
| 104 | if (Stade = StadeInit) or ((Simulation <> -1) and ((PSimulT.OptInit = 0) or ((PSimulT.OptInit = 1) and (TypeStade = 0)))) |
||
| 105 | then // Initialisation PV et P2 |
||
| 106 | begin
|
||
| 107 | // Initialisation du poids vif
|
||
| 108 | case TypeStade of |
||
| 109 | 0 : // Gestation |
||
| 110 | PVTot := PProfilT.Truies[Cycle].PdsSail ; |
||
| 111 | 1 : // Lactation |
||
| 112 | PVTot := PProfilT.Truies[Cycle].PdsApMB ; |
||
| 113 | end ;
|
||
| 114 | // Sensibilit? (Poids vif initial)
|
||
| 115 | if (Variable = 1) and (Stade = StadeInit) |
||
| 116 | then
|
||
| 117 | PVTot := PVTot * Variation ; |
||
| 118 | // Initialisation de l'?paisseur de lard
|
||
| 119 | case TypeStade of |
||
| 120 | 0 : // Gestation |
||
| 121 | P2 := PProfilT.Truies[Cycle].P2Sail ; |
||
| 122 | 1 : // Lactation |
||
| 123 | P2 := PProfilT.Truies[Cycle].P2MB ; |
||
| 124 | end ;
|
||
| 125 | // Sensibilit? (Epaisseur de lard initiale)
|
||
| 126 | if (Variable = 2) and (Stade = StadeInit) |
||
| 127 | then
|
||
| 128 | P2 := P2 * Variation ; |
||
| 129 | PVVTot := PVTot * PV2PVV ; |
||
| 130 | p := 2.28 + (0.178 * PVVTot) - (0.333 * P2) ; |
||
| 131 | l := -26.4 + (0.221 * PVVTot) + (1.331 * P2) ; |
||
| 132 | ENResCum := -1024 + (13.66 * PVVTot) + (45.94 * P2) ; |
||
| 133 | end ;
|
||
| 134 | // Initialisation de la s?quence alimentaire
|
||
| 135 | if (Simulation = -1) or (PSimulT.SeqAli[Cycle] = -1) |
||
| 136 | then // Profil |
||
| 137 | PSeqAliT := ListSeqAliT[FindIdxSeqAliT (FindNomSeqAliT (PProfilT.SeqAli))] |
||
| 138 | else // Simulation |
||
| 139 | PSeqAliT := ListSeqAliT[FindIdxSeqAliT (FindNomSeqAliT (PSimulT.SeqAli[Cycle]))] ; |
||
| 140 | case TypeStade of |
||
| 141 | 0 : // Gestation |
||
| 142 | for i := 1 to PSeqAliT.NbRuleGest do |
||
| 143 | RuleSeqAli[i] := PSeqAliT.RuleGest[i] ; |
||
| 144 | 1 : // Lactation |
||
| 145 | for i := 1 to PSeqAliT.NbRuleLact do |
||
| 146 | RuleSeqAli[i] := PSeqAliT.RuleLact[i] ; |
||
| 147 | else // ISSF |
||
| 148 | for i := 1 to PSeqAliT.NbRuleISSF do |
||
| 149 | RuleSeqAli[i] := PSeqAliT.RuleISSF[i] ; |
||
| 150 | end ;
|
||
| 151 | NumRuleSeqAli := 1 ;
|
||
| 152 | RuleSeqAliInit := 1 ;
|
||
| 153 | // Initialisation du plan de rationnement
|
||
| 154 | if (Simulation = -1) or (PSimulT.Ration[Cycle] = -1) |
||
| 155 | then // Profil |
||
| 156 | begin
|
||
| 157 | with RuleRation[1] do |
||
| 158 | begin
|
||
| 159 | ModeFin := -1 ;
|
||
| 160 | if TypeStade = 1 |
||
| 161 | then // Lactation |
||
| 162 | begin
|
||
| 163 | Equation := 3 ; // Curvilin?aire |
||
| 164 | a := PProfilT.Lact[Cycle] / 2 ; // Initial |
||
| 165 | c := PProfilT.Lact[Cycle] ; // Moyenne
|
||
| 166 | d := PProfilT.DureeLact ; // Dur?e
|
||
| 167 | end
|
||
| 168 | else // Gestation (ou ISSF) |
||
| 169 | begin
|
||
| 170 | Equation := 0 ; // Constant |
||
| 171 | a := PProfilT.Gest[Cycle] ; |
||
| 172 | end ;
|
||
| 173 | end ;
|
||
| 174 | Unite := PProfilT.Unite ; |
||
| 175 | end
|
||
| 176 | else // Simulation |
||
| 177 | begin
|
||
| 178 | PRationT := ListRationT[FindIdxRationT (FindNomRationT (PSimulT.Ration[Cycle]))] ; |
||
| 179 | case TypeStade of |
||
| 180 | 0 : // Gestation |
||
| 181 | for i := 1 to PRationT.NbRuleGest do |
||
| 182 | RuleRation[i] := PRationT.RuleGest[i] ; |
||
| 183 | 1 : // Lactation |
||
| 184 | for i := 1 to PRationT.NbRuleLact do |
||
| 185 | RuleRation[i] := PRationT.RuleLact[i] ; |
||
| 186 | else // ISSF |
||
| 187 | for i := 1 to PRationT.NbRuleISSF do |
||
| 188 | RuleRation[i] := PRationT.RuleISSF[i] ; |
||
| 189 | end ;
|
||
| 190 | case TypeStade of |
||
| 191 | 0 : // Gestation |
||
| 192 | Unite := PRationT.UniteGest ; |
||
| 193 | 1 : // Lactation |
||
| 194 | Unite := PRationT.UniteLact ; |
||
| 195 | else // ISSF |
||
| 196 | Unite := PRationT.UniteISSF ; |
||
| 197 | end ;
|
||
| 198 | end ;
|
||
| 199 | NumRuleRation := 1 ;
|
||
| 200 | RuleRationInit := 1 ;
|
||
| 201 | // Initialisation du logement
|
||
| 202 | if Simulation = -1 |
||
| 203 | then // Profil |
||
| 204 | PLogeT := ListLogeT[FindIdxLogeT (FindNomLogeT (PProfilT.Loge))] |
||
| 205 | else // Simulation |
||
| 206 | PLogeT := ListLogeT[FindIdxLogeT (FindNomLogeT (PSimulT.Logement))] ; |
||
| 207 | case TypeStade of |
||
| 208 | 0 : // Gestation |
||
| 209 | for i := 1 to PLogeT.NbRuleGest do |
||
| 210 | RuleLoge[i] := PLogeT.RuleGest[i] ; |
||
| 211 | 1 : // Lactation |
||
| 212 | for i := 1 to PLogeT.NbRuleLact do |
||
| 213 | RuleLoge[i] := PLogeT.RuleLact[i] ; |
||
| 214 | else // ISSF |
||
| 215 | for i := 1 to PLogeT.NbRuleISSF do |
||
| 216 | RuleLoge[i] := PLogeT.RuleISSF[i] ; |
||
| 217 | end ;
|
||
| 218 | NumRuleLoge := 1 ;
|
||
| 219 | RuleLogeInit := 1 ;
|
||
| 220 | // Dur?e du stade
|
||
| 221 | case TypeStade of |
||
| 222 | 0 : // Gestation |
||
| 223 | Duree := DureeGest ; |
||
| 224 | 1 : // Lactation |
||
| 225 | if Simulation = -1 |
||
| 226 | then // Profil |
||
| 227 | Duree := PProfilT.DureeLact |
||
| 228 | else // Simulation |
||
| 229 | Duree := PSimulT.DureeLact ; |
||
| 230 | else // ISSF |
||
| 231 | Duree := PProfilT.Truies[Cycle + 1].AgeSail - PProfilT.Truies[Cycle].AgeSail - DureeGest - PProfilT.DureeLact ;
|
||
| 232 | end ;
|
||
| 233 | // Coefficients de calibrage
|
||
| 234 | CoeffNR := Max (PProfilT.CoefNR[Cycle], 0) ;
|
||
| 235 | // Initialisation du stade
|
||
| 236 | PVNet := PVTot ; |
||
| 237 | case TypeStade of |
||
| 238 | 0 : // Gestation |
||
| 239 | begin
|
||
| 240 | NesTotaux := PProfilT.Porcelets[Cycle].NesTotaux ; |
||
| 241 | PdsNais := PProfilT.Porcelets[Cycle].PdsNais ; |
||
| 242 | PdsPt := NesTotaux * PdsNais ; |
||
| 243 | CorrPP := PdsPt * 1000 / Exp (9.095 - 17.69 * Exp (-0.0305 * DureeGest) + 0.0878 * NesTotaux) ; |
||
| 244 | CorrCU := (0.3 + (1.329 * PdsPt)) * 1000 / Exp (8.621 - 21.02 * Exp (-0.05302 * DureeGest) + 0.1114 * NesTotaux) ; |
||
| 245 | ENProdTot := Exp (11.72 - 8.62 * Exp (-0.01382 * 0) + 0.0932 * NesTotaux) / 1000 ; |
||
| 246 | ENProdTot := ENProdTot * CorrPP ; // Correction
|
||
| 247 | ENProdCum := ENProdTot ; |
||
| 248 | PProdTot := Exp (8.621 - 21.02 * Exp (-0.05302 * 0) + 0.1114 * NesTotaux) ; |
||
| 249 | PProdTot := PProdTot * CorrCu ; // Correction
|
||
| 250 | PProdCum := PProdTot ; |
||
| 251 | NRProdCum := 0 ;
|
||
| 252 | NRProdTot := Exp (8.09 - 8.71 * Exp (-0.01494 * 0) + 0.0872 * NesTotaux) / 6.25 ; |
||
| 253 | NRProdTot := NRProdTot * CorrPP ; // Correction
|
||
| 254 | NRProdCum := NRProdTot ; |
||
| 255 | EMEntInit := Power (PVVTot, 0.75) * EEGest * CoeffEntGest ;
|
||
| 256 | end ;
|
||
| 257 | 1 : // Lactation |
||
| 258 | begin
|
||
| 259 | NesTotaux := PProfilT.Porcelets[Cycle].NesTotaux ; |
||
| 260 | PdsNais := PProfilT.Porcelets[Cycle].PdsNais ; |
||
| 261 | Sevres := PProfilT.Porcelets[Cycle].Sevres ; |
||
| 262 | PdsSev := PProfilT.Porcelets[Cycle].PdsSev ; |
||
| 263 | if Simulation <> -1 |
||
| 264 | then // Ajustement du poids de sevrage ? la dur?e de lactation |
||
| 265 | begin
|
||
| 266 | GainProfil := (273 * PProfilT.DureeLact + (289 / 0.38) * Exp (-0.38 * PProfilT.DureeLact) - 289 / 0.38) / 1000 ; |
||
| 267 | GainSimul := (273 * PSimulT.DureeLact + (289 / 0.38) * Exp (-0.38 * PSimulT.DureeLact) - 289 / 0.38) / 1000 ; |
||
| 268 | PdsSev := PdsNais + GainSimul * (PdsSev - PdsNais) / GainProfil ; |
||
| 269 | end ;
|
||
| 270 | ENProdCum := 0 ;
|
||
| 271 | NRProdCum := 0 ;
|
||
| 272 | PProdCum := 0 ;
|
||
| 273 | // Gain de poids de la port?e
|
||
| 274 | GMQPt := (PdsSev - PdsNais) * Sevres / Duree ; |
||
| 275 | end ;
|
||
| 276 | else // ISSF |
||
| 277 | begin
|
||
| 278 | NesTotaux := 0 ;
|
||
| 279 | PdsNais := 0 ;
|
||
| 280 | EMEntInit := Power (PVVTot, 0.75) * EEGest * CoeffEntGest ;
|
||
| 281 | end ;
|
||
| 282 | end ;
|
||
| 283 | // Eviter le d?bordement du tableau de r?sultats
|
||
| 284 | if PResSimulT.NbJSim + Duree > MAX_LIG_TRUIE
|
||
| 285 | then
|
||
| 286 | Duree := MAX_LIG_TRUIE - PResSimulT.NbJSim ; |
||
| 287 | // Boucle jours
|
||
| 288 | for Jour := 1 to Duree do |
||
| 289 | begin
|
||
| 290 | Inc (PResSimulT.NbJSim) ; |
||
| 291 | PResSimulT.TabResult[78, PResSimulT.NbJSim] := Jour ;
|
||
| 292 | PResSimulT.TabResult[1, PResSimulT.NbJSim] := Age ;
|
||
| 293 | PResSimulT.TabResult[2, PResSimulT.NbJSim] := Cycle ;
|
||
| 294 | PResSimulT.TabResult[3, PResSimulT.NbJSim] := TypeStade ;
|
||
| 295 | PResSimulT.TabResult[58, PResSimulT.NbJSim] := PVTot ;
|
||
| 296 | PResSimulT.TabResult[60, PResSimulT.NbJSim] := PVNet ;
|
||
| 297 | PResSimulT.TabResult[67, PResSimulT.NbJSim] := P2 ;
|
||
| 298 | PResSimulT.TabResult[62, PResSimulT.NbJSim] := p ;
|
||
| 299 | PResSimulT.TabResult[65, PResSimulT.NbJSim] := l ;
|
||
| 300 | PVVNet := PVNet * PV2PVV ; |
||
| 301 | /////////////////////////
|
||
| 302 | // Apports journaliers //
|
||
| 303 | /////////////////////////
|
||
| 304 | // Aliment(s) distribu?(s)
|
||
| 305 | repeat
|
||
| 306 | ok := TRUE ; |
||
| 307 | with RuleSeqAli[NumRuleSeqAli] do |
||
| 308 | if ModeFin = 0 |
||
| 309 | then // Dur?e |
||
| 310 | if Jour - RuleSeqAliInit >= ValFin then ok := FALSE ; |
||
| 311 | if not (ok) |
||
| 312 | then // Changement de r?gle |
||
| 313 | begin
|
||
| 314 | Inc (NumRuleSeqAli) ; |
||
| 315 | RuleSeqAliInit := Jour ; |
||
| 316 | end ;
|
||
| 317 | until ok ;
|
||
| 318 | PResSimulT.TabResult[4, PResSimulT.NbJSim] := NumRuleSeqAli ;
|
||
| 319 | with RuleSeqAli[NumRuleSeqAli] do |
||
| 320 | begin
|
||
| 321 | PResSimulT.TabResult[7, PResSimulT.NbJSim] := NumAli1 ;
|
||
| 322 | PResSimulT.TabResult[8, PResSimulT.NbJSim] := NumAli2 ;
|
||
| 323 | // Composition aliment 1
|
||
| 324 | if NumAli1 = -1 |
||
| 325 | then
|
||
| 326 | begin
|
||
| 327 | RecCC1 := CCVide ; |
||
| 328 | for i := 0 to 12 do |
||
| 329 | TabAAtotal1[i] := 0 ;
|
||
| 330 | for i := 0 to 12 do |
||
| 331 | TabCUDAA1[i] := 0 ;
|
||
| 332 | end
|
||
| 333 | else
|
||
| 334 | begin
|
||
| 335 | PAliment := ListAliment[FindIdxAliment (FindNomAliment (NumAli1))] ; |
||
| 336 | RecCC1 := PAliment.CC ; |
||
| 337 | for i := 0 to 12 do |
||
| 338 | TabAAtotal1[i] := PAliment.AAtotal[i] ; |
||
| 339 | for i := 0 to 12 do |
||
| 340 | TabCUDAA1[i] := PAliment.CUDAA[i] ; |
||
| 341 | end ;
|
||
| 342 | // Composition aliment 2
|
||
| 343 | if NumAli2 = -1 |
||
| 344 | then
|
||
| 345 | begin
|
||
| 346 | RecCC2 := CCVide ; |
||
| 347 | for i := 0 to 12 do |
||
| 348 | TabAAtotal2[i] := 0 ;
|
||
| 349 | for i := 0 to 12 do |
||
| 350 | TabCUDAA2[i] := 0 ;
|
||
| 351 | end
|
||
| 352 | else
|
||
| 353 | begin
|
||
| 354 | PAliment := ListAliment[FindIdxAliment (FindNomAliment (NumAli2))] ; |
||
| 355 | RecCC2 := PAliment.CC ; |
||
| 356 | for i := 0 to 12 do |
||
| 357 | TabAAtotal2[i] := PAliment.AAtotal[i] ; |
||
| 358 | for i := 0 to 12 do |
||
| 359 | TabCUDAA2[i] := PAliment.CUDAA[i] ; |
||
| 360 | end ;
|
||
| 361 | // Calcul des % aliments
|
||
| 362 | if PctAli1Init = PctAli1Fin
|
||
| 363 | then
|
||
| 364 | PctAli1 := PctAli1Init |
||
| 365 | else // Transition |
||
| 366 | begin
|
||
| 367 | Ecart := PctAli1Fin - PctAli1Init ; |
||
| 368 | if ModeFin = 0 |
||
| 369 | then // Dur?e |
||
| 370 | PctAli1 := PctAli1Init + (Jour - RuleSeqAliInit) * Ecart / (ValFin - 1)
|
||
| 371 | else // Fin |
||
| 372 | PctAli1 := PctAli1Init + (Jour - RuleSeqAliInit) * Ecart / (Duree - RuleSeqAliInit) ; |
||
| 373 | end ;
|
||
| 374 | end ;
|
||
| 375 | PResSimulT.TabResult[9, PResSimulT.NbJSim] := PctAli1 ;
|
||
| 376 | PctAli2 := 100 - PctAli1 ;
|
||
| 377 | PResSimulT.TabResult[10, PResSimulT.NbJSim] := PctAli2 ;
|
||
| 378 | // Quantit?(s) distribu?e(s)
|
||
| 379 | repeat
|
||
| 380 | ok := TRUE ; |
||
| 381 | with RuleRation[NumRuleRation] do |
||
| 382 | if ModeFin = 0 |
||
| 383 | then // Dur?e |
||
| 384 | if Jour - RuleRationInit >= ValFin then ok := FALSE ; |
||
| 385 | if not (ok) |
||
| 386 | then // Changement de r?gle |
||
| 387 | begin
|
||
| 388 | Inc (NumRuleRation) ; |
||
| 389 | RuleRationInit := Jour ; |
||
| 390 | end ;
|
||
| 391 | until ok ;
|
||
| 392 | PResSimulT.TabResult[5, PResSimulT.NbJSim] := NumRuleRation ;
|
||
| 393 | with RuleRation[NumRuleRation] do |
||
| 394 | begin
|
||
| 395 | // Calcul des quantit?s
|
||
| 396 | case Equation of |
||
| 397 | 0 : // Constant |
||
| 398 | Quantite := a ; |
||
| 399 | 1 : // Lin?aire |
||
| 400 | Quantite := a + b * (Jour - RuleRationInit) ; |
||
| 401 | 2 : // Lin?aire-plateau |
||
| 402 | Quantite := LPvaleur (a, b, c, Jour - RuleRationInit + 1, d) ;
|
||
| 403 | 3 : // Curvilin?aire |
||
| 404 | Quantite := CLvaleur (a, c, Jour - RuleRationInit + 1, d) ;
|
||
| 405 | else
|
||
| 406 | Quantite := 0 ;
|
||
| 407 | end ;
|
||
| 408 | // Convertion de ED, EM, EN en quantit? si besoin
|
||
| 409 | case Unite of |
||
| 410 | 1 : // ED (MJ/j) |
||
| 411 | Ingere := Quantite |
||
| 412 | / (PctAli1 / 100 * RecCC1.ED_T * RecCC1.MS / 1000 |
||
| 413 | + PctAli2 / 100 * RecCC2.ED_T * RecCC2.MS / 1000) ; |
||
| 414 | 2 : // EM (MJ/j) |
||
| 415 | Ingere := Quantite |
||
| 416 | / (PctAli1 / 100 * RecCC1.EM_T * RecCC1.MS / 1000 |
||
| 417 | + PctAli2 / 100 * RecCC2.EM_T * RecCC2.MS / 1000) ; |
||
| 418 | 3 : // EN (MJ/j) |
||
| 419 | Ingere := Quantite |
||
| 420 | / (PctAli1 / 100 * RecCC1.EN_T * RecCC1.MS / 1000 |
||
| 421 | + PctAli2 / 100 * RecCC2.EN_T * RecCC2.MS / 1000) ; |
||
| 422 | 4 : // MS (kg/j) |
||
| 423 | Ingere := Quantite |
||
| 424 | / (PctAli1 / 100 * RecCC1.MS / 1000 |
||
| 425 | + PctAli2 / 100 * RecCC2.MS / 1000) ; |
||
| 426 | else // Quantit? (kg/j) |
||
| 427 | Ingere := Quantite ; |
||
| 428 | end ;
|
||
| 429 | end ;
|
||
| 430 | // Sensibilit? (Quantit? d'aliment)
|
||
| 431 | if (Variable = 5) |
||
| 432 | and ((Cycle >= CycleInitVar) and (Cycle <= CycleFinVar)) |
||
| 433 | and ((StadeVar = 2) or (TypeStade = StadeVar)) |
||
| 434 | then
|
||
| 435 | Ingere := Ingere * Variation ; |
||
| 436 | PResSimulT.TabResult[11, PResSimulT.NbJSim] := Ingere ;
|
||
| 437 | // Calcul des apports journaliers
|
||
| 438 | IngSec1 := Ingere * PctAli1 / 100 * RecCC1.MS / 1000 ; |
||
| 439 | IngSec2 := Ingere * PctAli2 / 100 * RecCC2.MS / 1000 ; |
||
| 440 | IngereSec := IngSec1 + IngSec2 ; |
||
| 441 | PResSimulT.TabResult[12, PResSimulT.NbJSim] := IngereSec ;
|
||
| 442 | // Energie brute ing?r?e
|
||
| 443 | PResSimulT.TabResult[13, PResSimulT.NbJSim] := IngSec1 * RecCC1.EB
|
||
| 444 | + IngSec2 * RecCC2.EB ; |
||
| 445 | // Energie digestible ing?r?e
|
||
| 446 | PResSimulT.TabResult[71, PResSimulT.NbJSim] := IngSec1 * RecCC1.ED_T
|
||
| 447 | + IngSec2 * RecCC2.ED_T ; |
||
| 448 | // Energie m?tabolisable ing?r?e
|
||
| 449 | EMIng := IngSec1 * RecCC1.EM_T + IngSec2 * RecCC2.EM_T ; |
||
| 450 | PResSimulT.TabResult[72, PResSimulT.NbJSim] := EMIng ;
|
||
| 451 | // Energie nette ing?r?e
|
||
| 452 | PResSimulT.TabResult[73, PResSimulT.NbJSim] := IngSec1 * RecCC1.EN_T
|
||
| 453 | + IngSec2 * RecCC2.EN_T ; |
||
| 454 | // Sensibilit? (Acides amin?s)
|
||
| 455 | if (Variable > 5) |
||
| 456 | and ((Cycle >= CycleInitVar) and (Cycle <= CycleFinVar)) |
||
| 457 | and ((StadeVar = 2) or (TypeStade = StadeVar)) |
||
| 458 | then
|
||
| 459 | begin
|
||
| 460 | TabAAtotal1[Variable - 5] := TabAAtotal1[Variable - 5] * Variation ; |
||
| 461 | TabAAtotal2[Variable - 5] := TabAAtotal2[Variable - 5] * Variation ; |
||
| 462 | if (Variable = 8) or (Variable = 12) |
||
| 463 | then // met+cys ou phe+tyr |
||
| 464 | begin
|
||
| 465 | TabAAtotal1[Variable - 6] := TabAAtotal1[Variable - 6] * Variation ; |
||
| 466 | TabAAtotal2[Variable - 6] := TabAAtotal2[Variable - 6] * Variation ; |
||
| 467 | end ;
|
||
| 468 | end ;
|
||
| 469 | // Acides amin?s totaux ing?r?s
|
||
| 470 | for i:= 0 to 12 do |
||
| 471 | PResSimulT.TabResult[14 + i, PResSimulT.NbJSim] := IngSec1 * TabAAtotal1[i]
|
||
| 472 | + IngSec2 * TabAAtotal2[i] ; |
||
| 473 | // Acides amin?s digestibles ing?r?s
|
||
| 474 | for i := 0 to 12 do |
||
| 475 | TabAAdig[i] := IngSec1 * TabAAtotal1[i] * TabCUDAA1[i] / 100
|
||
| 476 | + IngSec2 * TabAAtotal2[i] * TabCUDAA2[i] / 100 ;
|
||
| 477 | for i := 0 to 12 do |
||
| 478 | PResSimulT.TabResult[27 + i, PResSimulT.NbJSim] := TabAAdig[i] ;
|
||
| 479 | ////////////////
|
||
| 480 | // Simulation //
|
||
| 481 | ////////////////
|
||
| 482 | // Logement
|
||
| 483 | repeat
|
||
| 484 | ok := TRUE ; |
||
| 485 | with RuleLoge[NumRuleLoge] do |
||
| 486 | if ModeFin = 0 |
||
| 487 | then // Dur?e |
||
| 488 | if Jour - RuleLogeInit >= ValFin then ok := FALSE ; |
||
| 489 | if not (ok) |
||
| 490 | then // Changement de r?gle |
||
| 491 | begin
|
||
| 492 | Inc (NumRuleLoge) ; |
||
| 493 | RuleLogeInit := Jour ; |
||
| 494 | end ;
|
||
| 495 | until ok ;
|
||
| 496 | PResSimulT.TabResult[6, PResSimulT.NbJSim] := NumRuleLoge ;
|
||
| 497 | with RuleLoge[NumRuleLoge] do |
||
| 498 | begin
|
||
| 499 | if (Sol = 1) |
||
| 500 | then // Paill? |
||
| 501 | Temperature := Temp + 3
|
||
| 502 | else
|
||
| 503 | Temperature := Temp ; |
||
| 504 | // Sensibilit? (Temp?rature ambiante)
|
||
| 505 | if (Variable = 3) |
||
| 506 | and ((Cycle >= CycleInitVar) and (Cycle <= CycleFinVar)) |
||
| 507 | and ((StadeVar = 2) or (TypeStade = StadeVar)) |
||
| 508 | then
|
||
| 509 | Temperature := Round (Temperature * Variation) ; |
||
| 510 | Activite := Act ; |
||
| 511 | // Sensibilit? (Temps debout)
|
||
| 512 | if (Variable = 4) |
||
| 513 | and ((Cycle >= CycleInitVar) and (Cycle <= CycleFinVar)) |
||
| 514 | and ((StadeVar = 2) or (TypeStade = StadeVar)) |
||
| 515 | then
|
||
| 516 | Activite := Round (Activite * Variation) ; |
||
| 517 | end ;
|
||
| 518 | PResSimulT.TabResult[56, PResSimulT.NbJSim] := Temperature ;
|
||
| 519 | PResSimulT.TabResult[57, PResSimulT.NbJSim] := Activite ;
|
||
| 520 | case TypeStade of |
||
| 521 | 0 : // Gestation |
||
| 522 | begin
|
||
| 523 | // Energie m?tabolisable entretien
|
||
| 524 | EMEnt := Power (PVVTot, 0.75) * EEGest * CoeffEntGest ;
|
||
| 525 | PResSimulT.TabResult[74, PResSimulT.NbJSim] := EMEnt ;
|
||
| 526 | // Energie m?tabolisable thermor?gulation
|
||
| 527 | if (RuleLoge[NumRuleLoge].Typ = 1) |
||
| 528 | then // Collectif |
||
| 529 | if Temperature < TCICol
|
||
| 530 | then
|
||
| 531 | EMThe := (TCICol - Temperature) * Power (PVVTot, 0.75) * EThCol / 1000 |
||
| 532 | else
|
||
| 533 | EMThe := 0
|
||
| 534 | else
|
||
| 535 | if Temperature < TCIInd
|
||
| 536 | then
|
||
| 537 | EMThe := (TCIInd - Temperature) * Power (PVVTot, 0.75) * EThInd / 1000 |
||
| 538 | else
|
||
| 539 | EMThe := 0 ;
|
||
| 540 | PResSimulT.TabResult[75, PResSimulT.NbJSim] := EMThe ;
|
||
| 541 | // Energie m?tabolisable activit?
|
||
| 542 | EMAct := 0.3 * Power (PVVTot, 0.75) * (Activite - 240) / 1000 ; |
||
| 543 | PResSimulT.TabResult[76, PResSimulT.NbJSim] := EMAct ;
|
||
| 544 | // Energie contenus ut?rins
|
||
| 545 | ENProdTot := Exp (11.72 - 8.62 * Exp (-0.01382 * Jour) + 0.0932 * NesTotaux) / 1000 ; |
||
| 546 | ENProdTot := ENProdTot * CorrPP ; // Correction
|
||
| 547 | ENProd := ENProdTot - ENProdCum ; |
||
| 548 | ENProdCum := ENProdTot ; |
||
| 549 | EMProd := ENProd / KU ; |
||
| 550 | PResSimulT.TabResult[77, PResSimulT.NbJSim] := EMProd ;
|
||
| 551 | // Poids contenus ut?rins
|
||
| 552 | PProdTot := Exp (8.621 - 21.02 * Exp (-0.05302 * Jour) + 0.1114 * NesTotaux) ; |
||
| 553 | PProdTot := PProdTot * CorrCu ; // Correction
|
||
| 554 | PProd := PProdTot - PProdCum ; |
||
| 555 | PProdCum := PProdTot ; |
||
| 556 | // Azote contenus ut?rins
|
||
| 557 | NRProdTot := Exp (8.09 - 8.71 * Exp (-0.01494 * Jour) + 0.0872 * NesTotaux) / 6.25 ; |
||
| 558 | NRProdTot := NRProdTot * CorrPP ; // Correction
|
||
| 559 | NRProd := NRProdTot - NRProdCum ; |
||
| 560 | NRProdCum := NRProdTot ; |
||
| 561 | // R?tention azot?e potentielle
|
||
| 562 | Case Cycle of |
||
| 563 | 1 : // Port?e 1 |
||
| 564 | CorrNR := 0.5708 ;
|
||
| 565 | 2 : // Port?e 2 |
||
| 566 | CorrNR := 0.4345 ;
|
||
| 567 | else
|
||
| 568 | CorrNR := 0.3664 ;
|
||
| 569 | end ;
|
||
| 570 | if (Jour < 98) |
||
| 571 | then
|
||
| 572 | CorrJ := Jour |
||
| 573 | else
|
||
| 574 | CorrJ := Jour - 6 / 16 * (Jour - 98) ; |
||
| 575 | NRPot := NRProd + CoeffNR * (-0.43 + 45.92 * (CorrJ / 100) - 105.35 * Power ((CorrJ / 100), 2) + 64.388 * Power ((CorrJ / 100), 3) + CorrNR * (EMIng - EMThe - EMEntInit)) * 0.85 ; |
||
| 576 | PResSimulT.TabResult[40, PResSimulT.NbJSim] := NRPot * 6.25 ; |
||
| 577 | // R?tention possible par les apports d'acides amin?s
|
||
| 578 | TeneurLysPot := (0.060 * NRProd + 0.070 * (NRPot - NRProd)) / NRPot ; |
||
| 579 | NRAA[0] := (TabAAdig[0] - (AAm75[0] * Power (PVTot, 0.75) + IngereSec * AAendogene[0])) * kAA[0] / AAbody[0] / 6.25 ; |
||
| 580 | NRAA[1] := (TabAAdig[1] - (AAm75[1] * Power (PVTot, 0.75) + IngereSec * AAendogene[1])) * KLysGest / TeneurLysPot / 6.25 ; |
||
| 581 | for i := 2 to 12 do |
||
| 582 | begin
|
||
| 583 | if ProtIdGest[i] = 0 |
||
| 584 | then
|
||
| 585 | EqLysd := 0
|
||
| 586 | else
|
||
| 587 | EqLysd := TabAAdig[i] * 100 / ProtIdGest[i] ;
|
||
| 588 | NRAA[i] := (EqLysd - (AAm75[i] * Power (PVTot, 0.75) + IngereSec * AAendogene[i])) * KLysGest / TeneurLysPot / 6.25 ; |
||
| 589 | end ;
|
||
| 590 | // met+cys
|
||
| 591 | if ProtIdGest[2] + ProtIdGest[3] = 0 |
||
| 592 | then
|
||
| 593 | EqLysd := 0
|
||
| 594 | else
|
||
| 595 | EqLysd := (TabAAdig[2] + TabAAdig[3]) * 100 / (ProtIdGest[2] + ProtIdGest[3]) ; |
||
| 596 | NRAA[13] := (EqLysd - ((AAm75[2] + AAm75[3]) * Power (PVTot, 0.75) + IngereSec * (AAendogene[2] + AAendogene[3]))) * KLysGest / TeneurLysPot / 6.25 ; |
||
| 597 | // phe+tyr
|
||
| 598 | if ProtIdGest[6] + ProtIdGest[7] = 0 |
||
| 599 | then
|
||
| 600 | EqLysd := 0
|
||
| 601 | else
|
||
| 602 | EqLysd := (TabAAdig[6] + TabAAdig[7]) * 100 / (ProtIdGest[6] + ProtIdGest[7]) ; |
||
| 603 | NRAA[14] := (EqLysd - ((AAm75[6] + AAm75[7]) * Power (PVTot, 0.75) + IngereSec * (AAendogene[6] + AAendogene[7]))) * KLysGest / TeneurLysPot / 6.25 ; |
||
| 604 | for i := 0 to 14 do |
||
| 605 | if NRAA[i] > 0 |
||
| 606 | then
|
||
| 607 | PResSimulT.TabResult[41 + i, PResSimulT.NbJSim] := NRAA[i] * 6.25 ; |
||
| 608 | // Premier facteur limitant le d?pot de prot?ines
|
||
| 609 | NRPotAA := NRAA[0] ;
|
||
| 610 | for i := 1 to 14 do |
||
| 611 | if (i <> 3) and (i <> 7) // Ignorer 'cys' et 'tyr' |
||
| 612 | then
|
||
| 613 | if (NRAA[i] > 0) and (NRAA[i] < NRPotAA) |
||
| 614 | then
|
||
| 615 | NRPotAA := NRAA[i] ; |
||
| 616 | // R?tention azot?e
|
||
| 617 | NR := Min (NRPot, NRPotAA) ; |
||
| 618 | // R?tention azot?e maternelle
|
||
| 619 | NRRes := NR - NRProd ; |
||
| 620 | // Energie retenue dans les tissus maternels et production de chaleur
|
||
| 621 | ENProt := NRRes * 6.25 * 23.8 / 1000 ; |
||
| 622 | EMProt := ENProt / KProt ; |
||
| 623 | EMLip := EMIng - EMEnt - EMThe - EMAct - EMProd - EMProt ; |
||
| 624 | if EMLip > 0 |
||
| 625 | then // D?p?t |
||
| 626 | begin
|
||
| 627 | ENLip := EMLip * KLip ; |
||
| 628 | Chaleur := EMEnt + EMThe + EMAct + EMProd * (1 - KU) + EMProt * (1 - KProt) + EMLip * (1 - KLip) ; |
||
| 629 | end
|
||
| 630 | else // Mobilisation |
||
| 631 | begin
|
||
| 632 | ENLip := (EMIng - EMEnt - EMThe - EMAct - EMProd - EMProt) / KResGest ; |
||
| 633 | Chaleur := EMEnt + EMThe + EMAct + EMProd * (1 - KU) + EMProt * (1 - KProt) - EMLip * (1 - KResGest); |
||
| 634 | end;
|
||
| 635 | PResSimulT.TabResult[70, PResSimulT.NbJSim] := Chaleur ;
|
||
| 636 | ENRes := ENProt + ENLip ; |
||
| 637 | ENResCum := ENResCum + ENRes ; |
||
| 638 | // Excr?tion azot?e
|
||
| 639 | NExcF := (MATIng - MATdIng) / 6.25 ;
|
||
| 640 | NExcU := (MATdIng / 6.25) - NR ;
|
||
| 641 | // D?p?ts de prot?ines et lipides
|
||
| 642 | PDNet := NRRes * 6.25 ;
|
||
| 643 | PResSimulT.TabResult[63, PResSimulT.NbJSim] := PDNet ;
|
||
| 644 | PDTot := NR * 6.25 ;
|
||
| 645 | PResSimulT.TabResult[64, PResSimulT.NbJSim] := PDTot ;
|
||
| 646 | LD := ENLip / 39.7 * 1000 ; |
||
| 647 | PResSimulT.TabResult[66, PResSimulT.NbJSim] := LD ;
|
||
| 648 | p := p + PDNet / 1000 ;
|
||
| 649 | l := l + LD / 1000 ;
|
||
| 650 | // Poids vif
|
||
| 651 | R1 := -((p - 2.28) * 45.94 + (ENResCum + 1024) * 0.333) / APN ; |
||
| 652 | R2 := -((p - 2.28) * 1.3311 + (l + 24.5) * 0.333) / BPN ; |
||
| 653 | PVVNet := (R1 + R2) / 2 ;
|
||
| 654 | PVVTot := PVVNet + PProdTot / 1000 ;
|
||
| 655 | PVTot := PVVTot / PV2PVV ; |
||
| 656 | PResSimulT.TabResult[59, PResSimulT.NbJSim] := PVTot
|
||
| 657 | - PResSimulT.TabResult[58, PResSimulT.NbJSim] ;
|
||
| 658 | PVNet := PVVNet / PV2PVV ; |
||
| 659 | PResSimulT.TabResult[61, PResSimulT.NbJSim] := PVNet
|
||
| 660 | - PResSimulT.TabResult[60, PResSimulT.NbJSim] ;
|
||
| 661 | // Epaisseur de lard dorsale
|
||
| 662 | R1 := ((p - 2.28) * 13.66 - (ENResCum + 1024) * 0.178) / APN ; |
||
| 663 | R2 := ((p - 2.28) * 0.221 - (l + 26.4) * 0.178) / BPN ; |
||
| 664 | P2 := (R1 + R2) / 2 ;
|
||
| 665 | PResSimulT.TabResult[68, PResSimulT.NbJSim] := P2
|
||
| 666 | - PResSimulT.TabResult[67, PResSimulT.NbJSim] ;
|
||
| 667 | end ;
|
||
| 668 | 1 : // Lactation |
||
| 669 | begin
|
||
| 670 | // Energie m?tabolisable entretien
|
||
| 671 | EMEnt := Power (PVVNet, 0.75) * EELact * CoeffEntLact ;
|
||
| 672 | PResSimulT.TabResult[74, PResSimulT.NbJSim] := EMEnt ;
|
||
| 673 | // Energie m?tabolisable thermor?gulation
|
||
| 674 | EMThe := 0 ;
|
||
| 675 | PResSimulT.TabResult[75, PResSimulT.NbJSim] := EMThe ;
|
||
| 676 | // Energie m?tabolisable activit?
|
||
| 677 | EMAct := 0 ;
|
||
| 678 | PResSimulT.TabResult[76, PResSimulT.NbJSim] := EMAct ;
|
||
| 679 | // Production d'?nergie dans le lait
|
||
| 680 | RA := 0.0000023096 * Power (Duree, 4) - 0.00027619 * Power (Duree, 3) + 0.012889 * Power (Duree, 2) - 0.28116 * Duree + 4.799 ; |
||
| 681 | ENProd := (20.6 * GMQPt * 1000 - 376 * Sevres) * RA * Exp (-0.025 * Jour) * Exp (-Exp (0.5 - 0.1 * Jour)) / 1000 ; |
||
| 682 | ENProdCum := ENPRodCum + ENProd ; |
||
| 683 | EMProd := ENProd / KL ; |
||
| 684 | PResSimulT.TabResult[77, PResSimulT.NbJSim] := EMProd ;
|
||
| 685 | // Production d'azote dans le lait
|
||
| 686 | NRProd := (0.0257 * GMQPt * 1000 + 0.42 * Sevres) * RA * Exp (-0.025 * Jour) * Exp (- Exp (0.5 - 0.1 * Jour)) ; |
||
| 687 | NRProdCum := NRProdCum + NRProd ; |
||
| 688 | PProd := NRProd * 6.38 / 50 ; |
||
| 689 | PResSimulT.TabResult[69, PResSimulT.NbJSim] := PProd ;
|
||
| 690 | PProdCum := PProdCum + PProd ; |
||
| 691 | PResSimulT.TabResult[40, PResSimulT.NbJSim] := NRProd * 6.38 ; |
||
| 692 | // R?tention possible par les apports d'acides amin?s
|
||
| 693 | PResSimulT.TabResult[41, PResSimulT.NbJSim] := (TabAAdig[0] - (AAm75[0] * Power (PVVTot, 0.75) + IngereSec * AAendogene[0])) * kAA[0] / AAbody[0] ; |
||
| 694 | NRAA[1] := -14.2 + 1.335 * TabAAdig[1] - 0.629 * NRProd ; |
||
| 695 | PResSimulT.TabResult[42, PResSimulT.NbJSim] := (-14.2 + 1.335 * TabAAdig[1]) / 0.629 * 6.25 ; |
||
| 696 | for i := 2 to 12 do |
||
| 697 | begin
|
||
| 698 | EqLysd := TabAAdig[i] * 100 / ProtIdLact[i] ;
|
||
| 699 | NRAA[i] := -14.2 + 1.335 * EqLysd - 0.629 * NRProd ; |
||
| 700 | if EqLysd > 0 |
||
| 701 | then
|
||
| 702 | PResSimulT.TabResult[41 + i, PResSimulT.NbJSim] := (-14.2 + 1.335 * EqLysd) / 0.629 * 6.25 ; |
||
| 703 | end ;
|
||
| 704 | // met+cys
|
||
| 705 | EqLysd := (TabAAdig[2] + TabAAdig[3]) * 100 / (ProtIdLact[2] + ProtIdLact[3]) ; |
||
| 706 | NRAA[13] := -14.2 + 1.335 * EqLysd - 0.629 * NRProd ; |
||
| 707 | if EqLysd > 0 |
||
| 708 | then
|
||
| 709 | PResSimulT.TabResult[54, PResSimulT.NbJSim] := (-14.2 + 1.335 * EqLysd) / 0.629 * 6.25 ; |
||
| 710 | // phe+tyr
|
||
| 711 | EqLysd := (TabAAdig[6] + TabAAdig[7]) * 100 / (ProtIdLact[6] + ProtIdLact[7]) ; |
||
| 712 | NRAA[14] := -14.2 + 1.335 * EqLysd - 0.629 * NRProd ; |
||
| 713 | if EqLysd > 0 |
||
| 714 | then
|
||
| 715 | PResSimulT.TabResult[55, PResSimulT.NbJSim] := (-14.2 + 1.335 * EqLysd) / 0.629 * 6.25 ; |
||
| 716 | // Premier facteur limitant le d?pot de prot?ines
|
||
| 717 | NR := NRAA[1] ;
|
||
| 718 | for i := 2 to 14 do |
||
| 719 | if (i <> 3) and (i <> 7) // Ignorer 'cys' et 'tyr' |
||
| 720 | then
|
||
| 721 | if NRAA[i] < NR
|
||
| 722 | then
|
||
| 723 | NR := NRAA[i] ; |
||
| 724 | // R?tention azot?e maternelle
|
||
| 725 | if NR > 0 |
||
| 726 | then
|
||
| 727 | NR := 0 ;
|
||
| 728 | // Mobilisation des r?serves ?nerg?tiques
|
||
| 729 | EMRes := EMIng - EMEnt - EMThe - EMAct - EMProd ; |
||
| 730 | if EMRes > 0 |
||
| 731 | then // D?p?t |
||
| 732 | begin
|
||
| 733 | ENRes := EMRes * KLip ; |
||
| 734 | ENProt := 0 ;
|
||
| 735 | ENLip := ENRes ; |
||
| 736 | Chaleur := EMEnt + EMThe + EMAct + EMProd * (1 - KL) + EMRes * (1 - KLip) |
||
| 737 | end
|
||
| 738 | else // Mobilisation |
||
| 739 | begin
|
||
| 740 | ENRes := EMRes * KL / KResL ; |
||
| 741 | ENProt := Min (NR * 6.25 * 23.8 / 1000, ENRes * (1 - RatioLP)) ; |
||
| 742 | ENLip := ENRes - ENProt ; |
||
| 743 | Chaleur := EMEnt + EMThe + EMAct + (EMProd + EMRes) * (1 - KL) - EMRes * KL / KResL * (1 - KResL) ; |
||
| 744 | end ;
|
||
| 745 | ENResCum := ENResCum + ENRes ; |
||
| 746 | PResSimulT.TabResult[70, PResSimulT.NbJSim] := Chaleur ;
|
||
| 747 | NRRes := ENProt / 6.25 / 23.8 * 1000 ; |
||
| 748 | // Excr?tion azot?e
|
||
| 749 | NExcF := (MATIng - MATdIng) / 6.25 ;
|
||
| 750 | NExcU := (MATdIng / 6.25) - NR - NRProd ;
|
||
| 751 | // D?p?ts de prot?ines et lipides
|
||
| 752 | PDNet := NRRes * 6.25 ;
|
||
| 753 | PResSimulT.TabResult[63, PResSimulT.NbJSim] := PDNet ;
|
||
| 754 | PDTot := NR * 6.25 ;
|
||
| 755 | PResSimulT.TabResult[64, PResSimulT.NbJSim] := PDTot ;
|
||
| 756 | LD := ENLip / 39.7 * 1000 ; |
||
| 757 | PResSimulT.TabResult[66, PResSimulT.NbJSim] := LD ;
|
||
| 758 | p := p + PDNet / 1000 ;
|
||
| 759 | l := l + LD / 1000 ;
|
||
| 760 | // Poids vif
|
||
| 761 | R1 := -((p - 2.28) * 45.94 + (ENResCum + 1024) * 0.333) / APN ; |
||
| 762 | R2 := -((p - 2.28) * 1.3311 + (l + 24.5) * 0.333) / BPN ; |
||
| 763 | PVVNet := (R1 + R2) / 2 ;
|
||
| 764 | PVVTot := PVVNet ; |
||
| 765 | PVTot := PVVTot / PV2PVV ; |
||
| 766 | PResSimulT.TabResult[59, PResSimulT.NbJSim] := PVTot
|
||
| 767 | - PResSimulT.TabResult[58, PResSimulT.NbJSim] ;
|
||
| 768 | PVNet := PVVNet / PV2PVV ; |
||
| 769 | PResSimulT.TabResult[61, PResSimulT.NbJSim] := PVNet
|
||
| 770 | - PResSimulT.TabResult[60, PResSimulT.NbJSim] ;
|
||
| 771 | // Epaisseur de lard dorsale
|
||
| 772 | R1 := ((p - 2.28) * 13.66 - (ENResCum + 1024) * 0.178) / APN ; |
||
| 773 | R2 := ((p - 2.28) * 0.221 - (l + 26.4) * 0.178) / BPN ; |
||
| 774 | P2 := (R1 + R2) / 2 ;
|
||
| 775 | PResSimulT.TabResult[68, PResSimulT.NbJSim] := P2
|
||
| 776 | - PResSimulT.TabResult[67, PResSimulT.NbJSim] ;
|
||
| 777 | end ;
|
||
| 778 | else // ISSF |
||
| 779 | begin
|
||
| 780 | // Energie m?tabolisable entretien
|
||
| 781 | EMEnt := Power (PVVTot, 0.75) * EEGest * CoeffEntGest ;
|
||
| 782 | PResSimulT.TabResult[74, PResSimulT.NbJSim] := EMEnt ;
|
||
| 783 | // Energie m?tabolisable thermor?gulation
|
||
| 784 | case RuleLoge[NumRuleLoge].Typ of |
||
| 785 | 0 : // Logement individuel |
||
| 786 | if Temperature < TCIInd
|
||
| 787 | then
|
||
| 788 | EMThe := (TCIInd - Temperature) * Power (PVVTot, 0.75) * EThInd / 1000 |
||
| 789 | else
|
||
| 790 | EMThe := 0 ;
|
||
| 791 | 1 : // Logement Collectif |
||
| 792 | if Temperature < TCICol
|
||
| 793 | then
|
||
| 794 | EMThe := (TCICol - Temperature) * Power (PVVTot, 0.75) * EThCol / 1000 |
||
| 795 | else
|
||
| 796 | EMThe := 0 ;
|
||
| 797 | else
|
||
| 798 | EMThe := 0 ;
|
||
| 799 | end ;
|
||
| 800 | PResSimulT.TabResult[75, PResSimulT.NbJSim] := EMThe ;
|
||
| 801 | // Energie m?tabolisable acivit?
|
||
| 802 | EMAct := 0.3 * Power (PVVTot, 0.75) * (Activite - 240) / 1000 ; |
||
| 803 | PResSimulT.TabResult[76, PResSimulT.NbJSim] := EMAct ;
|
||
| 804 | // R?tention azot?e potentielle
|
||
| 805 | Case Cycle of |
||
| 806 | 1 : // Port?e 1 |
||
| 807 | CorrNR := 0.5708 ;
|
||
| 808 | 2 : // Port?e 2 |
||
| 809 | CorrNR := 0.4345 ;
|
||
| 810 | else
|
||
| 811 | CorrNR := 0.3664 ;
|
||
| 812 | end;
|
||
| 813 | if (Jour < 98) |
||
| 814 | then
|
||
| 815 | CorrJ := Jour |
||
| 816 | else
|
||
| 817 | CorrJ := Jour - 6 / 16 * (Jour - 98) ; |
||
| 818 | NRPot := (-0.43 + 45.92 * (CorrJ / 100) - 105.35 * Power ((CorrJ / 100), 2) + 64.388 * Power ((CorrJ / 100), 3) + CorrNR * (EMIng - EMThe - EMEntInit)) * 0.85 ; |
||
| 819 | PResSimulT.TabResult[40, PResSimulT.NbJSim] := NRPot * 6.25 ; |
||
| 820 | // R?tention possible par les apports d'acides amin?s
|
||
| 821 | NRAA[0] := (TabAAdig[0] - (AAm75[0] * Power (PVVTot, 0.75) + IngereSec * AAendogene[0])) * kAA[0] / AAbody[0] / 6.25 ; |
||
| 822 | NRAA[1] := (TabAAdig[1] - 0.036 * Power (PVVTot, 0.75)) * KLysGest / 0.065 / 6.25 ; |
||
| 823 | for i := 2 to 12 do |
||
| 824 | begin
|
||
| 825 | if ProtIdGest[i] = 0 |
||
| 826 | then
|
||
| 827 | EqLysd := 0
|
||
| 828 | else
|
||
| 829 | EqLysd := TabAAdig[i] * 100 / ProtIdGest[i] ;
|
||
| 830 | NRAA[i] := (EqLysd - 0.036 * Power (PVVTot, 0.75)) * KLysGest / 0.065 / 6.25 ; |
||
| 831 | end ;
|
||
| 832 | // met+cys
|
||
| 833 | if ProtIdGest[2] + ProtIdGest[3] = 0 |
||
| 834 | then
|
||
| 835 | EqLysd := 0
|
||
| 836 | else
|
||
| 837 | EqLysd := (TabAAdig[2] + TabAAdig[3]) * 100 / (ProtIdGest[2] + ProtIdGest[3]) ; |
||
| 838 | NRAA[13] := (EqLysd - 0.036 * Power (PVVTot, 0.75)) * KLysGest / 0.065 / 6.25 ; |
||
| 839 | // phe+tyr
|
||
| 840 | if ProtIdGest[6] + ProtIdGest[7] = 0 |
||
| 841 | then
|
||
| 842 | EqLysd := 0
|
||
| 843 | else
|
||
| 844 | EqLysd := (TabAAdig[6] + TabAAdig[7]) * 100 / (ProtIdGest[6] + ProtIdGest[7]) ; |
||
| 845 | NRAA[14] := (EqLysd - 0.036 * Power (PVVTot, 0.75)) * KLysGest / 0.065 / 6.25 ; |
||
| 846 | for i := 0 to 14 do |
||
| 847 | if NRAA[i] > 0 |
||
| 848 | then
|
||
| 849 | PResSimulT.TabResult[41 + i, PResSimulT.NbJSim] := NRAA[i] * 6.25 ; |
||
| 850 | // Premier facteur limitant le d?pot de prot?ines
|
||
| 851 | NRPotAA := NRAA[0] ;
|
||
| 852 | for i := 1 to 14 do |
||
| 853 | if (i <> 3) and (i <> 7) // Ignorer 'cys' et 'tyr' |
||
| 854 | then
|
||
| 855 | if (NRAA[i] > 0) and (NRAA[i] < NRPotAA) |
||
| 856 | then
|
||
| 857 | NRPotAA := NRAA[i] ; |
||
| 858 | // R?tention azot?e
|
||
| 859 | NR := Min (NRPot, NRPotAA) ; |
||
| 860 | // R?tention azot?e maternelle
|
||
| 861 | NRRes := NR ; |
||
| 862 | // Energie retenue dans les tissus maternels et production de chaleur
|
||
| 863 | ENProt := NRRes * 6.25 * 23.8 / 1000 ; |
||
| 864 | EMProt := ENProt / KProt ; |
||
| 865 | EMLip := EMIng - EMEnt - EMThe - EMAct - EMProt ; |
||
| 866 | if EMLip > 0 |
||
| 867 | then
|
||
| 868 | begin
|
||
| 869 | ENLip := EMLip * KLip ; |
||
| 870 | Chaleur := EMEnt + EMThe + EMAct + EMProt * (1 - KProt) + EMLip * (1 - KLip) ; |
||
| 871 | end
|
||
| 872 | else
|
||
| 873 | begin
|
||
| 874 | ENLip := (EMIng - EMEnt - EMThe - EMAct - EMProt) / KResGest ; |
||
| 875 | Chaleur := EMEnt + EMThe + EMAct + EMProt * (1 - KProt) - EMLip * (1 - KResGest); |
||
| 876 | end;
|
||
| 877 | PResSimulT.TabResult[70, PResSimulT.NbJSim] := Chaleur ;
|
||
| 878 | ENRes := ENProt + ENLip ; |
||
| 879 | ENResCum := ENResCum + ENRes ; |
||
| 880 | // Excr?tion azot?e
|
||
| 881 | NExcF := (MATIng - MATdIng) / 6.25 ;
|
||
| 882 | NExcU := (MATdIng / 6.25) - NR ;
|
||
| 883 | // D?p?ts de prot?ines et lipides
|
||
| 884 | PDNet := NRRes * 6.25 ;
|
||
| 885 | PResSimulT.TabResult[63, PResSimulT.NbJSim] := PDNet ;
|
||
| 886 | PDTot := NR * 6.25 ;
|
||
| 887 | PResSimulT.TabResult[64, PResSimulT.NbJSim] := PDTot ;
|
||
| 888 | LD := ENLip / 39.7 * 1000 ; |
||
| 889 | PResSimulT.TabResult[66, PResSimulT.NbJSim] := LD ;
|
||
| 890 | p := p + PDNet / 1000 ;
|
||
| 891 | l := l + LD / 1000 ;
|
||
| 892 | // Poids vif
|
||
| 893 | R1 := -((p - 2.28) * 45.94 + (ENResCum + 1024) * 0.333) / APN ; |
||
| 894 | R2 := -((p - 2.28) * 1.3311 + (l + 24.5) * 0.333) / BPN ; |
||
| 895 | PVVNet := (R1 + R2) / 2 ;
|
||
| 896 | PVVTot := PVVNet ; |
||
| 897 | PVTot := PVVTot / PV2PVV ; |
||
| 898 | PResSimulT.TabResult[59, PResSimulT.NbJSim] := PVTot
|
||
| 899 | - PResSimulT.TabResult[58, PResSimulT.NbJSim] ;
|
||
| 900 | PVNet := PVVNet / PV2PVV ; |
||
| 901 | PResSimulT.TabResult[61, PResSimulT.NbJSim] := PVNet
|
||
| 902 | - PResSimulT.TabResult[60, PResSimulT.NbJSim] ;
|
||
| 903 | // Epaisseur de lard dorsale
|
||
| 904 | R1 := ((p - 2.28) * 13.66 - (ENResCum + 1024) * 0.178) / APN ; |
||
| 905 | R2 := ((p - 2.28) * 0.221 - (l + 26.4) * 0.178) / BPN ; |
||
| 906 | P2 := (R1 + R2) / 2 ;
|
||
| 907 | PResSimulT.TabResult[68, PResSimulT.NbJSim] := P2
|
||
| 908 | - PResSimulT.TabResult[67, PResSimulT.NbJSim] ;
|
||
| 909 | end ;
|
||
| 910 | end ;
|
||
| 911 | Inc (Age) ; |
||
| 912 | end ; // Fin boucle jours |
||
| 913 | if TypeStade = 0 |
||
| 914 | then // Gestation |
||
| 915 | PVTot := PVNet ; // Mise bas
|
||
| 916 | end ; // Fin boucle stades |
||
| 917 | end ;
|
||
| 918 | |||
| 919 | end. |