Um modelo de turbulência baseado no conceito de vórtice

O fenômeno da turbulência está presente na maioria dos escoamentos observados na indústria e na natureza. Muitas são as considerações a respeito das dificuldades relacionadas à caracterização dos escoamentos turbulentos. Uma das muitas questões trata do procedimento de análise do problema através da descrição estatística dos campos por grandezas “médias”, o que leva ao problema de fechamento e à modelagem do tensor de Reynolds, normalmente com modelos baseados no conceito de viscosidade turbulenta. Os modelos de turbulência já existentes apresentam algumas deficiências na previsão do escoamento, além de outras limitações, o que justifica a busca por novas abordagens para o tratamento da turbulência. Neste trabalho, o problema de fechamento é tratado segundo a modelagem turbulenta baseada no conceito de viscosidade turbulenta. Um novo modelo de turbulência é proposto, que admite a existência de vórtices imersos no escoamento e aplica conceitos e definições relacionados à identificação de vórtices, com o uso do critério de identificação Q , que caracteriza a região do escoamento ocupada pelo vórtice. Propõe-se a investigação da aplicabilidade do critério Q em conjunto com o modelo k − ε , para o desenvolvimento de um novo modelo de turbulência chamado k − ε −Q . Validou-se a aplicabilidade do modelo através de um código numérico computacional para tratamento de escoamentos turbulentos. A solução numérica foi obtida através da discretização do domínio fluido, utilizando o método de volumes finitos e o método multigrid foi utilizado para resolver o sistema linear resultante. Como verificação, foi utilizado este modelo de turbulência para simular o escoamento em uma cavidade quadrada com tampa deslizante e o escoamento turbulento sobre um degrau. Os resultados obtidos foram confrontados com dados experimentais e demonstraram que o modelo aqui proposto se apresenta mais eficiente que o clássico modelo k − ε , no tratamento da turbulência nesses dois problemas clássicos.

Construção de um modelo experimental simulando condições de perfilagem de indução

A ferramenta de indução tem sido, desde 1949, largamente utilizada na prospecção de petróleo pois auxilia na localização de formações que contém hidrocarbonetos. Ainda existem, porém, várias restrições quanto ao uso deste equipamento. Camadas excessivamente finas, poços preenchidos com lamas muito condutivas e grandes diâmetros de invasão são exemplos de situações onde seu desempenho não é bom. Para que se possa minimizar estas limitações é necessário o conhecimento da resposta da sonda diante de condições conhecidas. Isto pode ser feito via modelamento analógico e/ou através de modelamento numérico. O presente trabalho consiste na montagem de um sistema de modelamento analógico capaz de simular o perfil de duas bobinas coaxiais frente a meios infinitos e camadas espessas invadidas. O meio geológico é representado através de solução salina de diferentes condutividades. O fator de escala linear é igual a 20 e o equipamento executa medidas na faixa dos 50 kHz aos 100 kHz. Este sistema será utilizado em futuros estudos envolvendo disposições de bobinas diferentes das normalmente utilizadas e em estudos de outros problemas encontrados na perfilagem de indução.

Estudo numérico de modelos de arrasto e do coeficiente de restituição no escoamento gás-sólido em leito fluidizado

Pós-graduação em Engenharia Mecânica - FEG

Simulação numérica de fluxo de águas subterrâneas na bacia do Rio Grande (BA)

The Urucuia Aquifer System represents a strategic water source in western Bahia. Its baseflow is responsible for the flow rate of the main tributaries of São Francisco river left bank in the dry season, including the Rio Grande, its main tributary in Bahia state. This river has a hydrological regime heavily affected by groundwater and is located in a region with conflicts about water resources. The aquifers geology is constituted by neocretacious sandstones of Urucuia Group subdivided in Posse Formation and Serra das Araras Formation. The embasement is formed by neoproterozoic rocks of Bambuí Group. This work focuses on an important tool application, the mathematical model, whose function is represent approximately and suitably the reality so that can assist in different scenarios simulations and make predictions. Many studies developed in this basin provided the conceptual model basis including a full free aquifer, lithological and hydraulical homogeneity in entire basin, null flux at plateau borders and aquifer base. The finite element method is the numerical method used and FEFLOW the computational algorithm. The simulated area was discretized in a single layer with 27.357,6 km² (314.432 elements and 320.452 nodes) totaling a 4.249,89 km³ volume. Were utilized 21 observation wells from CERB to calibrate the model. The terrain topography was obtained by SRTM data and the impermeable base was generated by interpolation of descriptive profiles from wells and electric vertical drilling from previous studies. Works in this area obtained mean recharge rates varying approximately from 20% to 25% of average precipitation, thus the values of model recharge zones varying in this range. Were distributed 4 hydraulic conductivity zones: (K1) west zone with K=6x10-5 m/s; (K2) center-east zone with K=3x10-4 m/s; (K3) far east zone with K=5x10-4 m/s; e (K4) east-north zone with K=1x10-5 m/s. Thereby was incorporated to the final conceptual model...

Comparação de metodologias na definição de fluxos acumulados a partir de modelos digitais de elevação do terreno aplicado a suscetibilidade de inundações

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Aportaciones al método de la celda en el diseño y análisis de un modelo de máquina rotativa trifásica de inducción magnética

Programa de doctorado: Tecnología industrial

Sistemas planos de estructuras (aproximación al modelo de análisis de piezas rectas de hormigón armado en estado último)

En esta tesis, se trata del establecimiento de las bases de análisis plástico de piezas de hormigón armado en servicio y estado último procurando sentar los soportes de funcionamiento que concuerdan con la realidad y normativa vigente. El intento de establecer estas bases, con el máximo rigor para su uso, ha llevado a consecuencias lejanas cualitativas y cuantitativas sobre los modelos elásticos y elastoplásticos en uso,-momento tope ect- que, incluyendo lo complejo del problema plástico que se trata, no adoptan formulaciones adecuadas, o los plantemientos analíticos, no son más coherentes con los resultados obtenidos, en el dimensionado, que otros de base más sencilla, limitándose a estrategias de análisis de sección basadas en condicones de deformación -como serían los métodos de Cross y Matrioficiales- , estableciendose un claro divorcio entre mecanismo de cálculo y geometría. Este trabajo, intenta abrir una puerta al análisis de la estructura bajo condiciones puramente de deformación y solo supone en este campo el primer paso de un largo cemino que iría de la pieza a la estructura y de ésta al edificio posteriormente, quedando así sentados los criterios básicos para el análisis plástico global de estructuras. Todo ello a partir de la imposibilidad de avanzar en el análisis de secciones, si no es penetrando en el proceso de plastificación de las mismas y el de sus armaduras, así como el de las zonas de afección de las diferentes plastificaciones locales. La inclusión de esta concepción de análisis, en lo referente a obtener gráficos reales de respuesta interna de la pieza, es posible gracias al uso, en este caso del miniordenador 9825-HP utilizado, en donde se aplica de forma reiterada los consabidos métodos de discretización elementales con los cuales se obtienen decorosas precisiones sobre los resultados. Este estudio, que en ocasiones, alude a normas, más discribir un comportamiento que sentar doctrina. Como se ha realizado solo par forjados se añade una colección de ábacos de utilización practica de éstos. Naturalmente las bases establecidas, se hacen extensivas a cualquier pieza de directriz recta de hormigón armado, en donde geometría y materiales -hormigón y acero- pueden ser variables y están contempladas dentro de programas usados. En el terreno de las conclusiones, los forjados no deben ser asimilados a vigas, como tradicionalmente venía haciendose, ya que los métodos de relajación empledos no coinciden con la realidad. Quedando estos límites plenamente establecidos en el desarrollo del trabajo. Se considera el tema relevante en lo reeferente a un análisis de la estructura tendente a un aprovechamento más adecuado de los materiales en un supuesto futuro de escasez de materia prima, quizás no lejano, y sobretodo en la posible aproximación al modelo matemático, que si bien normativa y bibliografía aludían constantemente, nunca precisaron las bases de su comportamiento.

Empuje sobre un muro de contención vertical, producido por la acción de una carga sobre un terreno elástico

En el presente trabajo se desarrolla un método de cálculo del empuje que produce sobre un muro una carga situada sobre un terreno adyacente a dicho muro. Por medio de hipótesis físicas sencillas, se llega a resultados que parecen diferir de los obtenidos experimentalmente en márgenes tolerables, dentro de la aproximación normalmente usada en los cálculos. Se incluye una generalización para dos muros iguales y paralelos, fundada en un artificio análogo al empleado en Aerodinamica, para estudiar la interferencia entre las paredes y el modelo en los ensayos en tunel aerodinámico. Finalmente se presenta una familia de curvas adimensionales que dan el empuje para sobrecargas geométricamente semejantes, situadas a diferentes distancias del muro. Familias de este tipo pueden ser de gran utilidad, ya que se ha comprobado que, para los pesos y presiones de inflado usuales, las huellas de las ruedas de diferentes aviones son geometricamente semejantes. Estas familias podrian servir para un ulterior intento de obtención de fórmulas aproximadas de más fácil aplicación.

Diseño de un modelo de riesgo integral de incendios forestales mediante técnicas multicriterio y su automatización en sistemas de información geográfica : una aplicación en la Comunidad Valenciana

Los incendios forestales son la principal causa de mortalidad de árboles en la Europa mediterránea y constituyen la amenaza más seria para los ecosistemas forestales españoles. En la Comunidad Valenciana, diariamente se despliega cerca de un centenar de vehículos de vigilancia, cuya distribución se apoya, fundamentalmente, en un índice de riesgo de incendios calculado en función de las condiciones meteorológicas. La tesis se centra en el diseño y validación de un nuevo índice de riesgo integrado de incendios, especialmente adaptado a la región mediterránea y que facilite el proceso de toma de decisiones en la distribución diaria de los medios de vigilancia contra incendios forestales. El índice adopta el enfoque de riesgo integrado introducido en la última década y que incluye dos componentes de riesgo: el peligro de ignición y la vulnerabilidad. El primero representa la probabilidad de que se inicie un fuego y el peligro potencial para que se propague, mientras que la vulnerabilidad tiene en cuenta las características del territorio y los efectos potenciales del fuego sobre el mismo. Para el cálculo del peligro potencial se han identificado indicadores relativos a los agentes naturales y humanos causantes de incendios, la ocurrencia histórica y el estado de los combustibles, extremo muy relacionado con la meteorología y las especies. En cuanto a la vulnerabilidad se han empleado indicadores representativos de los efectos potenciales del incendio (comportamiento del fuego, infraestructuras de defensa), como de las características del terreno (valor, capacidad de regeneración…). Todos estos indicadores constituyen una estructura jerárquica en la que, siguiendo las recomendaciones de la Comisión europea para índices de riesgo de incendios, se han incluido indicadores representativos del riesgo a corto plazo y a largo plazo. El cálculo del valor final del índice se ha llevado a cabo mediante la progresiva agregación de los componentes que forman cada uno de los niveles de la estructura jerárquica del índice y su integración final. Puesto que las técnicas de decisión multicriterio están especialmente orientadas a tratar con problemas basados en estructuras jerárquicas, se ha aplicado el método TOPSIS para obtener la integración final del modelo. Se ha introducido en el modelo la opinión de los expertos, mediante la ponderación de cada uno de los componentes del índice. Se ha utilizado el método AHP, para obtener las ponderaciones de cada experto y su integración en un único peso por cada indicador. Para la validación del índice se han empleado los modelos de Ecuaciones de Estimación Generalizadas, que tienen en cuenta posibles respuestas correlacionadas. Para llevarla a cabo se emplearon los datos de oficiales de incendios ocurridos durante el período 1994 al 2003, referenciados a una cuadrícula de 10x10 km empleando la ocurrencia de incendios y su superficie, como variables dependientes. Los resultados de la validación muestran un buen funcionamiento del subíndice de peligro de ocurrencia con un alto grado de correlación entre el subíndice y la ocurrencia, un buen ajuste del modelo logístico y un buen poder discriminante. Por su parte, el subíndice de vulnerabilidad no ha presentado una correlación significativa entre sus valores y la superficie de los incendios, lo que no descarta su validez, ya que algunos de sus componentes tienen un carácter subjetivo, independiente de la superficie incendiada. En general el índice presenta un buen funcionamiento para la distribución de los medios de vigilancia en función del peligro de inicio. No obstante, se identifican y discuten nuevas líneas de investigación que podrían conducir a una mejora del ajuste global del índice. En concreto se plantea la necesidad de estudiar más profundamente la aparente correlación que existe en la provincia de Valencia entre la superficie forestal que ocupa cada cuadrícula de 10 km del territorio y su riesgo de incendios y que parece que a menor superficie forestal, mayor riesgo de incendio. Otros aspectos a investigar son la sensibilidad de los pesos de cada componente o la introducción de factores relativos a los medios potenciales de extinción en el subíndice de vulnerabilidad. Summary Forest fires are the main cause of tree mortality in Mediterranean Europe and the most serious threat to the Spanisf forest. In the Spanish autonomous region of Valencia, forest administration deploys a mobile fleet of 100 surveillance vehicles in forest land whose allocation is based on meteorological index of wildlandfire risk. This thesis is focused on the design and validation of a new Integrated Wildland Fire Risk Index proposed to efficient allocation of vehicles and specially adapted to the Mediterranean conditions. Following the approaches of integrated risk developed last decade, the index includes two risk components: Wildland Fire Danger and Vulnerability. The former represents the probability a fire ignites and the potential hazard of fire propagation or spread danger, while vulnerability accounts for characteristics of the land and potential effects of fire. To calculate the Wildland Fire Danger, indicators of ignition and spread danger have been identified, including human and natural occurrence agents, fuel conditions, historical occurrence and spread rate. Regarding vulnerability se han empleado indicadores representativos de los efectos potenciales del incendio (comportamiento del fuego, infraestructurasd de defensa), como de las características del terreno (valor, capacidad de regeneración…). These indicators make up the hierarchical structure for the index, which, following the criteria of the European Commission both short and long-term indicators have been included. Integration consists of the progressive aggregation of the components that make up every level in risk the index and, after that, the integration of these levels to obtain a unique value for the index. As Munticriteria methods are oriented to deal with hierarchically structured problems and with situations in which conflicting goals prevail, TOPSIS method is used in the integration of components. Multicriteria methods were also used to incorporate expert opinion in weighting of indicators and to carry out the aggregation process into the final index. The Analytic Hierarchy Process method was used to aggregate experts' opinions on each component into a single value. Generalized Estimation Equations, which account for possible correlated responses, were used to validate the index. Historical records of daily occurrence for the period from 1994 to 2003, referred to a 10x10-km-grid cell, as well as the extent of the fires were the dependant variables. The results of validation showed good Wildland Fire Danger component performance, with high correlation degree between Danger and occurrence, a good fit of the logistic model used and a good discrimination power. The vulnerability component has not showed a significant correlation between their values and surface fires, which does not mean the index is not valid, because of the subjective character of some of its components, independent of the surface of the fires. Overall, the index could be used to optimize the preventing resources allocation. Nevertheless, new researching lines are identified and discussed to improve the overall performance of the index. More specifically the need of study the inverse relationship between the value of the wildfire Fire Danger component and the forested surface of each 10 - km cell is set out. Other points to be researched are the sensitivity of the index component´s weight and the possibility of taking into account indicators related to fire fighting resources to make up the vulnerability component.

Modelo proceso-respuesta de recesión de acantilados por variación del nivel del mar. Aplicación en la Costa de Holderness (Reino Unido).

El desarrollo de modelos predictivos de retroceso de acantilados costeros se encuentra limitado desde el punto de vista geomorfológico debido a la complejidad e interacción de los procesos de desarrollo espacio-temporal que tienen lugar en la zona costera. Los códigos de cálculo existentes en la actualidad incorporan modelos probabilísticos para resolver la estabilidad geomecánica de los materiales. Éstos pueden reproducir tasas de recesión históricas; sin embargo, el comportamiento bajo condiciones variables (cambio climático) no tiene por qué ser el mismo. Se presenta un modelo de proceso-respuesta de retroceso de acantilados costeros que incorpora el comportamiento geomecánico de costas compuestas por tills, incluye la acción protectora de los derrubios presentes al pie del acantilado frente a los procesos erosivos y permite evaluar el efecto del ascenso del nivel medio del mar asociado al cambio climático. Para ello, se acoplan la dinámica marina: nivel medio del mar, mareas y oleaje; con la evolución del terreno: erosión, desprendimiento rocoso y formación de talud de derrubios. La erosión del acantilado se calcula en función del nivel del mar, de las condiciones del oleaje incidente, de la pendiente de la plataforma rocosa y de la resistencia del material rocoso afectado en cada ciclo de marea. Al finalizar cada ciclo, se evalúa la ruptura del acantilado según un criterio geomecánico de roturas por vuelco y en caso de rotura, se genera un talud de derrubios al pie del mismo. El modelo ha sido validado frente a datos reales en la costa de Holderness, Yorkshire, Reino Unido. Los resultados obtenidos presentan un importante avance en los modelos de recesión costera, especialmente en su relación con las condiciones geomecánicas del medio y su respuesta ante condiciones variables del nivel medio del mar producidas por el cambio climático.

Un modelo predictivo de la evolución en costas rocosas.

El retroceso de las costas acantiladas es un fenómeno muy extendido sobre los litorales rocosos expuestos a la incidencia combinada de los procesos marinos y meteorológicos que se dan en la franja costera. Este fenómeno se revela violentamente como movimientos gravitacionales del terreno esporádicos, pudiendo causar pérdidas materiales y/o humanas. Aunque el conocimiento de estos riesgos de erosión resulta de vital importancia para la correcta gestión de la costa, su predicción es complicada. Los modelos de predicción publicados son escasos y con importantes inconvenientes: extrapolación, extienden información registros históricos; empíricos, sobre registros históricos estudian respuesta al cambio de un parámetro; estocásticos, determinan la cadencia y magnitud de los eventos futuros extrapolando las distribuciones de probabilidad extraídas de catálogos históricos; proceso-respuesta, de estabilidad y propagación del error inexplorada; en EDPs, computacionalmente costosos y poco exactos. En este trabajo se desarrolla un modelo combinado de proceso-respuesta basado en incorporar un balance de fuerzas de los mecanismos que actúan sobre el proceso erosivo en el frente del acantilado. El modelo simula la evolución espacio-temporal de un perfil-2D del acantilado, formado por materiales heterogéneos, acoplando la dinámica marina con la evolución del terreno en cada periodo de marea. Integra en cada periodo una función de erosión, dependiente de la pendiente de la zona afectada, que se desplaza sobre la onda de marea. Se ha estudiado el error de discretización del modelo y su propagación en el tiempo a partir de las soluciones exactas para los dos primeros periodos de marea para diferentes aproximaciones numéricas en la integración y de la pendiente. Los resultados obtenidos han permitido justificar las discretizaciones que minimizan el error y los métodos de aproximación más adecuados.

Modelado numérico de la fisuración por corrosión de armaduras

Se propone un modelo de elementos finitos, para el estudio de los efectos mecánicos de la corrosión del acero en el hormigón armado. Dicho modelo se basa en una imposición ordenada de la corrosión de los elementos de acero, mediante un incremento ficticio de temperatura que produce unos efectos análogos, a la vez que la fisuración del hormigón se recoge mediante un modelo conocido de fisura difusa. Se incluye un ejemplo resuelto que muestra el modo de aplicación de este modelo a problemas reales y el tipo de resultados que se pueden obtener, que se refieren especialmente al estudio de la propagación de la fisuración en función de la penetración de la corrosión en el acero.

Estudio numérico y experimental del efecto de tamaño en modo mixto de los paneles de yeso laminado y lana de roca

Los paneles sándwich de yeso laminado y lana de roca presentan una abundante patología de fisuración debida a flechas excesivas de forjados. Existe, por tanto, la necesidad de avanzar en la simulación y predicción de comportamiento bajo solicitaciones de tracción y cortante de ese tipo de paneles, a pesar de que en las aplicaciones habituales no tienen responsabilidad estructural. El comportamiento de este material puede ser considerado cuasi-frágil, y en base a ello en este trabajo ha sido estudiado haciendo uso de modelos de fisura cohesiva, cuya aplicación a otros materiales cuasifrágiles, como el hormigón, ha aportado resultados muy satisfactorios. En esta comunicación se presenta el trabajo realizado para estudiar el efecto del tamaño del elemento de yeso laminado y lana de roca en su comportamiento mecánico-resistente. Para ello se diseñó una campaña de ensayos en modo mixto sobre probetas de diferente tamaño. Se han realizado ensayos de flexión en tres puntos en modo mixto de unas probetas entalladas, geométricamente similares y de diferente tamaño, obteniéndose las curvas carga-desplazamiento y cargaabertura de la boca de la entalla. Para simular numéricamente el comportamiento en fractura del panel en modo mixto se ha utilizado un modelo de elementos finitos con fisura embebida basado en la fisura cohesiva en el que se introducen como entrada los parámetros obtenidos a partir de la experimentación de trabajos anteriores, obteniéndose un buen ajuste. En función de estos resultados se analiza el efecto del tamaño en los paneles. Sandwich panels of laminated gypsum and rockwool have an abundant pathology of cracking due to excessive slabs deflection. Therefore, it is necessary to progress in the simulation and prediction of behaviour under tensile and shear load of such panels, although in typical applications have no structural responsability. The behaviour of this material may be considered quasi-brittle and, based on this idea, in this work has been studied using a cohesive crack model that has been applied to other quasi-brittle materials, such as concrete, and has provided very satisfactory results. This communication presents the work carried out to study the size effect of the specimen of plasterboard and rockwool in its mechanical and resistant behaviour. The authors designed an experimental campaign under mixed mode composed by testing specimens of different sizes. Assymetrical three-point bending tests have been performed on notched specimens, geometrically similar and of different size, to obtain load-displacement and load-crack moutn opening displacement curves. To numerically simulate the mixed-mode fracture behaviour of the panels we have used a finite element model with embedded crack, based on the cohesive crack model, using as input the experimental parameters obtained in previous work, obtaining a good adjustment. Based on these results we analyze the size effect of the panels

Modelo de evaluación socioeconómica de los sistemas de bicicletas públicas

El objetivo de esta primera parte de este proyecto es realizar un modelo socioeconómico de evaluación para la instalación de un sistema de bicicleta pública con la ayuda del método de evaluación: el análisis multicriterio en el que se analizarán criterios de diferente tipología que se escogerán tras el estudio de unas encuestas realizadas a usuarios del sistema Vélo’v en Lyon y tras unas entrevistas realizadas a expertos y trabajadores en el campo de la bicicleta pública. Este modelo sirve como ayuda a la decisión en la fase de preparación de la creación de un SBP para asegurar el balance entre la rentabilidad financiera, la social y la medioambiental. Después de la toma de decisión, se procederá a la realización del proyecto y a su instalación y por último su funcionamiento. Los objetivos resultan: • Descripción de los sistemas de bicicleta pública y el marco de acción • Examinar los estudios realizados sobre el tema • Descripción de los trabajos realizados para este proyecto y posterior análisis para su aplicación • Determinar los requisitos de diseño y de dimensionamiento para la implantación de un SBP • Encontrar los criterios adecuados que se emplearan en el análisis multicriterio para la evaluación • Describir los indicadores para medir cada criterio • La ponderación de estos criterios según las encuestas realizadas a los usuarios • La realización de una alternativa 0 que describa la situación de la ciudad sin la realización del proyecto • Recomendación final del modelo a través de un cuadro de resultados. El objeto de la tercera parte del proyecto es la descripción técnica y gráfica de una caseta de servicio del sistema de préstamos de bicicleta pública en San Lorenzo de El Escorial. La Caseta proyectada tendrá funciones principales: un taller para el mantenimiento del servicio y un espacio de almacenamiento de bicicletas. La Caseta-Taller es una medida necesaria en la implantación de un servicio de bici pública puesto que el uso de las bicicletas necesita un mantenimiento. En un terreno accidentado y no siempre propicio para la circulación de bicicletas, estas sufren muchos degastes y necesitan reparaciones y mantenimiento continuo. Además, el servicio en San Lorenzo de El Escorial está sometido a flujos pendulares de desplazamientos puesto que la disponibilidad de bicicletas y de espacios libres en las estaciones es muy variable. Es necesario precisar de un espacio para el almacenamiento de bicicletas, ya sea para controlar la disponibilidad o bien para esperar a recibir el mantenimiento necesario. Aprovechando la construcción de un edificio nuevo se destinará un espacio en la Casa-Taller para el control y la gestión del sistema. En esta parte se detallará gráficamente la implantación de los puntos de préstamos en la vía urbana como complemento al dimensionamiento del anterior trabajo.

Estudio numérico y experimental del efecto de tamaño en modo mixto de los paneles de yeso laminado y lana de roca

Los paneles sándwich de yeso laminado y lana de roca presentan una abundante patología de fisuración debida a flechas excesivas de forjados. Existe, por tanto, la necesidad de avanzar en la simulación y predicción de comportamiento bajo solicitaciones de tracción y cortante de ese tipo de paneles, a pesar de que en las aplicaciones habituales no tienen responsabilidad estructural. El comportamiento de este material puede ser considerado cuasi-frágil, y en base a ello en este trabajo ha sido estudiado haciendo uso de modelos de fisura cohesiva, cuya aplicación a otros materiales cuasifrágiles, como el hormigón, ha aportado resultados muy satisfactorios. En esta comunicación se presenta el trabajo realizado para estudiar el efecto del tamaño del elemento de yeso laminado y lana de roca en su comportamiento mecánico-resistente. Para ello se diseñó una campaña de ensayos en modo mixto sobre probetas de diferente tamaño. Se han realizado ensayos de flexión en tres puntos en modo mixto de unas probetas entalladas, geométricamente similares y de diferente tamaño, obteniéndose las curvas carga-desplazamiento y cargaabertura de la boca de la entalla. Para simular numéricamente el comportamiento en fractura del panel en modo mixto se ha utilizado un modelo de elementos finitos con fisura embebida basado en la fisura cohesiva en el que se introducen como entrada los parámetros obtenidos a partir de la experimentación de trabajos anteriores, obteniéndose un buen ajuste. En función de estos resultados se analiza el efecto del tamaño en los paneles. Sandwich panels of laminated gypsum and rockwool have an abundant pathology of cracking due to excessive slabs deflection. Therefore, it is necessary to progress in the simulation and prediction of behaviour under tensile and shear load of such panels, although in typical applications have no structural responsability. The behaviour of this material may be considered quasi-brittle and, based on this idea, in this work has been studied using a cohesive crack model that has been applied to other quasi-brittle materials, such as concrete, and has provided very satisfactory results. This communication presents the work carried out to study the size effect of the specimen of plasterboard and rockwool in its mechanical and resistant behaviour. The authors designed an experimental campaign under mixed mode composed by testing specimens of different sizes. Assymetrical three-point bending tests have been performed on notched specimens, geometrically similar and of different size, to obtain load-displacement and load-crack moutn opening displacement curves. To numerically simulate the mixed-mode fracture behaviour of the panels we have used a finite element model with embedded crack, based on the cohesive crack model, using as input the experimental parameters obtained in previous work, obtaining a good adjustment. Based on these results we analyze the size effect of the panels.