[发明专利]一种原材料下料排样的求解方法

权利要求书

1.一种原材料下料排样的求解方法，其特征在于：包括如下步骤：

步骤S1、确定原材料的长度、每种下料工件的数量、每种下料工件的尺寸大小；

步骤S2、依据原材料的长度、每种下料工件的数量、每种下料工件的尺寸大小建立原始数学模型，对原始数学模型矩阵化得到简化后的数学模型；

步骤S3、建立初始单位矩阵，引入常量参数，利用初始单位矩阵和常量参数对简化后的数学模型进行初等行变化，得到新数学模型；

步骤S4、借助Gurobi优化器对新数学模型进行求解，得到新数学模型的初始可行解，初始可行解包括下料排样方案和决策变量，决策变量指不同排样方案对应的变量，即对应每一种排样方案所使用原材料的数量；

步骤S5、利用CG算法在初始可行解的基础上继续进行求解，得到最优松弛解,最优松弛解是指最优决策变量，若最优松弛解全为非负整数，则将最优决策变量求和的值与步骤S4中决策变量求和的值进行比较，选择求和最小的一方所产生的决策变量作为最终的决策变量，其所对应的下料排样方案为最终的原材料下料排样方案，若最优松弛解不全为非负整数，则进入步骤S6；

步骤S6、借助启发式算法对步骤S5获得的最优松弛解进行取整，得到非负整数解，将最优决策变量求和的值与步骤S4中决策变量求和的值进行比较，选择求和最小的一方所产生的决策变量作为最终的决策变量，其所对应的下料排样方案为最终的原材料下料排样方案；

所述步骤S2中的原始数学模型为：

第一目标函数：

第一约束条件：

，，，，（1）

其中，表示决策变量，具体表示第j种切割模式下所使用的原材料的数量；表示第j种切割模式中第i种待下料的产品的数量；表示第i种待下料的工件所需要的总数量；表示所需第i种待下料的工件的尺寸大小；L代表原材料的长度；表示非负整数；

所述步骤S2中简化后的数学模型具体为：

第二目标函数：

第二约束条件：

，，，

（2）

其中, 表示一个1行n列的矩阵；解向量矩阵表示是一个n行1列的矩阵，该矩阵的元素是决策变量，它的值被约束为大于或等于零且为整数；A表示排样方案的矩阵，是一个m行n列的矩阵，，当求解完成所有的切割模式后，由产生的切割模式组成矩阵A中所有的列，每一列代表一种切割模式，在迭代开始前，它是一个m行m列的单位矩阵；表示一个m行1列的列向量需求矩阵，该矩阵的元素值表示不同尺寸工件的需求数量，m表示所需求工件种类的数量；表示一个m行1列的矩阵，该矩阵是一个列向量变量的列矩阵，其中的元素值为非负整数，代表当前需要被求解的一种切割模式，用作为的索引，表示所有可能的切割模式集，表示当前迭代需要被求解的一种切割模式，其中元素表示在切割模式中，第i种工件的数量，；表示一个1行m列的行向量矩阵，该矩阵的元素代表m个不同工件的尺寸大小；表示在切割模式中，不同工件的尺寸乘以各自对应不同工件数量的乘积之和要小于等于原材料的长度；中的0 表示一个元素全为0的m行1列的列向量矩阵,其中，矩阵中所有的变量值被约束为大于或等于零，且必须小于或等于各类工件对应的需求数量，其值；

所述步骤S3中的新数学模型具体为：

假设决策变量的数量为m个，不同工件种类的数量为m个，每迭代一次，生成一个新列，也代表生成一种新的切割模式，当迭代求解生成第k种切割模式时，对应的第三目标函数和第三约束条件可写成如下：

第三目标函数：

第三约束条件：，

，

，

，and，，

，，

，（3）

在（3）式中，表示在当前迭代中，将切割模式对应的决策变量固定为一个常量参数，该常量参数等于前一次迭代求解得到的第k列对应决策变量的值；表示求解切割模式时，通过将列向量矩阵的第k个元素强制为0获得；同理，表示求解切割模式时，通过将矩阵的第k行的元素强制为0元素获得；表示列向量变量矩阵参与矩阵初等变化之后的矩阵表达式；表示决策变量矩阵中除第k列对应的决策变量为固定值外，其余决策变量全部被约束为大于或等于零；表示未知列向量变量矩阵参与矩阵A进行初等行变化之后，矩阵A中第k行第k列元素的表达式，其值等于1，其中是一个1行m列矩阵，可通过将矩阵的第k行构建的行向量矩阵；目标函数中是一个元素全为1的1行m列矩阵。