Neste trabalho estudamos a estabilidade assintótica de Lagrange para duas classes de sistemas, sob as hipóteses de detetabilidade fraca e de limitação do custo medio a longo prazo. Para sistemas lineares com saltos markovianos com rudo aditivo, a equivalência entre estabilidade e as condições mencionadas sera provada. Para sistemas dinâmicos generalizados, provaremos a estabilidade sob uma condição adicional / In this work we study Lagrange asymptotic stability for two classes of systems, under conditions of weak detectability and boundedness of the long run average cost. For Markov jump linear systems with additive noise, the equivalence between stability and the aforementioned conditions is proved. For generalized dynamical systems, we prove stability under an additional condition

Dados longitudinais na forma de contagens e na forma binária são muito comuns, os quais, frequentemente, podem ser analisados por distribuições de Poisson e de Bernoulli, respectivamente, pertencentes à família exponencial. Duas das principais limitações para modelar esse tipo de dados são: (1) a ocorrência de superdispersão, ou seja, quando a variabilidade dos dados não é adequadamente descrita pelos modelos, que muitas vezes apresentam uma relação pré-estabelecida entre a média e a variância, e (2) a correlação existente entre medidas realizadas repetidas vezes na mesma unidade experimental. Uma forma de acomodar a superdispersão é pela utilização das distribuições binomial negativa e beta binomial, ou seja, pela inclusão de um efeito aleatório com distribuição gama quando se considera dados provenientes de contagens e um efeito aleatório com distribuição beta quando se considera dados binários, ambos introduzidos de forma multiplicativa. Para acomodar a correlação entre as medidas realizadas no mesmo indivíduo podem-se incluir efeitos aleat órios com distribuição normal no preditor linear. Esses situações podem ocorrer separada ou simultaneamente. Molenberghs et al. (2010) propuseram modelos que generalizam os modelos lineares generalizados mistos Poisson-normal e Bernoulli-normal, incorporando aos mesmos a superdispersão. Esses modelos foram formulados e ajustados aos dados, usando-se o método da máxima verossimilhança. Entretanto, para um modelo de efeitos aleatórios, é natural pensar em uma abordagem Bayesiana. Neste trabalho, são apresentados modelos Bayesianos hierárquicos para dados longitudinais, na forma de contagens e binários que apresentam superdispersão. A análise Bayesiana hierárquica é baseada no método de Monte Carlo com Cadeias de Markov (MCMC) e para implementação computacional utilizou-se o software WinBUGS. A metodologia para dados na forma de contagens é usada para a análise de dados de um ensaio clínico em pacientes epilépticos e a metodologia para dados binários é usada para a análise de dados de um ensaio clínico para tratamento de dermatite. / Longitudinal count and binary data are very common, which often can be analyzed by Poisson and Bernoulli distributions, respectively, members of the exponential family. Two of the main limitations to model this data are: (1) the occurrence of overdispersion, i.e., the phenomenon whereby variability in the data is not adequately captured by the model, and (2) the accommodation of data hierarchies owing to, for example, repeatedly measuring the outcome on the same subject. One way of accommodating overdispersion is by using the negative-binomial and beta-binomial distributions, in other words, by the inclusion of a random, gamma-distributed eect when considering count data and a random, beta-distributed eect when considering binary data, both introduced by multiplication. To accommodate the correlation between measurements made in the same individual one can include normal random eects in the linear predictor. These situations can occur separately or simultaneously. Molenberghs et al. (2010) proposed models that simultaneously generalizes the generalized linear mixed models Poisson-normal and Bernoulli-normal, incorporating the overdispersion. These models were formulated and tted to the data using maximum likelihood estimation. However, these models lend themselves naturally to a Bayesian approach as well. In this paper, we present Bayesian hierarchical models for longitudinal count and binary data in the presence of overdispersion. A hierarchical Bayesian analysis is based in the Monte Carlo Markov Chain methods (MCMC) and the software WinBUGS is used for the computational implementation. The methodology for count data is used to analyse a dataset from a clinical trial in epileptic patients and the methodology for binary data is used to analyse a dataset from a clinical trial in toenail infection named onychomycosis.

[pt] Esta dissertação, além de apresentar uma abordagem de um caso prático real, fez reunir as técnicas estatísticas necessárias ao trato de experimentos envolvendo misturas. Foi visto que as metodologias adotadas em Projeto de Experimentos devem ser adaptadas para possibilitar o trato de problemas com misturas, já que há a necessidade de considerar a restrição básica desse tipo de experimento, o qual amarra a soma das proporções dos componentes, que deve ser sempre igual a 1, ou seja, 100%. O experimento do misto de retardo, objeto principal e motivador dessa dissertação, é um experimento com mistura, em que as proporções de todos os três componentes possuem restrições superiores e inferiores simultaneamente. Com essas restrições, o espaço fatorial restrito fica bem distorcido em relação ao simplex, havendo, portanto, a necessidade de geração de um design D-ótimo. Como houve a indicação de que a variância da resposta não é constante, no caso do misto de retardo, recorreu-se aos Modelos Lineares Generalizados, especificamente ao método da Quase- Verossimilhança. De posse do modelo adequado, pôde-se então determinar a proporção dos componentes do misto de retardo, tendo em vista o atendimento da especificação de projeto. / [en] This dissertation presents a real pactical case, and besides, it puts together the statistical techniques for the treatment of Mixture Experiments. It was presented, that the Design of Experiments techniques must be adapted in order to make possible the treatment of problems with mixtures, because the basic constraint in this type of experiment must be taken into account, that is, the sum of the proportions of all mixture components must be equal to 1 or 100%. The delay compound experiment, the main and motivating object in this dissertation, is a mixture experiment with simultaneous constraints in the proportions of all its three components. With these constraints, it is possible to observe a distortion in the restricted factorial design space in comparison to the simplex one. Therefore, it was necessary to generate a D-optimal design. When there was an indication that the response variance is not constant, in the case of the delay compound, the Generalized Linear Models, specifically the Quasi- Likelihood method was used to fit an adequate model. With the adequate model, it was possible to find the proportion of each component of the delay compound in order to attend the design specification.

Dentre as possíveis propostas encontradas na literatura estatística para analisar dados oriundos de estudos com observações correlacionadas, estão os modelos condicionais e os modelos marginais. Diversas técnicas têm sido propostas para a análise de diagnóstico nesses modelos. O objetivo deste trabalho é apresentar algumas das técnicas de diagnóstico disponíveis para os dois tipos de modelos e avaliá-las por meio de estudos de simulação. As técnicas apresentadas também foram aplicadas em um conjunto de dados reais. / Conditional and marginal models are among the possibilities in statistical literature to analyze data from studies with correlated observations. Several techniques have been proposed for diagnostic analysis in these models. The objective of this work is to present some of the diagnostic techniques available for both modeling approaches and to evaluate them by simulation studies. The presented techniques were also applied in a real dataset.

Made available in DSpace on 2011-05-04T12:42:00Z (GMT). No. of bitstreams: 0 Previous issue date: 2009 / (...) O objetivo desta tese foi modelar os fatores de risco associados à ocorrência de leptospirose urbana em diferentes contextos, com especial atenção para aspectos espaciais e temporais. Foram utilizadas técnicas de modelagem tais como, modelos generalizados aditivos e mistos. Também explorou-se técnicas de detecção de aglomerados espaço-temporais. (...) / Leptospirosis, a disease caused by pathogenic spirochete of the genus Leptospira, is one of the most widespread zoonoses in the world, considered a major public health problem associated with the lack of sanitation and poverty. It is endemic in Brazil, data from surveillance show that outbreaks of leptospirosis occur as cyclical annual epidemics during rainfalls. The aim of this thesis was modeling the risk factors associated with the occurrence of leptospirosis in di erent urban contexts, with particular attention to spatial and temporal aspects. We used some modeling techniques such as generalized additive and mixed models. Techniques for detection space-time clusters were also explored. This thesis has prioritized the use of free softwares - R, ubuntu linux operating system, LATEX , SatScan (this is not open source but free). This thesis was prepared in the form of three articles. In the rst article is presented a spatio-temporal analysis of leptospirosis cases occurrence in Rio de Janeiro between 1997 and 2002. Using the detection of space-time clusters - \outbreaks" method - were statistically signi cant only cluster ocorred in 1997 and 1998. Generalized Linear Mixed Models were used to evaluate the risk factors associated with the occurrence of cases that belonged to outbreaks in endemic cases. The cases belonging to the outbreaks are associated with the occurrence of rainfall over 4 mm (OR, 3.71; 95% CI, 1.83 - 7.51). There were no signi cant associations with socioeconomic covariates, in other words, being endemic or epidemic leptospirosis occurs in the same population. The second and third articles examined a seroprevalence survey and seroconversion cohort conducted in Pau da Lima community, Salvador, Bahia. In both Generalized Additive Models were used to t the exposure variables both in individuals and peridomicile context, as well as to estimate the spatial area of leptospirosis risk. The signi cant variables were: gender, age, presence of rats in the peridomicile, domicile near a trash collectin or an open sewer and domicile altitude above sea level. Studies show that individual and contextual variables explain much of the spatial variability of leptospirosis, but there are still factors that were not measured in the studies but which should be investigated. The maps of risk of seroprevalence and seroconversion show distinct regions where the spatial e ect is signi cantly di erent from the global average. It is still lack for a more robust integration between the professionals who develop and operate the GIS, epidemiologists and biostatistics. This integration represents an important advance enabling the development and use of these techniques in Public Health support. The study of prevalence and incidence of endemic areas, in the leptospirosis context, it is very complex and still grow up. The reunion of professional specialists from several areas of human knowledge (eg, clinicians, epidemiologists, geographers, biologists, statisticians, engineers, etc.), it is essential to advance the knowledge about the disease and their relationship to social inequality and environmental well to contribute to the creation of efficient and e ective measures to control endemic diseases.

Dentre as possíveis propostas encontradas na literatura estatística para analisar dados oriundos de estudos com observações correlacionadas, estão os modelos condicionais e os modelos marginais. Diversas técnicas têm sido propostas para a análise de diagnóstico nesses modelos. O objetivo deste trabalho é apresentar algumas das técnicas de diagnóstico disponíveis para os dois tipos de modelos e avaliá-las por meio de estudos de simulação. As técnicas apresentadas também foram aplicadas em um conjunto de dados reais. / Conditional and marginal models are among the possibilities in statistical literature to analyze data from studies with correlated observations. Several techniques have been proposed for diagnostic analysis in these models. The objective of this work is to present some of the diagnostic techniques available for both modeling approaches and to evaluate them by simulation studies. The presented techniques were also applied in a real dataset.

[ES] La Ley de Prevención de Riesgos Laborales de 8 de noviembre de 1995 (LPRL), en vigor desde el 10 de febrero de 1996, establece en su artículo 5: "tendrá por objeto la promoción de la mejora de las condiciones de trabajo dirigida a elevar el nivel de protección de la seguridad y la salud de los trabajadores en el trabajo." En esta Tesis se ha investigado la evolución de los índices de siniestralidad laboral y su relación con la evolución de diferentes variables explicativas relacionadas con el desarrollo normativo, el mercado de trabajo, la estructura productiva, las condiciones de empleo y las condiciones individuales, entre otras, para el caso de España y en el periodo 1995-2017, que abarca desde la promulgación de la LPRL hasta fechas recientes donde se disponía de los datos históricos necesarios. La investigación se ha centrado en los índices de salud más relevantes según su significado en términos de riesgo y/o sus componentes. El objetivo de la investigación ha sido el encontrar evidencias sobre relaciones causa-efecto entre índices y variables, a partir de las cuales extraer lecciones que facilitarán una mejor planificación de la acción preventiva. Para ello, se han propuesto varios modelos explicativos utilizando diferentes herramientas estadísticas, que han permitido formular de manera explícita y analizar la relación entre la evolución de los indicadores de salud ocupacional y la evolución de las principales variables explicativas. En términos generales puede concluirse que la implantación de dicha ley y normativa que la acompaña ha tenido un impacto positivo en las condiciones de trabajo y en consecuencia sobre el nivel de seguridad y salud de los trabajadores desde entonces y hasta la fecha. Sin embargo, se observan diferentes comportamientos cíclicos en la evolución de los indicadores, tales como los índices de incidencia, frecuencia y gravedad, que pone de manifiesto su dependencia de la naturaleza y comportamiento cíclico de algunas de las variables explicativas más importantes relacionadas con ciclos económicos, mercado de trabajo, estructura productiva, etc. Además, se observa como aspectos tales como la pertenencia a grupos de edad jóvenes o expertos, el nivel de estudios, determinadas categorías profesionales, y algunos sectores particulares tienen efectos significativos sobre los valores alcanzados por los índices de siniestralidad. En cambio, otros, como el trabajo a tiempo parcial o la contratación temporal no manifiestan tener tanta repercusión sobre los indicadores. / [CA] Partint de les dades corresponents als accidents ocorreguts en el període 1995-2017, es La Llei de Prevenció de Riscos Laborals de 8 de novembre de 1995 (*LPRL), en vigor des del 10 de febrer de 1996, estableix en el seu article 5: "tindrà per objecte la promoció de la millora de les condicions de treball dirigida a elevar el nivell de protecció de la seguretat i la salut dels treballadors en el treball." En aquesta Tesi s'ha investigat l'evolució dels índexs de sinistralitat laboral i la seua relació amb l'evolució de diferents variables explicatives relacionades amb el desenvolupament normatiu, el mercat de treball, l'estructura productiva, les condicions d'ocupació i les condicions individuals, entre altres, per al cas d'Espanya i en el període 1995-2017, que abasta des de la promulgació de la LPRL fins a dates recents on es disposava de les dades històriques necessàries. La investigació s'ha centrat en els índexs de salut més rellevants segons el seu significat en termes de risc i/o els seus components. L'objectiu de la investigació ha sigut el trobar evidències sobre relacions causa-efecte entre índexs i variables, a partir de les quals extraure lliçons que facilitaran una millor planificació de l'acció preventiva. Per a això, s'han proposat diversos models explicatius utilitzant diferents eines estadístiques, que han permés formular de manera explícita i analitzar la relació entre l'evolució dels indicadors de salut ocupacional i l'evolució de les principals variables explicatives. En termes generals pot concloure's que la implantació d'aquesta llei i normativa que l'acompanya ha tingut un impacte positiu en les condicions de treball i en conseqüència sobre el nivell de seguretat i salut dels treballadors des de llavors i fins hui. No obstant això, s'observen diferents comportaments cíclics en l'evolució dels indicadors, com ara els índexs d'incidència, freqüència i gravetat, que posa de manifest la seua dependència de la naturalesa i comportament cíclic d'algunes de les variables explicatives més importants relacionades amb cicles econòmics, mercat de treball, estructura productiva, etc. A més, s'observa com a aspectes com ara la pertinença a grups d'edat joves o experts, el nivell d'estudis, determinades categories professionals, i alguns sectors particulars tenen efectes significatius sobre els valors aconseguits pels índexs de sinistralitat. En canvi, uns altres, com el treball a temps parcial o la contractació temporal no manifesten tindre tanta repercussió sobre els indicadors. / [EN] The Occupational Risk Prevention Act of November 8, 1995 (ORPA), in force since February 10, 1996, establishes in its article 5: "will have as its objective the promotion of the improvement of working conditions aimed at raise the level of protection of the safety and health of workers at work. " This thesis has investigated the evolution of the occupational accident rates and their relationship with the evolution of different explanatory variables related to regulatory development, the labor market, the productive structure, employment conditions and individual conditions, among others, in the case of Spain and in the period 1995-2017, which ranges from the enactment of the LPRL to recent dates where the necessary historical data was available. Research has focused on the most relevant health indices according to their meaning in terms of risk and / or their components. The objective of the research has been to find evidence on cause-effect relationships between indices and variables, from which to extract lessons that will facilitate better planning of preventive action. To this end, several explanatory models have been proposed using different statistical tools, which have made it possible to explicitly formulate and analyze the relationship between the evolution of occupational health indicators and the evolution of the main explanatory variables. In general terms, it can be concluded that the implementation of said law and accompanying regulations has had a positive impact on working conditions and consequently on the level of health and safety of workers since then and to date. However, different cyclical behaviors are observed in the evolution of the indicators, such as incidence, frequency and severity indices, which highlights their dependence on the nature and cyclical behavior of some of the most important explanatory variables related to economic cycles, labor market, productive structure, etc. Furthermore, aspects such as belonging to young age groups or experts, educational level, certain professional categories, and some particular sectors are observed as having significant effects on the values reached by the accident rates. On the other hand, others, such as part-time work or temporary hiring, do not claim to have such an impact on the indicators. / Gallego Blasco, VS. (2021). Análisis de la incidencia de factores causales en la evolución de la siniestralidad laboral en España [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/168774 / TESIS

The Birnbaum-Saunders (BS) distribution is the most popular model used to describe lifetime process under fatigue. Throughout the years, this distribution has received a wide ranging of applications, demanding some more flexible extensions to solve more complex problems. One of the most well-known extensions of the BS distribution is the generalized Birnbaum- Saunders (GBS) family of distributions that includes the Birnbaum-Saunders special-case (BSSC) and the Birnbaum-Saunders generalized t (BSGT) models as special cases. Although the BS-SC distribution was previously developed in the literature, it was never deeply studied and hence, in this thesis, we provide a full Bayesian study and develop a tool to generate random numbers from this distribution. Further, we develop a very flexible regression model, that admits different degrees of skewness and kurtosis, based on the BSGT distribution using the generalized additive models for location, scale and shape (GAMLSS) framework. We also introduce a new extension of the BS distribution called the Birnbaum-Saunders power (BSP) family of distributions, which contains several special or limiting cases already published in the literature, including the GBS family. The main feature of the new family is that it can produce both unimodal and bimodal shapes depending on its parameter values. We also introduce this new family of distributions into the GAMLSS framework, in order to model any or all the parameters of the distribution using parametric linear and/or nonparametric smooth functions of explanatory variables. Throughout this thesis we present five different applications in real data sets in order to illustrate the developed theoretical results. / A distribuição Birnbaum-Saunders (BS) é o modelo mais popular utilizado para descrever processos de fadiga. Ao longo dos anos, essa distribuição vem recebendo aplicações nas mais diversas áreas, demandando assim algumas extensões mais flexíveis para resolver problemas mais complexos. Uma das extensões mais conhecidas na literatura é a família de distribuições Birnbaum-Saunders generalizada (GBS), que inclui as distribuições Birnbaum-Saunders casoespecial (BS-SC) e Birnbaum-Saunders t generalizada (BSGT) como modelos especiais. Embora a distribuição BS-SC tenha sido previamente desenvolvida na literatura, nunca foi estudada mais profundamente e, assim, nesta tese, um estudo bayesiano é desenvolvido acerca da mesma além de um novo gerador de números aleatórios dessa distribuição ser apresentado. Adicionalmente, um modelo de regressão baseado na distribuição BSGT é desenvolvido utilizando-se os modelos aditivos generalizados para locação, escala e forma (GAMLSS), os quais apresentam grande flexibilidade tanto para a assimetria como para a curtose. Uma nova extensão da distribuição BS também é apresentada, denominada família de distribuições Birnbaum-Saunders potência (BSP), que contém inúmeros casos especiais ou limites já publicados na literatura, incluindo a família GBS. A principal característica desta nova família é que ela é capaz de produzir formas tanto uni como bimodais dependendo do valor de seus parâmetros. Esta nova família também é introduzida na estrutura dos modelos GAMLSS para fornecer uma ferramenta capaz de modelar todos os parâmetros da distribuição como funções lineares e/ou não-lineares suavizadas de variáveis explicativas. Ao longo desta tese são apresentadas cinco diferentes aplicações em conjuntos de dados reais para ilustrar os resultados teóricos obtidos.

The Generalized Finite Element Method (GFEM) is a numerical method based on the Partition of Unity (PU) concept and inspired on both the Partition of Unity Method (PUM) and the hp-Cloud method. According to the GFEM, the PU is provided by first-degree Lagragian interpolation functions, defined over a mesh of elements similar to the Finite Element Method (FEM) meshes. In fact, the GFEM can be considered an extension of the FEM to which enrichment functions can be applied in specific regions of the problem domain to improve the solution. This technique has been successfully employed to solve problems presenting discontinuities and singularities, like those that arise in Fracture Mechanics. However, most publications on the method are related to linear analyses. The present thesis is a contribution to the few studies of nonlinear analyses of Solid Mechanics by means of the GFEM. One of its main topics is the derivation of a segment-to-segment generalized contact element based on the mortar method. Material and kinematic nonlinear phenomena are also considered in the numerical models. An Object-Oriented design was developed for the implementation of a GFEM nonlinear analyses framework written in Python programming language. The results validated the formulation and demonstrate the gains and possible drawbacks observed for the GFEM nonlinear approach. / O Método dos Elementos Finitos Generalizados (MEFG) é um método numérico baseado no conceito de partição da unidade (PU) e inspirado no Método da Partição da Unidade (MPU) e o método das Nuvens-hp. De acordo com o MEFG, a PU é obtida por meio de funções de interpolação Lagragianas de primeiro grau, definidas sobre uma rede de elementos similar àquela do Método dos Elementos Finitos (MEF). De fato, o MEFG pode ser considerado uma extensão do MEF para a qual se pode aplicar enriquecimentos em regiões específicas do domínio, buscando melhorias na solução. Esta técnica já foi aplicada com sucesso em problemas com descontinuidades e singularidades, como os originários da Mecânica da Fratura. Apesar disso, a maioria das publicações sobre o método está relacionada a análises lineares. A presente tese é uma contribuição aos poucos estudos relacionados a análises não-lineares de Mecânica dos Sólidos por meio do MEFG. Um de seus principais tópicos é o desenvolvimento de um elemento de contato generalizado do tipo segmento a segmento baseado no método mortar. Fenômenos não lineares devidos ao material e à cinemática também são considerados nos modelos numéricos. Um projeto de orientação a objetos para a implementação de uma plataforma de análises não-lineares foi desenvolvido, escrito em linguagem de programação Python. Os resultados validam a formulação e demonstram os ganhos e possíveis desvantagens da abordagem a problemas não lineares por meio do MEFG.

Este trabalho trata da combinação entre a formulação híbrida-mista de tensão (FHMT) (Freitas et al. (1996)), para a elasticidade plana, com o método dos elementos finitos generalizados (MEFG), Duarte et al. (2000). O MEFG se caracteriza como uma forma não-convencional do método dos elementos finitos (MEF) que resulta da incorporação a este de conceitos e técnicas dos métodos sem malha, como o enriquecimento nodal proposto do método das nuvens “hp". Como na FHMT são aproximados dois campos no domínio (tensão e deslocamento) e um no contorno (deslocamento), diferentes possibilidades de enriquecimento nodal são exploradas. Para a discretização do modelo híbrido-misto empregam-se elementos finitos quadrilaterais com funções de forma bilineares para o domínio e elementos lineares para o contorno. Essas funções são enriquecidas por funções polinomiais, trigonométricas, polinômios que proporcionam distribuição de tensões auto-equilibradas ou mesmo funções especiais relacionadas às soluções dos problemas de fratura. Uma extensão do teste numérico abordado em Zienkiewicz et al. (1986) é proposta como investigação inicial das condições necessárias para garantia de estabilidade da resposta numérica. O estudo da estabilidade é completado com a análise da condição de Babuška-Brezzi (inf-sup). Esta condição é aplicada nos elementos finitos quadrilaterais híbridos-mistos enriquecidos por meio de um teste numérico, denominado de inf-sup teste, desenvolvido com base no trabalho de Chapelle e Bathe (1993). Exemplos numéricos revelam que a FHMT é uma interessante alternativa para obtenção de boas estimativas para os campos de tensões e deslocamentos, usando-se enriquecimento sobre alguns nós de malhas pouco refinadas / This work presents a combination of hybrid-mixed stress model formulation (HMSMF) (Freitas et al. (1996)), to treat plane elasticity problems, with generalized finite element method (GFEM), (Duarte et al. (2000)). GFEM is characterized as a nonconventional formulation of the finite element method (FEM). GFEM is the result of the incorporation of concepts and techniques from meshless methods. One example of these techniques is the nodal enrichment that was formulated in the “hp" clouds method. Since two fields in domain (stress and displacement) and one in boundary (displacement) are approximated in the HMSMF, different possibilities of nodal enrichment are tested. For the discretization of the hybrid-mixed model quadrilateral finite elements with bilinear shape functions for the domain and linear elements for the boundary were employed. These functions are enriched with polynomial functions, trigonometric functions, polynomials that generate self-equilibrated stress distribution, or, even special functions connected with solutions of fracture problems. An extension of the numerical test cited in Zienkiewicz et al. (1986) is proposed as initial investigation of necessary conditions to assure the stability of the numerical answer. The stability study is completed with the analysis of the Babuška-Brezzi (inf-sup) condition. This last condition is applied to hybrid-mixed enrichment quadrilaterals finite elements by means of a numerical test, denominated inf-sup test, which was developed based on paper of Chapelle and Bathe (1993). Numerical examples reveal that HMSMF is an interesting alternative to obtain good estimates of the stress and displacement fields, using enrichment over some nodes of poor meshes