Abstract: Shape interpolation, as a bridge communicating static geometries and dynamic shape sequences, is a fundamental operation in digital geometry processing and computer animation. We propose a fast as-isometric-as-possible (AIAP) 3D mesh interpolation approach which casts the shape interpolation problem to finding an AIAP motion trajectory from the start shape to the end shape. This leads to a nonlinear optimization problem with all intermediate shapes as unknowns. The block-coordinate descent method is then employed to iteratively solve the optimization. In each iteration, we need to solve two linear equations whose dimensionality can further be reduced based on a decoupling strategy. Connection maps between orthogonal frames of adjacent edges are further introduced for producing an initial shape sequence in order to address the large-scale deformation problem. A propagation–optimization strategy is then presented to quickly reconstruct the orthogonal frames of all edges from connection maps as well as the orthogonal frame of a specified edge. Refinement of edge quality is available in our method due to the AIAP iterative procedure. In the end, a shape manipulation framework is established for shape sequence transfer and shape sequence editing.