The　superposition-finite　element　method　is　developed　to　analyze　the　mixed-mode　delamination　in　laminated　composites.　Both　a　coarse　global　mesh　and　an　overlaying　fine　local　mesh　are　integrated　into　the　finite　element　analysis　model.　The　whole　design　domain　is　discretized　by　a　uniform　global　mesh,　while　the　high　stress　regions　are　discretized　by　fine　local　meshes.　Local　mesh　is　built　independently　from　the　global　mesh,　which　greatly　simplifies　the　model　generation　procedures.　Strain　energy　release　rate　is　calculated　based　on　the　modified　virtual　crack　closure-integral　method,　which　is　used　to　describe　the　propagation　of　delamination　in　laminated　composites.　Mixed-mode　bending　tests　are　performed　for　unidirectional　and　cross-ply　carbon　fiber　reinforced　laminated　composites　to　characterize　the　interlaminar　fracture　behavior　under　mixed-mode　I/II　loading　conditions.　The　interlaminar　fracture　toughness　is　also　given　for　different　mixed-mode　ratio.　It　is　seen　that　fracture　resistance　for　laminated　composites　exhibit　R-curve　behavior　(increase　as　delamination　propagation).　However,　once　the　delamination　is　sufficiently　long,　the　interfacial　fracture　toughness　changes　slightly　as　the　delamination　extended　and　become　almost　independent　of　the　delamination　length.　The　results　of　finite　element　method　are　in　good　agreement　with　the　experimental　results　and　provides　a　basis　for　establishing　failure　criterion　used　in　damage　tolerance　analysis　of　composite　structures.