Tensile fracture simulations are performed on the (001) interfaces of Fe/Ni and Fe/Pd/Ni by molecular dynamics method, in order to discuss the interfacial strength by means of the eigenvalue and eigenvector of atomic elastic stiffness, Baij=Δσai/Δεj. Three models are considered, Fe/Ni stacked layers (Bimetal), 2 atomic interlayers of Pd between Fe and Ni phases (Pd 2layer) and same thickness Fe/Pd/Ni/Pd blocks (Trimetal). Before tensile simulation, distribution of ηa(1) < 0 atoms (ηa(1) is the 1st eigenvalue of the solution of BaijΔεj=ηaΔεi) and residual atomic stress in each atomic layers are discussed. Then tensile simulations are implemented; the peak strain and stress of Pd2layer are significantly reduced from that of Bimetal. Trimetal shows slightly lower maximum stress at slightly higher strain than Bimetal. Bimetal shows fracture in Ni phase near the interface, while Pd2layer and Trimetal do in Pd phase near the Ni/Pd and Fe/Pd interface, respectively. Void generation process or the initial stage of the fracture is then discussed with the distribution of ηa(1) < 0 atoms and corresponding eigenvector, {Δε1,Δε2,・・・,Δε6} ={Δεxx,Δεyy,・・・,Δγxy}. In all models the void is formed at the intersection between the interface and slip deformation, so that we can’t find clear mode for atomic plane opening but complicated modes for plastic deformation and inhomogeneous heterometal interface. However, large deformation modes or large negative ηa(1) definitely emerge at the opening void edges.