Case-based reasoning as a tool to improve the usability of numerical models

There may be advantages to be gained by combining Case-Based Reasoning (CBR) techniques with numerical models. In this paper we consider how CBR can be used as a flexible query engine to improve the usability of numerical models. Particularly they can help to solve inverse and mixed problems, and to solve constraint problems. We discuss this idea with reference to the illustrative example of a pneumatic conveyor. We describe a model of the problem of particle degradation in such a conveyor, and the problems faced by design engineers. The solution of these problems requires a system that allows iterative sharing of control between user, CBR system, and numerical model. This multi-initiative interaction is illustrated for the pneumatic conveyor by means of Unified Modeling Language (UML) collaboration and sequence diagrams. We show approaches to the solution of these problems via a CBR tool.

Thermal analysis during bone drilling using rigid polyurethane foams: numerical and experimental methodologies

Osteotomy or bone cutting is a common procedure in orthopaedic surgery, mainly in the treatment of fractures and reconstructive surgery. However, the excessive heat produced during the bone drilling process is a problem that counters the benefits of this type of surgery, because it can result in thermal osteonecrosis, bone reabsorption and damage the osseointegration of implants. The analysis of different drilling parameters and materials can allow to decrease the temperature during the bone drilling process and contribute to a greater success of this kind of surgical interventions. The main goal of this study was to build a numerical three-dimensional model to simulate the drilling process considering the type of bone, the influence of cooling and the bone density of the different composite materials with similar mechanical properties to the human bone and generally used in experimental biomechanics. The numerical methodology was coupled with an experimental methodology. The use of cooling proved to be essential to decrease the material damage during the drilling process. It was concluded that the materials with less porosity and density present less damage in drilling process. The developed numerical model proved to be a great tool in this kind of analysis. © 2016, The Brazilian Society of Mechanical Sciences and Engineering.

Three dimensional numerical simulation of thermohaline currents of the Persian Gulf

Observational data and a three dimensional numerical model (POM) are used to investigate the Persian Gulf outflow structure and its spreading pathway into the Oman Sea. The model is based on orthogonal curvilinear coordinate system in horizontal and train following coordinate (sigma coordinate) system in vertical. In the simulation, the horizontal diffusivity coefficients are calculated form Smogorinsky diffusivity formula and the eddy vertical diffusivities are obtained from a second turbulence closure model (namely Mellor-Yamada level 2.5 model of turbulence). The modeling area includes the east of the Persian Gulf, the Oman Sea and a part of the north-east of the Indian Ocean. In the model, the horizontal grid spacing was assumed to be about 3.5 km and the number of vertical levels was set to 32. The simulations show that the mean salinity of the PG outflow does not change substantially during the year and is about 39 psu, while its temperature exhibits seasonal variations. These lead to variations in outflow density in a way that is has its maximum density in late winter (March) and its minimum in mid-summer (August). At the entrance to the Oman Sea, the PG outflow turns to the right due to Coriolis Effect and falls down on the continental slope until it gains its equilibrium depth. The highest density of the outflow during March causes it to sink more into the deeper depths in contrast to that of August which the density is the lowest one. Hence, the neutral buoyancy depths of the outflow are about 500 m and 250 m for March and August respectively. Then, the outflow spreads in its equilibrium depths in the Oman Sea in vicinity of western and southern boundaries until it approach the Ras al Hamra Cape where the water depth suddenly begins to increase. Therefore, during March, the outflow that is deeper and wider relative to August, is more affected by the steep slope topography and as a result of vortex stretching mechanism and conservation of potential vorticity it separates from the lateral boundaries and finally forms an anti-cyclonic eddy in the Oman Sea. But during August the outflow moves as before in vicinity of lateral boundaries. In addition, the interaction of the PG outflow with tide in the Strait of Hormuz leads to intermittency in outflow movement into the Oman Sea and it could be the major reason for generations of Peddy (Peddies) in the Oman Sea.

Numerical study of circulation in the Caspian Sea and its effects on the southern coastal currents

In this study, numerical simulation of the Caspian Sea circulation was performed using COHERENS three-dimensional numerical model and field data. The COHERENS three-dimensional model and FVCOM were performed under the effect of the wind driven force, and then the simulation results obtained were compared. Simulation modeling was performed at the Caspian Sea. Its horizontal grid size is approximately equal to 5 Km and 30 sigma levels were considered. The numerical simulation results indicate that the winds' driven-forces and temperature gradient are the most important driving force factors of the Caspian circulation pattern. One of the effects of wind-driven currents was the upwelling phenomenon that was formed in the eastern shores of the Caspian Sea in the summer. The simulation results also indicate that this phenomenon occurred at a depth less than 40 meters, and the vertical velocity in July and August was 10 meters and 7 meters respectively. During the upwelling phenomenon period the temperatures on the east coast compared to the west coast were about 5°C lower. In autumn and winter, the warm waters moved from the south east coast to the north and the cold waters moved from the west coast of the central Caspian toward the south. In the subsurface and deep layers, these movements were much more structured and caused strengthening of the anti-clockwise circulation in the area, especially in the central area of Caspian. The obtained results of the two models COHERENS and FVCOM performed under wind driven-force show a high coordination of the two models, and so the wind current circulation pattern for both models is almost identical.

Numerical analysis of wooden slabs with perforations under fire conditions

Wood is considered an ideal solution for floors and roofs building construction, due the mechanical and thermal properties, associated with acoustic conditions. These constructions have good sound absorption, heat insulation and relevant architectonic characteristics. They are used in many civil applications: concert and conference halls, auditoriums, ceilings, walls… However, the high vulnerability of wooden elements submitted to fire conditions requires the evaluation of its structural behaviour with accuracy. The main objective of this work is to present a numerical model to assess the fire resistance of wooden cellular slabs with different perforations. Also the thermal behaviour of the wooden slabs will be compared considering different material insulation, with different sizes, inside the cavities. A transient thermal analysis with nonlinear material behaviour will be solved using ANSYS© program. This study allows to verify the fire resistance, the temperature evolution and the char-layer, throughout a wooden cellular slab with perforations and considering the insulation effect inside the cavities.

Effect of different feed-rate in bone drilling: experimental and numerical study

The behaviour of bone tissue during drilling has been subject of recent studies due to its great importance. Because of thermal nature of the bone drilling, high temperatures and thermal mechanical stresses are developed during drilling that affect the process quality. However, there is still a lack information with regard to the distribution of mechanical and thermal stresses during bone drilling. The present paper describes a sequentially coupled thermal-stress analysis to assess the mechanical and thermal stress distribution during bone drilling. A three-dimensional thermo-mechanical model was developed using the ANSYS/LSDYNA finite element code under different drilling conditions. The model incorporates the dynamic characteristics of drilling process, as well as the thermo-mechanical properties of the involved materials. Experimental tests with polyurethane foam materials were also carried out. It was concluded that the use of higher feed-rates lead to a decrease of normal stresses and strains in the foam materials. The experimental and numerical results were compared and showed good agreement. The proposed numerical model could be used to predict the better drilling parameters and minimize the bone injuries.

Development of a 1-D Catalyzed Diesel Particulate Filter Model for Simulation of the Performance and the Oxidation of Particulate Matter and Nitrogen Oxides Using Passive Oxidation and Active Regeneration Engine Experimental Data

Back-pressure on a diesel engine equipped with an aftertreatment system is a function of the pressure drop across the individual components of the aftertreatment system, typically, a diesel oxidation catalyst (DOC), catalyzed particulate filter (CPF) and selective catalytic reduction (SCR) catalyst. Pressure drop across the CPF is a function of the mass flow rate and the temperature of the exhaust flowing through it as well as the mass of particulate matter (PM) retained in the substrate wall and the cake layer that forms on the substrate wall. Therefore, in order to control the back-pressure on the engine at low levels and to minimize the fuel consumption, it is important to control the PM mass retained in the CPF. Chemical reactions involving the oxidation of PM under passive oxidation and active regeneration conditions can be utilized with computer numerical models in the engine control unit (ECU) to control the pressure drop across the CPF. Hence, understanding and predicting the filtration and oxidation of PM in the CPF and the effect of these processes on the pressure drop across the CPF are necessary for developing control strategies for the aftertreatment system to reduce back-pressure on the engine and in turn fuel consumption particularly from active regeneration. Numerical modeling of CPF's has been proven to reduce development time and the cost of aftertreatment systems used in production as well as to facilitate understanding of the internal processes occurring during different operating conditions that the particulate filter is subjected to. A numerical model of the CPF was developed in this research work which was calibrated to data from passive oxidation and active regeneration experiments in order to determine the kinetic parameters for oxidation of PM and nitrogen oxides along with the model filtration parameters. The research results include the comparison between the model and the experimental data for pressure drop, PM mass retained, filtration efficiencies, CPF outlet gas temperatures and species (NO2) concentrations out of the CPF. Comparisons of PM oxidation reaction rates obtained from the model calibration to the data from the experiments for ULSD, 10 and 20% biodiesel-blended fuels are presented.

Experimental and numerical investigation on short and long term performance of macro-synthetic fibre reinforced concrete materials

Fibre Reinforced Concretes are innovative composite materials whose applications are growing considerably nowadays. Being composite materials, their performance depends on the mechanical properties of both components, fibre and matrix and, above all, on the interface. The variables to account for the mechanical characterization of the material, could be proper of the material itself, i.e. fibre and concrete type, or external factors, i.e. environmental conditions. The first part of the research presented is focused on the experimental and numerical characterization of the interface properties and short term response of fibre reinforced concretes with macro-synthetic fibers. The experimental database produced represents the starting point for numerical models calibration and validation with two principal purposes: the calibration of a local constitutive law and calibration and validation of a model predictive of the whole material response. In the perspective of the design of sustainable admixtures, the optimization of the matrix of cement-based fibre reinforced composites is realized with partial substitution of the cement amount. In the second part of the research, the effect of time dependent phenomena on MSFRCs response is studied. An extended experimental campaign of creep tests is performed analysing the effect of time and temperature variations in different loading conditions. On the results achieved, a numerical model able to account for the viscoelastic nature of both concrete and reinforcement, together with the environmental conditions, is calibrated with the LDPM theory. Different type of regression models are also elaborated correlating the mechanical properties investigated, bond strength and residual flexural behaviour, regarding the short term analysis and creep coefficient on time, for the time dependent behaviour, with the variable investigated. The experimental studies carried out emphasize the several aspects influencing the material mechanical performance allowing also the identification of those properties that the numerical approach should consider in order to be reliable.

Numerical prediction of the lateral displacement capacity of flexural reinforced concrete shear walls

Previous earthquakes showed that shear wall damage could lead to catastrophic failures of the reinforced concrete building. The lateral load capacity of shear walls needs to be estimated to minimize associated losses during catastrophic events; hence it is necessary to develop and validate reliable and stable numerical methods able to converge to reasonable estimations with minimum computational effort. The beam-column 1-D line element with fiber-type cross-section model is a practical option that yields results in agreement with experimental data. However, shortcomings of using this model to predict the local damage response may come from the fact that the model requires fine calibration of material properties to overcome regularization and size effects. To reduce the mesh-dependency of the numerical model, a regularization method based on the concept of post-yield energy is applied in this work to both the concrete and the steel material constitutive laws to predict the nonlinear cyclic response and failure mechanism of concrete shear walls. Different categories of wall specimens known to produce a different response under in plane cyclic loading for their varied geometric and detailing characteristics are considered in this study, namely: 1) scaled wall specimens designed according to the European seismic design code and 2) unique full-scale wall specimens detailed according to the U.S. design code to develop a ductile behavior under cyclic loading. To test the boundaries of application of the proposed method, two full-scale walls with a mixed shear-flexure response and different values of applied axial load are also considered. The results of this study show that the use of regularized constitutive models considerably enhances the response predictions capabilities of the model with regards to global force-drift response and failure mode. The simulations presented in this thesis demonstrate the proposed model to be a valuable tool for researchers and engineers.

Root cause analysis applied to a finite element model's refinement of a negative stiffness structure

Negative Stiffness Structures are mechanical systems that require a decrease in the applied force to generate an increase in displacement. They are structures that possess special characteristics such as snap-through and bi-stability. All of these features make them particularly suitable for different applications, such as shock-absorption, vibration isolation and damping. From this point of view, they have risen awareness of their characteristics and, in order to match them to the application needed, a numerical simulation is of great interest. In this regard, this thesis is a continuation of previous studies in a circular negative stiffness structure and aims at refine the numerical model by presenting a new solution. To that end, an investigation procedure is needed. Amongst all of the methods available, root cause analysis was the chosen one to perform the investigation since it provides a clear view of the problem under analysis and a categorization of all the causes behind it. As a result of the cause-effect analysis, the main causes that have influence on the numerical results were obtained. Once all of the causes were listed, solutions to them were proposed and it led to a new numerical model. The numerical model proposed was of nonlinear type of analysis with hexagonal elements and a hyperelastic material model. The results were analyzed through force-displacement curves, allowing for the visualization of the structure’s energy recovery. When compared to the results obtained from the experimental part, it is evident that the trend is similar and the negative stiffness behaviour is present.

Numerical study on viscoelastic behavior of polypropylene fiber reinforced concrete subjected to temperature variations using LDPM theory

This thesis is focused on the viscoelastic behavior of macro-synthetic fiber-reinforced concrete (MSFRC) with polypropylene studied numerically when subjected to temperature variations (-30 oC to +60 oC). LDPM (lattice discrete particle model), a meso-scale model for heterogeneous composites, is used. To reproduce the MSFRC structural behavior, an extended version of LDPM that includes fiber effects through fiber-concrete interface micromechanics, called LDPM-F, is applied. Model calibration is performed based on three-point bending, cube, and cylinder test for plain concrete and MSFRC. This is followed by a comprehensive literature study on the variation of mechanical properties with temperature for individual fibers and plain concrete. This literature study and past experimental test results constitute inputs for final numerical simulations. The numerical response of MSFRC three-point bending test is replicated and compared with the previously conducted experimental test results; finally, the conclusions were drawn. LDPM numerical model is successfully calibrated using experimental responses on plain concrete. Fiber-concrete interface micro-mechanical parameters are subsequently fixed and LDPM-F models are calibrated based on MSFRC three-point bending test at room temperature. Number of fibers contributing crack bridging mechanism is computed and found to be in good agreement with experimental counts. Temperature variations model for individual constituents of MSFRC, fibers and plain concrete, are implemented in LDPM-F. The model is validated for MSFRC three-point bending stress-CMOD (crack mouth opening) response reproduced at -30 oC, -15 oC, 0 oC, +20 oC, +40 oC and +60 oC. It is found that the model can well describe the temperature variation behavior of MSFRC. At positive temperatures, simulated responses are in good agreement. Slight disagreement in negative regimes suggests an in-depth study on fiber-matrix interface bond behavior with varying temperatures.

Migração lateral da desembocadura do Rio Itapocú, SC, Brasil: evolução morfológica e condicionantes físicas

Desembocaduras são ambientes bastante dinâmicos e sujeitos à complexa interação entre fatores estabilizadores e desestabilizadores. Dependendo dessa interação, desembocaduras podem apresentar a tendência de migração ao longo de barreiras arenosas. Um dos mecanismos mais eficientes de transporte de sedimento paralelo à costa, e consequentemente migração de canais, são as correntes longitudinais geradas pelas ondas se aproximando obliquamente à costa. A motivação do presente trabalho é entender o comportamento morfodinâmico do sistema de desembocadura do rio Itapocú, localizado no centro-norte de Santa Catarina (SC), frente aos processos forçantes que atuam na sua migração ao longo da linha de costa. A morfologia dos pontais arenosos foi obtida a partir de levantamentos morfológicos com o uso de DGPS. Para analisar a refração de ondas foi utilizado o modelo numérico MIKE 21 SW, sendo considerados como condições de contorno os dados de ondas referentes ao ano de 2002 e os dados de ondas previstos referentes ao período de coleta. Os dados de saída do modelo foram utilizados para estimar a deriva litorânea potencial na região. Os resultados morfológicos obtidos demonstraram uma migração da desembocadura para o norte durante o período analisado, sendo mais intenso durante o inverno e o verão. Ondas incidentes do quadrante sul sofreram mais o fenômeno da refração e as ondas de leste apresentaram menor variação angular ao se aproximarem à costa. A deriva litorânea potencial anual para os dados de ondas de 2002 apresentou sentido norte-sul, com inversão de sentido durante o outono. Utilizando os dados de ondas previstas para o período dos levantamentos, a deriva litorânea potencial estimada apresentou sentido sul-norte, concordando com a migração observada. Na região próxima a desembocadura, nos pontais arenosos, a deriva potencial apresentou direção para o norte durante todas as estações. Os dados de descarga fluvial não apresentaram influência na migração do canal, porém apresentaram uma relação com a largura do mesmo sazonalmente.Os dados de morfologia juntamente com os dados de deriva litorânea referentes às ondas de 2004/2005 mostraram claramente a migração do canal para o norte sendo a deriva a principal contribuinte para a migração da desembocadura.

Evaluation of passive confinement in CFT columns

This paper presents the experimental results of 32 axially loaded concrete-filled steel tubular columns (CFT). The load was introduced only on the concrete core by means of two high strength steel cylinders placed at the column ends to evaluate the passive confinement provided by the steel tube. The columns were filled with structural concretes with compressive strengths of 30, 60, 80 and 100 MPa. The outer diameter (D) of the column was 114.3 mm, and the length/diameter (L/D) ratios considered were 3, 5, 7 and 10. The wall thicknesses of the tubes (t) were 3.35 mm and 6.0 mm, resulting in diameter/thickness (D/t) ratios of 34 and 19, respectively. The force vs. axial strain curves obtained from the tests showed, in general, a good post-peak behavior of the CFT columns, even for those columns filled with high strength concrete. Three analytical models of confinement for short concrete-filled columns found in the literature were used to predict the axial capacity of the columns tested. To apply these models to slender columns, a correction factor was introduced to penalize the calculated results, giving good agreement with the experimental values. Additional results of 63 CFT columns tested by other researchers were also compared to the predictions of the modified analytical models and presented satisfactory results. (C) 2009 Elsevier Ltd. All rights reserved.

Metodologia para análise da confiabilidade estrutural de escavações em rocha.

A methodology for rock-excavation structural-reliability analysis that uses Distinct Element Method numerical models is presented. The methodology solves the problem of the conventional numerical models that supply only punctual results and use fixed input parameters, without considering its statistical errors. The analysis of rock-excavation stability must consider uncertainties from geological variability, from uncertainty in the choice of mechanical behaviour hypothesis, and from uncertainties in parameters adopted in numerical model construction. These uncertainties can be analyzed in simple deterministic models, but a new methodology was developed for numerical models with results of several natures. The methodology is based on Monte Carlo simulations and uses principles of Paraconsistent Logic. It will be presented in the analysis of a final slope of a large-dimensioned surface mine.

Root Water Extraction under Combined Water and Osmotic Stress

Using a numerical implicit model for root water extraction by a single root in a symmetric radial flow problem, based on the Richards equation and the combined convection-dispersion equation, we investigated some aspects of the response of root water uptake to combined water and osmotic stress. The model implicitly incorporates the effect of simultaneous pressure head and osmotic head on root water uptake, and does not require additional assumptions (additive or multiplicative) to derive the combined effect of water and salt stress. Simulation results showed that relative transpiration equals relative matric flux potential, which is defined as the matric flux potential calculated with an osmotic pressure head-dependent lower bound of integration, divided by the matric flux potential at the onset of limiting hydraulic conditions. In the falling rate phase, the osmotic head near the root surface was shown to increase in time due to decreasing root water extraction rates, causing a more gradual decline of relative transpiration than with water stress alone. Results furthermore show that osmotic stress effects on uptake depend on pressure head or water content, allowing a refinement of the approach in which fixed reduction factors based on the electrical conductivity of the saturated soil solution extract are used. One of the consequences is that osmotic stress is predicted to occur in situations not predicted by the saturation extract analysis approach. It is also shown that this way of combining salinity and water as stressors yields results that are different from a purely multiplicative approach. An analytical steady state solution is presented to calculate the solute content at the root surface, and compared with the outputs of the numerical model. Using the analytical solution, a method has been developed to estimate relative transpiration as a function of system parameters, which are often already used in vadose zone models: potential transpiration rate, root length density, minimum root surface pressure head, and soil theta-h and K-h functions.