Sets C culturas /pastnat30, cfeno30, pastsem30, pastnat60, cfeno60, pastsem60, pastnat150, pastsem150, pastnat180, pastnat410/ A actividades de produção animal MPB /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 166.8, BovML 190.7/; Parameter K(U,C) Ajuda por hectare em MPB /U1.pastnat30 308, U1.cfeno30 219.9, U1.pastsem30 343.2, U1.pastnat60 276.6, U1.cfeno60 205.9, U1.pastsem60 304.5, U1.pastnat150 229.9, U1.pastsem150 246.6, U1.pastnat180 181.6, 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 257.7, U1.pastsem30 39.5, U1.pastnat60 0.0, U1.cfeno60 257.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 3000 U1.Pastsem30.pasts 300 250 1700 U1.Pastnat60.pastn 250 180 700 U1.cfeno60.feno 3000 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 -344 -535 -486 -584 bovML -360 -562 -510 -612 ; Table CMI(A,ZF) capacidade máx ing por unid pecuária R1 R2 R3 R4 bovM 707 1061 1282 1370 bovML 718 1078 1300 1390 ; 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 ajMPBn30 ajuda haMPB pastnat 30 ajMPBf30 ajuda haMPB cfeno 30 ajMPBs30 ajuda haMPB pastsem 30 ajMPBn60 ajuda haMPB pastnat 60 ajMPBf60 ajuda haMPB cfeno 60 ajMPBs60 ajuda haMPB pastsem 60 ajMPBn150 ajuda haMPB pastnat 150 ajMPBs150 ajuda haMPB pastnat 150 ajMPBn180 ajuda haMPB pastnat 180 ajMPBn410 ajuda haMPB 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 ; ajMPBn30 .. sum(U,X(U,'pastnat30')) =L= 30 ; ajMPBf30 .. sum(U,X(U,'cfeno30')) =L= 30 ; ajMPBs30 .. sum(U,X(U,'pastsem30')) =L= 30 ; ajMPBn60 .. sum(U,X(U,'pastnat60')) =L= 30 ; ajMPBf60 .. sum(U,X(U,'cfeno60')) =L= 20 ; ajMPBs60 .. sum(U,X(U,'pastsem60')) =L= 30 ; ajMPBn150 .. sum(U,X(U,'pastnat150')) =L= 90 ; ajMPBs150 .. sum(U,X(U,'pastsem150')) =L= 12 ; ajMPBn180 .. sum(U,X(U,'pastnat180')) =L= 30 ; ajMPBn410 .. sum(U,X(U,'pastnat410')) =L= 230 ; CNland .. Sum(A,Q(A)*1.18) - 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 MPB/all/ ; Solve MPB using LP maximizing Z;