In this paper,we consider the maximal positive definite solution of the nonlinear matrix equation.By using the idea of Algorithm 2.1 in ZHANG(2013),a new inversion-free method with a stepsize parameter is proposed to ...In this paper,we consider the maximal positive definite solution of the nonlinear matrix equation.By using the idea of Algorithm 2.1 in ZHANG(2013),a new inversion-free method with a stepsize parameter is proposed to obtain the maximal positive definite solution of nonlinear matrix equation X+A^(*)X|^(-α)A=Q with the case 0<α≤1.Based on this method,a new iterative algorithm is developed,and its convergence proof is given.Finally,two numerical examples are provided to show the effectiveness of the proposed method.展开更多
This study introduces a comprehensive theoretical framework for accurately calculating the electronic band-structure of strained long-wavelength InAs/GaSb type-Ⅱsuperlattices.Utilizing an eight-band k·p Hamilto⁃...This study introduces a comprehensive theoretical framework for accurately calculating the electronic band-structure of strained long-wavelength InAs/GaSb type-Ⅱsuperlattices.Utilizing an eight-band k·p Hamilto⁃nian in conjunction with a scattering matrix method,the model effectively incorporates quantum confinement,strain effects,and interface states.This robust and numerically stable approach achieves exceptional agreement with experimental data,offering a reliable tool for analyzing and engineering the band structure of complex multi⁃layer systems.展开更多
Transportation structures such as composite pavements and railway foundations typically consist of multi-layered media designed to withstand high bearing capacity.A theoretical understanding of load transfer mechanism...Transportation structures such as composite pavements and railway foundations typically consist of multi-layered media designed to withstand high bearing capacity.A theoretical understanding of load transfer mechanisms in these multi-layer composites is essential,as it offers intuitive insights into parametric influences and facilitates enhanced structural performance.This paper employs an improved transfer matrix method to address the limitations of existing theoretical approaches for analyzing multi-layer composite structures.By establishing a twodimensional composite pavement model,it investigates load transfer characteristics and validates the accuracy through finite element simulation.The proposed method offers a straightforward analytical approach for examining internal interactions between structural layers.Case studies indicate that the concrete surface layer is the main load-bearing layer for most vertical normal and shear stresses.The soil base layer reduces the overall mechanical response of the substructure,while horizontal actions increase the risk of interfacial slip and cracking.Structural optimization analysis demonstrates that increasing the thickness of the concrete surface layer,enhancing the thickness and stiffness of the soil base layer,or incorporating gradient layers can significantly mitigate these risks of interfacial slip and cracking.The findings of this study can guide the optimization design,parameter analysis,and damage prevention of multi-layer composite structures.展开更多
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 giv...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.展开更多
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 r...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.展开更多
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 con...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.展开更多
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 inf...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.展开更多
The first order differential matrix equations of the host shell and constrained layer for a sandwich rotational shell are derived based on the thin shell theory.Employing the layer wise principle and first order shear...The first order differential matrix equations of the host shell and constrained layer for a sandwich rotational shell are derived based on the thin shell theory.Employing the layer wise principle and first order shear deformation theory, only considering the shearing deformation of the viscoelastic layer, the integrated first order differential matrix equation of a passive constrained layer damping rotational shell is established by combining with the normal equilibrium equation of the viscoelastic layer.A highly precise transfer matrix method is developed by extended homogeneous capacity precision integration technology.The numerical results show that present method is accurate and effective.展开更多
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 stud...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.展开更多
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...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.展开更多
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...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.展开更多
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 w...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.展开更多
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 st...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.展开更多
The bi-conjugate gradients(Bi-CG)and bi-conjugate residual(Bi-CR)methods are powerful tools for solving nonsymmetric linear systems Ax=b.By using Kronecker product and vectorization operator,this paper develops the Bi...The bi-conjugate gradients(Bi-CG)and bi-conjugate residual(Bi-CR)methods are powerful tools for solving nonsymmetric linear systems Ax=b.By using Kronecker product and vectorization operator,this paper develops the Bi-CG and Bi-CR methods for the solution of the generalized Sylvester-transpose matrix equationp i=1(Ai X Bi+Ci XTDi)=E(including Lyapunov,Sylvester and Sylvester-transpose matrix equations as special cases).Numerical results validate that the proposed algorithms are much more efcient than some existing algorithms.展开更多
A large database is desired for machine learning(ML) technology to make accurate predictions of materials physicochemical properties based on their molecular structure.When a large database is not available,the develo...A large database is desired for machine learning(ML) technology to make accurate predictions of materials physicochemical properties based on their molecular structure.When a large database is not available,the development of proper featurization method based on physicochemical nature of target proprieties can improve the predictive power of ML models with a smaller database.In this work,we show that two new featurization methods,volume occupation spatial matrix and heat contribution spatial matrix,can improve the accuracy in predicting energetic materials' crystal density(ρ_(crystal)) and solid phase enthalpy of formation(H_(f,solid)) using a database containing 451 energetic molecules.Their mean absolute errors are reduced from 0.048 g/cm~3 and 24.67 kcal/mol to 0.035 g/cm~3 and 9.66 kcal/mol,respectively.By leave-one-out-cross-validation,the newly developed ML models can be used to determine the performance of most kinds of energetic materials except cubanes.Our ML models are applied to predict ρ_(crystal) and H_(f,solid) of CHON-based molecules of the 150 million sized PubChem database,and screened out 56 candidates with competitive detonation performance and reasonable chemical structures.With further improvement in future,spatial matrices have the potential of becoming multifunctional ML simulation tools that could provide even better predictions in wider fields of materials science.展开更多
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 ad...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.展开更多
Objective To analyze the 2022 multiple-country monkeypox outbreak and assess its importation risk into China.Methods Data was from United States Centers for Disease Control and Prevention.We described the global epide...Objective To analyze the 2022 multiple-country monkeypox outbreak and assess its importation risk into China.Methods Data was from United States Centers for Disease Control and Prevention.We described the global epidemic and calculated concentration index to measure economic-related inequality.Importation risk into China was evaluated and ranked by using risk matrix method and Borda count method,respectively.Results As of July 29,2022,of 79 countries or territories,39(49.37%,39/79),17(21.52%,17/79),6(7.59%,6/79),12(15.19%,12/79),and 5(6.33%,5/79)country or territories identified cases<10,10-,51-,101-,and>1,000.There were economic-related health disparities exist in the distribution of cases(the concentration index=0.42,P=0.027),and the inequality disadvantageous to the rich(pro-poor).There were 12(15.38%,12/78),15(19.23%,15/78),6(7.69%,6/78),and 45(57.69%,45/78)countries or territories with extremely high,high,moderate,and low importation risk.United States and France ranked first with the highest Borda points of 156,and counts of zero.Conclusion Of 78 countries or territories,the key attention need be paid to the United States and France,relatively.As the epidemic progresses,preparing prevention and control measures to further reduce importation risk was crucial.展开更多
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 highperfo...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).展开更多
To reduce the uncertainty and reworks in complex projects,a novel mechanism is systematically developed in this paper based on two classical design structure matrix(DSM)clustering methods:Loop searching method(LSM)and...To reduce the uncertainty and reworks in complex projects,a novel mechanism is systematically developed in this paper based on two classical design structure matrix(DSM)clustering methods:Loop searching method(LSM)and function searching method(FSM).Specifically,the optimal working areas for the two clustering methods are first obtained quantitatively in terms of non-zero fraction(NZF)and singular value modularity index(SMI),in which the whole working area is divided into six sub-zones.Then,a judgement procedure is proposed for conveniently choosing the optimal DSM clustering method,which makes it easy to determine which DSM clustering method performs better for a given case.Subsequently,a conceptual model is constructed to assist project managers in effectively analyzing the network of projects and greatly reducing reworks in complex projects by defining preventive actions.Finally,the aircraft design process is presented to show how the proposed judgement mechanism can be utilized to reduce the reworks in actual projects.展开更多
The following is proved: 1) The linear independence of assumed stress modes is the necessary and sufficient condition for the nonsingular flexibility matrix; 2) The equivalent assumed stress modes lead to the identica...The following is proved: 1) The linear independence of assumed stress modes is the necessary and sufficient condition for the nonsingular flexibility matrix; 2) The equivalent assumed stress modes lead to the identical hybrid element. The Hilbert stress subspace of the assumed stress modes is established. So, it is easy to derive the equivalent orthogonal normal stress modes by Schmidt's method. Because of the resulting diagonal flexibility matrix, the identical hybrid element is free from the complex matrix inversion so that the hybrid efficiency, is improved greatly. The numerical examples show that the method is effective.展开更多
基金Supported in part by Natural Science Foundation of Guangxi(2023GXNSFAA026246)in part by the Central Government's Guide to Local Science and Technology Development Fund(GuikeZY23055044)in part by the National Natural Science Foundation of China(62363003)。
文摘In this paper,we consider the maximal positive definite solution of the nonlinear matrix equation.By using the idea of Algorithm 2.1 in ZHANG(2013),a new inversion-free method with a stepsize parameter is proposed to obtain the maximal positive definite solution of nonlinear matrix equation X+A^(*)X|^(-α)A=Q with the case 0<α≤1.Based on this method,a new iterative algorithm is developed,and its convergence proof is given.Finally,two numerical examples are provided to show the effectiveness of the proposed method.
文摘This study introduces a comprehensive theoretical framework for accurately calculating the electronic band-structure of strained long-wavelength InAs/GaSb type-Ⅱsuperlattices.Utilizing an eight-band k·p Hamilto⁃nian in conjunction with a scattering matrix method,the model effectively incorporates quantum confinement,strain effects,and interface states.This robust and numerically stable approach achieves exceptional agreement with experimental data,offering a reliable tool for analyzing and engineering the band structure of complex multi⁃layer systems.
基金supported by Fundamental Research Funds for the Central Universities(No.lzujbky-2024-05)Innovation Foundation of Provincial Education Department of Gansu(2024B-005)+2 种基金Scientific Department of Gansu(24CXGA083,24CXGA024,JK2024-28,JK2024-32 and 23CXJA0007)Industrial Support Plan Project of Provincial Education Department of Gansu(2025CYZC-003 and CYZC-2024-10)the Hunan Natural Science Foundation Science and Education Joint Fund Project(2022JJ60109).
文摘Transportation structures such as composite pavements and railway foundations typically consist of multi-layered media designed to withstand high bearing capacity.A theoretical understanding of load transfer mechanisms in these multi-layer composites is essential,as it offers intuitive insights into parametric influences and facilitates enhanced structural performance.This paper employs an improved transfer matrix method to address the limitations of existing theoretical approaches for analyzing multi-layer composite structures.By establishing a twodimensional composite pavement model,it investigates load transfer characteristics and validates the accuracy through finite element simulation.The proposed method offers a straightforward analytical approach for examining internal interactions between structural layers.Case studies indicate that the concrete surface layer is the main load-bearing layer for most vertical normal and shear stresses.The soil base layer reduces the overall mechanical response of the substructure,while horizontal actions increase the risk of interfacial slip and cracking.Structural optimization analysis demonstrates that increasing the thickness of the concrete surface layer,enhancing the thickness and stiffness of the soil base layer,or incorporating gradient layers can significantly mitigate these risks of interfacial slip and cracking.The findings of this study can guide the optimization design,parameter analysis,and damage prevention of multi-layer composite structures.
基金Supported by Natural Science Fundations of China and Shanghai.
文摘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.
基金AHKJT of China under Grant Nos.1708085QE121 and 1808085ME147AHEDU of China under Grant No.TSKJ2017B13
文摘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.
基金supported by the Research Fund for the Doctoral Program of Higher Education of China(Grants 20113219110025,20133219110037)the National Natural Science Foundation of China(Grants 11102089,61304137)the Program for New Century Excellent Talents in University(NCET-10-0075)
文摘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.
基金Projects(51478049,51778068)supported by the National Natural Science Foundation of ChinaProject(14JJ2075,2019JJ40301)supported by the Hunan Natural Science Foundation of China+1 种基金Project(17A010)supported by the Scientific Research Fund of Hunan Provincial Education Department of ChinaProject(2017GK4034)supported by the Major Technological Achievements Transformation Program of Hunan Strategic Emerging Industries of China
文摘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.
基金supported by the National Natural Science Foundation of China (No.10662003)Educational Commission of Guangxi Province of China (No.200807MS109)
文摘The first order differential matrix equations of the host shell and constrained layer for a sandwich rotational shell are derived based on the thin shell theory.Employing the layer wise principle and first order shear deformation theory, only considering the shearing deformation of the viscoelastic layer, the integrated first order differential matrix equation of a passive constrained layer damping rotational shell is established by combining with the normal equilibrium equation of the viscoelastic layer.A highly precise transfer matrix method is developed by extended homogeneous capacity precision integration technology.The numerical results show that present method is accurate and effective.
基金supported by the Natural Science Foundation of China Government (10902051)the Natural Science Foundation of Jiangsu Province (BK2008046)the German Science Foundation
文摘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.
文摘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.
基金Project(2007CB613706) supported by the National Basic Research Program of ChinaProject(00900054R4001) supported by Innovation Project for Talents of BJUTProject(00900054K4004) supported by the Science Foundation for Youths of BJUT
文摘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.
文摘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.
基金Project supported by the Program for New Century Excellent Talents in Universities(NCET)by the Ministry of Education of China(No.NCET-04-0373)
文摘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.
文摘The bi-conjugate gradients(Bi-CG)and bi-conjugate residual(Bi-CR)methods are powerful tools for solving nonsymmetric linear systems Ax=b.By using Kronecker product and vectorization operator,this paper develops the Bi-CG and Bi-CR methods for the solution of the generalized Sylvester-transpose matrix equationp i=1(Ai X Bi+Ci XTDi)=E(including Lyapunov,Sylvester and Sylvester-transpose matrix equations as special cases).Numerical results validate that the proposed algorithms are much more efcient than some existing algorithms.
基金support from the Ministry of Education(MOE) Singapore Tier 1 (RG8/20)。
文摘A large database is desired for machine learning(ML) technology to make accurate predictions of materials physicochemical properties based on their molecular structure.When a large database is not available,the development of proper featurization method based on physicochemical nature of target proprieties can improve the predictive power of ML models with a smaller database.In this work,we show that two new featurization methods,volume occupation spatial matrix and heat contribution spatial matrix,can improve the accuracy in predicting energetic materials' crystal density(ρ_(crystal)) and solid phase enthalpy of formation(H_(f,solid)) using a database containing 451 energetic molecules.Their mean absolute errors are reduced from 0.048 g/cm~3 and 24.67 kcal/mol to 0.035 g/cm~3 and 9.66 kcal/mol,respectively.By leave-one-out-cross-validation,the newly developed ML models can be used to determine the performance of most kinds of energetic materials except cubanes.Our ML models are applied to predict ρ_(crystal) and H_(f,solid) of CHON-based molecules of the 150 million sized PubChem database,and screened out 56 candidates with competitive detonation performance and reasonable chemical structures.With further improvement in future,spatial matrices have the potential of becoming multifunctional ML simulation tools that could provide even better predictions in wider fields of materials science.
基金Supported by National Natural Science Foundation of China (No.51275348)College Students Innovation and Entrepreneurship Training Program of Tianjin University (No.201210056339)
文摘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.
基金partly supported by the National Natural Science Foundation of China[grant numbers 72122001,71934002]the National Key Research and Development Project of China[grant number 2021ZD0114101,2021ZD0114104,2021ZD0114105]+2 种基金National Statistical Science Research Project[grant numbers 2021LY038]the Fundamental Research Funds for the Central Universities supported by Global Center for Infectious Disease and Policy Research&Global Health and Infectious Diseases Group,of Peking University[grant number 202204]National Science and Technology Project on Development Assistance for Technology,Developing China-ASEAN Public Health Research and Development Collaborating Center[grant number KY202101004]
文摘Objective To analyze the 2022 multiple-country monkeypox outbreak and assess its importation risk into China.Methods Data was from United States Centers for Disease Control and Prevention.We described the global epidemic and calculated concentration index to measure economic-related inequality.Importation risk into China was evaluated and ranked by using risk matrix method and Borda count method,respectively.Results As of July 29,2022,of 79 countries or territories,39(49.37%,39/79),17(21.52%,17/79),6(7.59%,6/79),12(15.19%,12/79),and 5(6.33%,5/79)country or territories identified cases<10,10-,51-,101-,and>1,000.There were economic-related health disparities exist in the distribution of cases(the concentration index=0.42,P=0.027),and the inequality disadvantageous to the rich(pro-poor).There were 12(15.38%,12/78),15(19.23%,15/78),6(7.69%,6/78),and 45(57.69%,45/78)countries or territories with extremely high,high,moderate,and low importation risk.United States and France ranked first with the highest Borda points of 156,and counts of zero.Conclusion Of 78 countries or territories,the key attention need be paid to the United States and France,relatively.As the epidemic progresses,preparing prevention and control measures to further reduce importation risk was crucial.
基金Financial supports from the Natural Science Foundation of China(41761134089,41977218)Six Talent Peaks Project of Jiangsu Province(RJFW-003)the Fundamental Research Funds for the Central Universities(14380103)are gratefully acknowledged.
文摘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).
基金supported by the National Natural Science Foundation of China (Nos. 71471087, 71071076, 61673209)the Funding for Outstanding Doctoral Dissertation in Nanjing University of Aeronautics and Astronautics (No. BCXJ17-11)the Research and Innovation Program for Graduate Education of Jiangsu Province (No. KYZZ160145)
文摘To reduce the uncertainty and reworks in complex projects,a novel mechanism is systematically developed in this paper based on two classical design structure matrix(DSM)clustering methods:Loop searching method(LSM)and function searching method(FSM).Specifically,the optimal working areas for the two clustering methods are first obtained quantitatively in terms of non-zero fraction(NZF)and singular value modularity index(SMI),in which the whole working area is divided into six sub-zones.Then,a judgement procedure is proposed for conveniently choosing the optimal DSM clustering method,which makes it easy to determine which DSM clustering method performs better for a given case.Subsequently,a conceptual model is constructed to assist project managers in effectively analyzing the network of projects and greatly reducing reworks in complex projects by defining preventive actions.Finally,the aircraft design process is presented to show how the proposed judgement mechanism can be utilized to reduce the reworks in actual projects.
文摘The following is proved: 1) The linear independence of assumed stress modes is the necessary and sufficient condition for the nonsingular flexibility matrix; 2) The equivalent assumed stress modes lead to the identical hybrid element. The Hilbert stress subspace of the assumed stress modes is established. So, it is easy to derive the equivalent orthogonal normal stress modes by Schmidt's method. Because of the resulting diagonal flexibility matrix, the identical hybrid element is free from the complex matrix inversion so that the hybrid efficiency, is improved greatly. The numerical examples show that the method is effective.
