0

The goal of the docplex model below is to choose the trucks with total minimum cost to fulfill all orders. I treat the truck type as 1 truck in the following code, but in fact, a truck type can have infinite number of truck, in that case, how should I modified the constraints and objective functions to fit that?

import pandas as pd
import docplex.mp
from docplex.mp.model import Model
df_order = pd.DataFrame([[1,5000],[2,1000],[3,2000]], columns =['OrderID','Weight']) 
orders = df_order['OrderID'].values.tolist()
df_order.set_index('OrderID', inplace=True)

df_truck_types = pd.DataFrame([[1,'TYPE1',20000,'SP TRANSPORTS',40000],[2,'TYPE2',20000,'SP TRANSPORTS',40000],[3,'TYPE3',10000,'SP TRANSPORTS',30000]], columns =['TruckTypeID','VehicleType','Capacity','Transporter','TruckCost']) 
truck_types = df_truck_types['TruckTypeID'].values.tolist()
df_truck_types.set_index('TruckTypeID', inplace=True)

def get_order(o):
    return df_order.loc[o]

def get_truck_type(f):
    return df_truck_types.loc[f]

# Decision Variable and Objective Function
mdl = Model(name='itc_load_planning')

Assignment = mdl.binary_var_matrix(orders, truck_types, name='Assignment')
#for objective 1
TruckCost=mdl.sum(Assignment[o,f] * get_truck_type(f).TruckCost for o in orders for f in truck_types)

mdl.minimize(TruckCost)

# C1
for o in orders:
    mdl.add_constraint(mdl.sum(Assignment[o,f] for f in truck_types)==1, "C1")
    
# solve model
print("===solution===")
si = mdl.solve()
mdl.print_solution()
william007
  • 17,375
  • 25
  • 118
  • 194

1 Answers1

1

This is not an OPL model but a docplex python model.

To deal with trucks instead of truck types I would change

Assignment = mdl.binary_var_matrix(orders, truck_types, name='Assignment')

into

Assignment = mdl.binary_var_matrix(orders, trucks, name='Assignment')

And a truck in trucks would have a capacity so you will have to add a constraint about the capacity limit in a truck

In OPL you can write this very quickly:

tuple order
{
  key int id;
  int weight;
}

{order} orders={<1,5000>,<2,1000>,<3,2000>};

tuple trucktype
{
  key int TruckTypeID;
  string VehicleType;
  int Capacity;
  string Transporter;
  int TruckCost;
};

{trucktype} trucktypes=
{
  <1,"TYPE1",20000,"SP TRANSPORTS",40000>,
  <2,"TYPE2",20000,"SP TRANSPORTS",40000>,
  <3,"TYPE3",10000,"SP TRANSPORTS",30000>
};

int nbTruckPerType=3;

tuple truck
{
  trucktype TruckType;
  int rank;
}

{truck} trucks={<i,j> | i in trucktypes,j in 1..nbTruckPerType};

dvar boolean assign[orders][trucks];
dvar boolean istruckused[trucks];


minimize sum(t in trucks) istruckused[t]*t.TruckType.TruckCost;
subject to
{
  forall(t in trucks) ctUseTruck:istruckused[t]==(1<=sum(o in orders) assign[o][t]);
  
  forall(o in orders) ctOrders:sum(t in trucks) assign[o][t]==1;
  
  forall(t in trucks) ctCapacity:sum(o in orders) assign[o][t]*o.weight<=t.TruckType.Capacity;
}
Alex Fleischer
  • 9,276
  • 2
  • 12
  • 15