The existing three-parameter single-step time integration methods, such as the Generalized-a method, improve numerical dissipation by modifying equilibrium equation at time points, which cause them to lose accuracy du...The existing three-parameter single-step time integration methods, such as the Generalized-a method, improve numerical dissipation by modifying equilibrium equation at time points, which cause them to lose accuracy due to the interpolation of load vectors. Moreover, these three-parameter methods do not present an available formulation applied to a general secondorder non linear differential equatio n. To solve these problems, this paper proposes an innovative three-parameter single-step method by introducing an additional variable into update equations. Although the present method is spectrally identical to the Generalized-cx method for undamped systems, it possesses higher accuracy since it strictly satisfies the equilibrium equation at time points, and can be readily used to solve nonlinear equations. By the analysis of accuracy, stability, numerical dissipation and dispersion, the optimal second-order implicit and explicit schemes are generated, which can maximize low-frequency accuracy when high-frequency dissipation is specified. To check the performance of the proposed method, several numerical experiments are conducted and the proposed method is compared with a few up-to-date methods.展开更多
Based on the weighted residual method,a single-step time integration algorithm with higher-order accuracy and unconditional stability has been proposed,which is superior to the second-order accurate algorithms in trac...Based on the weighted residual method,a single-step time integration algorithm with higher-order accuracy and unconditional stability has been proposed,which is superior to the second-order accurate algorithms in tracking long-term dynamics.For improving such a higher-order accurate algorithm,this paper proposes a two sub-step higher-order algorithm with unconditional stability and controllable dissipation.In the proposed algorithm,a time step interval[t_(k),t_(k)+h]where h stands for the size of a time step is divided into two sub-steps[t_(k),t_(k)+γh]and[t_(k)+γh,t_(k)+h].A non-dissipative fourth-order algorithm is used in the rst sub-step to ensure low-frequency accuracy and a dissipative third-order algorithm is employed in the second sub-step to lter out the contribution of high-frequency modes.Besides,two approaches are used to design the algorithm parameterγ.The rst approach determinesγby maximizing low-frequency accuracy and the other determinesγfor quickly damping out highfrequency modes.The present algorithm usesρ_(∞)to exactly control the degree of numerical dissipation,and it is third-order accurate when 0≤ρ_(∞)<1 and fourth-order accurate whenρ_(∞)=1.Furthermore,the proposed algorithm is self-starting and easy to implement.Some illustrative linear and nonlinear examples are solved to check the performances of the proposed two sub-step higher-order algorithm.展开更多
As thermal protection substrates for wearable electronics,functional soft composites made of polymer materials embedded with phase change materials and metal layers demonstrate unique capabilities for the thermal prot...As thermal protection substrates for wearable electronics,functional soft composites made of polymer materials embedded with phase change materials and metal layers demonstrate unique capabilities for the thermal protection of human skin.Here,we develop an analytical transient phase change heat transfer model to investigate the thermal performance of a wearable electronic device with a thermal protection substrate.The model is validated by experiments and the finite element analysis(FEA).The effects of the substrate structure size and heat source power input on the temperature management efficiency are investigated systematically and comprehensively.The results show that the objective of thermal management for wearable electronics is achieved by the following thermal protection mechanism.The metal thin film helps to dissipate heat along the in-plane direction by reconfiguring the direction of heat flow,while the phase change material assimilates excessive heat.These results will not only promote the fundamental understanding of the thermal properties of wearable electronics incorporating thermal protection substrates,but also facilitate the rational design of thermal protection substrates for wearable electronics.展开更多
Highly accurate closed-form eigensolutions for flutter of three-dimensional(3D)panel with arbitrary combinations of simply supported(S),glide(G),clamped(C)and free(F)boundary conditions(BCs),such as cantilever panels,...Highly accurate closed-form eigensolutions for flutter of three-dimensional(3D)panel with arbitrary combinations of simply supported(S),glide(G),clamped(C)and free(F)boundary conditions(BCs),such as cantilever panels,are achieved according to the linear thin plate theory and the first-order piston theory as well as the complex modal analysis,and all solutions are in a simple and explicit form.The iterative Separation-of-Variable(iSOV)method proposed by the pre-sent authors is employed to obtain the highly accurate eigensolutions.The flutter mechanism is studied with the benefit of eigenvalue properties from mathematical senses.The effects of boundary conditions,chord-thickness ratios,aerodynamic damping,aspect ratios and in-plane loads on flut-ter properties are examined.The results are compared with those of Kantorovich method and Galerkin method,and also coincide well with analytical solutions in literature,verifying the accu-racy of the present closed-form results.It is revealed that,(A)the flutter characteristics are domi-nated by the cross section properties of panels in the direction of stream flow;(B)two types of flutter,called coupled-mode flutter and zero-frequency flutter which includes zero-frequency single-mode flutter and buckling,are observed;(C)boundary conditions and in-plane loads can affect both flutter boundary and flutter type;(D)the flutter behavior of 3D panel is similar to that of the two-dimensional(2D)panel if the aspect ratio is up to a certain value;(E)four to six modes should be used in the Galerkin method for accurate eigensolutions,and the results converge to that of Kantorovich method which uses the same mode functions in the direction perpendicular to the stream flow.The present analysis method can be used as a reference for other stability issues characterized by complex eigenvalues,and the highly closed-form solutions are useful in parameter designs and can also be taken as benchmarks for the validation of numerical methods.展开更多
Based on the eigenvector expansion idea, the Multiscale Eigenelement Method(MEM)was proposed by the author and co-workers. MEM satisfies two equivalent conditions, one condition is the equivalence of strain energy, an...Based on the eigenvector expansion idea, the Multiscale Eigenelement Method(MEM)was proposed by the author and co-workers. MEM satisfies two equivalent conditions, one condition is the equivalence of strain energy, and the other is the deformation similarity. These two equivalent conditions character the structure-preserving property of a multiscale analysis method. The equivalence of strain energy is necessary for achieving accurate macro behaviors such as lower order frequencies, while the deformation similarity is essential for predicting accurate micro behaviors such as stresses. The MEM has become a powerful multiscale method for the analysis of composite structures because of its high accuracy and efficiency. In this paper, the research advances of MEM are reviewed and all types of eigenelement methods are compared, focusing on superiorities and deficiencies from practical viewpoint. It is concluded that the eigenelement methods with smooth shape functions are more suitable for the analysis of macro behaviors such as lower order frequencies, and the eigenelement methods with piecewise shape functions are suitable for the analysis of both macro and micro behaviors.展开更多
基金the National NaturalScience Foundation of China (11672019, 11372021. and 37686003).
文摘The existing three-parameter single-step time integration methods, such as the Generalized-a method, improve numerical dissipation by modifying equilibrium equation at time points, which cause them to lose accuracy due to the interpolation of load vectors. Moreover, these three-parameter methods do not present an available formulation applied to a general secondorder non linear differential equatio n. To solve these problems, this paper proposes an innovative three-parameter single-step method by introducing an additional variable into update equations. Although the present method is spectrally identical to the Generalized-cx method for undamped systems, it possesses higher accuracy since it strictly satisfies the equilibrium equation at time points, and can be readily used to solve nonlinear equations. By the analysis of accuracy, stability, numerical dissipation and dispersion, the optimal second-order implicit and explicit schemes are generated, which can maximize low-frequency accuracy when high-frequency dissipation is specified. To check the performance of the proposed method, several numerical experiments are conducted and the proposed method is compared with a few up-to-date methods.
基金supported by the National Natural Science Foundation of China(Grant Numbers 11872090,11672019,11472035).
文摘Based on the weighted residual method,a single-step time integration algorithm with higher-order accuracy and unconditional stability has been proposed,which is superior to the second-order accurate algorithms in tracking long-term dynamics.For improving such a higher-order accurate algorithm,this paper proposes a two sub-step higher-order algorithm with unconditional stability and controllable dissipation.In the proposed algorithm,a time step interval[t_(k),t_(k)+h]where h stands for the size of a time step is divided into two sub-steps[t_(k),t_(k)+γh]and[t_(k)+γh,t_(k)+h].A non-dissipative fourth-order algorithm is used in the rst sub-step to ensure low-frequency accuracy and a dissipative third-order algorithm is employed in the second sub-step to lter out the contribution of high-frequency modes.Besides,two approaches are used to design the algorithm parameterγ.The rst approach determinesγby maximizing low-frequency accuracy and the other determinesγfor quickly damping out highfrequency modes.The present algorithm usesρ_(∞)to exactly control the degree of numerical dissipation,and it is third-order accurate when 0≤ρ_(∞)<1 and fourth-order accurate whenρ_(∞)=1.Furthermore,the proposed algorithm is self-starting and easy to implement.Some illustrative linear and nonlinear examples are solved to check the performances of the proposed two sub-step higher-order algorithm.
基金Project supported by the National Natural Science Foundation of China(No.11772030)the Aeronautical Science Foundation of China(No.2018ZC51030)the Opening fund of State Key Laboratory of Structural Analysis for Industrial Equipment of Dalian University of Technology(No.GZ19117)。
文摘As thermal protection substrates for wearable electronics,functional soft composites made of polymer materials embedded with phase change materials and metal layers demonstrate unique capabilities for the thermal protection of human skin.Here,we develop an analytical transient phase change heat transfer model to investigate the thermal performance of a wearable electronic device with a thermal protection substrate.The model is validated by experiments and the finite element analysis(FEA).The effects of the substrate structure size and heat source power input on the temperature management efficiency are investigated systematically and comprehensively.The results show that the objective of thermal management for wearable electronics is achieved by the following thermal protection mechanism.The metal thin film helps to dissipate heat along the in-plane direction by reconfiguring the direction of heat flow,while the phase change material assimilates excessive heat.These results will not only promote the fundamental understanding of the thermal properties of wearable electronics incorporating thermal protection substrates,but also facilitate the rational design of thermal protection substrates for wearable electronics.
基金supported by the National Natural Science Foundation of China(Nos.11872090,11672019,11472035)。
文摘Highly accurate closed-form eigensolutions for flutter of three-dimensional(3D)panel with arbitrary combinations of simply supported(S),glide(G),clamped(C)and free(F)boundary conditions(BCs),such as cantilever panels,are achieved according to the linear thin plate theory and the first-order piston theory as well as the complex modal analysis,and all solutions are in a simple and explicit form.The iterative Separation-of-Variable(iSOV)method proposed by the pre-sent authors is employed to obtain the highly accurate eigensolutions.The flutter mechanism is studied with the benefit of eigenvalue properties from mathematical senses.The effects of boundary conditions,chord-thickness ratios,aerodynamic damping,aspect ratios and in-plane loads on flut-ter properties are examined.The results are compared with those of Kantorovich method and Galerkin method,and also coincide well with analytical solutions in literature,verifying the accu-racy of the present closed-form results.It is revealed that,(A)the flutter characteristics are domi-nated by the cross section properties of panels in the direction of stream flow;(B)two types of flutter,called coupled-mode flutter and zero-frequency flutter which includes zero-frequency single-mode flutter and buckling,are observed;(C)boundary conditions and in-plane loads can affect both flutter boundary and flutter type;(D)the flutter behavior of 3D panel is similar to that of the two-dimensional(2D)panel if the aspect ratio is up to a certain value;(E)four to six modes should be used in the Galerkin method for accurate eigensolutions,and the results converge to that of Kantorovich method which uses the same mode functions in the direction perpendicular to the stream flow.The present analysis method can be used as a reference for other stability issues characterized by complex eigenvalues,and the highly closed-form solutions are useful in parameter designs and can also be taken as benchmarks for the validation of numerical methods.
基金supported by the National Natural Science Foundation of China (Nos. 11672019, 11372021 and 37686003)the Academic Excellence Foundation of BUAA for PhD Students (No. 2017038)
文摘Based on the eigenvector expansion idea, the Multiscale Eigenelement Method(MEM)was proposed by the author and co-workers. MEM satisfies two equivalent conditions, one condition is the equivalence of strain energy, and the other is the deformation similarity. These two equivalent conditions character the structure-preserving property of a multiscale analysis method. The equivalence of strain energy is necessary for achieving accurate macro behaviors such as lower order frequencies, while the deformation similarity is essential for predicting accurate micro behaviors such as stresses. The MEM has become a powerful multiscale method for the analysis of composite structures because of its high accuracy and efficiency. In this paper, the research advances of MEM are reviewed and all types of eigenelement methods are compared, focusing on superiorities and deficiencies from practical viewpoint. It is concluded that the eigenelement methods with smooth shape functions are more suitable for the analysis of macro behaviors such as lower order frequencies, and the eigenelement methods with piecewise shape functions are suitable for the analysis of both macro and micro behaviors.