No 6 (2024)
- Year: 2024
- Articles: 10
- URL: https://ered.pstu.ru/index.php/mechanics/issue/view/434
- DOI: https://doi.org/10.15593/perm.mech/2024.6
EFFECT OF ELECTROLYTIC HYDROGENATION ON MECHANICAL AND ACOUSTIC PROPERTIES OF STAINLESS STEEL
Abstract
The work studies the effect of hydrogen embrittlement on AISI 420 martensitic stainless steel using the acoustic method. We evaluate an estimate of changes in mechanical characteristics and informative parameters of acoustic testing by loading the samples pre-saturated with hydrogen using the electrolytic method. It was found that the parameters of Rayleigh waves depend on the degree of hydrogen embrittlement both in the unloaded and loaded states, starting from the yield point and up to the fracture. It was shown that the dependences of the propagation velocity of Rayleigh waves can be used in problems of an early assessment of hydrogen embrittlement in a material. The results of studying the influence of hydrogen concentration on the Rayleigh wave parameters serve as the basis for constructing a comprehensive acoustic method to determine the physical and mechanical characteristics of structural materials. These data make it possible to estimate the critical concentration of hydrogen, at which 80 % of the initial mechanical properties is preserved. In addition, they allow forming a criterion for the limit state of stainless steel under static loading using parameters of acoustic waves. The Rayleigh wave propagation velocities were calculated for materials in the limit state. Dependences of coefficients of the strength degradation and hydrogen embrittlement on the acoustic properties degradation coefficient were found. These relationships can be approximated by a sigmoid Boltzmann function with high correlation coefficients.
PNRPU Mechanics Bulletin. 2024;(6):5-14
INVESTIGATION OF PROPERTIES OF THE COMPOSITE MATERIAL BASED ON THE EPOXY AND SILICA SAND POLYDISPERSE MIXTURE
Abstract
The paper studies physical and mechanical properties of an epoxy and the silica sand polydisperse mixture. As a dispersion medium we used the epoxy L with hardener EPH 161, and as a disperse material we used the silica sand mixtures of two varieties, here named fine and coarse ones (two characteristic particle sizes). The particles microphotographs are given. The particle size distribution histograms are constructed. Student's t-test was calculated to determine the difference between the mean values of two independent samples. The statistical difference between the two samples was shown to be statistically distinguishable. The mixture porosity and permeability depending on the fine sand proportion were measured. It is shown that each component’s porosity has almost the same values, and the porosity has its minimum if the fine sand content is 40 %. The dependence of permeability on the fine sand proportion is plotted. It is calculated how the effective pore radius affects the proportion of the fine particles. Mechanical tests of the polymerized samples were performed using the universal testing machine ZWIC Z-250. Non-monotonic dependence of ultimate tensile and compressive strength of the polymerized samples depending on the fine particle proportion was revealed. Pairwise linear correlation coefficients of mechanical characteristics with the medium porosity have been calculated. It is shown that there is an inverse relationship between tensile and compressive strains
PNRPU Mechanics Bulletin. 2024;(6):15-24
INFLUENCE OF THE TYPE OF KINETIC FUNCTIONS ON CONVERSION IN A THIN LAYER PLACED ON A SUBSTRATE UNDER THE CONTROL OF A MOVING HEAT SOURCE
Abstract
Advanced combined laser technologies assume a conjunction of surface treatment capabilities with material synthesis. In this case, preliminary temperature estimations based only on thermophysical models are unreasonable. It is necessary to consider chemical reactions that can take place both in solid and liquid phases. For this purpose, we present a brief review of reactions involving solids and types of kinetic functions that are used to process experiments in chemistry and chemical kinetics. A variety of kinetic functions provides a wide range of possibilities for analyzing experiments, which contributes to a deeper understanding of the processes involved. Further, on the basis of a one-dimensional reduced model with a total chemical reaction and taking into account the accumulation of the liquid phase, we show that the type of kinetic function, which reflects the reaction mechanism at the microlevel, affects both the maximum temperature and the degree of the reaction’s completion. Reduced models allow us to investigate each factor under study more effectively, which can subsequently be useful for the process optimizations. For the Al+TiO2 system, we present dependences of characteristics of the quasistationary stage on parameters reflecting the role of the synthesis conditions. In this system, the transformations can be described with two total reactions: the first gives the hardening phase, the second leads to a change in the matrix composition. The dependence of the maximum temperature on the power density of the heat source appears to be linear. After processing, the coating contains both reaction products and initial substances. The specific composition is determined by the processing conditions.
PNRPU Mechanics Bulletin. 2024;(6):25-37
LOCATION OF FORCE IMPACT BY INDICATOR PIEZOELECTRIC MDS COATING BASED ON INFORMATIVE VOLTAGE PULSES
Abstract
An electromechanical mathematical model is developed to locate a single force effect of pressing a rigid ball particle into the sensor surface of an indicator piezoelectric MDS coating with an integrated double spiral of electrodes. We consider formation of informative voltage pulses on arcs of electrodes of a double spiral in a circular perturbation zone of an indicator coating due to piezoelectric deformation between the electrodes when pressing a rigid ball particle. Shape and duration of informative voltage pulses as functions of time are found by solving the corresponding differential equation taking into account piezoelectric and geometric characteristics, value, nonuniformity, duration and distance from the epicentre of the electrode/piezoelectric/electrode disturbance arcs acting on each of the piezocells within the circular disturbance zone of the indicator coating. The number of pulses generated and recorded is proportional to the ratio of the perturbation zone radius and the sensor spiral pitch of the indicator coating. Here the perturbation zone radius is proportional to the diagnosed force. A simple, convenient and applicable solution of location is obtained, i.e. finding polar coordinates of the epicentre of the force effect, while the accuracy of the location improves when the number of the spiral turns increases, i.e. the pitch of the sensor spiral decreases, relative to the radius of the perturbation zone.
PNRPU Mechanics Bulletin. 2024;(6):38-46
NUMERICAL SIMULATION OF OSCILLATORY PROCESSES IN A MODEL CHANNEL WITH TWO PARTITIONS FILLED WITH HYDROGEN
Abstract
This article discusses the propagation of a sound wave in a channel with a rectangular cross-section with two fiberglass partitions. The time dependences of pressure and displacements for points located in different parts of the model channel are obtained. Hydrogen and air are considered as the working media. An algorithm has been developed for solving the boundary value problem of the numerical simulation of oscillatory processes in the model’s channel with a rectangular cross-section, taking into account the bidirectional interaction (2-way Fluid-Structure Interaction or 2FSI reduction) between deformable partitions and the flow of the working environment, using the ANSYS engineering analysis system. A mathematical model corresponding to the boundary value problem has been developed. Initial and boundary conditions are set and control points along the length of the model channel and on fiberglass partitions are selected. The main results of the study are presented in the form of dependences of the amplitude of displacement of the partition and pressure in the working media (air or hydrogen) on time at characteristic points of the channel. It was found that transients in air and hydrogen differ in amplitude and frequency of oscillation; the dependence of the influence of the number of partitions on fluctuations in the gaseous medium was found. The considered system can become a model for working out the conditions for the occurrence of self-oscillations for different working environments. Predicting the behavior of the pipeline structure, in a model representation, under the influence of different gases during transportation will allow us to assess the influence of the transported medium on the acoustic characteristics of the system. The study is highly important for finding areas of noise reduction during the transportation of any gas mixtures through pipelines.
PNRPU Mechanics Bulletin. 2024;(6):47-59
EFFECT OF VARIATION OF COEFFICIENT OF THERMAL EXPANSION ON THERMOMECHANICAL DEFORMATION OF CYLINDRICAL IN-PLANE AUXETICS BASED ON RE-ENTRANT CELLS
Abstract
Active development of mechanical metamaterials has currently led to the widespread application of auxetic structures in various applications with differing loading conditions. This research explores thermomechanical behaviour of novel in-plane cylindrical auxetic lattice structures by studying the correlation between their deformational characteristics and the coefficient of thermal expansion (CTE) of the material. Unlike traditional auxetic cylinders, the plane of auxeticity in the developed models is oriented perpendicular to the cylinder axis, which defines their specific behavior. To understand the behavior of auxetic lattices under combined thermal and mechanical loads, computational experiments were conducted based on the finite element method (FEM). Deformations of both rectangular and cylindrical lattice structures were investigated. The relationships between transverse deformation and CTE were obtained and compared. The influence of CTE on the structural Poisson's ratio of rectangular auxetic lattices was assessed, which is a key parameter characterizing the auxetic behavior of the structure under thermomechanical loading. The feasibility of modeling the mechanical behavior of auxetic cylinders using an orthotropic mechanical model with effective material properties was verified. The constants for defining such a material model were obtained by simulating a numerical experiment on tensile and shear testing of the rectangular auxetic lattices along coordinate axes. Using an example with artificial material properties, it was demonstrated that the predominant deformation mechanism, caused by the opposing effects of mechanical and thermal loads, can be controlled by selecting a material with an appropriate CTE. This allows for regulating the structural response to changes in temperature and mechanical load. Such results can be utilized for creating cylindrical auxetic lattice structures subjected to mechanical and thermal deformations in applications requiring controlled thermomechanical responses.
PNRPU Mechanics Bulletin. 2024;(6):60-71
On the Approximation of Derivatives in Hexahedral 8-Node Finite Elements
Abstract
In numerical solutions of elasticity and plasticity problems, finite elements with a reduced integration are often used, especially in solving dynamic problems. In this case, for 8-node 3D elements, one point of the numerical integration is used instead of 8. In this case, it is actually assumed that the strains and stresses are constant within an element. In this case, the traditional technique of constructing the stiffness matrix for an element of a standard shape in the form of a cube with a subsequent mapping of the actual finite elements of an arbitrary shape and size onto the standard one is not necessary. Instead, the stiffness matrix can be constructed directly for a finite element of an arbitrary shape. In this case, it is expressed through the coefficients of grid operators approximating the first partial derivatives of the displacement field in the finite element. The paper considers a new approach to approximating derivatives when constructing the stiffness matrix for a 3D 8-node finite element with one integration point. The theoretical basis for this approach is the further development of the class of rare mesh FEM schemes. The obtained formulas allow one to construct incompatible FEM schemes with improved properties. The paper discusses the problems of hourglass instability, shear and volume locking. A new effective approach to solving the hourglass problem is proposed. The possibility of applying new derivative approximation formulas to finite elements of a degenerate form with a number of nodes less than eight is also discussed. It is shown that they remain applicable in a standard way in this case too. The results of the study are confirmed by the presented numerical solution results of the model static elasticity problems.
PNRPU Mechanics Bulletin. 2024;(6):72-80
Features of Non-Isothermal Grain Boundary Diffusion in Ti3Al
Abstract
A two-dimensional diffusion-kinetic model is proposed to research the effect of oxygen diffusion along grain boundaries on the oxidation dynamics of the intermetallic alloy Ti3Al. The contribution of grain boundaries is evaluated by comparing the dynamics of processes in a structure with clearly defined grains and boundaries and in a material with effective properties, where the diffusion coefficient was calculated depending on the fraction of the boundary phase. Oxygen diffusion occurs in a mixed kinetic mode typical of additively fabricated nanoscale structures. The structure with an explicit consideration of grains and boundaries in the model has its symmetry. Rectangular grains are located relative to each other similar to "brickwork" so that they form triple junctions. The material with effective properties is a continuous rectangular region, in which the fraction of the boundary phase is taken into account through the diffusion coefficient. A constant oxygen source is specified on the surface. The problem is solved numerically in dimensionless variables. An implicit difference scheme of splitting by coordinates is used to solve the diffusion equation. To solve the kinetic equations, a method similar to the explicit Euler method is used with the organization of the iterative process. The results are compared for the isothermal mode and for the conditions of linear heating with subsequent cooling. The study is carried out for the initial stage of oxidation of the nanosized intermetallic alloy Ti3Al. The contribution of grain boundaries to the oxidation dynamics is estimated in the range of the fraction of the boundary phase from 0.1 to 0.5, which changes due to the variation of the grain sizes relative to the boundary width. The results obtained are in a qualitative agreement with the literature data.
PNRPU Mechanics Bulletin. 2024;(6):81–93
A NEW VERSION OF THE NONLINEAR MODEL FOR SHEAR FLOW OF TIXOTROPIC VISCOELASTIC MEDIA WITH A POWER DEPENDENCE OF RELAXATION TIME ON CURRENT STRUCTUREDNESS
Abstract
We formulate and analyze a new version of the nonlinear constitutive equation for shear flow of tixotropic viscoelastic media accounting for interaction of deformation and structure evolution (kinetics of its formation and destruction) proposed earlier. Viscosity and shear modulus appearing in the Maxwell-type model are chosen to be power functions of current structuredness instead of exponential functions in the first version of the model. The new model is also governed by an increasing material function and six positive parameters, but the set of two nonlinear differential equations for two unknown functions of time, namely, stress and structuredness (or relative cross-links density), to which the model is reduced in one-dimensional case, differ from the first version. We analyze mathematical properties of the model assuming the material function and six material parameters are arbitrary and show that the new model inherits all valuable basic properties that have been discovered analyzing the first version. We confirm the existence and uniqueness of its equilibrium point and that equilibrium stress and structuredness depend monotonically on every material parameter and on shear rate. We derive general equations for the flow and viscosity curves generated by the model and confirm that the first one increases and the second one decreases while the shear rate grows. Thus, the model describes the basic phenomena observed for the simple shear flow of shear thinning fluids, but a number of qualitative properties of the flow and viscosity curves differ from the first version. We found two applicability indicators for each version of the model which are easy to check in tests. Thus, the new version of the model is a useful additional tool aimed at simulating shear flows of tixotropic viscoelastic media.
PNRPU Mechanics Bulletin. 2024;(6):94–112
THE MARTENSITE STABILIZATION EFFECT IN TITANIUM NICKELIDE AFTER PRELIMINARY DEFORMATION BY COOLING UNDER STRESS IN AN INCOMPLETE TRANSFORMATION INTERVAL
Abstract
When using thermosensitive sensors and servos, based on shape memory alloys (SMA), it is necessary to take into account the influence of shifting factors capable the temperatures of martensitic transformations (MP). One of these factors is the martensite stabilization effect (MSE), which manifests itself as an increase in reverse MP temperatures after preliminary deformation. Previous MSE studies allowed, based on experimental data, to hypothesize that this effect is due to damage to the boundaries orientation of the martensite domains. This suggests that minimizing or eliminating MSE is possible while minimizing damage to these boundaries. One of the ways to achieve this is to cool the SMA from the austenitic state under constant load in an incomplete temperature range of transformation. Numerical experiments performed within the framework of a microstructural model that takes into account MSE made it possible to determine the dependence of the temperature shift of the reverse transformation on the value of the fraction of direct transformation achieved during cooling. The results showed that the temperature shift remains zero or negligible until a certain critical fraction of the transformation is reached.
PNRPU Mechanics Bulletin. 2024;(6):113–122