Search results

When establishing a mathematical model to simulate solid mechanics, considering realistic geometries, special tools are needed to translate measured data, possibly with noise, into idealized geometrical entities. As an engineering application, wheel-rail contact interactions are fundamental in the dynamic modeling of railway vehicles. Many approaches used to solve the contact problem require a continuous parametric description of the geometries involved. However, measured wheel and rail profiles are often available as sets of discrete points. A reconstruction method is needed to transform discrete data into a continuous geometry.

The authors present an approximation method based on optimization to solve the problem of fitting a set of points with an arc spline. It consists of an initial guess based on a curvature function estimated from the data, followed by a least-squares optimization to improve the approximation. The authors also present a segmentation scheme that allows the method to increment the number of segments of the spline, trying to keep it at a minimal value, to satisfy a given error tolerance.

The paper provides a better understanding of arc splines and how they can be deformed. Examples with parametric curves and slightly noisy data from realistic wheel and rail profiles show that the approach is successful.

The developed methods have theoretical value. Furthermore, they have practical value since the approximation approach is better suited to deal with the reconstruction of wheel/rail profiles than interpolation, which most methods use to some degree.

The purpose of this study is to investigate the post-buckling analysis of functionally graded columns by using three analytical, approximate and numerical methods. A pre-defined function as an initial assumption for the post-buckling path is introduced to solve the differential equation. The finite difference method is used to approximate the lateral deflection of the column based on the differential equation. Moreover, the finite element method is used to derive the tangent stiffness matrix of the column.

The non-linear buckling analysis of functionally graded materials is carried out by using three analytical, finite difference and finite element methods. The elastic deformation and Euler-Bernoulli beam theory are considered to establish the constitutive and kinematics relations, respectively. The governing differential equation of the post-buckling problem is derived through the energy method and the calculus variation.

An incremental iterative solution and the perturbation of the displacement vector at the critical buckling point are performed to determine the post-buckling path. The convergence of the finite element results and the effects of geometric and material characteristics on the post-buckling path are investigated.

The key point of the research is to compare three methods and to detect error sources by considering the derivation process of relations. This comparison shows that a non-incremental solution in the analytical and finite difference methods and an initial assumption in the analytical method lead to an error in results. However, the post-buckling path in the finite element method is traced by the updated tangent stiffness matrix in each load step without any initial limitation.

Currently, several studies have been published using sensorized insoles for estimating ground reaction force using plantar pressure. However, information on design parameters, manufacturing techniques and guidelines for developing insoles is scarce, often leaving gaps that do not allow reproducing the insole. This study aims to empirically investigate the main parameters of constructing a sensorized insole for application in human gait.

Two devices were built to evaluate the force sensors. The first focuses on the construction of the sensors with different settings: the density of the sensor’s conductive trails (thickness and distance of the trails) and the inertia of the sensors (use of spacers to prevent unwanted readings). The second device focuses on the data capture and processing system: resolution of the analog–digital converter, acquisition rate and sensor activation level.

The resolution increase of the analog–digital converter and acquisition rate do not contribute to noise increase. Reducing the sensors’ coverage area can increase sensorized insole capacity. The inertia of the sensors can be adjusted using spacers without changing the electrical circuit and acquisition system.

Most sensorized insoles use commercial sensors. For this reason, it is not possible a full customization. This paper maps the main variables to manufacture custom sensors and data acquisition systems. This work also presents a case study where it is possible to see the influence of the parameters in the correlation between the sensorized insole and an instrumented treadmill with a force platform.

The ability to characterise highly conducting objects, that may also be highly magnetic, by the complex symmetric rank–2 magnetic polarizability tensor (MPT) is important for metal detection applications including discriminating between threat and non-threat objects in security screening, identifying unexploded anti-personnel landmines and ordnance and identifying metals of high commercial value in scrap sorting. Many everyday non-threat items have both a large electrical conductivity and a magnetic behaviour, which, for sufficiently weak fields and the frequencies of interest, can be modelled by a high relative magnetic permeability. This paper aims to discuss the aforementioned idea.

The numerical simulation of the MPT for everyday non-threat highly conducting magnetic objects over a broad range of frequencies is challenging due to the resulting thin skin depths. The authors address this by employing higher order edge finite element discretisations based on unstructured meshes of tetrahedral elements with the addition of thin layers of prismatic elements. Furthermore, computer aided design (CAD) geometrical models of the non-threat and threat object are often not available and, instead, the authors extract the geometrical features of an object from an imaging procedure.

The authors obtain accurate numerical MPT characterisations that are in close agreement with experimental measurements for realistic physical objects. The assessment of uncertainty shows the impact of geometrical and material parameter uncertainties on the computational results.

The authors present novel computations and measurements of MPT characterisations of realistic objects made of magnetic materials. A novel assessment of uncertainty in the numerical predictions of MPT characterisations for uncertain geometry and material parameters is included.

1) The OE-MLPG penalty meshfree method is developed to solve cracked structure.(2) Smartening the numerical meshfree method by combining the particle swarm optimization (PSO) optimization algorithms and Voronoi computational geometric algorithm. (3). Selection of base functions, finding optimal penalty factor and distribution of appropriate nodal points to the accuracy of calculation in the meshless local Petrov–Galekrin (MLPG) meshless method.

Using appropriate shape functions and distribution of nodal points in local domains and sub-domains and choosing an approximation or interpolation method has an effective role in the application of meshless methods for the analysis of computational fracture mechanics problems, especially problems with geometric discontinuity and cracks. In this research, computational geometry technique, based on the Voronoi diagram (VD) and Delaunay triangulation and PSO algorithm, are used to distribute nodal points in the sub-domain of analysis (crack line and around it on the crack plane).

By doing this process, the problems caused by too closeness of nodal points in computationally sensitive areas that exist in general methods of nodal point distribution are also solved. Comparing the effect of the number of sentences of basic functions and their order in the definition of shape functions, performing the mono-objective PSO algorithm to find the penalty factor, the coefficient, convergence, arrangement of nodal points during the three stages of VD implementation and the accuracy of the answers found indicates, the efficiency of V-E-MLPG method with Ns = 7 and ß = 0.0037–0.0075 to estimation of 3D-stress intensity factors (3D-SIFs) in computational fracture mechanics.

The present manuscript is a continuation of the studies (Ref. [33]) carried out by the authors, about; feasibility assessment, improvement and solution of challenges, introduction of more capacities and capabilities of the numerical MLPG method have been used. In order to validate the modeling and accuracy of calculations, the results have been compared with the findings of reference article [34] and [35].

This study aims to assess the robustness and accuracy of the face-centred finite volume (FCFV) method for the simulation of compressible laminar flows in different regimes, using numerical benchmarks.

The work presents a detailed comparison with reference solutions published in the literature –when available– and numerical results computed using a commercial cell-centred finite volume software.

The FCFV scheme provides first-order accurate approximations of the viscous stress tensor and the heat flux, insensitively to cell distortion or stretching. The strategy demonstrates its efficiency in inviscid and viscous flows, for a wide range of Mach numbers, also in the incompressible limit. In purely inviscid flows, non-oscillatory approximations are obtained in the presence of shock waves. In the incompressible limit, accurate solutions are computed without pressure correction algorithms. The method shows its superior performance for viscous high Mach number flows, achieving physically admissible solutions without carbuncle effect and predictions of quantities of interest with errors below 5%.

The FCFV method accurately evaluates, for a wide range of compressible laminar flows, quantities of engineering interest, such as drag, lift and heat transfer coefficients, on unstructured meshes featuring distorted and highly stretched cells, with an aspect ratio up to ten thousand. The method is suitable to simulate industrial flows on complex geometries, relaxing the requirements on mesh quality introduced by existing finite volume solvers and alleviating the need for time-consuming manual procedures for mesh generation to be performed by specialised technicians.

The aim of this paper is to provide a simple yet accurate and efficient geometric method for thermal homogenization of impregnated and non-impregnated coil winding technologies based on the concept of thermal resistance.

For regular windings, the periodic microscopic cell in the winding space is identified. Also, for irregular windings, the average microscopic cell of the winding is determined. An approximation is used to calculate the thermal resistance of the winding cell. Based on this approximation, the winding insulation is considered as a circular ring around the wire. Mathematical equations are obtained to calculate the equivalent thermal resistance of the cell. The equivalent thermal conductivity of the winding is calculated using equivalent thermal resistance of the cell. Winding thermal homogenization is completed by determining the equivalent thermal properties of the cell.

The thermal pattern of different windings is simulated and compared with the results of different homogenization methods. The results show that the proposed method is applicable for a wide range of windings in terms of winding scheme, packing factor and winding insulation. Also, the results show that the proposed method is more accurate than other winding homogenization methods in calculating the equivalent thermal conductivity of the winding.

In this paper, the change of electrical resistance of the winding with temperature and thermal contact between the sub-components are ignored. Also, liquid insulators, such as oils, and rectangular wires were not investigated. Research in these topics is considered as future work.

Unlike other homogenization methods, the proposed method can be applied to non-impregnated and irregular windings. Also, compared to other homogenization methods, the proposed method has a simpler formulation that makes it easier to program and implement. All of these indicate the efficiency of the proposed method in the thermal analysis of the winding.

The conventional design method relies on a priori knowledge, which limits the rapid and efficient development of electronic packaging structures. The purpose of this study is to propose a hybrid method of data-driven inverse design, which couples adaptive surrogate model technology with optimization algorithm to to enable an efficient and accurate inverse design of electronic packaging structures.

The multisurrogate accumulative local error-based ensemble forward prediction model is proposed to predict the performance properties of the packaging structure. As the forward prediction model is adaptive, it can identify respond to sensitive regions of design space and sample more design points in those regions, getting the trade-off between accuracy and computation resources. In addition, the forward prediction model uses the average ensemble method to mitigate the accuracy degradation caused by poor individual surrogate performance. The Particle Swarm Optimization algorithm is then coupled with the forward prediction model for the inverse design of the electronic packaging structure.

Benchmark testing demonstrated the superior approximate performance of the proposed ensemble model. Two engineering cases have shown that using the proposed method for inverse design has significant computational savings while ensuring design accuracy. In addition, the proposed method is capable of outputting multiple structure parameters according to the expected performance and can design the packaging structure based on its extreme performance.

Because of its data-driven nature, the inverse design method proposed also has potential applications in other scientific fields related to optimization and inverse design.

The purpose of this paper is to improve autonomous flight performance of a fixed-wing unmanned aerial vehicle (UAV) via simultaneous morphing wingtip and control system design conducted with optimization, computational fluid dynamics (CFD) and machine learning approaches.

The main wing of the UAV is redesigned with morphing wingtips capable of dihedral angle alteration by means of folding. Aircraft dynamic model is derived as equations depending only on wingtip dihedral angle via Nonlinear Least Squares regression machine learning algorithm. Data for the regression analyses are obtained by numerical (i.e. CFD) and analytical approaches. Simultaneous perturbation stochastic approximation (SPSA) is incorporated into the design process to determine the optimal wingtip dihedral angle and proportional-integral-derivative (PID) coefficients of the control system that maximizes autonomous flight performance. The performance is defined in terms of trajectory tracking quality parameters of rise time, settling time and overshoot. Obtained optimal design parameters are applied in flight simulations to test both longitudinal and lateral reference trajectory tracking.

Longitudinal and lateral autonomous flight performances of the UAV are improved by redesigning the main wing with morphing wingtips and simultaneous estimation of PID coefficients and wingtip dihedral angle with SPSA optimization.

This paper originally discusses the simultaneous design of innovative morphing wingtip and UAV flight control system for autonomous flight performance improvement. The proposed simultaneous design idea is conducted with the SPSA optimization and a machine learning algorithm as a novel approach.

The lack of geometric and dimensional accuracy of parts produced by additive manufacturing (AM) is directly related to the machine, material and process used. This paper aims to propose a method for the analysis and compensation of machine-related geometric errors applicable to any AM machine, regardless of the manufacturing process and technology used.

For this purpose, an error calculation model inspired by those used in computerized numerical control machines and coordinate measuring machines was developed. The error functions of the model were determined from the position deviations of a set of virtual points that are not sensitive to material and process errors. These points were obtained from the measurement of an ad hoc designed and manufactured master artefact. To validate the model, off-line compensation was applied to both the original designed artefact and an example part.

The geometric deviations in both cases were significantly smaller than those found before applying the geometric compensation. Dimensional enhancements were also achieved on the example part by using a correction parameter available in the three-dimensional printing software, whose value was adjusted from the measurement of the geometrically compensated master artefact.

The errors that persist in the part derive from both material and process. Compensation for these type of errors requires a detailed analysis of the influencing parameters, which will be the subject of future research.

The use of the virtual-point-based error model increases the quality of additively manufactured parts and can be used in any AM system.