Publicações

N/A

N/A

N/A

N/A

N/A

N/A

N/A

Diferentes modelos e métodos foram desenvolvidos, tendo como base a abordagem Lagrangiana, e sãogeralmente conhecidos genericamente como modelos SPH (Smoothed Particle Hydrodynamics). Neste artigoapresentam-se as potencialidades de um método Lagrangiano na modelação de escoamentos com superfícielivre, o modelo SPHysics, através da sua aplicação a diferentes casos de interacção entre ondas e estruturascosteiras.

Wave overtopping is a violent natural event thatinvolves highly complex phenomenon such as large deformationof free surface, turbulence and eddy vortices, strong interactionbetween the wave and the structure. Models based on SmoothedParticle Hydrodynamics (SPH), that used a mesh-free technique,are an option to address wave overtopping and other phenomenainvolved on the interaction between waves and coastal structures.In the present paper, SPHysics model is validated and applied forwave propagation and wave overtopping of an impermeableseawall. Validation and convergence study is carried outconsidering several parameters such as the initial particle densityand the eXSPH parameter. Free surface elevation in several gaugesand overtopping discharge over the structure are analyzed andcompared to experimental data and other numerical results. Avery satisfactory agreement is obtained with experimentalmeasurements. Finally, the numerical model is applied formodelling wave propagation with breaking and overtopping ofan impermeable sea wall coastal defence structure, a commonstructure employed at the Portuguese coast. Numerical resultsare compared with experimental data from model scale testscarried at the National Civil Engineering Laboratory (LNEC).Good agreement is obtained for both free surface elevation andovertopping discharge over the structure.

Considera-se um conjunto de 10 anos de dados de parâmetros de agitação marítima provenientes da bóia de Faro, do Instituto Hidrográfico. Mau grado, esta colecção apresenta um razoável número de falhas ao longo desses dez anos, com variadas extensões, que se pretendem colmatar. Adoptaram-se duas formas distintas de tratamento do problema, uma mediante utilização de uma metodologia de colmatação com base em valores adjacentes dasérie existente (pequenas falhas), e outra mediante substituição dos valores em falta com observações colhidas em bóia-ondógrafo próxima (bóia-ondógrafo de Cádiz, dos Puertos del Estado, à cota -450 m Z.H.), após adequada validação dessas observações.Nesta comunicação dá-se conta da análise e processamento das séries temporais iniciais edescrevem-se os métodos utilizados para a colmatação dos registos em falha. Assim, faz-se uma descrição estatística dos dados brutos, constituída pelas séries temporais e histogramas dos parâmetros mais relevantes, procede-se à identificação dos registos em falha e descrevem-se e comparam-se métodos de colmatação de falhas das séries temporais.Para o presente conjunto de dados, conclui-se que a metodologia mais apropriada para acolmatação de falhas superiores a um dia é a que resulta da utilização de dados de uma bóia próxima.

N/A

En esta ponencia se ha presentado la solución de proyecto para la corrección delmuelle de atraque de los ferries del puerto de Vila do Porto. La estructura, que tienecomo principal objetivo disminuir la agitación y el rebase en la dársena sin ocupardemasiado espacio, está compuesta de una plataforma horizontal asentada enpilares rectangulares por encima de un talud de escollera. Se presentan los resultadosde ensayos en modelo físico 3D y los obtenidos con dos modelos numéricos RANS en lo que se refiere a presiones en la plataforma y reducción de la reflexión: el modeloCOBRAS-UC y el FLUENT.El modelado físico permitió verificar el comportamiento de la estructura: la agitación yel rebase son ahora menores en toda la dársena, o sea, entre el muelle comercial y lazona de atraque del Ferry, y las presiones máximas y la estabilidad han sidoverificadas. En cuanto al modelado numérico, el modelo FLUENT (2D), que seencuentra todavía en validación para aplicaciones de estudios de interacción olaestructura,se ha comparado con los resultados del modelo COBRAS-UC y se haverificado que presenta problemas de convergencia cuando se usa un esquema deintegración en el tiempo de segundo orden y difusión numérica de la fracción devolumen cuando se usa un esquema de primer orden (aunque resuelva los problemasde divergencia). El FLUENT 3D presenta, además, un elevado tiempo de cálculo,particularmente en 3D.

Neste trabalho estuda-se a interacção de ondas regulares com cilindros verticais fixos aofundo. Na simulação numérica deste problema tridimensional viscoso e com superfície livreutiliza-se o código FLUENT que resolve as equações de Navier-Stokes em valor médio(equações de Reynolds) e utiliza uma técnica VoF (Volume of Fluid) na identificação dasuperfície livre.Apresentam-se os parâmetros relativos à discretização e valida-se o modelo numérico porcomparação com resultados analítico-experimentais e obtidos com o código WAMITamplamente validado.As diferenças entre resultados são correlacionadas com os efeitos não lineares ligados àamplitude relativa (irrotacionalidade do escoamento) e à relação entre o diâmetro do cilindro e ocomprimento da onda.

Neste trabalho testa-se a capacidade do código numérico COBRAS-UC,baseado nas equação RANS e utilizando uma técnica VOF, na simulação dos efeitos nãolineares resultantes da interacção de dois grupos de ondas. O caso de teste consiste napropagação, em sentidos opostos (interacção fraca) de dois grupos de ondas emcondições de água profunda. Estudos analíticos, experimentais e numéricos anteriorespermitiram identificar a variação da velocidade de fase como a única consequência dainteracção não linear. A comparação de anteriores resultados com os obtidos no âmbitodeste trabalho permitem concluir sobre a precisão e potencialidades do códigoCOBRAS-UC na simulação de fenómenos não lineares exigentes do ponto numérico.

A análise da eficiência das estruturas costeiras em termos de galgamento, é classicamenterealizada através de estudos em modelo físico. No entanto, a modelação numérica podecomplementar estes estudos experimentais.A aplicação de dois tipos de modelos numéricos na simulação do galgamento de estruturacosteiras impermeáveis, cada modelo com as suas vantagens e limitações, é analisada nesteartigo: um modelo Eulariano, COBRAS-UC, e um modelo Lagrangiano, SPHysics.Os resultados de elevação de superfície livre e de caudal médio galgado obtido pelos doismodelos são comparados com dados de ensaios obtidos em modelo físico à escala 1:10, parao estudo do galgamento numa estrutura de protecção marginal impermeável.

N/A

The effect of aspect ratio on steady laminarfluid flow past an elliptical cylinder is investigatednumerically at two low Reynolds numbers, 20 and40. The two-dimensional Navier-Stokes equationsare solved using an original fully coupled resolutionmethod, without any transformation of continuityequation. The wake

N/A

Os efeitos das alterações climáticas aumentaram a necessidade de investigar a capacidade da vegetação em atenuar a força das ondas, nas zonas costeiras. Apoiando-se em modelação física, analítica e numérica da propagação de ondas sobre vegetação, este estudo teve como objetivos: i) analisar a diferença na dissipação de ondas sobre campos de vegetação rígida e flexível, ii) caracterizar o movimento de vegetação flexível, iii) estimar o coeficiente de arrasto de vegetação rígida e vegetação flexível, e iv) analisar a influência da flexibilidade (movimentoda vegetação) e da força de inércia no coeficiente de arrasto.

Na presente dissertação, desenvolveu-se uma metodologia para caracterização do movimento de um navio amarrado no interior de portos na presença de um campo de ondas. Esta metodologia consiste em utilizar vários modelos numéricos para caracterizar convenientemente: a) a propagação das ondas desde o largo até ao local de acostagem; b) a modificação das características das ondas no interior do porto; c) a ação das ondas sobre o navio; d) a resposta do navio; e) as forças resultantes no sistema de marração. Os modelos utilizados são:· SWAN: modelo não linear espectral para simular a propagação de ondas em zonas costeiras extensas;· DREAMS: modelo linear para simular a propagação de ondas em zonas costeiras abrigadas;· BOUSS-WMH: modelo não linear para simular a propagação de ondas regulares e irregulares em zonas costeiras abrigadas;· GMALHA: gerador de malhas de elementos finitos triangulares;· WAMIT®: que resolve, no domínio da frequência, os problemas de radiação e de difração na interação de um corpo flutuante livre com as ondas;· NPP: gerador de painéis planos sobre o casco do navio;· HYDRO: que transforma os resultados da interação do corpo flutuante livre com as ondas do domínio da frequência para o domínio do tempo;· BAS: que resolve as equações de movimento de um navio amarrado.O desempenho do modelo BOUSS-WMH, em situações de grande não linearidade, foi analisado através da simulação de quatro casos de teste também simulados em laboratório. O BOUSS-WMH provou ser um modelo robusto capaz de representar a deformação da superfície livre de forma adequada, embora, em zonas com forte não linearidade, tenha apresentado algumas limitações, nomeadamente na formaçãode harmónicas de ordem superior à terceira. Verificou-se que os resultados numéricos tendiam a sobrestimar os valores da altura de onda em zonas pouco profundas.Para minorar estas limitações e tornar o modelo mais robusto, fiável e adequado para ser aplicado em casos reais, foram adicionados e implementados novos termos às equações do modelo BOUSS-WMH, nomeadamente, a reflexão parcial, o atrito de fundo e a rebentação (fenómenos físicos inerentes às águas pouco profundas e zonas confinadas).A reflexão parcial foi implementada a partir da adição de um termo de viscosidade na equação de conservação da quantidade de movimento. O modelo com reflexão parcial foi aplicado a dois casos de teste, sendo um deles um caso real (Porto de Vila do Porto, ilha de Santa Maria, Açores), e os resultados foram comparados com medições em ensaios em modelo físico. Verificou-se que a implementação da reflexão parcial diminui em 44% o erro médio quadrático face à simulação com fronteiras totalmente refletoras.O atrito de fundo foi implementado a partir da adição de um termo representativo da tensão tangencial de fundo à equação de conservação da quantidade de movimento. O modelo com atrito de fundo foi aplicado a dois casos de teste, sendo um deles um caso real (Praia de Faro). Verificou-se que a introdução do atrito de fundo dissipa energia como pretendido e reduz as instabilidades numéricas.A rebentação foi implementada com a introdução de um termo de viscosidade turbulenta. Dois critérios distintos foram introduzidos para definir o início e o final da rebentação: a) FSA, baseado na aceleração da elevação da superfície livre; b) RTFN, baseado no número Froude relativo da cava. O modelo com rebentação foi aplicado a três casos de teste estudados em laboratório. A comparação dos resultados do modelo com as medições dos ensaios físicos permite dizer que, de um modo geral, o modelo numérico simula muito bem a localização da rebentação, bem como a dissipação de energia a ela associada. O critério de rebentação baseado no método de RTFN tende a sobrestimar a altura das ondas, mas apresenta índices de concordância maiores do que os obtidos com o critério FSA. A vantagem deste método reside no menor número de parâmetros a calibrar. Os resultados alcançados foram muito bons, dada a complexidade do fenómeno de rebentação das ondas e a simplicidade dos termos utilizados na sua modelação.O gerador GMALHA também foi alvo de melhoramentos nos algoritmos de geração, de forma a possibilitar a existência de ilhas no interior do domínio de cálculo.O modelo numérico WAMIT® foi desenvolvido para aplicações em mar aberto sem muitas estruturas refletoras próximas do navio. Todas as estruturas devem ser discretizadas, o que aumenta exponencialmente o esforço computacional exigido, e são obrigatoriamente totalmente reflectoras. Ora esta é uma limitação importante quando se está no interior de um porto, primeiro devido às dimensões e depois devido às fronteiras parcialmente refletoras. Assim, o modelo WAMIT® foi modificado de modo a incluir as características das ondas modeladas previamente com um dos modelos de propagação. Foram desenvolvidas duas novas rotinas HASK_D e HASK_B, respetivamente para o modelo linear DREAMS e para o modelo não linear BOUSS-WMH. Estas rotinas e as modificações ao modelo WAMIT® foram validadas com a aplicação a casos simples onde o modelo WAMIT® original também podia ser aplicado. A presença de um campo de ondas difratadas pela presença de uma estrutura fixa foi simulada com sucesso. Existem diferenças pouco significativas, principalmente nas altas frequências. Uma explicação para estas diferenças relaciona-se com o facto de que as ondas com períodos baixos apresentam maiores instabilidades numéricas à medida que diminui o número de pontos por comprimento de onda na malha de elementos finitos.Com o intuito de facilitar o tratamento da informação, que, num estudo deste género, tem um volume e complexidade consideráveis, foram desenvolvidas interfaces para cada um dos modelos e estes foram acoplados numa ferramenta integrada

Esta dissertação divide-se em duas partes principais, uma dedicada ao desenvolvimento e validação de um novo modelo numérico para a propagação de ondas fracamente não lineares com a inclusão da dissipação de energia por rebentação, e uma dedicada ao desenvolvimento de um sistema integrado de modelação da agitação marítima em zonas costeiras e portuárias, GUIOMAR, para a utilização simplificada de modelos de propagação de ondas e para o apoio à tomada de decisão em estudos de engenharia costeira (quer correntes, quer situações de emergência).Na primeira parte, as equações de Beji e Nadaoka (1997a,b), são extendidas de forma a incluírem dissipação de energia por rebentação, através da consideração das tensões de Reynolds na equação da conservação da quantidade de movimento e da dedução de uma formulação alternativa da equação não linear de declive suave dependente do tempo. O modelo desenvolvido neste trabalho é aplicado a três testes de ensaios em modelo físico, que abrangem praias de declive constante e um perfil barra-fossa. Os testes efectuados permitiram seleccionar o critério de rebentação e definir regras para a determinação dos valores para as condições críticas de início e de fim da rebentação, calibrar o valor a atribuir ao parâmetro dos termos dissipativos e estabelecer regras para a sua definição, assim como para avaliar o desempenho do modelo. Foram obtidos bons resultados, embora o empolamento da onda tenha sido subestimado. Verificou-se uma estreita relação entre o valor a utilizar para a condição crítica de início da rebentação e a declividade da onda à entrada do domínio de cálculo, e entre o valor óptimo para o parâmetro dos termos dissipativos e o número de Ursell à entrada do domínio.Na segunda parte, apresenta-se a estrutura do sistema GUIOMAR, assim como as interfaces suas constituintes. O sistema GUIOMAR permite a utilização simplificada de modelos de agitação marítima, incluindo a construção, manipulação e visualização de dados e resultados. Estes factores contribuem para a sua afirmação como uma ferramenta de apoio à tomada de decisão em estudos de engenharia costeira e portuária. Foca-se também a avaliação das variabilidades de origem geográfica nos resultados dos modelos, nomeadamente as associadas à escala e resolução dos dados batimétricos utilizados, aos métodos de interpolação mais comuns e à subjectividade inerente à definição das linhas de costa. Da análise do comportamento dos modelos de agitação marítima constatou-se a necessidade de ter em atenção no respectivo estudo, aos dados de batimetria utilizados e às interpolações a que se recorre. São delineadas sugestões/recomendações para um futuro utilizador do sistema GUIOMAR aquando da utilização de métodos de interpolação.

Portugal is one of many countries in the world that suffers from coastal erosion.Conventional ways of protecting a coastline appear to entail some disadvantages. Aninnovative and interesting way of protecting a local coastal zone by means of multifunctionalartificial reefs avoids some of them. A multi-functional artificial reef is asubmerged breakwater which, besides the helping to protect the local coastline, canoffer other purposes; in particular, it may enhance the surfing possibilities and theenvironmental value of the local area. The structure has several positive side-effects:first, it provides an unimpaired visual amenity; second, it offers tourist and economicbenefits by improving the surfing conditions; third, it can enhance the environmentalvalue of the area where it is built, and fourth, with a proper design the down drifterosion can be minimal.An optimal reef design for the Leirosa beach, located to the south of Figueira da Foz,midway along the Portuguese West Atlantic coast, is investigated. In order to achievethis optimal design several steps were conducted. A preliminary design, as a result of atheoretical study, achieved step by step, is proposed for a multi-functional artificial reefmaking use of the theory and of the state of the art. The proposed reef geometry is usedas initial design in the numerical and physical tests, which were executed to analyze thecapacity of a multi-functional reef breakwater to protect the local coastline of Leirosa,Portugal, and to increase the local surfing possibilities. Further, a physical and anumerical 2DV study have been executed to get insight in the influence of the length ofthe reef and in the submergence on the Iribarren number, which is an indication for thebreaker type. In this study, also the optimum height for the reef has been defined.Subsequently, in order to be able to investigate the three design parameters: the reefangle, the geometry of the reef (with or without a platform) and the horizontaldimensions in 2DH for the optimal geometry, the Boussinesq-type COULWAVE modelwas validated; finally, a 2DH numerical study has been conducted with COULWAVEto define the best values for these parameters.In terms of surfability, and for the conditions of the local coastline of Leirosa, thefollowing values for the main parameters were found: a reef angle of 66°, a structureheight of 3.20 m, a reef geometry composed by a delta without platform, a reefsubmergence of 1.50 m and a structure seaward slope of 1:10.The calibration of the numerical model COULWAVE has shown that it can be usedvery well to simulate the hydrodynamics around an MFAR. However, improvementsare recommended on the actual version of COULWAVE. These improvements includemainly the introduction of bottom porous structures. Even though the mean velocityfields give an indication about the sediment transport, they are not sufficient to get a deep insight about the influence of an MFAR on the morphology around the structure.For a proper design of the reef, particularly for the protection of a local coastline, amorphological study has to be executed. / Portugal é um dos muitos países no mundo que sofre os efeitos adversos da erosãocosteira. É sabido que as formas convencionais de protecção costeira apresentamalgumas desvantagens. Uma forma particularmente interessante de proteger oucontribuir para a protecção local de um troço da zona costeira, minimizando algumasdas desvantagens das estruturas tradicionais, é através da construção de recifesartificiais multi-funcionais. Um recife artificial multi-funcional é um quebra-marsubmerso construído com objectivos de ajudar a proteger localmente a zona costeira eoferecer outras valências, nomeadamente, aumentar as possibilidades locais para aprática de surf e contribuir para a valorização ambiental da região.Um recife artificial com estas características apresenta diversos aspectos positivos. Comefeito, esta estrutura: (1) não tem qualquer impacto visual negativo; (2) oferece osbenefícios económicos resultantes da atracção turística, ao melhorar as condições locaispara práticas de surf e/ou outros desportos náuticos (mergulho e pesca); (3) permiteaumentar o valor ambiental da região, e (4) com um projecto adequado permite reduzirao mínimo os efeitos adversos da erosão local ou prolongar mesmo a praia por efeito detômbolo.Neste trabalho são investigados os parâmetros característicos de um recife artificialmulti-funcional, tendo como caso de estudo um troço da zona costeira da Leirosa,situado a cerca de 15 km a sul da Figueira da Foz, Portugal. Para obter as característicasóptimas da estrutura a construir foram realizados vários estudos ao longo de diversasetapas. Como resultado do estado da arte e de um estudo teórico, realizado passo apasso, é proposta uma geometria preliminar de recife artificial multi-funcional. Estageometria preliminar do recife é em seguida utilizada como projecto inicial de diversostestes em modelo numérico e em modelo físico para analisar a sua capacidade naprotecção da praia e do cordão dunar da Leirosa e, em simultâneo, analisar aspossibilidades locais para a prática de surf.Além dos testes realizados em modelo físico, foi efectuado um estudo numéricoutilizando o modelo 2DV COBRAS-UC para obter informações sobre a influência docomprimento do recife e da submersão do mesmo no número Iribarren, sendo este umparâmetro que fornece uma boa indicação do tipo de rebentação a esperar. Neste estudofoi igualmente investigada a altura ideal para o recife. Posteriormente, procedeu-se àvalidação de um modelo numérico 2DH do tipo Boussinesq, o modelo COULWAVE,tendo como objectivo analisar a influência de três parâmetros de projecto com particulardesenvolvimento no plano horizontal: o ângulo do recife, a geometria do recife (com ousem plataforma) e as suas dimensões horizontais. Uma vez validado, foi este modelo numérico utilizado para analisar a influência daqueles parâmetros e identificar os seusvalores mais adequados para a definição de uma geometria optimizada do recife.Em termos de características para a prática de surf, e para as condições locais do litoralda Leirosa, foram identificados os seguintes valores para os principais parâmetros: umaestrutura de recife com ângulo de 66º, altura de 3.20 m e uma geometria composta porum delta sem plataforma, uma submersão do recife de 1.50 m e um declive de barlamarda estrutura de 1:10.A calibração do modelo numérico COULWAVE mostrou que este pode ser usado comsucesso para simular a hidrodinâmica em torno de um recife artificial multi-funcional.No entanto, recomendam-se melhorias na actual versão deste modelo. Tais melhoriasincluem principalmente a possibilidade de modelação com fundos porosos. Embora oscampos da velocidade média forneçam uma primeira indicação sobre o transporte desedimentos, esta informação não é por si só suficientemente profunda para obterconclusões definitivas quanto à influência de uma estrutura com estas característicassobre a morfologia local. Para uma concepção mais adequada da estrutura final dorecife, em particular para efeitos de protecção costeira, deverá ser elaborado um estudomorfodinâmico.

N/A

N/A

N/A

N/A

N/A

N/A

N/A

N/A

N/A

N/A

N/A

N/A

N/A

N/A

N/A

N/A

The project LIFEGARACHICO (LIFE20 CCA/ES/001641) proposes the creation of an effective flexible adaptation framework for the coastal unicipalities of Macaronesia, making specific local interventions to increase their resilience against extreme coastal events resulting from climate change. Social surveys were applied in Praia da Vitória, Terceira, Azores, to understand the civilians perception of coastal flooding events.

The project LIFE-GARACHICO (LIFE20 CCA/ES/001641) proposes the creation of an effective flexible adaptation framework for the coastal municipalities of Macaronesia, making specific local interventions to increase their resilience against extreme coastal events resulting from climate change. Social surveys were applied in both localities in order to understand the civilians perception of coastal flooding events

Realização de ensaios em modelo físico reduzido 3D para a avaliação da estabilidade do quebra-mar destacado do Porto Inlgês, Maio, Cabo-Verde

N/A

Adverse sea conditions can cause emergency situations associated to wave overtopping, which endanger the safety of people and goods, with negative impacts for society, the economy and the environment. Therefore, a methodology to assess the overtopping risk in port and coastal areas is essential for a proper planning and management of these areas.The Portuguese National Laboratory for Civil Engineering (LNEC) has been developing the HIDRALERTA system, an integrated decision-support tool for port and coastal management, which focus in preventing and supporting the management of emergency situations and the long-term planning of interventions in the study areas. It enables the user to calculate the risk for various port and coastal activities, starting with the characterization of sea-waves, wind field and tide levels at the study regions.HIDRALERTA was already tested at different places in Portugal: Costa da Caparica beach, in Almada, and port of Praia da Vitória, in the Azores. However, in both cases, the evaluation of wave overtopping / flooding consequences was carried out using a quite simple approach. In fact, such consequences were assessed only at a global and qualitative level, based on information provided by the responsible authorities.This paper aims at testing a quantitative methodology to evaluate overtopping consequences: the Analytic Network Process (ANP). This methodology was proposed to solve complex decision-making problems and its objective is to calculate the priorities among decision elements which will define their final weights on the consequences. The ANP approach can produce interdependencies between criteria and compute the respective weight of each criterion. Such a network model with dependence and feedback improves the priorities derived from judgements and makes prediction more accurate. So, the ANP allows comparison among clusters of elements. Furthermore, in this methodology we take into account the impacts of alternatives on the importance of criteria, by normalizing the comparison matrix. Thus, the given alternatives can influence the ranking of criteria.Application of the ANP methodology to the case study of Praia da Vitória port is presented and a consequences level map is obtained.

Apresenta-se um novo sistema de alerta para navios amarrados em portos que se encontra em desenvolvimento no Laboratório Nacional de Engenharia Civil. O sistema SWAMS_ALERTA é um sistema de previsão e alerta baseado na avaliação do risco associado ao comportamento de navios amarrados em zonas portuárias.Os riscos associados a navios amarrados decorrem de movimentos excessivos, quer no plano horizontal, quer no plano vertical. Estes movimentos podem condicionar as atividades portuárias, nomeadamente cargas e descargas, mas em casos extremos podem levar à ocorrência de situações de emergência, tais como rotura de cabos, rotura de cabeços de amarração ou até colisões com o cais. As consequências deste tipo de situações envolvem sempre grandes prejuízos materiais e por vezes humanos.O sistema SWAMS_ALERTA utiliza as medições e previsões de agitação marítima para a determinação dos seus efeitos em termos de movimentos nos seis graus de liberdade e forças nas amarras e defensas quando o navio se encontra estacionado no cais. Para tal, o sistema recorre a uma série de modelos numéricos, interligados entre si. A comparação destes valores com valores máximos admissíveis pré-estabelecidos permite a avaliação, em tempo real, de situações de emergência e a emissão de alertas dirigidos às entidades portuárias.O sistema é constituído por 4 módulos:I

N/A

N/A

N/A