Numerical Methods for a Class of Quadratic Matrix Equations

Quadratic matrix equations arise in many elds of scienti c computing and engineering applications.In this paper,we consider a class of quadratic matrix equations.Under a certain condition,we rst prove the existence of minimal nonnegative solution for this quadratic matrix equation,and then propose some numerical methods for solving it.Convergence analysis and numerical examples are given to verify the theories and the numerical methods of this paper.

Transfer matrix method for free and forced vibrations of multi-level functionally graded material stepped beams with different boundary conditions

Functionally graded materials(FGMs)are a novel class of composite materials that have attracted significant attention in the field of engineering due to their unique mechanical properties.This study aims to explore the dynamic behaviors of an FGM stepped beam with different boundary conditions based on an efficient solving method.Under the assumptions of the Euler-Bernoulli beam theory,the governing differential equations of an individual FGM beam are derived with Hamilton’s principle and decoupled via the separation-of-variable approach.Then,the free and forced vibrations of the FGM stepped beam are solved with the transfer matrix method(TMM).Two models,i.e.,a three-level FGM stepped beam and a five-level FGM stepped beam,are considered,and their natural frequencies and mode shapes are presented.To demonstrate the validity of the method in this paper,the simulation results by ABAQUS are also given.On this basis,the detailed parametric analyses on the frequencies and dynamic responses of the three-level FGM stepped beam are carried out.The results show the accuracy and efficiency of the TMM.

Global Epidemic of Ebola Virus Disease and the Importation Risk into China: An Assessment Based on the Risk Matrix Method

Objective To analyze the global epidemic status of the Ebola virus disease(EVD) and assess the importation risk into China.Methods Data from World Health Organization reports were used. We described the global epidemic status of EVD from 1976–2021, and assessed and ranked the importation risk of EVD from the diseaseoutbreaking countries into China using the risk matrix and Borda count methods, respectively.Results From 1976–2021, EVD mainly occurred in western and central Africa, with the highest cumulative number of cases(14,124 cases) in Sierra Leone, and the highest cumulative fatality rate(85%) in the Congo. Outbreaks of EVD have occurred in the Democratic Republic of the Congo and Guinea since 2018. The importation risk into China varies across countries with outbreaks of disease.The Democratic Republic of the Congo had an extremely high risk(23 Borda points), followed by Guinea and Liberia. Countries with a moderate importation risk were Nigeria, Uganda, Congo, Sierra Leone,Mali, and Gabon, while countries with a low importation risk included Sudan, Senegal, and Co te d’Ivoire.Conclusion China is under the risk of EVD importation with the globalization and severe epidemic status of EVD. Key attention need to be paid to the Democratic Republic of the Congo, Guinea, and Liberia. Therefore, it is necessary to prevent and prepare in advance for importation risk in China.

An Extended Numerical Method by Stancu Polynomials for Solution of Integro-Differential Equations Arising in Oscillating Magnetic Fields

In this study, the Bernstein collocation method has been expanded to Stancu collocation method for numerical solution of the charged particle motion for certain configurations of oscillating magnetic fields modelled by a class of linear integro-differential equations. As the method has been improved, the Stancu polynomials that are generalization of the Bernstein polynomials have been used. The method has been tested on a physical problem how the method can be applied. Moreover, numerical results of the method have been compared with the numerical results of the other methods to indicate the efficiency of the method.

A Modified Eigenvector Method for Blind Deconvolution of MIMO Systems Using the Matrix Pseudo-Inversion Lemma

Recently we have developed an eigenvector method (EVM) which can achieve the blind deconvolution (BD) for MIMO systems. One of attractive features of the proposed algorithm is that the BD can be achieved by calculating the eigenvectors of a matrix relevant to it. However, the performance accuracy of the EVM depends highly on computational results of the eigenvectors. In this paper, by modifying the EVM, we propose an algorithm which can achieve the BD without calculating the eigenvectors. Then the pseudo-inverse which is needed to carry out the BD is calculated by our proposed matrix pseudo-inversion lemma. Moreover, using a combination of the conventional EVM and the modified EVM, we will show its performances comparing with each EVM. Simulation results will be presented for showing the effectiveness of the proposed methods.

Revealing the Hidden Mathematical Beauties of the Cayley-Hamilton Method

The inversion of a non-singular square matrix applying a Computer Algebra System (CAS) is straightforward. The CASs make the numeric computation efficient but mock the mathematical characteristics. The algorithms conducive to the output are sealed and inaccessible. In practice, other than the CPU timing, the applied inversion method is irrelevant. This research-oriented article discusses one such process, the Cayley-Hamilton (C.H.) [1]. Pursuing the process symbolically reveals its unpublished hidden mathematical characteristics even in the original article [1]. This article expands the general vision of the original named method without altering its practical applications. We have used the famous CAS Mathematica [2]. We have briefed the theory behind the method and applied it to different-sized symbolic and numeric matrices. The results are compared to the named CAS’s sealed, packaged library commands. The codes are given, and the algorithms are unsealed.

Application of the N + 2 Transversal Network Method to the Study of a Coupled Resonator Filter

This paper presents a new approach to synthesize admittance function polynomials and coupling matrices for coupled resonator filters. The N + 2 transversal network method is applied to study a coupled resonator filter. This method allowed us to determine the polynomials of the reflection and transmission coefficients. A study is made for a 4 poles filter with 2 transmission zeros between the N + 2 transversal network method and the one found in the literature. A MATLAB code was designed for the numerical simulation of these coefficients for the 6, 8, and 10 pole filter with 4 transmission zeros.

Another SSOR Iteration Method

Kellogg gave a version of the Peaceman-Radford method. In this paper, we introduce a SSOR iteration method which uses Kellogg’s method. The new algorithm has some advantages over the traditional SSOR algorithm. A Cyclic Reduction algorithm is introduced via a decoupling in Kellogg’s method.

GENERALIZED MATRIX MULTISPLITTING RELAXATION METHODS AND THEIR CONVERGENCE

In this paper, we set up a general framework of parallel matrix mullisplitting relaxation methods for solving large scale system of linear equations. We investigate the convergence properties of this framework and give several sufficient conditions ensuring it to converge as well as diverge. At last, we conclude a necessary and sufficient condition for the convergence of this framework when the coefficient matrix is an L-matrix.

Transverse free vibration analysis of a tapered Timoshenko beam on visco-Pasternak foundations using the interpolating matrix method

The characteristics of transverse free vibration of a tapered Timoshenko beam under an axially conservative compression resting on visco-Pasternak foundations are investigated by the interpolating matrix method. The research is executed in view of a three-parameter foundation which includes the eff ects of the Winkler coeffi cient, Pasternak coeffi cient and damping coeffi cient of the elastic medium. The governing equations of free vibration of a non-prismatic Timoshenko beam under an axially conservative force resting on visco-Pasternak foundations are transformed into ordinary diff erential equations with variable coeffi cients in light of the bending rotation angle and transverse displacement. All the natural frequencies orders together with the corresponding mode shapes of the beam are calculated at the same time, and a good convergence and accuracy of the proposed method is verifi ed through two numerical examples. The infl uences of foundation mechanical characteristics together with rotary inertia and shear deformation on natural frequencies of the beam with diff erent taper ratios are analyzed. A comprehensive parametric numerical study is carried out emphasizing the primary parameters that describe the dynamic property of the beam.

Transfer matrix method for determination of the natural vibration characteristics of elastically coupled launch vehicle boosters

The analysis of natural vibration characteristics has become one of important steps of the manufacture and dynamic design in the aerospace industry. This paper presents a new scenario called virtual cutting in the context of the transfer matrix method of linear multibody systems closed- loop topology for computing the free vibration characteristics of elastically coupled flexible launch vehicle boosters. In this approach, the coupled system is idealized as a triple-beam system-like structure coupled by linear translational springs, where a non-uniform free-free Euler-Bemoulli beam is used. A large thrust-to-weight ratio leads to large axial accelera- tions that result in an axial inertia load distribution from nose to tail. Consequently, it causes the development of significant compressive forces along the length of the launch vehicle. Therefore, it is important to take into account this effect in the transverse vibration model. This scenario does not need the global dynamics equations of a system, and it has high computational efficiency and low memory requirements. The validity of the presented scenario is achieved through com- parison to other approaches published in the literature.

Discrete time transfer matrix method for dynamics of multibody system with flexible beams moving in space

In this paper, by defining new state vectors and developing new transfer matrices of various elements mov- ing in space, the discrete time transfer matrix method of multi-rigid-flexible-body system is expanded to study the dynamics of multibody system with flexible beams moving in space. Formulations and numerical example of a rigid- flexible-body three pendulums system moving in space are given to validate the method. Using the new method to study the dynamics of multi-rigid-flexible-body system mov- ing in space, the global dynamics equations of system are not needed, the orders of involved matrices of the system are very low and the computational speed is high, irrespec- tive of the size of the system. The new method is simple, straightforward, practical, and provides a powerful tool for multi-rigid-flexible-body system dynamics.

Carbon matrix effects on the micro-structure and performance of Pt nanowire cathode prepared by decal transfer method

High performance cathode for polymer electrolyte membrane fuel cell was prepared by depositing Pt nanowires in a carbon matrix coated on a substrate, and using decal transfer method to fabricate the membrane electrode assembly. The effects of carbon and ionomer contents on the electrode micro-structure and fuel cell performance are investigated by physical characterization and single cell testing. The Pt nanowires are gradient distributed across the cathode thickness, and more Pt exists near the membrane. Both the carbon and ionomer contents can affect the Pt nanowires distribution and aggregation. In addition, the carbon loading dominates the transport distance of gas and proton, and the ionomer content affects the triple phase boundaries and porosity in the cathode. The optimal structure of Pt nanowire cathode is obtained at 0.10 mg·cm^-2 carbon loading and 10 wt% ionomer.

Determination of initial cable force of cantilever casting concrete arch bridge using stress balance and influence matrix methods

Cantilever casting concrete arch bridge using form traveller has a broad application prospect.However,it is difficult to obtain reasonable initial cable force in construction stage.In this study,stress balance and influence matrix methods were developed to determine the initial cable force of cantilever casting concrete arch bridge.The stress balance equation and influence matrix of arch rib critical section were established,and the buckle cable force range was determined by the allowable stress of arch rib critical section.Then a group of buckle cable forces were selected and substituted into the stress balance equation,and the reasonable initial buckle cable force was determined through iteration.Based on the principle of force balance,the initial anchor cable force was determined.In an engineering application example,it is shown that the stress balance and influence matrix methods for the determination of initial cable force are feasible and reliable.The initial cable forces of arch rib segments only need to be adjusted once in the corresponding construction process,which improves the working efficiency and reduces the construction risk.It is found that the methods have great advantages for determining initial cable force in cantilever casting construction process of concrete arch bridge.

SiC nanoparticles reinforced magnesium matrix composites fabricated by ultrasonic method

SiC nanoparticles reinforced magnesium matrix composites were fabricated by ultrasonic method.The AZ91 alloy and SiC nanoparticles with the average diameter of 50 nm were used as the matrix alloy and the reinforcement,respectively.The addition of nanoparticles was 0.1%,0.3%,and 0.5%(mass fraction) of the composites.The results of microstructural evaluation and mechanical properties indicate that the nanoparticles can be dispersed into magnesium alloys efficiently and uniformly with the aid of ultrasonic vibration.As compared with the matrix alloys,the grains of composites were refined and the mechanical properties of composites were improved significantly.The SEM and DSC analyses show that the SiC nanoparticles can act as the heterogeneous nucleation of α-Mg.Also,the strengthening mechanism responsible for the composites reinforced with SiC nanoparticles was discussed.

Precise integration method without inverse matrix calculation for structural dynamic equations

The precise integration method proposed for linear time-invariant homogeneous dynamic systems can provide accurate numerical results that approach an exact solution at integration points. However, difficulties arise when the algorithm is used for non-homogeneous dynamic systems due to the inverse matrix calculation required. In this paper, the structural dynamic equalibrium equations are converted into a special form, the inverse matrix calculation is replaced by the Crout decomposition method to solve the dynamic equilibrium equations, and the precise integration method without the inverse matrix calculation is obtained. The new algorithm enhances the present precise integration method by improving both the computational accuracy and efficiency. Two numerical examples are given to demonstrate the validity and efficiency of the proposed algorithm.

Method of reverberation ray matrix for static analysis of planar framed structures composed of anisotropic Timoshenko beam members

Based on the method of reverberation ray matrix（MRRM）, a reverberation matrix for planar framed structures composed of anisotropic Timoshenko（T） beam members containing completely hinged joints is developed for static analysis of such structures.In the MRRM for dynamic analysis, amplitudes of arriving and departing waves for joints are chosen as unknown quantities. However, for the present case of static analysis, displacements and rotational angles at the ends of each beam member are directly considered as unknown quantities. The expressions for stiffness matrices for anisotropic beam members are developed. A corresponding reverberation matrix is derived analytically for exact and unified determination on the displacements and internal forces at both ends of each member and arbitrary cross sectional locations in the structure. Numerical examples are given and compared with the finite element method（FEM） results to validate the present model. The characteristic parameter analysis is performed to demonstrate accuracy of the present model with the T beam theory in contrast with errors in the usual model based on the Euler-Bernoulli（EB） beam theory. The resulting reverberation matrix can be used for exact calculation of anisotropic framed structures as well as for parameter analysis of geometrical and material properties of the framed structures.

Dynamics Analysis of a Parallel Mill-turn Tool Spindle Head Driven by Dual-linear Motors Using Extended Transfer Matrix Method

The hybrid dynamics of multi-rigid-body and multi-flexible-body system becomes the mainstream of multi-body dynamics.Currently there lacks a compact approach to model the hybrid dynamics,especially in modern machine tool application,due to the difficulty of solving the hybrid equations or the limitation of current software when dealing with the hybrid dynamics.The extended transfer matrix method（E-TMM）,which extends elements in three-dimensional space with higher matrixes,is proposed to simplify the modeling process of the hybrid dynamics.The E-TMM modeling approaches of 3 basic elements including 3D vibrant rigid body,joint and flexible body are studied in details.A parallel mill-turn tool spindle head unit driven by dual-linear motors is chosen as a plant to demonstrate the E-TMM modeling process.By using E-TMM,the spindle head unit is simplified as a topological network consisting of the three types of element,i.e.,3D vibrant rigid body,joint and flexible body,including 11 rigid bodies,14 joints and 1 3D-Timoshenko beam.Then the dynamic model of the system can be easily obtained by deducing the element-network by means of state vector transformation.The dynamic characteristics of the spindle head,such as natural frequencies,dynamic flexibility,etc.can be predicted by solving the obtained model.Experiment verification indicates that the E-TMM is valid with enough accuracy in the dynamic analysis of the parallel mill-turn tool spindle head.The E-TMM is capable of modeling the dynamics of machine tool structure with no requirements of deducing and solving the sophisticated differential equations.Moreover,the E-TMM provides a simple and elegant tool for hybrid dynamic analysis in future dynamic design of machine tools.

Accelerated Matrix Recovery via Random Projection Based on Inexact Augmented Lagrange Multiplier Method

In this paper, a unified matrix recovery model was proposed for diverse corrupted matrices. Resulting from the separable structure of the proposed model, the convex optimization problem can be solved efficiently by adopting an inexact augmented Lagrange multiplier (IALM) method. Additionally, a random projection accelerated technique (IALM+RP) was adopted to improve the success rate. From the preliminary numerical comparisons, it was indicated that for the standard robust principal component analysis (PCA) problem, IALM+RP was at least two to six times faster than IALM with an insignificant reduction in accuracy; and for the outlier pursuit (OP) problem, IALM+RP was at least 6.9 times faster, even up to 8.3 times faster when the size of matrix was 2 000×2 000.

MatDEM-fast matrix computing of the discrete element method

Discrete element method can effectively simulate the discontinuity,inhomogeneity and large deformation and failure of rock and soil.Based on the innovative matrix computing of the discrete element method,the highperformance discrete element software MatDEM may handle millions of elements in one computer,and enables the discrete element simulation at the engineering scale.It supports heat calculation,multi-field and fluidsolid coupling numerical simulations.Furthermore,the software integrates pre-processing,solver,postprocessing,and powerful secondary development,allowing recompiling new discrete element software.The basic principles of the DEM,the implement and development of the MatDEM software,and its applications are introduced in this paper.The software and sample source code are available online(http://matdem.com).