本发明公开了一种板壳起皱失稳极限图的数值模拟求解及绘制方法,涉及金属板材塑性加工技术领域,包括以下步骤:获取金属板壳试件成形的条件要素;根据金属板壳试件成形的条件要素建立金属板壳试件成形的数值起皱失稳模型;根据所建立的数值起皱失稳模型建立金属板壳试件起皱单元簇起皱极限应变求解模型;根据金属板壳试件起皱失稳区域状态类型判定板壳起皱失稳极限图表征形式;根据所建立的起皱单元簇起皱极限应变求解模型和所判定的板壳起皱失稳极限图表征形式建立板壳起皱失稳极限图。本发明准确预测板壳材料起皱失稳,本发明操作简单、节能环保、效率高,计算精准,计算过程直观,结果可靠性高,可广泛应用于工程实践中。