• Structural Engineering and Mechanics
• /
• v.39 no.6
• /
• pp.861-875
• /
• 2011

The aim of the present article is to study the micropolar thermoelastic interactions in an infinite Kelvin-Voigt type viscoelastic thermally conducting plate. The coupled dynamic thermoelasticity and generalized theories of thermoelasticity, namely, Lord and Shulman's and Green and Lindsay's are employed by assuming the mechanical behaviour as dynamic to study the problem. The model has been simplified by using Helmholtz decomposition technique and the resulting equations have been solved by using variable separable method to obtain the secular equations in isolated mathematical conditions for homogeneous isotropic micropolar thermo-viscoelastic plate for symmetric and skew-symmetric wave modes. The dispersion curves, attenuation coefficients, amplitudes of stresses and temperature distribution for symmetric and skew-symmetric modes are computed numerically and presented graphically for a magnesium crystal.

• Journal of the Korean Geotechnical Society
• /
• v.38 no.9
• /
• pp.45-52
• /
• 2022

Because discontinuity in the rock mass and contact of soil-structure interaction exhibits coupled thermal-hydromechanical (THM) behavior, it is necessary to develop an interface element based on the full governing equations. In this study, we derive force equilibrium, fluid continuity, and energy equilibrium equations for the interface element. Additionally, we present a stiffness matrix of the elastoplastic mechanical model for the interface element. The developed interface element uses six nodes for displacement and four nodes for water pressure and temperature in a two-dimensional analysis. The fully coupled THM analysis for fluid injection into a fault can model the complicated evolution of injection pressure due to decreasing effective stress in the fault and thermal contraction of the surrounding rock mass. However, the result of hydromechanical analysis ignoring thermal phenomena overestimates hydromechanical variables.

• Tribology and Lubricants
• /
• v.40 no.1
• /
• pp.17-23
• /
• 2024

This study analyzes the thermal effects on the performance of an air foil thrust bearing (AFTB) using COMSOL Multiphysics to approximate actual bearing behavior under real conditions. An AFTB is a sliding-thrust bearing that uses air as a lubricant to support the axial load. The AFTB consists of top and bump foils and supports the rotating disk through the hydrodynamic pressure generated by the wedge effect from the inclined surface of the top foil and the elastic deformation of the bump foils, similar to a spring. The use of air as a lubricant has some advantages such as low friction loss and less heat generation, enabling air bearings to be widely used in high-speed rotating systems. However, even in AFTB, the effects of energy loss due to viscosity at high speeds, interface frictional heat, and thermal deformation of the foil caused by temperature increase cannot be ignored. Foil deformation derived from the thermal effect influences the minimum decay in film thickness and enhances the film pressure. For these reasons, performance analyses of isothermal AFTBs have shown few discrepancies with real bearing behavior. To account for this phenomenon, a thermal-fluid-structure analysis is conducted to describe the combined mechanics. Results show that the load capacity under the thermal effect is slightly higher than that obtained from isothermal analysis. In addition, the push and pull effects on the top foil and bump foil-free edges can be simulated. The differences between the isothermal and thermal behaviors are discussed.

• Coupled systems mechanics
• /
• v.7 no.4
• /
• pp.373-393
• /
• 2018

This paper is concerned with the wave propagation behavior of rotating functionally graded temperature-dependent nanoscale beams subjected to thermal loading based on nonlocal strain gradient stress field. Uniform, linear and nonlinear temperature distributions across the thickness are investigated. Thermo-elastic properties of FG beam change gradually according to the Mori-Tanaka distribution model in the spatial coordinate. The nanobeam is modeled via a higher-order shear deformable refined beam theory which has a trigonometric shear stress function. The governing equations are derived by Hamilton's principle as a function of axial force due to centrifugal stiffening and displacement. By applying an analytical solution and solving an eigenvalue problem, the dispersion relations of rotating FG nanobeam are obtained. Numerical results illustrate that various parameters including temperature change, angular velocity, nonlocality parameter, wave number and gradient index have significant effect on the wave dispersion characteristics of the understudy nanobeam. The outcome of this study can provide beneficial information for the next generation researches and exact design of nano-machines including nanoscale molecular bearings and nanogears, etc.

• Tunnel and Underground Space
• /
• v.31 no.6
• /
• pp.593-609
• /
• 2021

We proposed a numerical method for the thermo-mechanical behavior of rock fracture using a grain-based distinct element model (GBDEM) and simulated thermally induced fracture slip. The present study is the benchmark simulation performed as part of DECOVALEX-2023 Task G, which aims to develop a numerical method to estimate the coupled thermo-hydro-mechanical processes within the crystalline rock fracture network. We represented the rock sample as an assembly of Voronoi grains and calculated the interaction of the grains (blocks) and their interfaces (contacts) using a distinct element code, 3DEC. Based on an equivalent continuum approach, the micro-parameters of grains and contacts were determined to reproduce rock as an elastic material. Then, the behavior of the fracture embedded in the rock was characterized by the contacts with Coulomb shear strength and tensile strength. In the benchmark simulation, we quantitatively examined the effects of the boundary stress and thermal stress due to heat conduction on fracture behavior, focusing on the mechanism of thermally induced fracture slip. The simulation results showed that the developed numerical model reasonably reproduced the thermal expansion and thermal stress increment, the fracture stress and displacement and the effect of boundary condition. We expect the numerical model to be enhanced by continuing collaboration and interaction with other research teams of DECOVALEX-2023 Task G and validated in further study experiments.

• Advances in nano research
• /
• v.12 no.1
• /
• pp.81-99
• /
• 2022

The aim of current work is to evaluate thermo-electrical characteristics of graphene nanoplatelets Reinforced Composite (GNPRC) coupled with magneto-electro-elastic (MEE) face sheet. In this regard, a cylindrical smart nanocomposite made of GNPRC with an external MEE layer is considered. The bonding between the layers are assumed to be perfect. Because of the layer nature of the structure, the material characteristics of the whole structure is regarded as graded. Both mechanical and thermal boundary conditions are applied to this structure. The main objective of this work is to determine critical temperature and critical voltage as a function of thermal condition, support type, GNP weight fraction, and MEE thickness. The governing equation of the multilayer nanocomposites cylindrical shell is derived. The generalized differential quadrature method (GDQM) is employed to numerically solve the differential equations. This method is integrated with Deep Learning Network (DNN) with ADADELTA optimizer to determine the critical conditions of the current sandwich structure. This the first time that effects of several conditions including surrounding temperature, MEE layer thickness, and pattern of the layers of the GNPRC is investigated on two main parameters critical temperature and critical voltage of the nanostructure. Furthermore, Maxwell equation is derived for modeling of the MEE. The outcome reveals that MEE layer, temperature change, GNP weight function, and GNP distribution patterns GNP weight function have significant influence on the critical temperature and voltage of cylindrical shell made from GNP nanocomposites core with MEE face sheet on outer of the shell.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.20 no.1
• /
• pp.95-105
• /
• 2012

The selective catalytic reduction (SCR) system is a highly-effective aftertreatment device for NOx reduction of diesel engines. Generally, the ammonia ($NH_3$) was generated from reaction mechanism of SCR in the SCR system using the liquid urea as the reluctant. Therefore, the precise urea dosing control is a very important key for NOx and $NH_3$ slip reduction in the SCR system. This paper investigated NOx and $NH_3$ emission characteristics of urea-SCR dosing system based on model-based control algorithm in order to reduce NOx. In the map-based control algorithm, target amount of urea solution was determined by mass flow rate of exhaust gas obtained from engine rpm, torque and $O_2$ for feed-back control NOx concentration should be measured by NOx sensor. Moreover, this algorithm can not estimate $NH_3$ absorbed on the catalyst. Hence, the urea injection can be too rich or too lean. In this study, the model-based control algorithm was developed and evaluated on the numerical model describing physical and chemical phenomena in SCR system. One channel thermo-fluid model coupled with finely tuned chemical reaction model was applied to this control algorithm. The vehicle test was carried out by using map-based and model-based control algorithms in the NEDC mode in order to evaluate the performance of the model based control algorithm.

• Tunnel and Underground Space
• /
• v.30 no.4
• /
• pp.359-381
• /
• 2020

Within the framework of DECOVALEX-2019 Task D, full-scale engineered barriers experiment (FEBEX) at Grimsel Test Site was numerically simulated to investigate an applicability of implemented Barcelona basic model (BBM) into TOUGH2-MP/FLAC3D simulator, which was developed for the prediction of the coupled thermo-hydro-mechanical behavior of bentonite buffer. And the calculated heater power, temperature, relative humidity, total stress, saturation, water content and dry density were compared with in situ data monitored in the various sections. In general, the calculated heater power and temperature provided a fairly good agreement with experimental observations, however, the difference between power of heater #1 and that of heater #2 could not captured in the numerical analysis. It is necessary to consider lamprophyre with low thermal conductivity around heater #1 and non-simplified installation progresses of bentonite blocks in the tunnel for better modeling results. The evolutions and distributions of relative humidity were well reproduced, but hydraulic model needs to be modified because the re-saturation process was relatively fast near the heaters. In case of stress evolutions due to the thermal and hydraulic expansions, the computed stress was in good agreement with the data. But, the stress is slightly higher than the measured in situ data at the early stage of the operation, because gap between rock mass and bentonite blocks have not been considered in the numerical simulations. The calculated distribution of saturation, water content, and dry density along the radial distance showed good agreement with the observations after the first and final dismantling. The calculated dry density near the center of the FEBEX tunnel and heaters were overestimated compared with the observations. As a result, the saturation and water content were underestimated with the measurements. Therefore, numerical model of permeability is needed to modify for the production of better numerical results. It will be possible to produce the better analysis results and more realistically predict the coupled THM behavior in the bentonite blocks by performing the additional studies and modifying the numerical model based on the results of this study.

• Tunnel and Underground Space
• /
• v.31 no.4
• /
• pp.270-288
• /
• 2021

We proposed a numerical method to simulate the hydro-mechanical behavior of rock fracture using a grain-based distinct element model (GBDEM) in the paper. As a part of DECOVALEX-2023 Task G, we verified the method via benchmarks with analytical solutions. DECOVALEX-2023 Task G aims to develop a numerical method to estimate the coupled thermo-hydro-mechanical processes within the crystalline rock fracture network. We represented the rock sample as a group of tetrahedral grains and calculated the interaction of the grains and their interfaces using 3DEC. The micro-parameters of the grains and interfaces were determined by a new methodology based on an equivalent continuum approach. In benchmark modeling, a single fracture embedded in the rock was examined for the effects of fracture inclination and roughness, the boundary stress condition and the applied pressure. The simulation results showed that the developed numerical model reasonably reproduced the fracture slip induced by boundary stress condition, the fracture opening induced by fluid injection, the stress distribution variation with fracture inclination, and the fracture roughness effect. In addition, the fracture displacements associated with the opening and slip showed good agreement with the analytical solutions. We expect the numerical model to be enhanced by continuing collaboration and interaction with other research teams of DECOVALEX-2023 Task G and validated in further study experiments.

• Proceedings of the Korean Powder Metallurgy Institute Conference
• /
• 1997.10a
• /
• pp.7-7
• /
• 1997

During sintering of very porous green bodies, as obtained by compaction of hard powders - such as tungsten carbide or ceramics - or by injection moulding, important shrinkage occurs. Due to heterogeneous green density field, gravity effects, friction on the support, thermal gradients, etc., this shrinkage is often non-uniform, which' may induce significant shape changes. As the ratio of compact dimension to powder size is very high, the mechanics of continuum is relevant to model such phenomena. Thus numerical techniques, such as the finite element method can be used to simulate the sintering process and predict the final shape of the sintered part. Such type of simulation has much been developed in the last decade firstly for hot isostatic pressing and next for die compaction. Finite element modelling has been recently applied to free sintering. The simulation of sintering should be based on constitutive equations describing the thermo-mechanical behaviour of the material under any state of stress and any temperature which may arise within the sintering body. These equations can be drawn either from experimental data or from micromechanical models. The experiments usually consist in free sintering and sinter-forging tests. Indeed applying more complex loading conditions at high temperature under controlled atmosphere is delicate. Micromechanical models describe the constitutive behaviour of aggregates of spheres from the deformation of two-sphere contact either by viscous flow or grain boundary diffusion. Such models are not able to describe complex microstructure and mechanisms as observed in real materials but they can give some basic information on the formulation of constitutive equations. Practically both experimental and theoretical approaches can be coupled to identify the constitutive equations. Such procedure has been performed for modelling the sintering of compacts obtained by die pressing of a mixture of tungsten carbide and cobalt powders. The constitutive behaviour of this material during sintering has been described by a linear viscous constitutive model, whose functions have been fitted from results of free sintering and sinter-forging experiments. This model has next been introduced in ABAQUS finite element code to simulate the sintering of heterogeneous green compacts of various geometries at constant temperature. Examples of simulations are shown and compared with experiments.