本发明公开一种基于启发式算法和遗传算法相结合的多维度约束的智能装载方法，面向楔形货格装载货物的环节，先采用启发式算法产生高质量初始解，而后采用遗传算法对初始解进行迭代寻优，将启发式算法的简洁高效和遗传算法的优化性能充分结合，从而大幅度提升货格的空间使用效率。本发明同时针对具体问题，给出了实用的程序开发方案，提高了装载运输的效率。
1.一种多维度约束下的货格智能装载方法，针对矩形货格建立空间直角坐标系统，由背侧、左侧和底板三面相交形成的顶点为坐标原点，竖向为Z轴，背侧的宽度方向为X轴，左侧的宽度方向为Y轴；其特征是：所述货格在竖向为多层结构，每一层呈横排设置有多个单元格，处在两端的端部单元格由背板和外侧板形成楔形端头，楔形端头的俯视形状划分为中部矩形区域以及与中部矩形区域邻接的直角三角区域；每一层中两只端部单元格之间的各中部单元格均为矩形格，各矩形格及中部矩形区域均具有所述矩形货格的特征，所述矩形货格的特征是指前侧为敞口，背侧、左侧、右侧、顶板及底板均为封闭；所述装载方法是按如下步骤进行：步骤1、获取物品信息和货格信息
所述物品信息包括：物品总数、物品种类、各类物品总数，各物品编号、各类物品的重量、各类物品的体积、各类物品的耐压值、各类物品的取用频次；所述物品的耐压值是指物品能够承受的最大压力值，所述物品的取用频次区分为低频次和高频次；
所述货格信息包括：货格中单元格总数，各单元格的承重量、各单元格的编号，以及各单元格的形状；
形状和尺寸约束规则：装载在货格中的物品完全被包容在货格内；
重量分布均匀约束规则：单元格间重量差异不超过重量差异设定值，保持货格稳定；
承重约束规则：各单元格内的物品重量不超过单元格承载重量设定值；
承压约束规则：单元格中上部物品对下部物品的压力不超过下部物品的承压限值；
取用频次约束规则：取用频次更高的物品在货格内更靠敞口侧放置，以利取用；
步骤3、根据设定的各约束规则，利用已获取的物品信息和货格信息，通过启发式算法和改进遗传算法按如下步骤得到最优装载方案：步骤3.1、通过计算获得物品的优先级
P＝α1×P1+α2×P2 (1)，P1表示物品A的体积比值，P2表示物品A的耐压比值，并有：P1＝VA/V，VA为物品A的体积，V为待装载物品中体积最大的物品的体积；
P2＝PA/P，PA为物品A的耐压值，P为待装载物品中耐压值最大的物品的耐压值；
利用式(1)分别计算获得每个物品的优先级，按照优先级从大到小的顺序对所有物品进行排列，获得待装载物品列表；
在所述待装载物品列表中，依次选取体积小且取用频次为低频次的物品安置于端部单元格中的直角三角区域，保留端部单元格中的中部矩形区域为可用空间，将已得到安置的物品从所述待装载物品列表中删除，获得待装载物品更新列表；
针对已放置有物品，但尚有剩余空间的矩形格，按照矩形货格的特征对剩余空间进行划分，形成体积更小的各独立可用空间；
将所有独立可用空间按体积从大到小排列，得到货格的可用空间集；
所述按照矩形货格的特征对剩余空间进行划分是指：物品正上方的剩余空间为一个独立可用空间；其它位置上的剩余空间划分为多个矩形分区，相邻两个矩形分区合并形成合并区，并取其中合并体积最大的合并区为一个独立可用空间，未合并的各矩形分区分别为各自的独立可用空间；
从待装载物品更新列表中的第一个物品开始，按如下方式依次设计每个物品的装载方案：步骤3.4a、针对当前物品，从可用空间集中的第一个可用空间开始检测：若当前可用空间能够容纳当前物品，且在容纳当前物品后，当前可用空间所属的单元格满足承重约束规则，则进行步骤3.4b；否则，按步骤3.4a相同的方式检测可用空间集中的下一个可用空间；
步骤3.4b、根据取用频次约束规则生成当前物品放置方案：根据取用情况将物品的取用频次定义为高频次或低频次两种情况；
若物品取用频次为高频次，以矩形货格直角坐标系中Y轴最外侧顶点为外侧参考点进行放置，若所述外侧参考点已经放置物品，则优先考虑由外侧参考点沿X轴正方向移动的空间区域，然后考虑由外侧参考点沿Y轴负方向移动的空间区域，最后考虑由外侧参考点沿Z轴正方向移动的空间区域；
若物品取用频次为低频次，以矩形货格直角坐标系原点为内侧参考点进行放置，若内侧参考点已经放置物品，则优先考虑由内侧参考点沿X轴正方向移动的空间区域，然后考虑由内侧参考点沿Y轴正方向移动的空间区域，最后考虑由内侧参考点沿Z轴正方向移动的空间区域；
步骤3.4c、采集获得物品位置信息，所述物品位置信息包括：物品所在单元格的编号；
步骤3.4d、将已经获得物品位置信息的物品从待装载物品更新列表中删除，生成新的待装载物品列表，以所述新的待装载物品列表替代所述待装载物品更新列表；计算货格的当前剩余空间集，以所述当前剩余空间集代替所述可用空间集；计算各单元格内所有物品的总重量；
步骤3.4e、重复步骤3.4a到步骤3.4d，直至满足下述条件之一停止启发式算法的迭代，获得初始解：条件一、待装载物品更新列表为空；
步骤3.5、采用遗传算法对初始解进行迭代寻优步骤3.5a、获取遗传算法参数，包括：种群规模N，变异概率Pm，交叉概率Pc，迭代次数M；
步骤3.5b、调整遗传算法参数，获得初始种群：调整体积权重α1、耐压值权重α2，运行步骤3.1至步骤3.4，获得N组不同的启发式方案作为遗传算法的初始种群；
将每个物品的物品位置信息作为染色体信息，根据交叉和变异概率选择染色体，然后进行交叉和变异操作；
当满足下列条件之一时，终止迭代寻优的过程，获得最优装载方案：条件一、迭代次数达到给定的迭代次数M次；
2.根据权利要求1所述的多维度约束下的货格智能装载方法，其特征是：针对所获得的装载方案，通过计算获得各评价指标，所述评价指标包括：各单元格空间利用率Si、货格空间利用率S、各单元格承重率Ri、货格承重率R，以及货格重量分配率D；
TVi为单元格i内放置的物品的总体积，Vi为单元格i的容积；
CTV为货格内放置物品的总体积，CV为货格总容积；
TWi为单元格i内放置的物品的总重量，Wi为单元格i的最大承重；
TWmax为货格中的装载重量最大的单元格的装载重量，TWmin为货格中装载重量最小的单元格的装载重量，TWavg为货格中所有单元格的平均装载重量；
TWi为单元格i内放置的物品的总重量，i＝1,2,…,CN；CN为货格中单元格的总数；
根据所述评价指标判断所获得的装载方案是否满足设定要求。
3.根据权利要求2所述的多维度约束下的货格智能装载方法，其特征是：调整所述体积权重α1和耐压值权重α2得到不同的装载方案，针对所获得的装载方案利用计算获得的各评价指标决策实施的装载方案。
多维度约束下的楔形货格智能装载方法
技术领域
[0001] 本发明属于物流装载技术领域，具体地说是一种基于启发式算法和遗传算法相结合的多维度约束的智能装载算法。
背景技术
[0002] 货物的配载是物流配送过程中的主要环节。合理的配载配装可以给物流配送中心带来更大的经济效益及社会效益。然而到目前为止，国内大部分配送中心在实际操作中仍然凭借经验进行配载和配装，这种原始的配载和配装既不能充分利用空间，也增加了运输成本。因此，充分利用装载空间，降低运输成本，将是亟需解决的一大难题。
发明内容
[0003] 本发明是为了克服现有技术存在的不足之处，提供一种基于启发式算法和遗传算法相结合的多维度约束的智能装载算法，以期能最大的利用楔形容器空间，提高装载空间的利用率，从而降低运输成本。
[0004] 本发明为解决技术问题采用如下技术方案：
[0005] 本发明多维度约束下的货格智能装载方法，其特点是：所述货格在竖向为多层结构，每一层呈横排设置有多个单元格，处在两端的端部单元格由背板和外侧板形成楔形端头，楔形端头的俯视形状划分为中部矩形区域以及与中部矩形区域邻接的直角三角区域；每一层中两只端部单元格之间的各中部单元格均为矩形格，各矩形格及中部矩形区域均具有矩形货格的特征，所述矩形货格的特征是指前侧为敞口，背侧、左侧、右侧、顶板及底板均为封闭；针对所述矩形货格建立空间直角坐标系统，由背侧、左侧和底板三面相交形成的顶点为坐标原点，竖向为Z轴，背侧的宽度方向为X轴，左侧的宽度方向为Y轴；所述装载方法是按如下步骤进行：
[0006] 步骤1、获取物品信息和货格信息
[0007] 所述物品信息包括：物品总数、物品种类、各类物品总数，各物品编号、各类物品的重量、各类物品的体积、各类物品的耐压值、各类物品的取用频次；所述物品的耐压值是指物品能够承受的最大压力值，所述物品的取用频次区分为低频次和高频次；
[0008] 所述货格信息包括：货格中单元格总数，各单元格的承重量、各单元格的编号，以及各单元格的形状；
[0009] 步骤2、设定如下各约束规则：
[0010] 形状和尺寸约束规则：装载在货格中的物品完全被包容在货格内；
[0011] 重量分布均匀约束规则：单元格间重量差异不超过重量差异设定值，保持货格稳定；
[0012] 承重约束规则：各单元格内的物品重量不超过单元格承载重量设定值；
[0013] 承压约束规则：单元格中上部物品对下部物品的压力不超过下部物品的承压限值；
[0014] 取用频次约束规则：取用频次更高的物品在货格内更靠敞口侧放置，以利取用；
[0015] 步骤3、根据设定的各约束规则，利用已获取的物品信息和货格信息，通过启发式算法和改进遗传算法按如下步骤得到最优装载方案：
[0016] 步骤3.1、通过计算获得物品的优先级
[0017] 利用式(1)计算获得物品A的优先级P：
[0018] P＝α1×P1+α2×P2 (1)，
[0019] P1表示物品A的体积比值，P2表示物品A的耐压比值，并有：
[0020] P1＝VA/V，VA为物品A的体积，V为待装载物品中体积最大的物品的体积；
[0021] P2＝PA/P，PA为物品A的耐压值，P为待装载物品中耐压值最大的物品的耐压值；
[0022] α1为体积权重，α2为耐压值权重；
[0023] 利用式(1)分别计算获得每个物品的优先级，按照优先级从大到小的顺序对所有物品进行排列，获得待装载物品列表；
[0024] 步骤3.2、填充端部单元格中的直角三角区域
[0025] 在所述待装载物品列表中，依次选取体积小且取用频次为低频次的物品安置于端部单元格中的直角三角区域，保留端部单元格中的中部矩形区域为可用空间，将已得到安置的物品从所述待装载物品列表中删除，获得待装载物品更新列表；
[0026] 步骤3.3、按如下方式获得货格的可用空间集
[0027] 所有空置的矩形格均为独立可用空间；
[0028] 针对已放置有物品，但尚有剩余空间的矩形格，按照矩形货格的特征对剩余空间进行划分，形成体积更小的各独立可用空间；
[0029] 将所有独立可用空间按体积从大到小排列，得到货格的可用空间集；
[0030] 所述按照矩形货格的特征对剩余空间进行划分是指：物品正上方的剩余空间为一个独立可用空间；其它位置上的剩余空间划分为多个矩形分区，相邻两个矩形分区合并形成合并区，并取其中合并体积最大的合并区为一个独立可用空间，未合并的各矩形分区分别为各自的独立可用空间；
[0031] 步骤3.4、放置物品
[0032] 从待装载物品更新列表中的第一个物品开始，按如下方式依次设计每个物品的装载方案：
[0033] 步骤3.4a、针对当前物品，从可用空间集中的第一个可用空间开始检测：
[0034] 若当前可用空间能够容纳当前物品，且在容纳当前物品后，当前可用空间所属的单元格满足承重约束规则，则进行步骤3.4b；否则，按步骤3.4a相同的方式检测可用空间集中的下一个可用空间；
[0035] 步骤3.4b、根据取用频次约束规则生成当前物品放置方案：
[0036] 根据取用情况将物品的取用频次定义为高频次或低频次两种情况；
[0037] 若物品取用频次为高频次，以矩形货格直角坐标系中Y轴最外侧顶点为外侧参考点进行放置，若所述外侧参考点已经放置物品，则优先考虑由外侧参考点沿X轴正方向移动的空间区域，然后考虑由外侧参考点沿Y轴负方向移动的空间区域，最后考虑由外侧参考点沿Z轴正方向移动的空间区域；
[0038] 若物品取用频次为低频次，以矩形货格直角坐标系原点为内侧参考点进行放置，若内侧参考点已经放置物品，则优先考虑由内侧参考点沿X轴正方向移动的空间区域，然后考虑由内侧参考点沿Y轴正方向移动的空间区域，最后考虑由内侧参考点沿Z轴正方向移动的空间区域。
[0039] 步骤3.4c、采集获得物品位置信息，所述物品位置信息包括：
[0040] 物品所在单元格的编号；
[0041] 物品的中心点在对应单元格中的坐标；
[0042] 物品映射在X、Y、Z轴上的区间；
[0043] 同一单元格中处在物品正上方及正下方的物品编号；
[0044] 步骤3.4d、将已经获得物品位置信息的物品从待装载物品更新列表中删除，生成新的待装载物品列表，以所述新的待装载物品列表替代所述待装载物品更新列表；计算货格的当前剩余空间集，以所述当前剩余空间集代替所述可用空间集；计算各单元格内所有物品的总重量；
[0045] 步骤3.4e、重复步骤3.4a到步骤3.4d，直至满足下述条件之一停止启发式算法的迭代，获得初始解：
[0046] 条件一、待装载物品更新列表为空；
[0047] 条件二、可用空间集为空；
[0048] 步骤3.5、采用遗传算法对初始解进行迭代寻优
[0049] 步骤3.5a、获取遗传算法参数，包括：种群规模N，变异概率Pm，交叉概率Pc，迭代次数M；
[0050] 步骤3.5b、调整遗传算法参数，获得初始种群：
[0051] 调整体积权重α1、耐压值权重α2，运行步骤3.1至步骤3.4，获得N组不同的启发式方案作为遗传算法的初始种群；
[0052] 步骤3.5c、交叉变异操作
[0053] 将每个物品的物品位置信息作为染色体信息，根据交叉和变异概率选择染色体，然后进行交叉和变异操作；
[0054] 步骤3.5d、迭代
[0055] 更新种群信息以及变量信息，进入下一次迭代；
[0056] 步骤3.5e、终止条件
[0057] 当满足下列条件之一时，终止迭代寻优的过程，获得最优装载方案：
[0058] 条件一、迭代次数达到给定的迭代次数M次；
[0059] 条件二、获得的方案满足要求。
[0060] 本发明多维度约束下的货格智能装载方法的特点也在于：针对所获得的装载方案，通过计算获得各评价指标，所述评价指标包括：各单元格空间利用率Si、货格空间利用率S、各单元格承重率Ri、货格承重率R，以及货格重量分配率D；
[0061]
[0062]
[0063] 其中：
[0064] TVi为单元格i内放置的物品的总体积，Vi为单元格i的容积；
[0065] CTV为货格内放置物品的总体积，CV为货格总容积；
[0066] TWi为单元格i内放置的物品的总重量，Wi为单元格i的最大承重；
[0067] CTW为货格内物品的总重量，CW为货格总承重；
[0068] TWmax为货格中的装载重量最大的单元格的装载重量，TWmin为货格中装载重量最小的单元格的装载重量，TWavg为货格中所有单元格的平均装载重量；
[0069] TWi为单元格i内放置的物品的总重量，i＝1,2,…,CN；CN为货格中单元格的总数；
[0070] 根据所述评价指标判断所获得的装载方案是否满足设定要求。
[0071] 本发明多维度约束下的货格智能装载方法的特点也在于：调整所述体积权重α1和耐压值权重α2得到不同的装载方案，针对所获得的装载方案利用计算获得的各评价指标决策实施的装载方案。
[0072] 与已有技术相比，本发明有益效果体现在：
[0073] 1、本发明设定一系列的决策变量用来精准描述各维度约束，在算法输入阶段确保所有约束能够完全体现；通过将约束规则和算法运行高效融合，保证算法输出方案的可行性，并对非常规情况安排特定处理；最终获得方案及其各项性能指标输出，使装载方案能够直观地展示，同时全方位地显示装载方案的性能。
[0074] 2、本发明考虑了多维度约束：货格的形状、尺寸限制；货物摆放重量应分布均匀、总重量不可超出货格的最大承重；考虑物品的耐压性、取用频次。在算法设计中，将约束规则和算法运行高效融合，保证算法输出方案的可行性，并对非常规情况安排了特殊处理规则。
[0075] 3.本发明通过启发式算法和遗传算法相结合的方式，先由启发式算法产生高质量初始解，再由遗传算法对初始解进行迭代，寻求最优解，实现最大利用楔形容器空间，提高装载空间的利用率，从而降低运输成本。
附图说明
[0076] 图1为本发明多维度约束下的智能装载方法的总体流程；
[0077] 图2a为本发明中货格中单层结构示意图；
[0078] 图2b为本发明中货格中楔形端头俯视图；
[0079] 图3为本发明中物品周围的各子空间示意图。
具体实施方式
[0080] 本实施例中多维度约束下的货格智能装载方法，其货格的结构形式是：货格在竖向为多层结构，每一层呈横排设置有多个单元格，处在两端的端部单元格由背板和外侧板形成楔形端头，楔形端头的俯视形状划分为中部矩形区域以及与中部矩形区域邻接的直角三角区域；每一层中两只端部单元格之间的各中部单元格均为矩形格，各矩形格及中部矩形区域均具有矩形货格的特征，所述矩形货格的特征是指前侧为敞口，背侧、左侧、右侧、顶板及底板均为封闭；针对所述矩形货格建立空间直角坐标系统，由背侧、左侧和底板三面相交形成的顶点为坐标原点，竖向为Z轴，背侧的宽度方向为X轴，左侧的宽度方向为Y轴。
[0081] 参见图1，本实施例中多维度约束下的货格智能装载方法按如下步骤进行：
[0082] 步骤1、获取物品信息和货格信息
[0083] 所述物品信息包括：物品总数、物品种类、各类物品总数，各物品编号、各类物品的重量、各类物品的体积、各类物品的耐压值、各类物品的取用频次；所述物品的耐压值是指物品能够承受的最大压力值，所述物品的取用频次区分为低频次和高频次；
[0084] 所述货格信息包括：货格中单元格总数，各单元格的承重量、各单元格的编号，以及各单元格的形状。
[0085] 步骤2、设定如下各约束规则：
[0086] 形状和尺寸约束规则：装载在货格中的物品完全被包容在货格内；
[0087] 重量分布均匀约束规则：单元格间重量差异不超过重量差异设定值，保持货格稳定；
[0088] 承重约束规则：各单元格内的物品重量不超过单元格承载重量设定值；
[0089] 承压约束规则：单元格中上部物品对下部物品的压力不超过下部物品的承压限值；
[0090] 取用频次约束规则：取用频次更高的物品在货格内更靠敞口侧放置，以利取用。
[0091] 步骤3、根据设定的各约束规则，利用已获取的物品信息和货格信息，通过启发式算法和改进遗传算法按如下步骤得到最优装载方案：
[0092] 步骤3.1、通过计算获得物品的优先级
[0093] 利用式(1)计算获得物品A的优先级P：
[0094] P＝α1×P1+α2×P2 (1)，
[0095] P1表示物品A的体积比值，P2表示物品A的耐压比值，并有：
[0096] P1＝VA/V，VA为物品A的体积，V为待装载物品中体积最大的物品的体积；
[0097] P2＝PA/P，PA为物品A的耐压值，P为待装载物品中耐压值最大的物品的耐压值；
[0098] α1为体积权重，α2为耐压值权重；针对不同的产品，对于体积权重和耐压值权重可以有不同的考虑，通过调整体积权重α1和耐压值权重α2予实现。
[0099] 利用式(1)分别计算获得每个物品的优先级，按照优先级从大到小的顺序对所有物品进行排列，获得待装载物品列表。
[0100] 步骤3.2、填充端部单元格中的直角三角区域
[0101] 在所述待装载物品列表中，依次选取体积小且取用频次为低频次的物品安置于端部单元格中的直角三角区域，保留端部单元格中的中部矩形区域为可用空间，将已得到安置的物品从所述待装载物品列表中删除，获得待装载物品更新列表。
[0102] 步骤3.3、按如下方式获得货格的可用空间集
[0103] 所有空置的矩形格均为独立可用空间；
[0104] 针对已放置有物品，但尚有剩余空间的矩形格，按照矩形货格的特征对剩余空间进行划分，形成体积更小的各独立可用空间；除了用于填充货格中直角三角区域的物品外，企业也可以根据自身需求放置有固定摆放位置或者具有特殊要求的物品，比如：具有排它性的物品可以单独放置。
[0105] 将所有独立可用空间按体积从大到小排列，得到货格的可用空间集；
[0106] 所述按照矩形货格的特征对剩余空间进行划分是指：物品正上方的剩余空间为一个独立可用空间；其它位置上的剩余空间划分为多个矩形分区，相邻两个矩形分区合并形成合并区，并取其中合并体积最大的合并区为一个独立可用空间，未合并的各矩形分区分别为各自的独立可用空间。
[0107] 步骤3.4、放置物品
[0108] 从待装载物品更新列表中的第一个物品开始，按如下方式依次设计每个物品的装载方案：
[0109] 步骤3.4a、针对当前物品，从可用空间集中的第一个可用空间开始检测：
[0110] 若当前可用空间能够容纳当前物品，且在容纳当前物品后，当前可用空间所属的单元格满足承重约束规则，则进行步骤b；否则，按步骤a相同的方式检测可用空间集中的下一个可用空间。
[0111] 步骤3.4b、根据取用频次约束规则生成当前物品放置方案：
[0112] 根据取用情况将物品的取用频次定义为高频次或低频次两种情况；
[0113] 若物品取用频次为高频次，以矩形货格直角坐标系中Y轴最外侧顶点为外侧参考点进行放置，若所述外侧参考点已经放置物品，则优先考虑由外侧参考点沿X轴正方向移动的空间区域，然后考虑由外侧参考点沿Y轴负方向移动的空间区域，最后考虑由外侧参考点沿Z轴正方向移动的空间区域；
[0114] 若物品取用频次为低频次，以矩形货格直角坐标系原点为内侧参考点进行放置，若内侧参考点已经放置物品，则优先考虑由内侧参考点沿X轴正方向移动的空间区域，然后考虑由内侧参考点沿Y轴正方向移动的空间区域，最后考虑由内侧参考点沿Z轴正方向移动的空间区域。取用频次的高或低可以由企业根据产品的情况自行设定，比如：对于汽修企业，其检修工具的取用频次为高；但是对于运输企业，检修工具只有在出现故障时才需要取用，因此取用频次为低。
[0115] 步骤3.4c、采集获得物品位置信息，所述物品位置信息包括：
[0116] 物品所在单元格的编号；
[0117] 物品的中心点在对应单元格中的坐标；
[0118] 物品映射在X、Y、Z轴上的区间；
[0119] 同一单元格中处在物品正上方及正下方的物品编号。
[0120] 步骤3.4d、将已经获得物品位置信息的物品从待装载物品更新列表中删除，生成新的待装载物品列表，以所述新的待装载物品列表替代所述待装载物品更新列表；计算货格的当前剩余空间集，以所述当前剩余空间集代替所述可用空间集；计算各单元格内所有物品的总重量。
[0121] 步骤3.4e、重复步骤3.4a到步骤3.4d，直至满足下述条件之一停止启发式算法的迭代，获得初始解：
[0122] 条件一、待装载物品更新列表为空；
[0123] 条件二、可用空间集为空。
[0124] 步骤3.5、采用遗传算法对初始解进行迭代寻优
[0125] 步骤3.5a、获取遗传算法参数，包括：种群规模N，变异概率Pm，交叉概率Pc，迭代次数M。
[0126] 步骤3.5b、调整遗传算法参数，获得初始种群：
[0127] 调整体积权重α1、耐压值权重α2，运行步骤3.1至步骤3.4，获得N组不同的启发式方案作为遗传算法的初始种群。
[0128] 步骤3.5c、交叉变异操作
[0129] 将每个物品的物品位置信息作为染色体信息，根据交叉和变异概率选择染色体，然后进行交叉和变异操作。
[0130] 步骤3.5d、迭代
[0131] 更新种群信息以及变量信息，进入下一次迭代。
[0132] 步骤3.5e、终止条件
[0133] 当满足下列条件之一时，终止迭代寻优的过程，获得最优装载方案：
[0134] 条件一、迭代次数达到给定的迭代次数M次；
[0135] 条件二、获得的方案满足要求。
[0136] 具体实施中，针对所获得的装载方案，通过计算获得各评价指标，所述评价指标包括：各单元格空间利用率Si、货格空间利用率S、各单元格承重率Ri、货格承重率R，以及货格重量分配率D；
[0137]
[0138]
[0139] 其中：
[0140] TVi为单元格i内放置的物品的总体积，Vi为单元格i的容积；
[0141] CTV为货格内放置物品的总体积，CV为货格总容积；
[0142] TWi为单元格i内放置的物品的总重量，Wi为单元格i的最大承重；
[0143] CTW为货格内物品的总重量，CW为货格总承重；
[0144] TWmax为货格中的装载重量最大的单元格的装载重量，TWmin为货格中装载重量最小的单元格的装载重量，TWavg为货格中所有单元格的平均装载重量；
[0145] TWi为单元格i内放置的物品的总重量，i＝1,2,…,CN；CN为货格中单元格的总数；
[0146] 各评价指标从不同角度反映了方案的性能，其中各单元格空间利用率Si和货格空间利用率S反映了货格容积的利用率，其指标值越接近100％，则表明货格的空间利用越好；各单元格承重率Ri和货格承重率R反映了货格承重能力的利用率，其指标值越接近100％，表明货格承重的利用越好；货格重量分配率D反映了货格之间重量分布的均匀程度，其指标值越接近0，表明越均匀，货格总体的稳定性越好。
[0147] 根据所述评价指标判断所获得的装载方案是否满足设定要求。
[0148] 具体实施中，调整体积权重α1和耐压值权重α2得到不同的装载方案，针对所获得的装载方案利用计算获得的各评价指标决策实施的装载方案。
[0149] 图2a所示为货格单层示意图，其中间可有任意数量矩形格；处在两端的端部单元格21为楔形端头，两只端部单元格21之间的各中部单元格22均为矩形格。
[0150] 图2b所示为楔形端头俯视图，楔形端头的俯视形状划分为中部矩形区域23以及与中部矩形区域邻接的直角三角区域24。
[0151] 图3所示为本发明中位于已放置物品5的周围的各子空间示意图，包括正上方独立可用空间1、正前方独立可用空间2、右侧独立可用空间3，以及对角位置独立可用空间4。
