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 calculatin...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.展开更多
Dry friction damping structures are widely-used in aero-engines to mitigate vibration.The nonlinear nature of friction and the two-dimensional in-plane motion on the contact interface bring challenges to accurately an...Dry friction damping structures are widely-used in aero-engines to mitigate vibration.The nonlinear nature of friction and the two-dimensional in-plane motion on the contact interface bring challenges to accurately and efficiently predict the forced response of frictionally damped structures.The state-of-the-art Multi-Harmonic Balance Method(MHBM)on quasi-3D contact model in engineering cannot precisely capture the kinematics on the friction interface although the efficiency is high.The full-3D contact model can describe the constitutive relationship of the interface in a more accurate manner;however,the efficiency and convergence are not guaranteed for large-scale models.In this paper,a semi-analytical MHBM on full-3D contact model is proposed.The original Trajectory Tracking Method(TTM)for evaluating the contact force is reformulated to make the calculation more concise and the derivation of the Analytical Jacobian Matrix(AJM)feasible.Based on the chain rule of derivation,the AJM which is the core to upgrade the performance is deduced.Through a shrouded blade finite element model,the accuracy and efficiency of the proposed method are compared with both the MHBM on full-3D contact model with numerical Jacobian matrix and the MHBM on quasi-3D contact model with AJM.The results show that the AJM improves significantly the efficiency of the MHBM on full-3D contact model.The time cost of the proposed method is in the same order of magnitude as that of the MHBM on quasi-3D contact model.We also confirm that the full-3D contact model is necessary for the dynamic analyses of shrouded blades.If one uses the quasi-3D model,the estimation relative error of damping can even reach 31.8%in some cases.In addition,the AJM also brings benefits for stability analysis.It is highly recommended that engineers use the MHBM on full-3D contact model for the dynamic analysis and design of shrouded blades.展开更多
The primary impediments impeding the implementation of high-order methods in simulating viscous flow over complex configurations are robustness and convergence.These challenges impose significant constraints on comput...The primary impediments impeding the implementation of high-order methods in simulating viscous flow over complex configurations are robustness and convergence.These challenges impose significant constraints on computational efficiency,particularly in the domain of engineering applications.To address these concerns,this paper proposes a robust implicit high-order discontinuous Galerkin(DG)method for solving compressible Navier-Stokes(NS)equations on arbitrary grids.The method achieves a favorable equilibrium between computational stability and efficiency.To solve the linear system,an exact Jacobian matrix solving strategy is employed for preconditioning and matrix-vector generation in the generalized minimal residual(GMRES)method.This approach mitigates numerical errors in Jacobian solution during implicit calculations and facilitates the implementation of an adaptive Courant-Friedrichs-Lewy(CFL)number increasing strategy,with the aim of improving convergence and robustness.To further enhance the applicability of the proposed method for intricate grid distortions,all simulations are performed in the reference domain.This practice significantly improves the reversibility of the mass matrix in implicit calculations.A comprehensive analysis of various parameters influencing computational stability and efficiency is conducted,including CFL number,Krylov subspace size,and GMRES convergence criteria.The computed results from a series of numerical test cases demonstrate the promising results achieved by combining the DG method,GMRES solver,exact Jacobian matrix,adaptive CFL number,and reference domain calculations in terms of robustness,convergence,and accuracy.These analysis results can serve as a reference for implicit computation in high-order calculations.展开更多
Due to the diversity of work requirements and environment,the number of degrees of freedom(DOFs)and the complexity of structure of industrial robots are constantly increasing.It is difficult to establish the accurate ...Due to the diversity of work requirements and environment,the number of degrees of freedom(DOFs)and the complexity of structure of industrial robots are constantly increasing.It is difficult to establish the accurate dynamical model of industrial robots,which greatly hinders the realization of a stable,fast and accurate trajectory tracking control.Therefore,the general expression of the constraint relation in the explicit dynamic equation of the multi-DOF industrial robot is derived on the basis of solving the Jacobian matrix and Hessian matrix by using the kinematic influence coefficients method.Moreover,an explicit dynamic equation with general constraint relation expression is established based on the Udwadia-Kalaba theory.The problem of increasing the time of establishing constraint relationship when the multi-DOF industrial robots complete complex task constraints is solved.With the SCARA robot as the research object,the simulation results show that the proposed method can provide a new idea for industrial robot system modeling with complex constraints.展开更多
文摘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.
基金financially supported by the National Natural Science Foundation of China(Nos.52175071,91860205)the Major Projects of Aero-engines and Gas turbines(No.J2019-IV-023-0091)。
文摘Dry friction damping structures are widely-used in aero-engines to mitigate vibration.The nonlinear nature of friction and the two-dimensional in-plane motion on the contact interface bring challenges to accurately and efficiently predict the forced response of frictionally damped structures.The state-of-the-art Multi-Harmonic Balance Method(MHBM)on quasi-3D contact model in engineering cannot precisely capture the kinematics on the friction interface although the efficiency is high.The full-3D contact model can describe the constitutive relationship of the interface in a more accurate manner;however,the efficiency and convergence are not guaranteed for large-scale models.In this paper,a semi-analytical MHBM on full-3D contact model is proposed.The original Trajectory Tracking Method(TTM)for evaluating the contact force is reformulated to make the calculation more concise and the derivation of the Analytical Jacobian Matrix(AJM)feasible.Based on the chain rule of derivation,the AJM which is the core to upgrade the performance is deduced.Through a shrouded blade finite element model,the accuracy and efficiency of the proposed method are compared with both the MHBM on full-3D contact model with numerical Jacobian matrix and the MHBM on quasi-3D contact model with AJM.The results show that the AJM improves significantly the efficiency of the MHBM on full-3D contact model.The time cost of the proposed method is in the same order of magnitude as that of the MHBM on quasi-3D contact model.We also confirm that the full-3D contact model is necessary for the dynamic analyses of shrouded blades.If one uses the quasi-3D model,the estimation relative error of damping can even reach 31.8%in some cases.In addition,the AJM also brings benefits for stability analysis.It is highly recommended that engineers use the MHBM on full-3D contact model for the dynamic analysis and design of shrouded blades.
基金supported by the National Natural Science Foundation of China(Grant No.12102247)the Technology Development Program(Grant No.JCKY2022110C119).
文摘The primary impediments impeding the implementation of high-order methods in simulating viscous flow over complex configurations are robustness and convergence.These challenges impose significant constraints on computational efficiency,particularly in the domain of engineering applications.To address these concerns,this paper proposes a robust implicit high-order discontinuous Galerkin(DG)method for solving compressible Navier-Stokes(NS)equations on arbitrary grids.The method achieves a favorable equilibrium between computational stability and efficiency.To solve the linear system,an exact Jacobian matrix solving strategy is employed for preconditioning and matrix-vector generation in the generalized minimal residual(GMRES)method.This approach mitigates numerical errors in Jacobian solution during implicit calculations and facilitates the implementation of an adaptive Courant-Friedrichs-Lewy(CFL)number increasing strategy,with the aim of improving convergence and robustness.To further enhance the applicability of the proposed method for intricate grid distortions,all simulations are performed in the reference domain.This practice significantly improves the reversibility of the mass matrix in implicit calculations.A comprehensive analysis of various parameters influencing computational stability and efficiency is conducted,including CFL number,Krylov subspace size,and GMRES convergence criteria.The computed results from a series of numerical test cases demonstrate the promising results achieved by combining the DG method,GMRES solver,exact Jacobian matrix,adaptive CFL number,and reference domain calculations in terms of robustness,convergence,and accuracy.These analysis results can serve as a reference for implicit computation in high-order calculations.
基金the Beijing Municipal Scienceand Technology Project (No.KM202111417006)the Academic Research Projects of Beijing Union University (Nos.ZK10202305 and ZK80202004)the Beijing Municipal Science and Technology Project (No.KM202111417005)。
文摘Due to the diversity of work requirements and environment,the number of degrees of freedom(DOFs)and the complexity of structure of industrial robots are constantly increasing.It is difficult to establish the accurate dynamical model of industrial robots,which greatly hinders the realization of a stable,fast and accurate trajectory tracking control.Therefore,the general expression of the constraint relation in the explicit dynamic equation of the multi-DOF industrial robot is derived on the basis of solving the Jacobian matrix and Hessian matrix by using the kinematic influence coefficients method.Moreover,an explicit dynamic equation with general constraint relation expression is established based on the Udwadia-Kalaba theory.The problem of increasing the time of establishing constraint relationship when the multi-DOF industrial robots complete complex task constraints is solved.With the SCARA robot as the research object,the simulation results show that the proposed method can provide a new idea for industrial robot system modeling with complex constraints.
