This study was carried out on the background of Sutong Bridge project based on fracture mechanics, aiming at analyzing the growth mechanism of fatigue cracks of a bridge under the load of vehicles. Stress intensity factor (SIF) can be calculated by various methods. Three steel plates with different kinds of cracks were taken as the samples in this study. With the combination of finite element analysis software ABAQUS and the J integral method, SIF values of the samples were calculated. After that, the extended finite element method in the simulation of fatigue crack growth was introduced, and the simulation of crack growth paths under different external loads was analyzed. At last, we took a partial model from the Sutong Bridge and supposed its two dangerous parts already had fine cracks; then simulative vehicle load was added onto the U-rib to predict crack growth paths using the extended finite element method.
Experimental design and computational model for predicting debonding initiation and propagation are of interest of scientists and engineers. The design and model are expected to explain the phenomenon for a wide range of loading rates. In this work, a method to measure and quantify debonding strength at various loading rates is proposed. The method is experimentally verified using data obtained from a static test and a pulse-type dynamic test. The proposed method involves the cohesive zone model, which can uniquely be characterized with a few parameters. Since those parameters are difficult to be measured directly, indirect inference is deployed where the parameters are inferred by minimizing discrepancy of mechanical response of a numerical model and that of the experiments. The main finding suggests that the design is easy to be used for the debonding characterization and the numerical model can accurately predict the debonding for the both loading cases. The cohesive strength of the stress-wave case is significantly higher than that of the static case; meanwhile, the cohesive energy is twice larger.
The present paper is devoted to the discussion and review of the non-destructive testing methods mainly based on vibration and wave propagation. In the first part, the experimental methods of actuating and analyzing the signal (vibration) are discussed. The piezoelectric elements, fiber optic sensors and Laser Scanning Doppler Vibrometer (SLDV) method are described. Effective detecting of the flaws needs very accurate theoretical models. Thus, the numerical methods, e.g. finite element, spectral element method and numerical models of the flaws in isotropic and composite materials are presented. Moreover, the detection of the damage in structures, which are subjected to cyclic or static loads, is based on the analyzing of the change in natural frequency of the whole structure, the change of internal impedance of the material and the change in guided waves propagating through the investigated structure. All these cases are characterized in detail. At the end of this paper, several applications of the structural health monitoring systems in machine design and operation are presented.
The paper presents a model of a rapping system of an electrostatic precipitator. The rapping system consists of a set of collecting electrodes hanging on a suspension bar and braced together in a brushing bar. The suspension and brushing bars are modeled using the rigid finite element method, while the collecting plates are modeled using the hybrid method. The method combines the rigid finite element method with the classical finite element method. As a result, the mass matrix is diagonal. Some results of numerical simulations concerning free vibrations of the collecting plates and the influence of the number of elements, into which the plate is divided, on the vibrations of the rapping system are presented.
The finite element method (FEM) is one of the most frequently used numerical methods for finding the approximate discrete point solution of partial differential equations (PDE). In this method, linear or nonlinear systems of equations, comprised after numerical discretization, are solved to obtain the numerical solution of PDE. The conjugate gradient algorithms are efficient iterative solvers for the large sparse linear systems. In this paper the performance of different conjugate gradient algorithms: conjugate gradient algorithm (CG), biconjugate gradient algorithm (BICG), biconjugate gradient stabilized algorithm (BICGSTAB), conjugate gradient squared algorithm (CGS) and biconjugate gradient stabilized algorithm with l GMRES restarts (BICGSTAB(l)) is compared when solving the steady-state axisymmetric heat conduction problem. Different values of l parameter are studied. The engineering problem for which this comparison is made is the two-dimensional, axisymmetric heat conduction in a finned circular tube.
Various methods for steady-state and transient analysis of temperature distribution and efficiency of continuous-plate fins are presented. For a constant heat transfer coefficient over the fin surface, the plate fin can be divided into imaginary rectangular or hexangular fins. At first approximate methods for determining the steady-state fin efficiency like the method of equivalent circular fin and the sector method are discussed. When the fin geometry is complex, thus transient temperature distribution and fin efficiency can be determined using numerical methods. A numerical method for transient analysis of fins with complex geometry is developed. Transient temperature distributions in continuous fins attached to oval tubes is computed using the finite volume - finite element methods. The developed method can be used in the transient analysis of compact heat exchangers to calculate correctly the heat flow rate transferred from the finned tubes to the fluid.
The paper presents the dynamic behaviour of three-layer annular plates with damaged facings. The plate is composed of thin laminated, fibre-reinforced composite facings and thicker, foam core. Failure of the plate facings is modelled as fibre or matrix cracks. The plate loaded in the plane of facings with quickly increasing radially compressed forces loses its dynamic stability. Evaluation of the critical state of the plate with failures was carried out using both analytical and numerical solutions. The comparison of results between plates with material properties treated as isotropic, quasi-isotropic and composite has been conducted. Presented tables and figures create the image of dynamic responses of examined composite plates with structure failures.
Based on the example of the pilot area in Kiev the influence of the increased static load on the superstructure of the stress-strain state of the slope was studied. The efficiency of the proposed methodology when considering the work of "home-slope-retaining structure" depending on natural and anthropogenic factors was demonstrated.
Some materials-related microstructural problems calculated using the phase-field method are presented. It is well known that the phase field method requires mesh resolution of a diffuse interface. This makes the use of mesh adaptivity essential especially for fast evolving interfaces and other transient problems. Complex problems in 3D are also computationally challenging so that parallel computations are considered necessary. In this paper, a parallel adaptive finite element scheme is proposed. The scheme keeps the level of node and edge for 2D and level of node and face for 3D instead of the complete history of refinements to facilitate derefinement. The information is local and exchange of information is minimized and also less memory is used. The parallel adaptive algorithms that run on distributed memory machines are implemented in the numerical simulation of dendritic growth and capillary-driven flows.
Numerical methods are mostly used to predict the acoustic pressure inside duct systems. In this paper, the development of a numerical method based on the convected Helmholtz equation to compute the acoustic pressure inside an axisymmetric duct is presented. A validation of the proposed method was done by a comparison with the analytical formulation for simple cases of hard wall and lined ducts. The effect of the flow on the acoustic pressure inside these ducts was then evaluated by computing this field with different Mach numbers.
The rigid finite element method (RFEM) has been used mainly for modelling systems with beam-like links. This paper deals with modelling of a single set of electrodes consisting of an upper beam with electrodes, which are shells with complicated shapes, and an anvil beam. Discretisation of the whole system, both the beams and the electrodes, is carried out by means of the rigid finite element method. The results of calculations concerned with free vibrations of the plates are compared with those obtained from a commercial package of the finite element method (FEM), while forced vibrations of the set of electrodes are compared with those obtained by means of the hybrid finite element method (HFEM) and experimental measurements obtained on a special test stand.
The following paper presents the solution to the problem of searching the best shape - structural form of the bottoms and optimal dimensions of the main cylinder of the carding machine with consideration to the criterion of minimal deflection amplitude. The ANSYS package of the Finite Element Method has been used for the analysis. Polak-Ribery conjugate gradient method has been applied for searching the optimal solution, basing on the parametric model of the cylinder written with the use of Ansys Parametric Design Language. As a result of the performed analyses, reduction of maximum deflection value at approximately 80 percent has been obtained. Optimal cylinder dimensions enable application of a new textile technology - microfibre carding and improvement in the quality of traditional carding technology of woollen and wool-like fibres.
The article presents a numerical model of the concrete heat accumulator for solar heating systems. Model uses control volume finite element method with an explicit solution method for time integration. The use of an explicit method is an essential advantage in the simulation of time-dependent changes in temperature of the air at the accumulator inlet. The study compares the results of numerical model calculations of the accumulator heating with experimental measurements and with computational fluid dynamics modeling. The comparison shows a good correlation between the results of calculation using the model and the results of measurements.
Miniature heat exchangers are used to provide higher cooling capacity for new technologies. This means a reduction in their size and cost but the identical power. The paper presents the method for determination of boiling heat transfer coefficient for a rectangular minichannel of 0.1 mm depth, 40 mm width and 360 mm length with asymmetric heating. Experimental research has focused on the transition from single phase forced convection to nucleate boiling, i.e., the zone of boiling incipience. The ‘boiling front’ location has been determined from the temperature distribution of the heated wall obtained from liquid crystal thermography. The experiment has been carried out with R-123, mass flux 220 kg/(m2s), pressure at the channel inlet 340 kPa. Local values of heat transfer coefficient were calculated on the basis of empirical data from the experiment following the solution of the two-dimensional inverse heat transfer problem. This problem has been solved with the use of the finite element method in combination with Trefftz functions. Temperature approximates (linear combinations of Trefftz functions) strictly fulfill the governing equations. In presented method the inverse problem is solved in the same way as the direct problem. The results confirmed that considerable heat transfer enhancement takes place at boiling incipience in the minichannel flow boiling. Moreover, under subcooling boiling, local heat coefficients exhibit relatively low values.
Th This paper presents a simple method of evaluating the load resistance and stiffness of corrugatedsheets locally loaded with suspended technical fixtures. As a part of this research, parametricnumerical analyses of corrugated sheets of different span, and with differently located concentratedforces, were carried out. Stress distributions in the individual folds in the elastic range and inthe ultimate limit state were identified. On their basis, equivalent concentrated load factors forthe individual folds in the sheet were determined. The load factors enable analyses of the loadresistance and stiffness of corrugated sheets loaded with concentrated forces, which can be helpfulin design practice.
A three Dimensional finite element model (FEM) incorporating the anisotropic properties and temperature profile of hot mix asphalt (HMA) pavement was developed to predict the structural responses of HMA pavement subject to heavy loads typically encountered in the field. In this study, ABAQUS was adopted to model the stress and strain relationships within the pavement structure. The results of the model were verified using data collected from the Korean Highway Corporation Test Road (KHCTR). The results demonstrated that both the base course and surface course layers follow the anisotropic behavior and the incorporation of the temperature profile throughout the pavement has a substantial effect on the pavement response predictions that impact pavement design. The results also showed that the anisotropy level of HMA and base material can be reduced to as low as 80% and 15% as a result of repeated loading, respectively.
The present article investigates the dynamic behavior of a fully assembled turbogenerator system influenced by misalignment. In the past, most of the researchers have neglected the foundation flexibility in the turbogenerator systems in their study, to overcome this modelling error a more realistic model of a turbogenerator system has been attempted by considering flexible shafts, flexible coupling, flexible bearings and flexible foundation. Equations of motion for fully assembled turbogenerator system including flexible foundations have been derived by using finite element method. The methodology developed based on least squares technique requires forced response information to quantify the bearing–coupling–foundation dynamic parameters of the system associated with different faults along with residual unbalances. The proposed methodology is tested for the various level of measurement noise and modelling error in the system parameters, i.e., 5% deviation in E (modulus of elasticity) and ρ (density), respectively, for robustness of the algorithm. In a practical sense, the condition analyzed in the present article relates to the identification of misalignment and other dynamic parameters viz. bearing and residual unbalance in a rotor integrated with flexible foundation.
A two-scale numerical homogenization approach was used for granular materials. At small-scale level, granular micro-structure was simulated using the discrete element method. At macroscopic level, the finite element method was applied. An up-scaling technique took into account a discrete model at each Gauss integration point of the FEM mesh to derive numerically an overall constitutive response of the material. In this process, a tangent operator was generated with the stress increment corresponding to the given strain increment at the Gauss point. In order to detect a loss of the solution uniqueness, a determinant of the acoustic tensor associated with the tangent operator was calculated. Some elementary geotechnical tests were numerically calculated using a combined DEM-FEM technique.
The main idea of this work is to demonstrate an application of the generalized perturbation-based Stochastic Finite Element Method for a determination of the reliability indicators concerning elastic stability for a certain spectrum of the civil engineering structures. The reliability indicator is provided after the Eurocode according to the First Order Reliability Method, and computed using the higher order Taylor expansions with random coefficients. Computational implementation provided by the hybrid usage of the FEM system ROBOT and the computer algebra system MAPLE enables for reliability analysis of the critical forces in the most popular civil engineering structures like simple Euler beam, 2 and 3D single and multi-span steel frames, as well as polyethylene underground cylindrical shell. A contrast of the perturbation-based numerical approach with the Monte-Carlo simulation technique for the entire variability of the input random dispersion included into the Euler problem demonstrates the probabilistic efficiency of the perturbation method proposed.
The paper presents a tool for accurate evaluation of high field concentrations near singular lines, such as contours of cracks, notches and grains intersections, in 3D problems solved the BEM. Two types of boundary elements, accounting for singularities, are considered: (i) edge elements, which adjoin a singular line, and (ii) intermediate elements, which while not adjoining the line, are still under strong influence of the singularity. An efficient method to evaluate the influence coefficients and the field intensity factors is suggested for the both types of the elements. The method avoids time expensive numerical evaluation of singular and hypersingular integrals over the element surface by reduction to 1D integrals. The method being general, its details are explained by considering a representative examples for elasticity problems for a piece-wise homogeneous medium with cracks, inclusions and pores. Numerical examples for plane elements illustrate the exposition. The method can be extended for curvilinear elements.
The Bulletin of the Polish Academy of Sciences: Technical Sciences (Bull.Pol. Ac.: Tech.) is published bimonthly by the Division IV Engineering Sciences of the Polish Academy of Sciences, since the beginning of the existence of the PAS in 1952. The journal is peer‐reviewed and is published both in printed and electronic form. It is established for the publication of original high quality papers from multidisciplinary Engineering sciences with the following topics preferred: Artificial and Computational Intelligence, Biomedical Engineering and Biotechnology, Civil Engineering, Control, Informatics and Robotics, Electronics, Telecommunication and Optoelectronics, Mechanical and Aeronautical Engineering, Thermodynamics, Material Science and Nanotechnology, Power Systems and Power Electronics. Journal Metrics: JCR Impact Factor 2018: 1.361, 5 Year Impact Factor: 1.323, SCImago Journal Rank (SJR) 2017: 0.319, Source Normalized Impact per Paper (SNIP) 2017: 1.005, CiteScore 2017: 1.27, The Polish Ministry of Science and Higher Education 2017: 25 points. Abbreviations/Acronym: Journal citation: Bull. Pol. Ac.: Tech., ISO: Bull. Pol. Acad. Sci.-Tech. Sci., JCR Abbrev: B POL ACAD SCI-TECH Acronym in the Editorial System: BPASTS.
The paper presents a multi-scale mathematical model dedicated to a comprehensive simulation of resistance heating combined with the melting and controlled cooling of steel samples. Experiments in order to verify the formulated numerical model were performed using a Gleeble 3800 thermo-mechanical simulator. The model for the macro scale was based upon the solution of Fourier-Kirchhoff equation as regards predicting the distribution of temperature fields within the volume of the sample. The macro scale solution is complemented by a functional model generating voluminal heat sources, resulting from the electric current flowing through the sample. The model for the micro-scale, concerning the grain growth simulation, is based upon the probabilistic Monte Carlo algorithm, and on the minimization of the system energy. The model takes into account the forming mushy zone, where grains degrade at the melting stage – it is a unique feature of the micro-solution. The solution domains are coupled by the interpolation of node temperatures of the finite element mesh (the macro model) onto the Monte Carlo cells (micro model). The paper is complemented with examples of resistance heating results and macro- and micro-structural tests, along with test computations concerning the estimation of the range of zones with diverse dynamics of grain growth.
The aim of the study is to identify the relevant aspects of numerical analysis of impact of projectiles with soft cores into a package composed of thin flexible plies located on the plastic backing. In order to illustrate the problem, normal impact of 7.62 mm TT projectile into an unclamped package comprising 36 plies of Dyneema SB71 supported on the plastic backing was selected. The problem was solved with the use of the finite element method (FEM) with the explicit integration scheme (central difference method) of motion equations in the matrix form. Based on the conducted numerical computations, it was revealed that obtaining the extreme deformations of a projectile soft core and the backing material in Lagrangian description requires employment of adaptive methods. The proposed R-adaptive method performs its role but must be used carefully due to the mass loss which may appear during calculations.
The field of mechanical manufacturing is becoming more and more demanding on machining accuracy. It is essential to monitor and compensate the deformation of structural parts of a heavy-duty machine tool. The deformation of the base of a heavy-duty machine tool is an important factor that affects machining accuracy. The base is statically indeterminate and complex in load. It is difficult to reconstruct deformation by traditional methods. A reconstruction algorithm for determining bending deformation of the base of a heavy-duty machine tool using inverse Finite Element Method (iFEM) is presented. The base is equivalent to a multi-span beam which is divided into beam elements with support points as nodes. The deflection polynomial order of each element is analysed. According to the boundary conditions, the deformation compatibility conditions and the strain data measured by Fiber Bragg Grating (FBG), the deflection polynomial coefficients of a beam element are determined. Using the coordinate transformation, the deflection equation of the base is obtained. Both numerical verification and experiment were carried out. The deflection obtained by the reconstruction algorithm using iFEM and the actual deflection measured by laser displacement sensors were compared. The accuracy of the reconstruction algorithm is verified.