scholarly journals Learning data-driven reduced elastic and inelastic models of spot-welded patches

2021 ◽  
Vol 22 ◽  
pp. 32
Author(s):  
Agathe Reille ◽  
Victor Champaney ◽  
Fatima Daim ◽  
Yves Tourbier ◽  
Nicolas Hascoet ◽  
...  

Solving mechanical problems in large structures with rich localized behaviors remains a challenging issue despite the enormous advances in numerical procedures and computational performance. In particular, these localized behaviors need for extremely fine descriptions, and this has an associated impact in the number of degrees of freedom from one side, and the decrease of the time step employed in usual explicit time integrations, whose stability scales with the size of the smallest element involved in the mesh. In the present work we propose a data-driven technique for learning the rich behavior of a local patch and integrate it into a standard coarser description at the structure level. Thus, localized behaviors impact the global structural response without needing an explicit description of that fine scale behaviors.

Author(s):  
Jeremy J. Laflin ◽  
Kurt S. Anderson ◽  
Michael Hans

Since computational performance is critically important for simulations to be used as an effective tool to study and design dynamic systems, the computing performance gains offered by Graphics Processing Units (GPUs) cannot be ignored. Since the GPU is designed to execute a very large number of simultaneous tasks (nominally Single Instruction Multi-Data (SIMD)), recursive algorithms in general, such as the DCA, are not well suited to be executed on GPU-type architecture. This is because each level of recursion is dependent on the previous level. However, there are some ways that the GPU can be leveraged to increase computational performance when using the DCA to form and solve the equations of motion for articulated multibody systems with a very large number of degrees-of-freedom. Computational performance of dynamic simulations is highly dependent on the nature of the underlying formulation and the number of generalized coordinates used to characterize the system. Therefore, algorithms that scale in a more desirable (lower order) fashion with the number of degrees-of-freedom are generally preferred when dealing with large (N > 10) systems. However, the utility of using simulations as a scientific tool is directly related to actual compute time. The DCA, and other top performing methods, have demonstrated the desirable property of the required compute time scaling linearly with (O(n)) with the number of degrees-of-freedom (n) and sublinearly (O(logn) performance when implemented in parallel. However for the DCA, total compute time could be further reduced by exploiting the large number of independent operations involved in the first few levels of recursion. A simple chain-type pendulum example is used to explore the feasibility of using the GPU to execute the assembly and disassembly operations for the levels of recursion that contain enough bodies for this process to be computationally advantageous. A multi-core CPU is used to perform the operations in parallel using Open MP for the remaining levels. The number of levels of recursion that utilizes the GPU is varied from zero to all levels. The data corresponding to zero utilization of the GPU provides the reference compute-time in which the assembly and disassembly operations necessary at each level are performed in parallel using Open MP. The computational time required to simulate the system for one time-step where the GPU is utilized for various levels of recursion is compared to the reference compute time also varying the number of bodies in the system. A decrease in the compute-time when using the GPU is demonstrated relative to the reference compute-time even for systems of moderate size n < 1000 for arrangements using the GPU. This is a lower number of bodies than was expected for this test case and confirms that the GPU can bring significant increases in computational efficiency for large systems, while preserving the attractive sub-linear scalability (w.r.t. compute time) of the DCA.


2021 ◽  
Vol 241 ◽  
pp. 112405
Author(s):  
Debarshi Sen ◽  
James Long ◽  
Hao Sun ◽  
Xander Campman ◽  
Oral Buyukozturk

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.


2018 ◽  
Author(s):  
Jihyeon Kim ◽  
Narakorn Srinil

Subsea jumper is the steel pipe structure to connect wellhead and subsea facilities such as manifolds or processing units in order to transport the produced multiphase flows. Generally, the jumper consists of a goalpost with two loop structures and a straight pipe between them, carrying the multiphase oil and gas from the producing well. Due to the jumper pipe characteristic geometry and multi-fluid properties, slug flows may take place, creating problematic fluctuating forces causing the jumper oscillations. Severe dynamic fluctuations cause the risk of pipe deformations and resonances resulting from the hydrodynamic momentum/pressure forces which can lead to unstable operating pressure and decreased production rate. Despite the necessity to design subsea jumper with precise prediction on the process condition and the awareness of slug flow risks, it is challenging to experimentally evaluate, identify and improve the modified design in terms of the facility scale, time and cost efficiency. With increasing high computational performance, numerical analysis provides an alternative approach to simulate multiphase flow-induced force effects on the jumper. The present paper discusses the modelling of 3-D flow simulations in a subsea jumper for understanding the development process of internal slug flows causing hydrodynamic forces acting on the pipe walls and bends. Based on the fluctuating pressure calculated by the fluid solver, dynamic responses of the jumper pipe are assessed by a one-way interaction approach to evaluate deformation and stress. A potential resonance is discussed with the jumper modal analysis. Results from the structural response analyses show dominant multi-modal frequencies due to intermittent slug flow frequencies. Numerical results and observed behaviors may be useful for a comparison with other simulation and experiment.


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.


2018 ◽  
Vol 2018 ◽  
pp. 1-17 ◽  
Author(s):  
Qi Chen ◽  
Xugang Wang ◽  
Jing Yang

An indirect Gauss pseudospectral method based path-following guidance law is presented in this paper. A virtual target moving along the desired path with explicitly specified speed is introduced to formulate the guidance problem. By establishing a virtual target-fixed coordinate system, the path-following guidance is transformed into a terminal guidance with impact angle constraints, which is then solved by using indirect Gauss pseudospectral method. Meanwhile, the acceleration dynamics are modeled as the first-order lag to the command. Using the receding horizon technique a closed-loop guidance law, which considers generalized weighting functions (even discontinuous) of both the states and the control cost, is derived. The accuracy and effectiveness of the proposed guidance law are validated by numerical comparisons. A STM32 Nucleo board based on the ARM Cortex-M7 processor is used to evaluate the real-time computational performance of the proposed indirect Gauss pseudospectral method. Simulations for various types of desired paths are presented to show that the proposed guidance law has better performance when compared with the existing results for pure pursuit, a nonlinear guidance law, and trajectory shaping path-following guidance and provides more degrees of freedom in path-following guidance design applications.


2020 ◽  
Author(s):  
◽  
Uriel Jacket Tresor Demby's

In the context of articulated robotic manipulators, the Forward Kinematics (FK) is a highly non-linear function that maps joint configurations of the robot to poses of its endeffector. Furthermore, while in the most useful cases these functions are neither injective (one-to-one) nor surjective (onto), depending on the robot configuration -- i.e. the sequence of prismatic versus revolute joints, and the number of Degrees of Freedom (DoF) -- the associated Inverse Kinematics (IK) problem may be practically or even theoretically impossible to be solved analytically. Therefore, in the past decades, several approximate methods have been developed for many instances of IK problems. The approximate methods can be divided into two distinct categories: data-driven and numerical approaches. In the first case, data-driven approaches have been successfully used for small workspace domains (e.g., task-driven applications), but not fully explored for large ones, i.e. in task-independent applications where a more general IK is required. Similarly, and despite many successful implementations over the years, numerical solutions may fail if an improper matrix inverse is employed (e.g., Moore-Penrose generalized inverse). In this research, we propose a systematic, robust and accurate numerical solution for the IK problem using the Unit-Consistent (UC) and the Mixed (MX) Inverse methods to invert the Jacobians derived from the Denavit-Hartenberg (D-H) representation of the FK for any robot. As we demonstrate, this approach is robust to whether the system is underdetermined (less than 6 DoF) or overdetermined (more than 6 DoF). We compare the proposed numerical solution to data driven solutions using different robots -- with DoF varying from 3 to 7. We conclude that numerical solutions are easier to implement, faster, and more accurate than most data-driven approaches in the literature, specially for large workspaces as in task-independent applications. We particularly compared the proposed numerical approach against two data-driven approaches: Multi-Layer Perceptron (MLP) and Adaptive Neuro-Fuzzy Inference System (ANFIS), while exploring various architectures of these Neural Networks (NN): i.e. number of inputs, number of outputs, depth, and number of nodes in the hidden layers.


2008 ◽  
Vol 17 (02) ◽  
pp. 265-273 ◽  
Author(s):  
JAMES A. ISENBERG

The analysis of a general multibody physical system governed by Einstein's equations is quite difficult, even if numerical methods (on a computer) are used. Some of the difficulties — many coupled degrees of freedom, dynamic instability — are associated with the presence of gravitational waves. We have developed a number of "waveless approximation theories" (WAT's) which repress the gravitational radiation and thereby simplify the analysis. The matter, according to these theories, evolves dynamically. The gravitational field, however, is determined at each time step by a set of elliptic equations with matter sources. There is reason to believe that for many physical systems, the WAT-generated system evolution is a very accurate approximation to that generated by the full Einstein theory.


Author(s):  
Mohammad Reza Tabeshpour ◽  
Reza Hedayatpour

Having deep view in structural response of tension leg platform is important issue not only for response analysis but also for engineering design. Coupling between surge and heave motions of tension leg platform is such a problem. Here, tension leg platform motions are considered only in surge and heave degrees of freedom without pitch effect. The coupled term of heave is a nonlinear differential equation. Because the focus of this article is on this term, therefore, Duffing equation of motion in the surge direction is linearized. The wave forces are calculated using Airy’s wave theory and Morison’s equation, ignoring the diffraction effects. Current force also can be very important in dynamic analysis of tension leg platform. Because it affects the term of heave that is coupled with surge. It is shown that the effect of surge motion coupling on heave motion is very important in large displacement of surge motion in many sea states. The main result is that the coupling effects appeared in some frequencies such as heave and surge frequency, twice the frequency of wave, twice the natural surge frequency, and summation and difference of frequency of wave and surge frequency.


Sign in / Sign up

Export Citation Format

Share Document