运用非连续变形分析法与三维有限元法相结合的方法，以GPS资料作为位移速率和震源机制的约束条件，通过数值模拟研究了青藏高原及其东侧邻区构造地块的运动、变形、相互作用及其与近30年来发生于该区的大地震之间的关系。研究中引入了以应力与摩擦强度的比值定义的断层“失稳危险度”，通过数值模拟计算得到了研究区地壳块体边界断层的失稳危险度分布。结果表明，失稳危险度高的地段与近期该区发生的M_S≥7.0地震所在的位置基本一致，其中龙门山断裂带上包括汶川和芦山大地震的发震断层均为失稳危险度最高值地区。计算得到的应变率强度分布图显示，青藏高原东部边缘整条地带均为应变率强度的陡变带，特别是以龙门山断裂带上的陡变最为明显，其西侧应变率强度为东侧的近4倍，而且，这个带位于宽度相同、走向与龙门山断裂带走向相一致的高应变能密度带中，表明这两次大地震前，作为其发震断层的龙门山断裂带已积累了相当高的应变能，失稳危险度高，处于力学上的不稳定状态。模拟计算得到在上地壳层中，2001年昆仑山口西M_S8.1地震引起汶川、芦山地震发震断层的库仑破裂应力增加约0.016 MPa，相当于龙门山断裂带约两年的应力积累，也就是说，使汶川、芦山地震发震断层的失稳破裂提前了约两年。 此外，关于2008年汶川M_S8.0地震的模拟计算表明，汶川地震的发生也使包括芦山地震发震断层的龙门山断裂带西南段和东昆仑断裂带东南端的库仑破裂应力增大，应变能积累增强，这说明汶川M_S8.0地震的发生对已处于失稳危险度较高状态的2013年芦山地震和2017年九寨沟地震发震断层的提前失稳破裂起到了促进作用。