Efficient Solution of Maggi’s Equations

Author(s):  
Javier Garci´a de Jalo´n ◽  
Alfonso Callejo ◽  
Andre´s F. Hidalgo ◽  
Mari´a D. Gutie´rrez

According to a recent paper by Laulusa and Bauchau [1], Maggi’s formulation is a simple and stable way to solve the dynamic equations of constrained multibody systems. Among the difficulties of Maggi’s formulation, Laulusa and Bauchau quoted the need for an appropriate choice (and change, when necessary) of independent coordinates, as well as the high cost of computing and updating the basis of the tangent null space of constraint equations. In this paper, index-1 Lagrange’s equations are first considered, including the not-so-rare case of having a singular mass matrix and redundant constraints. The existence and uniqueness of solution for acceleration vector and Lagrange multipliers vector is studied in a very simple way. Then, following Von Schwerin [2], Maggi’s formulation is described as the most efficient way (globally speaking) to solve these index-1 equations. Next, an improved double-step method, which implements the matrix transformations of Maggi’s formulation in an efficient way, is described. Finally, two large real-life examples are presented.

Author(s):  
Javier García de Jalón ◽  
Alfonso Callejo ◽  
Andrés F. Hidalgo

According to a recent paper (Laulusa and Bauchau, 2008, “Review of Classical Approaches for Constraint Enforcement in Multibody Systems,” ASME J. Comput. Nonlinear Dyn., 3(1), 011004), Maggi’s formulation is a simple and stable way to solve the dynamic equations of constrained multibody systems. Among the difficulties of Maggi’s formulation, Laulusa and Bauchau quoted the need for an appropriate choice (and change, when necessary) of independent coordinates, as well as the high cost of computing and updating the basis of the tangent null space of constraint equations. In this paper, index-1 Lagrange’s equations are first considered, including the not-so-rare case of having a singular mass matrix and redundant constraints. The existence and uniqueness of solution for acceleration vector and Lagrange multipliers vector is studied in a very simple way. Then, following Von Schwerin (Von Schwerin, Multibody System Simulation. Numerical Methods, Algorithms and Software, Springer, New York, 1999), Maggi’s formulation is described as the most efficient way (in general) to solve these index-1 equations. Next, an improved double-step method, which implements the matrix transformations of Maggi’s formulation in an efficient way, is described. Finally, two large real-life examples are presented.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Hai-Fei Zhu ◽  
Xiao-Wei Sun ◽  
Ting Song ◽  
Xiao-Dong Wen ◽  
Xi-Xuan Liu ◽  
...  

AbstractIn view of the influence of variability of low-frequency noise frequency on noise prevention in real life, we present a novel two-dimensional tunable phononic crystal plate which is consisted of lead columns deposited in a silicone rubber plate with periodic holes and calculate its bandgap characteristics by finite element method. The low-frequency bandgap mechanism of the designed model is discussed simultaneously. Accordingly, the influence of geometric parameters of the phononic crystal plate on the bandgap characteristics is analyzed and the bandgap adjustability under prestretch strain is further studied. Results show that the new designed phononic crystal plate has lower bandgap starting frequency and wider bandwidth than the traditional single-sided structure, which is due to the coupling between the resonance mode of the scatterer and the long traveling wave in the matrix with the introduction of periodic holes. Applying prestretch strain to the matrix can realize active realtime control of low-frequency bandgap under slight deformation and broaden the low-frequency bandgap, which can be explained as the multiple bands tend to be flattened due to the localization degree of unit cell vibration increases with the rise of prestrain. The presented structure improves the realtime adjustability of sound isolation and vibration reduction frequency for phononic crystal in complex acoustic vibration environments.


Entropy ◽  
2021 ◽  
Vol 23 (4) ◽  
pp. 384
Author(s):  
Rocío Hernández-Sanjaime ◽  
Martín González ◽  
Antonio Peñalver ◽  
Jose J. López-Espín

The presence of unaccounted heterogeneity in simultaneous equation models (SEMs) is frequently problematic in many real-life applications. Under the usual assumption of homogeneity, the model can be seriously misspecified, and it can potentially induce an important bias in the parameter estimates. This paper focuses on SEMs in which data are heterogeneous and tend to form clustering structures in the endogenous-variable dataset. Because the identification of different clusters is not straightforward, a two-step strategy that first forms groups among the endogenous observations and then uses the standard simultaneous equation scheme is provided. Methodologically, the proposed approach is based on a variational Bayes learning algorithm and does not need to be executed for varying numbers of groups in order to identify the one that adequately fits the data. We describe the statistical theory, evaluate the performance of the suggested algorithm by using simulated data, and apply the two-step method to a macroeconomic problem.


2020 ◽  
Vol 20 (1) ◽  
Author(s):  
Vincent Vandewalle ◽  
Alexandre Caron ◽  
Coralie Delettrez ◽  
Renaud Périchon ◽  
Sylvia Pelayo ◽  
...  

Abstract Background Usability testing of medical devices are mandatory for market access. The testings’ goal is to identify usability problems that could cause harm to the user or limit the device’s effectiveness. In practice, human factor engineers study participants under actual conditions of use and list the problems encountered. This results in a binary discovery matrix in which each row corresponds to a participant, and each column corresponds to a usability problem. One of the main challenges in usability testing is estimating the total number of problems, in order to assess the completeness of the discovery process. Today’s margin-based methods fit the column sums to a binomial model of problem detection. However, the discovery matrix actually observed is truncated because of undiscovered problems, which corresponds to fitting the marginal sums without the zeros. Margin-based methods fail to overcome the bias related to truncation of the matrix. The objective of the present study was to develop and test a matrix-based method for estimating the total number of usability problems. Methods The matrix-based model was based on the full discovery matrix (including unobserved columns) and not solely on a summary of the data (e.g. the margins). This model also circumvents a drawback of margin-based methods by simultaneously estimating the model’s parameters and the total number of problems. Furthermore, the matrix-based method takes account of a heterogeneous probability of detection, which reflects a real-life setting. As suggested in the usability literature, we assumed that the probability of detection had a logit-normal distribution. Results We assessed the matrix-based method’s performance in a range of settings reflecting real-life usability testing and with heterogeneous probabilities of problem detection. In our simulations, the matrix-based method improved the estimation of the number of problems (in terms of bias, consistency, and coverage probability) in a wide range of settings. We also applied our method to five real datasets from usability testing. Conclusions Estimation models (and particularly matrix-based models) are of value in estimating and monitoring the detection process during usability testing. Matrix-based models have a solid mathematical grounding and, with a view to facilitating the decision-making process for both regulators and device manufacturers, should be incorporated into current standards.


Author(s):  
Michel Mandjes ◽  
Birgit Sollie

AbstractThis paper considers a continuous-time quasi birth-death (qbd) process, which informally can be seen as a birth-death process of which the parameters are modulated by an external continuous-time Markov chain. The aim is to numerically approximate the time-dependent distribution of the resulting bivariate Markov process in an accurate and efficient way. An approach based on the Erlangization principle is proposed and formally justified. Its performance is investigated and compared with two existing approaches: one based on numerical evaluation of the matrix exponential underlying the qbd process, and one based on the uniformization technique. It is shown that in many settings the approach based on Erlangization is faster than the other approaches, while still being highly accurate. In the last part of the paper, we demonstrate the use of the developed technique in the context of the evaluation of the likelihood pertaining to a time series, which can then be optimized over its parameters to obtain the maximum likelihood estimator. More specifically, through a series of examples with simulated and real-life data, we show how it can be deployed in model selection problems that involve the choice between a qbd and its non-modulated counterpart.


Author(s):  
Martin M. Tong

Numerical solution of the dynamics equations of a flexible multibody system as represented by Hamilton’s canonical equations requires that its generalized velocities q˙ be solved from the generalized momenta p. The relation between them is p = J(q)q˙, where J is the system mass matrix and q is the generalized coordinates. This paper presents the dynamics equations for a generic flexible multibody system as represented by p˙ and gives emphasis to a systematic way of constructing the matrix J for solving q˙. The mass matrix is shown to be separable into four submatrices Jrr, Jrf, Jfr and Jff relating the joint momenta and flexible body mementa to the joint coordinate rates and the flexible body deformation coordinate rates. Explicit formulas are given for these submatrices. The equations of motion presented here lend insight to the structure of the flexible multibody dynamics equations. They are also a versatile alternative to the acceleration-based dynamics equations for modeling mechanical systems.


Mathematics ◽  
2018 ◽  
Vol 6 (11) ◽  
pp. 268 ◽  
Author(s):  
Kuddusi Kayaduman ◽  
Fevzi Yaşar

In 1978, the domain of the Nörlund matrix on the classical sequence spaces lp and l∞ was introduced by Wang, where 1 ≤ p < ∞. Tuğ and Başar studied the matrix domain of Nörlund mean on the sequence spaces f0 and f in 2016. Additionally, Tuğ defined and investigated a new sequence space as the domain of the Nörlund matrix on the space of bounded variation sequences in 2017. In this article, we defined new space and and examined the domain of the Nörlund mean on the bs and cs, which are bounded and convergent series, respectively. We also examined their inclusion relations. We defined the norms over them and investigated whether these new spaces provide conditions of Banach space. Finally, we determined their α­, β­, γ­duals, and characterized their matrix transformations on this space and into this space.


2021 ◽  
Vol 71 (6) ◽  
pp. 1375-1400
Author(s):  
Feyzi Başar ◽  
Hadi Roopaei

Abstract Let F denote the factorable matrix and X ∈ {ℓp , c 0, c, ℓ ∞}. In this study, we introduce the domains X(F) of the factorable matrix in the spaces X. Also, we give the bases and determine the alpha-, beta- and gamma-duals of the spaces X(F). We obtain the necessary and sufficient conditions on an infinite matrix belonging to the classes (ℓ p (F), ℓ ∞), (ℓ p (F), f) and (X, Y(F)) of matrix transformations, where Y denotes any given sequence space. Furthermore, we give the necessary and sufficient conditions for factorizing an operator based on the matrix F and derive two factorizations for the Cesàro and Hilbert matrices based on the Gamma matrix. Additionally, we investigate the norm of operators on the domain of the matrix F. Finally, we find the norm of Hilbert operators on some sequence spaces and deal with the lower bound of operators on the domain of the factorable matrix.


2021 ◽  
Vol 0 (0) ◽  
Author(s):  
Kirill Svit ◽  
Konstantin Zhuravlev ◽  
Sergey Kireev ◽  
Karl K. Sabelfeld

Abstract A stochastic model of nanocrystals clusters formation is developed and applied to simulate an aggregation of cadmium sulfide nanocrystals upon evaporation of the Langmuir–Blodgett matrix. Simulations are compared with our experimental results. The stochastic model suggested governs mobilities both of individual nanocrystals and its clusters (arrays). We give a comprehensive analysis of the patterns simulated by the model, and study an influence of the surrounding medium (solvent) on the aggregation processes. In our model, monomers have a finite probability of separation from the cluster which depends on the temperature and binding energy between nanocrystals, and can also be redistributed in the composition of the cluster, leading to its compaction. The simulation results obtained in this work are compared with the experimental data on the aggregation of CdS nanocrystals upon evaporation of the Langmuir–Blodgett matrix. This system is a typical example from real life and is noteworthy in that the morphology of nanocrystals after evaporation of the matrix cannot be described exactly by a model based only on the motion of individual nanocrystals or by a cluster-cluster aggregation model.


1995 ◽  
Vol 398 ◽  
Author(s):  
A.D. Sequeira ◽  
H.A. Calderon ◽  
G. Kostorz

ABSTRACTThe influence of coherency strains produced by the γ-γ’ lattice mismatch, δ, on the decomposition process of Ni-Al-Mo alloys with a bimodal size distribution is presented. Samples with δ ranging from positive to negative, were investigated in a double-step aging procedure. The evolution of the microstructure and the kinetics of coarsening were studied using transmission electron microscopy (TEM). The lattice mismatch between the matrix and the different classes of precipitates was determined by high-resolution high-temperature x-ray diffraction. It is shown that the strain fields produced by the lattice mismatch can influence dramatically the decomposition of metallic alloys. It is suggested that the reduction of the coarsening rate of the large precipitates, the fast coarsening rate of the small precipitates and the distortions detected in the matrix are all direct consequences of the elastic fields produced by the γ-γ’ lattice mismatch.


Sign in / Sign up

Export Citation Format

Share Document