Les modèles linéaires généralisés sont une généralisation des modèles de régression linéaire, et sont très utilisés dans le domaine du vivant. Le modèle de régression logistique, l'un des modèles de cette classe, très souvent utilisé dans les études biomédicales demeure le modèle de régression le plus approprié quand il s'agit de modéliser une variable discrète de nature binaire. Dans cette thèse, nous nous intéressons au problème de l'inférence statistique dans le modèle de régression logistique, en présence d'individus immunes dans la population d'étude.Dans un premier temps, nous considérons le problème de l'estimation dans le modèle de régression logistique en présence d'individus immunes, qui entre dans le cadre des modèles de régression à excès de zéros (ou zéro-inflatés). Un individu est dit immune s'il n'est pas exposé à l'événement d'intérêt. Le statut d'immunité est inconnu sauf si l'événement d'intérêt a été observé. Nous développons une méthode d'estimation par maximum de vraisemblance en proposant une modélisation conjointe de l'immunité et des risques d'infection. Nous établissons d'abord l'identifiabilité du modèle proposé. Puis, nous montrons l'existence de l'estimateur du maximum de vraisemblance des paramètres de ce modèle. Nous montrons ensuite,la consistance de cet estimateur, et nous établissons sa normalité asymptotique. Enfin, nous étudions, au moyen de simulations, leur comportement sur des échantillons de taille finie.Dans un deuxième temps, nous nous intéressons à la construction de bandes de confiance simultanées pour la probabilité d'infection, dans le modèle de régression logistique avec fraction immune. Nous proposons trois méthodes de constructions de bandes de confiance pour la fonction de régression. La première méthode (méthodede Scheffé) utilise la propriété de normalité asymptotique de l'estimateur du maximum de vraisemblance, et une approximation par une loi du khi deux pour approcher le quantile nécessaire à la construction des bandes. La deuxième méthode utilise également la propriété de normalité asymptotique de l'estimateur du maximum de vraisemblance et est basée sur une égalité classique de (Landau & Sheep 1970). La troisième méthode (méthode bootstrap) repose sur des simulations, pour estimer le quantile approprié de la loi du supremum d'un processus gaussien. Enfin, nous évaluons, au moyen de simulations, leurs propriétés sur des échantillons de taille finie.Enfin, nous appliquons les résultats de modélisation à des données réelles surla dengue. Il s'agit d'une maladie vectorielle tropicale à transmission strictement inter-humaine. Les résultats montrent que les probabilités d'infection estimées à partir de notre approche de modélisation sont plus élevées que celles obtenues à partir d'un modèle de régression logistique standard qui ne tient pas compte d'une possible immunité. En particulier, les estimations fournies par notre approche suggèrent que le sous-poids constitue un facteur de risque majeur de l'infection par la dengue, indépendamment de l'âge.

Cette thèse aborde deux sujets peu traités dans la littérature concernant le théorie des valeurs extrêmes : celui des observations en présence de covariables et celui des mesures de dépendance pour des paires d'observations. Dans la première partie de cette thèse, nous avons considéré le cas où la variable d'intérêt est observée simultanément avec une covariable, pouvant être fixe ou aléatoire. Dans ce contexte, l'indice de queue dépend de la covariable et nous avons proposé des estimateurs de ce paramètre dont nous avons étudié les propriétés asymptotiques. Leurs comportements à distance finie ont été validés par simulations. Puis, dans la deuxième partie, nous nous sommes intéressés aux extrêmes multivariés et plus particulièrement à mesurer la dépendance entre les extrêmes. Dans une situation proche de l'indépendance asymptotique, il est très difficile de mesurer cette dépendance et de nouveaux modèles doivent être introduits. Dans ce contexte, nous avons adapté un outil de géostatistique, le madogramme, et nous avons étudié ses propriétés asymptotiques. Ses performances sur simulations et données réelles ont également été exhibées. Cette thèse offre de nombreuses perspectives, tant sur le plan pratique que théorique dont une liste non exhaustive est présentée en conclusion de la thèse.

Dans ce travail nous démontrons un théorème de simplification uniforme concernant les équations différentielles ordinaires du second ordre singulièrement perturbées au voisinage d'un point dégénéré, appelé point tournant. Il s'agit d'une version analytique d'un résultat formel dû à Hanson et Russell, qui généralise un théorème connu de Sibuya. Pour traiter ce problème, nous utilisons les développements asymptotiques combinés Gevrey introduits par Fruchard et Schäfke. Dans une première partie nous rappelons les définitions et théorèmes principaux de cette récente théorie. Nous établissons trois résultats généraux que nous utilisons ensuite dans la seconde partie de ce manuscrit pour démontrer le théorème principal de réduction analytique annoncé. Enfin nous considérons des équations différentielles ordinaires d'ordre supérieur à deux, singulièrement perturbées à point tournant, et nous démontrons un théorème de réduction analytique.

L'objectif principal du travail réalisé au cours de la thèse consistera à proposer une démarche théorique rigoureuse visant à intégrer les éléments microstructuraux et les gradients de déformation dans une approche micromécanique. La thèse comportera deux volets : - Dans une première partie, on s'attachera à établir un cadre théorique rigoureux permettant d'intégrer les effets du gradient de la déformation et les longueurs caractéristiques de la microstructure sur le comportement effectif des matériaux composites - L'approche par développement asymptotique conduit à la résolution d'une succession problèmes d'élasticité tridimensionnelle posée sur une cellule élémentaire du milieu périodique. La résolution de ces problèmes d'élasticité, et par conséquent la détermination des propriétés effectives du composite, nécessite la mise en œuvre d'une méthode de résolution numérique. Dans cette seconde partie du travail, il s'agira de proposer une méthode de résolution basée sur la transformée de Fourier rapide (Méthode FFT)

Cette thèse est dédiée à la modélisation et à l'analyse mathématique de modèles asymptotiques utilisés en océanographie décrivant la propagation des ondes internes à l'interface entre deux couches de fluides de densités différentes, soumis à la seule force de gravité.L'objectif de cette thèse est de construire et justifier de nouveaux modèles asymptotiques prenant en compte la variation de la topographie. Pour ce faire, on pose plusieurs hypothèses de petitesse sur la profondeur de l'eau et sur les déformations à l'interface et au fond. On s'intéresse plus particulièrement à deux régimes de variations topographiques, celui de moyenne amplitude et celui de lentes variations de grande amplitude.La première partie de cette thèse consiste à justifier rigoureusement et étudier mathématiquement (existence, unicité, stabilité et convergence de la solution) deux classes de modèles asymptotiques. Une classe de modèles couplés et une classe de modèles scalaires. Cette dernière classe est caractérisée par la description de la propagation unidirectionnelle des ondes internes.Dans la deuxième partie on propose un schéma numérique pour résoudre le modèle asymptotique couplé dérivé dans la première partie dans le cadre d'un font plat. Ce modèle existant dans la littérature a été reformulé d'une façon plus appropriée pour la résolution numérique en gardant le même ordre de précision que l'original et en améliorant ses propriétés de dispersion. Enfin nous présentons plusieurs simulations numériques pour valider notre schéma. / This thesis is dedicated to the modeling and the mathematical analysis of asymptotic models used in oceanography describing the propagation of internal waves at the interface between two layers of fluids of different densities, under the only influence of gravity.We aim here at constructing and justifying new asymptotic models taking into account variable topography. To this end, we assume several smallness assumptions on the depth of the water and on the deformations at the interface and at the bottom. We are interested in two topographic regimes, one for variations of medium amplitude and one for slow variations with large amplitude.In the first part of this thesis we rigorously justify and mathematically study (existence, uniqueness, stability and convergence of the solution) two classes of asymptotic models. A class of coupled models and a class of scalar models. The latter class is characterized by the description of the propagation of unidirectional internal waves. In the second part we propose a numerical resolution for the coupled asymptotic model derived in the first part restricted to the flat bottom case. This existing model in the literature has been rewritten under a new formulation more suitable for numericalresolution with the same order of precision as the standard one but with improved frequency dispersion. Finally, we present several numerical simulations to validate our scheme.

On s'intéresse à la géométrie globale et asymptotique de certaines variétés riemanniennes non compactes. Dans une première partie, on étudie la topologie et la géométrie à l'infini des variétés riemanniennes à courbure (de Ricci) positive ayant un rapport asymptotique de courbure fini. On caractérise le cas non effondré via la notion de cône asymptotique et on donne des conditions suffisantes sur le groupe fondamental pour garantir un non effondrement. La seconde partie est dédiée à l'étude des solutions de Type III du flot de Ricci à courbure positive et aux solitons gradients de Ricci expansifs (points fixes de Type III) présentant une décroissance quadratique de la courbure. On montre l'existence et l'unicité des cônes asymptotiques de tels points fixes. On donne également des conditions suffisantes de nature algébrique et géométrique pour garantir une géométrie de révolution de tels solitons. Dans une troisième partie, on caractérise la géométrie des solitons gradients de Ricci stables à courbure positive et à croissance volumique linéaire. Puis, on s'intéresse au non effondrement des variétés riemanniennes de dimension trois à courbure de Ricci positive ayant un rapport asymptotique de courbure fini. / We study the global and asymptotic geometry of non-compact Riemannian manifolds. First, we study the topology and geometry at infinity of Riemannian manifolds with nonnegative (Ricci) curvature and finite asymptotic curvature ratio. We focus on the non-collapsed case with the help of asymptotic cones and we give sufficient conditions on the fundamental group to guarantee non-collapsing. The second part is dedicated to the study of (non-negatively curved) Type III Ricci flow solutions. We mainly analyze the asymptotic geometry of Type III self-similar solutions (expanding gradient Ricci soliton) with finite asymptotic curvature ratio. We prove the existence and uniqueness of their asymptotic cones. We also give algebraic and geometric sufficient conditions to guarantee rotational symmetry of such metrics. In the last part, we characterize the geometry of steady gradient Ricci solitons with nonnegative sectional curvature and linear volume growth. Finally, we study the non-collapsing of three dimensional Riemannian manifold with nonnegative Ricci curvature and finite asymptotic curvature ratio.

Dans cette thèse on s'intéresse à un modèle mathématique décrivant l'adaptation du développement des populations de moustiques face à l'usage intensif des insecticides durant la nuit (moustiquaires imprégnées, répulsifs en spray, répulsifs avec diffuseur électrique, ...).Le modèle proposé dans cette thèse est structuré en âge et dépend du temps/moment où le moustique pique pour prendre son repas. Ceci nous conduità des modèles du type ultra parabolique. Le terme de renouvellement de lapopulation de moustiques est non-local, comme pour tous les problèmes démographiques, mais comporte ici un noyau qui permet à la nouvelle générationd'adapter son temps de piqure (repas). Ceci est dû à la sélection de certainsmoustiques qui piquent plus tôt ou plus tard que les autres moustiques, suite àla pression imposée par l'usage intensif des pesticides à l'intérieur des habitats et en particulier durant la nuit. Les conditions aux bords par rapport au moment de piqure (repas) seront périodiques car selon les espèces, les moustiques prennent toujours leurs repas au même moment de la journée.Les principaux résultats peuvent être classés dans 4 parties.Dans la première partie on présente un modèle structuré en âge décrivant laplasticité du moustique dans un environnement non contrôlé. On montre quele problème est bien posé via la théorie des semi-groupes. Le comportementasymptotique est décrit grâce à l'étude du spectre de l'opérateur A générateurdu C0 semi-groupe. On prouve également l'existence ou la non existence dessolutions stationnaires (sous certaines hypothèses).Dans la deuxième partie on s'intéresse à un problème de contrôle optimalde la population de moustiques. Le contrôle correspond à la proportion demoustiques éliminée et dépend du temps, de l'âge des moustiques et du tempsoù le moustique pique pour se nourrir. On démontre d’abord l’existence desolutions grâce à un argument de point fixe puis on établit des résultats decomparaisons pour notre problème. On établit ensuite l'existence d'un contrôleoptimal puis on dérive le système d'optimalité.Dans la troisième partie on s'intéresse à la question de contrôlabilité exacte locale pour le problème décrivant la capacité des moustiques à adapter leurdynamique face à l'usage intensif des insecticides. On établit une nouvelleinégalité de type Carleman pour le modèle structuré en âge avec diffusionet une condition au bord de renouvellement non-locale et des conditions auxbords périodiques par rapport au temps de piqure des moustiques.Dans la quatrième partie on s'intéresse au comportement en temps longd'un modèle non linéaire décrivant l'adaptation de la population des moustiques à l'usage intensif des insecticides. Quand le contrôle est petit (usage limité des insecticides) alors la population mature de moustiques devient grandeavec le temps et quand le contrôle est grand (usage intensif des insecticides)la population mature de moustiques devient petite avec le temps. Dans le casintermédiaire on obtient un modèle avec retard en temps pour la populationmature de moustiques qui peut être gouvernée par une sur-équation et unesous-équation. Finalement on montre que la sous-équation admet des ondesvoyageuses et la population mature de moustiques sera donc comprise entreces ondes voyageuses et les sur-solutions. / This dissertation is concerned with an age structured problem modelling mosquito plasticity. The main results can be divided into four parts.The first part presents an age structured problem modelling mosquito plasticity in a natural environment. We first investigate the analytical asymptotic solution through studying the spectrum of an operator A which is the infinitesimal generator of a C0-semigroup. Additionally, we get the existence and nonexistence of nonnegative steady solutions under some conditions.In the second part, we study the optimal control of an age structured problem. Firstly, we prove the existence of solutions and the comparison principle for a generalized system. Then, we prove the existence of the optimal control for the best harvesting. Finally, we establish necessary optimality conditions.In the third part, we investigate the local exact controllability of an age structured problem modelling the ability of malaria vectors to shift their biting time to avoid the stressful environmental conditions generated by the use of indoor residual spraying (IRs) and insecticide-treated nets (ITNs). We establish a new Carleman's inequality for our age diffusive model with non local birth processus and periodic biting-time boundary conditions.In the fourth part, we model a mosquito plasticity problem and investigate the large time behavior of matured population under different control strategies. Firstly, we prove that when the control is small, then the matured population will become large for large time and when the control is large, then the matured population will become small for large time. In the intermediate case, we derive a time-delayed model for the matured population which can be governed by a sub-equation and a super-equation. Finally, we prove the existence of traveling fronts for the sub-equation and use it to prove that the matured population will finally be between the positive states of the sub-equation and super-equation.

Cette thèse porte sur l’analyse mathématique de l’interaction d’un fluide non-visqueux avec des structures immergées. Plus précisément, elle est structurée autour de deux axes principaux. L’un d’eux est l’analyse asymptotique du mouvement d’une particule infinitésimale en milieu liquide. L’autre concerne l’interaction entre des vagues et une structure immergée. La première partie de la thèse repose sur l’analyse mathématique d’un système d’équations différentielles ordinaires non-linéaires d’ordre 2 modélisant le mouvement d’un solide infiniment petit dans un fluide incompressible en 2D. Les inconnues du modèle décrivent la position du solide, c’est-à-dire la position du centre de masse et son angle de rotation. Les équations proviennent de la deuxième loi de Newton avec un prototype de force de type Kutta-Joukowski. Plus précisément, nous étudions la dynamique de ce système lorsque l’inertie du solide tend vers 0. Les principaux outils utilisés sont des développements asymptotiques multiéchelles en temps. Pour la dynamique de la position du centre de masse, l’étude met en évidence des analogies avec le mouvement d’une particule chargée dans un champ électromagnétique et la théorie du centre-guide. En l’occurrence, le mouvement du centreguide est donné par une équation de point-vortex. La dynamique de l’angle est quant à elle donnée par une équation de pendule non-linéaire lentement modulée. Des régimes très différents se distinguent selon les données initiales. Pour de petites vitesses angulaires initiales la méthode de Poincaré-Lindstedt fait apparaitre une modulation des oscillations rapides, alors que pour de grandes vitesses angulaires initiales, un movement giratoire bien plus irrégulier est observé. C’est une conséquence particulière et assez spectaculaire de l’enchevêtrement des trajectoires homocliniques. La deuxième partie de la thèse porte sur le problème des vagues dans le cas où le domaine occupé par le fluide est à surface libre et avec un fond plat sur lequel un objet solide se translate horizontalement sous l’effet des forces de pression du fluide. Nous avons étudié deux systèmes asymptotiques qui décrivent le cas d’un fluide parfait incompressible en faible profondeur. Ceux-ci correspondent respectivement aux équations de Saint-Venant et de Boussinesq. Grâce à leur caractère bien-posé en temps long, les modèles traités permettent de prendre en compte certains effets de la mécanique du solide, comme les forces de friction, ainsi que les effets non-hydrostatiques. Notre analyse théorique a été complétée par des études numériques. Nous avons développé un schéma de différences finies d’ordre élevé et nous l’avons adapté à ce problème couplé afin de mettre en évidence les effets d’un solide (dont le mouvement est limité à des translations sur le fond) sur les vagues qui passent au dessus de lui. A la suite de ces travaux, nous avons souligné l’influence des forces de friction sur ce genre de systèmes couplés ainsi que sur le déferlement des vagues. Quant à l’amortissement dû aux effets hydrodynamiques, une vague ressemblance avec le phénomène de l’eau morte est mise en évidence. / This PhD thesis concerns the mathematical analysis of the interaction of an inviscid fluid with immersed structures. More precisely it revolves around two main problems: one of them is the asymptotic analysis of an infinitesimal immersed particle, the other one being the interaction of water waves with a submerged solid object. Concerning the first problem, we studied a system of second order non-linear ODEs, serving as a toy model for the motion of a rigid body immersed in a two-dimensional perfect fluid. The unknowns of the model describe the position of the object, that is the position of its center of mass and the angle of rotation; the equations arise from Newton’s second law with the consideration of a Kutta-Joukowski type lift force. It concerns the detailed analysis of the dynamic of this system when the solid inertia tends to 0. For the evolution of the position of the solid’s center of mass, the study highlights similarities with the motion of a charged particle in an electromagnetic field and the wellknown “guiding center approximation”; it turns out that the motion of the corresponding guiding center is given by a point-vortex equation. As for the angular equation, its evolution is given by a slowly-in-time modulated non-linear pendulum equation. Based on the initial values of the system one can distinguish qualitatively different regimes: for small angular velocities, by the Poincaré-Lindstedt method one observes a modulation in the fast time-scale oscillatory terms, for larger angular velocities however erratic rotational motion is observed, a consequence of Melnikov’s observations on the presence of a homoclinic tangle. About the other problem, the Cauchy problem for the water waves equations is considered in a fluid domain which has a free surface on the upper vertical limit and a flat bottom on which a solid object moves horizontally, its motion determined by the pressure forces exerted by the fluid. Two shallow water asymptotic regimes are detailed, well-posedness results are obtained for both the Saint-Venant and the Boussinesq system coupled with Newton’s equation characterizing the solid motion. Using the particular structure of the coupling terms one is able to go beyond the standard scale for the existence time of solutions to the Boussinesq system with a moving bottom. An extended numerical study has also been carried out for the latter system. A high order finite difference scheme is developed, extending the convergence ratio of previous, staggered grid based models. The discretized solid mechanics are adapted to represent important features of the original model, such as the dissipation due to the friction term. We observed qualitative differences for the transformation of a passing wave over a moving solid object as compared to an immobile one. The movement of the solid not only influences wave attenuation but it affects the shoaling process as well as the wave breaking. The importance of the coefficient of friction is also highlighted, influencing qualitative and quantitative properties of the coupled system. Furthermore, we showed the hydrodynamic damping effects of the waves on the solid motion, reminiscent of the so-called dead water phenomenon.

Cette thèse s'inscrit dans le cadre de la modélisation de systèmes biologiques complexes. En premier lieu (travail en collaboration avec l'équipe de biologistes de Louis Casteilla), nous introduisons un modèle individus-centré pour étudier l'émergence de structures cellulaires de forme lobulaire dans un réseau organisé de fibres. Une étude paramétrique sur les résultats numériques ainsi que des méthodes de traitement d'images sur les données biologiques nous permettent de montrer que l'émergence de structures biologiquement cohérentes peut être reproduite par un modèle basé essentiellement sur des règles mécaniques entre les cellules et le réseau de fibres. L'originalité de ce modèle réside dans la modélisation de structures géométriques complexes (réseaux de fibres) par un ensemble d'unités élémentaires connectées interagissant à l'aide de fonctionnelles simples. Cette nouveauté a donné lieu à d'autres types de travaux en cours présentés comme perspectives directes de ce travail. Les deuxième et troisième parties portent sur la dérivation d'un modèle cinétique puis macroscopique pour étudier la matrice extra-cellulaire dans sa globalité. L'originalité de notre résultat réside dans l'obtention d'un système d'équations fermé décrivant l'évolution de la distribution des fibres individuelles et des liens de fibres. La limite hydrodynamique de l'équation cinétique est obtenue à l'aide de techniques non conventionnelles dû au un manque d'équation de conservation pour le système étudié. Dans le cas d'une densité homogène de fibres, nous prouvons l'existence de solutions au modèle macroscopique, et les simulations numériques montrent une bonne correspondance entre le modèle macroscopique et son homologue microscopique. Finalement (travail en collaboration avec S. Motsch), nous nous intéressons à l'influence d’interactions de type répulsion cellules-cellules dans un modèle de croissance tumorale. Nous montrons que le modèle macroscopique dérivé des équations microscopiques fait apparaître une instabilité et proposons une version modifiée de l'équation macroscopique que nous sommes capables de relier à la dynamique cellulaire. Les simulations numériques montrent la bonne correspondance entre les deux modèles. / In a first part (work in collaboration with the team of biologists of L. Casteilla) we propose an Individual based model for studying the emergence of lobule-like structures of cells in an organized fiber network. A parametric analysis on the numerical results as well as image processing methods on the biological images enable us to show that biologically-relevant structures can be reproduced by a model mostly based on cell-fiber mechanical interactions. The originality of this model mainly relies in the modelling of complex geometrical structures such as fiber networks as sets of connected elementary units interacting through simple functionals. This novelty has been used to build other types of models presented as direct perspective of this work. The second and third parts lie in the derivation of kinetic and macroscopic models for an interconnected fiber network, closely linked to the microscopic one. The originality of this work lies in the obtained closed system of two evolution equations: one for the distribution of individual fibers and one for the fiber links. In the case of homogeneous fiber density, we show existence of stationary solutions to the macroscopic equation, and numerical simulations show the good correspondence between the microscopic and macroscopic models. Finally (work in collaboration with S. Motsch), we are interested in the role of cell-cell interactions in the invasion properties (speed, geometry...) of a growing mass of cells. We show that the macroscopic model derived from the microscopic one features instabilities, and we propose a modified macroscopic model that we are able to link to the particle dynamics. The numerical simulations show the relevance of the macroscopic model to describe the microscopic dynamics at large scale.

L'objectif de cette thèse est de développer des techniques de modélisation et de simulation multi-échelle avancées et efficaces pour étudier les phénomènes d'instabilité dans trois structures d'ingénierie courantes: membrane, film/substrat et structures sandwich, en combinant la technique des coefficients de Fourier lentement variables (TSVFC) et la méthode numérique asymptotique (ANM). À cette fin, basée sur les équations de la plaque de Von Karman, la TSVFC été utilisée pour développer un modèle de Fourier à bidimensionnel (2D) qui a également été implémenté dans ABAQUS via sa sous-routine UEL. Ensuite, un 2D modèle de Fourier est construit pour le film/substrat. En outre, en utilisant leurs caractéristiques de déformation, un 1D modèle de Fourier est développé en utilisant à la fois le TSVFC et le CUF. Par la suite, sur la base d'une cinématique Zig-Zag d'ordre supérieur, un 2D modèle de Fourier est déduit pour une plaque sandwich. Les équations directrices pour les modèles ci-dessus sont discrétisées par la méthode des éléments finis, et les systèmes non linéaires résultants sont résolus par le solveur non linéaire efficace et robuste ANM. Ces modèles sont ensuite adoptés pour étudier les instabilités dans ces structures. Les résultats montrent que les modèles établis peuvent simuler avec précision et efficacité divers phénomènes d'instabilité. En outre, on constate que l'instabilité membranaire est sensible aux conditions aux limites et qu'il existe un paramètre sans dimension presque constant près du point de bifurcation pour différents cas de charge et paramètres géométriques, ce qui peut être utile pour prédire rapidement l'apparition des rides / The main aim of this thesis is to develop advanced and efficient multi-scale modeling and simulation techniques to study instability phenomena in three common engineering structures, i.e., membrane, film/substrate and sandwich structures, by combining the Technique of Slowly Variable Fourier Coefficients (TSVFC) and the Asymptotic Numerical Method (ANM). Towards this end, based on the Von Karman plate equations, the TSVFC has been firstly used to develop a two-dimensional (2D) Fourier double-scale model for membrane, which has also been implemented into ABAQUS via its subroutine UEL. Then a 2D Fourier model is constructed for film/substrate. Further, making use of deformation features of the film/substrate, a 1D Fourier model is developed by using both the TSVFC and the Carrera’s Unified Formulation (CUF). Subsequently, based on high-order kinematics belonging to Zig-Zag theory, a 2D Fourier model is deduced for sandwich plate. The governing equations for the above models are discretized by the Finite Element Method, and the resulting nonlinear systems are solved by the efficient and robust nonlinear solver ANM. These models are then adopted to study instabilities in these structures. Results show that the established models could accurately and efficiently simulate various instability phenomena. Besides, it’s found that the membrane instability is very sensitive to boundary conditions, and there exists a dimensionless parameter that is almost constant near bifurcation point for various loading cases and geometric parameters, which may be helpful for fast predicting the occurrence of wrinkles