Stepwise analysis of pantographic beams subjected to impulsive loads

2020 ◽  
Vol 26 (1) ◽  
pp. 62-79
Author(s):  
Emilio Turco

Materials based on pantographic unit cells have very interesting mechanical peculiarities. For these reasons they are largely studied from a theoretical, experimental, and numerical point of view. Numerical simulations furnish an important contribution for the the design and optimization of such materials and, more generally, for metamaterials. Here, we consider the influence of inertial forces, removing the hypothesis of quasistatic loading. By using an intrinsically discrete model, inspired by Hencky’s ideas, already tested in a series of published works, here we add the contribution of inertial forces and, in the framework of stepwise schemes, we re-experience an adaptive integration scheme capable of reconstructing the best structural response corresponding to a prefixed time step. Several numerical simulations, although preparatory, inspire some remarks on materials based on pantographic cells and outline the way for future challenges.

2021 ◽  
pp. 108128652110108
Author(s):  
Emilio Turco ◽  
Emilio Barchiesi ◽  
Francesco dell’Isola

This contribution presents the results of a campaign of numerical simulations aimed at better understanding the propagation of longitudinal waves in pantographic beams within the large-deformation regime. Initially, we recall the key features of a Lagrangian discrete spring model, which was introduced in previous works and that was tested extensively as capable of accurately forecasting the mechanical response of structures based on the pantographic motif, both in statics and dynamics. Successively, a stepwise integration scheme used to solve equations of motions is briefly discussed. The key content of the present contribution concerns the thorough presentation of some selected numerical simulations, which focus in particular on the propagation of stretch profiles induced by impulsive loads. The study takes into account different tests, by varying the number of unit cells, i.e., the total length of the system, spring stiffnesses, the shape of the impulse, as well as its properties such as duration and peak amplitude, and boundary conditions. Some conjectures about the form of traveling waves are formulated, to be confirmed by both further numerical simulations and analytical investigations.


Author(s):  
Xian Luo ◽  
Martin R. Maxey ◽  
George E. Karniadakis

We re-formulate and demonstrate a new method for particulate flows, the so-called “Smoothed Profile” method (SPM) first proposed in [1]. The method uses a fixed computational mesh, which does not conform to the geometry of the particles. The particles are represented by certain smoothed indicator profiles to construct a smooth body force density term added into the Navier-Stokes equations. The SPM imposes accurately and efficiently the rigid-body constraint inside the particles. In particular, while the original method employs a fully-explicit time-integration scheme, we develop a high-order semi-implicit splitting scheme, which we implement in the context of spectral/hp element discretization. We show that the modeling error of SPM has a non-monotonic dependence on the time step size Δt. The optimum time step size balances the thickness of the Stokes layer and that of the profile interface. Subsequently, we present several numerical simulations, including flow past three-dimensional complex-shaped particles and two interacting microspheres, which are compared against full direct numerical simulations and the force coupling method (FCM).


Author(s):  
M. El Badaoui ◽  
V. Cahouet ◽  
F. Guillet ◽  
J. Daniere ◽  
P. Velex

Abstract The early detection of failures in geared systems is an important industrial problem which has still to be addressed from both an experimental and theoretical viewpoint. The proposed paper combines some extensive numerical simulations of a single stage geared unit with localized tooth faults and the use of several detection techniques whose performances are compared and critically assessed. A model aimed at simulating the contributions of local tooth defects such as spalling to the gear dynamic behavior is set up. The pinion and the gear of a pair are modeled as to two rigid cylinders with all six degrees of freedom connected by a series of springs which represent gear body and gear tooth compliances on the base plane. Classical shaft finite elements including torsional, flexural and axial displacements can be superimposed on the gear element together with some lumped stiffnesses, masses, inertias, ... which account for the load machines, bearings and couplings. Tooth defects are modeled by a distribution of normal deviations over a zone which can be located anywhere on the active tooth flanks. Among the numerous available signal processing techniques used in vibration monitoring, cepstrum analysis is sensitive, reliable and it can be adapted to a complex geared system with several meshes. From an analytical analysis of the equations of motion, two complementary detection techniques based upon the acceleration power cesptrum are proposed. The equations of motion and the contact problem between mating flanks are simultaneously solved by coupling an implicit time-step integration scheme and a unilateral normal contact algorithm. The results of the numerical simulations are used as a data base for the proposed detection techniques. The combined influence of the defect location, depth and extent is analyzed for two examples of spur and helical gears with various profile modifications and the effectiveness of the two complementary detection methods is discussed before some conclusions are drawn.


1999 ◽  
Vol 123 (3) ◽  
pp. 422-430 ◽  
Author(s):  
M. El Badaoui ◽  
V. Cahouet ◽  
F. Guillet ◽  
J. Danie`re ◽  
P. Velex

The early detection of failures in geared systems is an important industrial problem which has still to be addressed from both an experimental and theoretical viewpoint. The proposed paper combines some extensive numerical simulations of a single stage geared unit with localized tooth faults and the use of several detection techniques whose performances are compared and critically assessed. A model aimed at simulating the contributions of local tooth defects such as spalling to the gear dynamic behavior is set up. The pinion and the gear of a pair are assimilated to two rigid cylinders with all six degrees of freedom connected by a series of springs which represent gear body and gear tooth compliances on the base plane. Classical shaft finite elements including torsional, flexural and axial displacements can be superimposed to the gear element together with some lumped stiffnesses, masses, inertias, … which account for the load machines, bearings and couplings. Tooth defects are modeled by a distribution of normal deviations over a zone which can be located anywhere on the active tooth flanks. Among the numerous available signal processing techniques used in vibration monitoring, cepstrum analysis is sensitive, reliable and it can be adapted to complex geared system with several meshes. From an analytical analysis of the equations of motion, two complementary detection techniques based upon acceleration power cepstrum are proposed. The equations of motion and the contact problem between mating flanks are simultaneously solved by coupling an implicit time-step integration scheme and a unilateral normal contact algorithm. The results of the numerical simulations are used as a data base for the proposed detection techniques. The combined influence of the defect location, depth and extent is analyzed for two examples of spur and helical gears with various profile modifications and the effectiveness of the two complementary detection methods is discussed before some conclusions are drawn.


Algorithms ◽  
2021 ◽  
Vol 14 (5) ◽  
pp. 129
Author(s):  
Yuan Li ◽  
Ni Zhang ◽  
Yuejiao Gong ◽  
Wentao Mao ◽  
Shiguang Zhang

Compared with continuous elements, discontinuous elements advance in processing the discontinuity of physical variables at corner points and discretized models with complex boundaries. However, the computational accuracy of discontinuous elements is sensitive to the positions of element nodes. To reduce the side effect of the node position on the results, this paper proposes employing partially discontinuous elements to compute the time-domain boundary integral equation of 3D elastodynamics. Using the partially discontinuous element, the nodes located at the corner points will be shrunk into the element, whereas the nodes at the non-corner points remain unchanged. As such, a discrete model that is continuous on surfaces and discontinuous between adjacent surfaces can be generated. First, we present a numerical integration scheme of the partially discontinuous element. For the singular integral, an improved element subdivision method is proposed to reduce the side effect of the time step on the integral accuracy. Then, the effectiveness of the proposed method is verified by two numerical examples. Meanwhile, we study the influence of the positions of the nodes on the stability and accuracy of the computation results by cases. Finally, the recommended value range of the inward shrink ratio of the element nodes is provided.


2013 ◽  
Vol 2013 ◽  
pp. 1-21 ◽  
Author(s):  
Rita Greco ◽  
Francesco Trentadue

Response sensitivity evaluation is an important element in reliability evaluation and design optimization of structural systems. It has been widely studied under static and dynamic forcing conditions with deterministic input data. In this paper, structural response and reliability sensitivities are determined by means of the time domain covariance analysis in both classically and nonclassically damped linear structural systems. A time integration scheme is proposed for covariance sensitivity. A modulated, filtered, white noise input process is adopted to model the stochastic nonstationary loads. The method allows for the evaluation of sensitivity statistics of different quantities of dynamic response with respect to structural parameters. Finally, numerical examples are presented regarding a multistorey shear frame building.


2013 ◽  
Vol 135 (6) ◽  
Author(s):  
R. Fargère ◽  
P. Velex

A global model of mechanical transmissions is introduced which deals with most of the possible interactions between gears, shafts, and hydrodynamic journal bearings. A specific element for wide-faced gears with nonlinear time-varying mesh stiffness and tooth shape deviations is combined with shaft finite elements, whereas the bearing contributions are introduced based on the direct solution of Reynolds' equation. Because of the large bearing clearances, particular attention has been paid to the definition of the degrees-of-freedom and their datum. Solutions are derived by combining a time step integration scheme, a Newton–Raphson method, and a normal contact algorithm in such a way that the contact conditions in the bearings and on the gear teeth are simultaneously dealt with. A series of comparisons with the experimental results obtained on a test rig are given which prove that the proposed model is sound. Finally, a number of results are presented which show that parameters often discarded in global models such as the location of the oil inlet area, the oil temperature in the bearings, the clearance/elastic couplings interactions, etc. can be influential on static and dynamic tooth loading.


Author(s):  
SD Yu ◽  
BC Wen

This article presents a simple procedure for predicting time-domain vibrational behaviors of a multiple degrees of freedom mechanical system with dry friction. The system equations of motion are discretized by means of the implicit Bozzak–Newmark integration scheme. At each time step, the discontinuous frictional force problem involving both the equality and inequality constraints is successfully reduced to a quadratic mathematical problem or the linear complementary problem with the introduction of non-negative and complementary variable pairs (supremum velocities and slack forces). The so-obtained complementary equations in the complementary pairs can be solved efficiently using the Lemke algorithm. Results for several single degree of freedom and multiple degrees of freedom problems with one-dimensional frictional constraints and the classical Coulomb frictional model are obtained using the proposed procedure and compared with those obtained using other approaches. The proposed procedure is found to be accurate, efficient, and robust in solving non-smooth vibration problems of multiple degrees of freedom systems with dry friction. The proposed procedure can also be applied to systems with two-dimensional frictional constraints and more sophisticated frictional models.


Water ◽  
2021 ◽  
Vol 13 (6) ◽  
pp. 785
Author(s):  
Arman Rokhzadi ◽  
Musandji Fuamba

This paper studies the air pressurization problem caused by a partially pressurized transient flow in a reservoir-pipe system. The purpose of this study is to analyze the performance of the rigid column model in predicting the attenuation of the air pressure distribution. In this regard, an analytic formula for the amplitude and frequency will be derived, in which the influential parameters, particularly, the driving pressure and the air and water lengths, on the damping can be seen. The direct effect of the driving pressure and inverse effect of the product of the air and water lengths on the damping will be numerically examined. In addition, these numerical observations will be examined by solving different test cases and by comparing to available experimental data to show that the rigid column model is able to predict the damping. However, due to simplified assumptions associated with the rigid column model, the energy dissipation, as well as the damping, is underestimated. In this regard, using the backward Euler implicit time integration scheme, instead of the classical fourth order explicit Runge–Kutta scheme, will be proposed so that the numerical dissipation of the backward Euler implicit scheme represents the physical dissipation. In addition, a formula will be derived to calculate the appropriate time step size, by which the dissipation of the heat transfer can be compensated.


Author(s):  
Fabien Bigot ◽  
François-Xavier Sireta ◽  
Eric Baudin ◽  
Quentin Derbanne ◽  
Etienne Tiphine ◽  
...  

Ship transport is growing up rapidly, leading to ships size increase, and particularly for container ships. The last generation of Container Ship is now called Ultra Large Container Ship (ULCS). Due to their increasing sizes they are more flexible and more prone to wave induced vibrations of their hull girder: springing and whipping. The subsequent increase of the structure fatigue damage needs to be evaluated at the design stage, thus pushing the development of hydro-elastic simulation models. Spectral fatigue analysis including the first order springing can be done at a reasonable computational cost since the coupling between the sea-keeping and the Finite Element Method (FEM) structural analysis is performed in frequency domain. On the opposite, the simulation of non-linear phenomena (Non linear springing, whipping) has to be done in time domain, which dramatically increases the computation cost. In the context of ULCS, because of hull girder torsion and structural discontinuities, the hot spot stress time series that are required for fatigue analysis cannot be simply obtained from the hull girder loads in way of the detail. On the other hand, the computation cost to perform a FEM analysis at each time step is too high, so alternative solutions are necessary. In this paper a new solution is proposed, that is derived from a method for the efficient conversion of full scale strain measurements into internal loads. In this context, the process is reversed so that the stresses in the structural details are derived from the internal loads computed by the sea-keeping program. First, a base of distortion modes is built using a structural model of the ship. An original method to build this base using the structural response to wave loading is proposed. Then a conversion matrix is used to project the computed internal loads values on the distortion modes base, and the hot spot stresses are obtained by recombination of their modal values. The Moore-Penrose pseudo-inverse is used to minimize the error. In a first step, the conversion procedure is established and validated using the frequency domain hydro-structure model of a ULCS. Then the method is applied to a non-linear time domain simulation for which the structural response has actually been computed at each time step in order to have a reference stress signal, in order to prove its efficiency.


Sign in / Sign up

Export Citation Format

Share Document