为充分利用建筑结构在地震作用下的局部非线性特征，采用数值子结构方法将原本复杂的非线性结构分析转化为主结构等效线弹性分析和屈服构件隔离子结构非线性分析，可有效兼顾计算效率和模拟精度。然而，该方法中主结构串行计算，其应用于超大规模结构的非线性分析中仍有一定局限性。为此，将数值子结构方法与传统的区域分解法相结合，提出了数值子结构并行计算方法。将整体结构划分为若干弹性子域，各子域同时进行等效线弹性分析，子域内屈服构件由隔离子结构进行非线性分析。在子域回代求解时，设计了降阶的牛顿迭代算法，以子域内屈服构件单元的位移场作为基本未知量，可进一步降低子域迭代分析的矩阵运算规模。以一15层3跨平面钢框架结构的地震弹塑性全过程分析为例进行数值模拟，结果表明，所提出的方法准确、可靠，当应用于具有局部材料非线性的大规模结构分析时，该方法可以大幅降低矩阵运算规模，显著提高计算效率，将大问题简化为若干小的子问题，可实现更大、更复杂结构的非线性分析。