Sets C culturas /pastnat30, cfeno30, pastsem30, pastnat60, cfeno60, pastsem60, pastnat150, pastsem150, pastnat180, pastnat410/ A actividades de produção animal PRODI /BovM, BovML/ U unidades de utilização /U1/ F alimentos /pastn, feno, pasts/ FP(F) alimentos produzidos na empresa /pastn, feno, pasts/ Tr transferência de alimentos /T1-2, T2-3, T3-4, T4-1/ ZF Períodos de alimentação /R1*R4/ ZT trabalho alocado a actividade pecuária /vaqueiroM, vaqueiroML/ UC(U,C) crossing possibilities UUtilização com culturas /U1.pastnat30, U1.cfeno30, U1.pastsem30, U1.pastnat60, U1.cfeno60, U1.pastsem60, U1.pastnat150, U1.pastsem150, U1.pastnat180, U1.pastnat410/ CF(C,F) crossing possibilities culturas alimentos da empresa /pastnat30.pastn, cfeno30.feno, pastsem30.pasts, pastnat60.pastn, cfeno60.feno, pastsem60.pasts, pastnat150.pastn, pastsem150.pasts, pastnat180.pastn, pastnat410.pastn/ Parameter R(A) Rendimento das actividades pecuárias /BovM 143, BovML 164 /; Parameter K(U,C) Ajuda por hectare em PRODI /U1.pastnat30 247.5, U1.cfeno30 186.8, U1.pastsem30 269.6, U1.pastnat60 228, U1.cfeno60 179.4, U1.pastsem60 245.7, U1.pastnat150 198.8, U1.pastsem150 209.8, U1.pastnat180 169.5, U1.pastnat410 150/ ; Scalar CV custo hora contrat vaqueiro /10.6/; Parameter ChaC(U,C) custo VARIÁVEL por ha cultura /U1.pastnat30 0.0, U1.cfeno30 305.7, U1.pastsem30 39.5, U1.pastnat60 0.0, U1.cfeno60 305.7, U1.pastsem60 39.5, U1.pastnat150 0.0, U1.pastsem150 39.5, U1.pastnat180 0.0, U1.pastnat410 0.0/; Table AF(U,C,F,ZF) produção alimento por cultura e por hectare e período R1 R2 R3 R4 U1.Pastnat30.pastn 250 180 700 U1.cfeno30.feno 4400 U1.Pastsem30.pasts 300 250 1700 U1.Pastnat60.pastn 250 180 700 U1.cfeno60.feno 4400 U1.Pastsem60.pasts 300 250 1700 U1.Pastnat150.pastn 250 180 700 U1.Pastsem150.pasts 300 250 1700 U1.Pastnat180.pastn 250 180 700 U1.Pastnat410.pastn 250 180 700 ; Table COF(U,A,F,ZF) Alim disponível por classe anim e período de alimentação R1 R2 R3 R4 U1.BovM.pastn -1 -1 -1 -1 U1.BovM.feno -1 -1 -1 -1 U1.BovM.pasts -1 -1 -1 -1 U1.BovML.pastn -1 -1 -1 -1 U1.BovML.feno -1 -1 -1 -1 U1.BovML.pasts -1 -1 -1 -1 ; Table TRF(F,TR,ZF) transferências de alimentos entre períodos R1 R2 R3 R4 pastn.T1-2 -1 1 pastn.T2-3 -1 1 pastn.T3-4 -1 0.85 feno.T3-4 -1 1 feno.T4-1 1 -1 feno.T1-2 -1 1 pasts.T1-2 -1 1 pasts.T2-3 -1 1 pasts.T3-4 -1 0.85 ; Table UFL(A,F,ZF) UFL por kgMS (past e feno) e kg alimento (ACC e TLuz) R1 R2 R3 R4 BovM.pastn 0.8 0.75 0.5 0.35 BovM.pasts 0.85 0.85 0.6 0.40 BovM.feno 0.6 0.6 0.6 0.6 BovML.pastn 0.8 0.75 0.5 0.35 BovML.pasts 0.85 0.85 0.6 0.40 BovML.feno 0.6 0.6 0.6 0.6 ; Table DR(A,F,ZF) Coeficiente matéria seca R1 R2 R3 R4 BovM.pastn -1 -1 -1 -1 BovM.feno -1 -1 -1 -1 BovM.pasts -1 -1 -1 -1 BovML.pastn -1 -1 -1 -1 BovML.feno -1 -1 -1 -1 BovML.pasts -1 -1 -1 -1 ; Table NUFL(A,ZF) necessidades UFL por unidade pecuária R1 R2 R3 R4 bovM -333 -520 -472 -567 bovML -348 -543 -494 -593 ; Table CMI(A,ZF) capacidade máx ing por unid pecuária R1 R2 R3 R4 bovM 710 1065 1287 1376 bovML 721 1081 1306 1396 ; Parameter HT(A,ZT) horas trabalho por unid pecuária /bovM.vaqueiroM 13, bovML.vaqueiroML 13/ ; Variables Z resultado de exploração X(U,C) ha por cultura Q(A) nº unidades pecuárias Transf(F,TR) transferência de alimentos Consm(A,F,ZF) consumo de alimentos Contvaq(ZT) contrat vaqueiro ; Positive variables X,Q,Transf,Consm,Contvaq ; Equations income função objectivo land1 restrição área feno land2 restrição área pastsem land3 restrição área total ajPRODIn30 ajuda haPRODI pastnat 30 ajPRODIf30 ajuda haPRODI cfeno 30 ajPRODIs30 ajuda haPRODI pastsem 30 ajPRODIn60 ajuda haPRODI pastnat 60 ajPRODIf60 ajuda haPRODI cfeno 60 ajPRODIs60 ajuda haPRODI pastsem 60 ajPRODIn150 ajuda haPRODI pastnat 150 ajPRODIs150 ajuda haPRODI pastnat 150 ajPRODIn180 ajuda haPRODI pastnat 180 ajPRODIn410 ajuda haPRODI pastnat 410 CNland Blabor eq trab para bovinos Feedbal(F,ZF) restrição de alimentos UFLBal(A,ZF) UFL balanço de consumo alimentos Maxing(A,ZF) cap máx ing ; INCOME .. Sum(A,Q(A)*R(A))+ sum((U,C),K(U,C)*X(U,C))-(Contvaq('vaqueiroM')*CV)- (Contvaq('vaqueiroML')*CV)- Sum((U,C),ChaC(U,C)*X(U,C)) =E= Z ; Land1 .. Sum((U)$UC(U,'cfeno30'),X(U,'cfeno30'))+ Sum((U)$UC(U,'cfeno60'),X(U,'cfeno60')) =L= 50 ; Land2 .. Sum((U)$UC(U,'pastsem60'),X(U,'pastsem60')) +Sum((U)$UC(U,'pastsem60'),X(U,'pastsem60')) +Sum((U)$UC(U,'pastsem150'),X(U,'pastsem150')) =L= 72 ; Land3 .. Sum(U,X(U,'pastnat30'))+Sum(U,X(U,'cfeno30'))+Sum(U,X(U,'pastsem30'))+ Sum(U,X(U,'pastnat60'))+Sum(U,X(U,'cfeno60'))+Sum(U,X(U,'pastsem60')) +Sum(U,X(U,'pastnat150'))+Sum(U,X(U,'pastsem150'))+Sum(U,X(U,'pastnat180')) +Sum(U,X(U,'pastnat410')) =L= 410 ; ajPRODIn30 .. sum(U,X(U,'pastnat30')) =L= 30 ; ajPRODIf30 .. sum(U,X(U,'cfeno30')) =L= 30 ; ajPRODIs30 .. sum(U,X(U,'pastsem30')) =L= 30 ; ajPRODIn60 .. sum(U,X(U,'pastnat60')) =L= 30 ; ajPRODIf60 .. sum(U,X(U,'cfeno60')) =L= 20 ; ajPRODIs60 .. sum(U,X(U,'pastsem60')) =L= 30 ; ajPRODIn150 .. sum(U,X(U,'pastnat150')) =L= 90 ; ajPRODIs150 .. sum(U,X(U,'pastsem150')) =L= 12 ; ajPRODIn180 .. sum(U,X(U,'pastnat180')) =L= 30 ; ajPRODIn410 .. sum(U,X(U,'pastnat410')) =L= 230 ; CNland .. Sum(A,Q(A)*1.16) - Sum(U,X(U,'pastnat30'))-Sum(U,X(U,'cfeno30'))-Sum(U,X(U,'pastsem30'))- Sum(U,X(U,'pastnat60'))-Sum(U,X(U,'cfeno60'))-Sum(U,X(U,'pastsem60')) -Sum(U,X(U,'pastnat150'))-Sum(U,X(U,'pastsem150'))-Sum(U,X(U,'pastnat180')) =E= 0 ; BLabor .. Sum((A),HT(A,'vaqueiroM')*Q('bovM'))+ Sum((A),HT(A,'vaqueiroML')*Q('bovML')) =L= (contvaq('vaqueiroM')+ contvaq('vaqueiroML')) ; Feedbal(F,ZF) .. Sum((U,C),AF(U,C,F,ZF)*X(U,C)$UC(U,C))+ Sum(TR,TRF(F,TR,ZF)*Transf(F,TR))+ Sum((U,A),COF(U,A,F,ZF)*Consm(A,F,ZF)$A(A)) =G= 0 ; UFLBAL(A,ZF) .. NUFL(A,ZF)*Q(A)+Sum(FP,UFL(A,FP,ZF)*consm(A,FP,ZF)$A(A)) =G= 0 ; Maxing(A,ZF) .. (CMI(A,ZF)*Q(A))+Sum(FP,DR(A,FP,ZF)*Consm(A,FP,ZF)$A(A)) =G= 0 ; model PRODI/all/ ; Solve PRODI using LP maximizing Z;