Multi-pulse　homoclinic　orbits　and　chaotic　dynamics　of　a　supercritical　composite　panel　with　free　layer　damping　treatment　in　subsonic　flow　are　investigated　considering　one　to　two　internal　resonance.　Inviscid　potential　flow　theory　is　employed　to　exhibit　the　aerodynamic　pressure　and　Kelvin＇s　model　is　used　to　describe　the　viscoelastic　property　of　the　free　damping　layer.　By　Hamilton＇s　principle,　the　governing　equation　of　the　composite　panel　in　the　subcritical　regime　is　derived.　In　the　supercritical　regime,　the　buckling　configuration　is　solved　analytically　and　the　PDE　is　obtained　by　introducing　a　displacement　transformation　for　nontrivial　equilibrium　configuration.　Then　the　governing　equation　in　the　first　supercritical　region　is　transformed　into　a　discretized　nonlinear　gyroscopic　system　via　assumed　modes　and　then　Galerkin＇s　method.　The　method　of　multiple　scales　and　canonical　transformation　are　applied　to　reduce　the　equations　of　motions　to　the　near-integrable　Hamiltonian　standard　forms.　The　Energy-Phase　method　is　employed　to　demonstrate　the　existence　of　chaotic　dynamics　by　identifying　the　existence　of　multipulse　jumping　orbits　in　the　perturbed　phase　space.　The　global　solutions　are　finally　interpreted　in　terms　of　the　physical　motion　of　the　gyroscopic　continua　and　the　dynamical　mechanism　of　chaotic　pattern　conversion　between　the　forward　traveling　wave　motion　and　the　complex　bidirectional　traveling　wave　motion　are　discussed.　(C)　2016　Elsevier　Ltd.　All　rights　reserved.