Chemically stable metal-organic frameworks (MOFs) featuring interconnected hierarchical pores have proven to be promising for a remarkable variety of applications. Nevertheless, the framework's susceptibility to capillary-force-induced pore collapse, especially during water evacuation, has often limited practical applications. Methodologies capable of predicting the relative magnitudes of these forces as functions of the pore size, chemical composition of the pore walls, and fluid loading would be valuable for resolution of the pore collapse problem. Here, we report that a molecular simulation approach centered on evacuation-induced nanocavitation within fluids occupying MOF pores can yield the desired physical-force information. The computations can spatially pinpoint evacuation elements responsible for collapse and the chemical basis for mitigation of the collapse of modified pores. Experimental isotherms and difference-electron density measurements of the MOF NU-1000 and four chemical variants validate the computational approach and corroborate predictions regarding relative stability, anomalous sequence of pore-filling, and chemical basis for mitigation of destructive forces.