Based on the principle of piecewise linearization, the incremental forms of microstructure evolution models were integrated into the thermo-mechanical coupled finite element(FE) model to simulate nonlinear microstru...Based on the principle of piecewise linearization, the incremental forms of microstructure evolution models were integrated into the thermo-mechanical coupled finite element(FE) model to simulate nonlinear microstructure evolution during multi-pass hot deformation. This is an unsteady-state deformation where dynamic recrystallization(DRX), meta-dynamic recrystallization(MDRX), static recrystallization(SRX) and grain growth(GG) take place during hot deformation or deformation interval. The distributions of deformation and microstructure for cylindrical AZ31 sample during single-pass and double-pass hot compressions were quantitatively calculated and compared with the metallographic observation. It is shown that both the deformation and microstructure are non-uniformly distributed due to the presence of friction between the die and the flat end of sample. The average grain size and its standard deviation under the double-pass hot compression are slightly smaller than those under single-pass compression.The simulated average grain sizes agree well with the experiments, which validates that the developed FE model on the basis of incremental forms of microstructure evolution models is reasonable.展开更多
By the improvement of Riks’and Crisfield’s arc-length method,the adaptiveparameter incremental method is preasted for predicting the snapping response ofstructures. Its justification is fulfilled. Finally,the effect...By the improvement of Riks’and Crisfield’s arc-length method,the adaptiveparameter incremental method is preasted for predicting the snapping response ofstructures. Its justification is fulfilled. Finally,the effectiveness of this method isdemonstrated by solving the snapping response of spherical caps subjected to centrallydistributed pressures.展开更多
We have deduced incremental harmonic balance an iteration scheme in the (IHB) method using the harmonic balance plus the Newton-Raphson method. Since the convergence of the iteration is dependent upon the initial va...We have deduced incremental harmonic balance an iteration scheme in the (IHB) method using the harmonic balance plus the Newton-Raphson method. Since the convergence of the iteration is dependent upon the initial values in the iteration, the convergent region is greatly restricted for some cases. In this contribution, in order to enlarge the convergent region of the IHB method, we constructed the zeroth-order deformation equation using the homotopy analysis method, in which the IHB method is employed to solve the deformation equation with an embedding parameter as the active increment. Taking the Duffing and the van der Pol equations as examples, we obtained the highly accurate solutions. Importantly, the presented approach renders a convenient way to control and adjust the convergence.展开更多
Dielectric elastomer(DE) is suitable in soft transducers for broad applications,among which many are subjected to dynamic loadings, either mechanical or electrical or both. The tuning behaviors of these DE devices cal...Dielectric elastomer(DE) is suitable in soft transducers for broad applications,among which many are subjected to dynamic loadings, either mechanical or electrical or both. The tuning behaviors of these DE devices call for an efficient and reliable method to analyze the dynamic response of DE. This remains to be a challenge since the resultant vibration equation of DE, for example, the vibration of a DE balloon considered here is highly nonlinear with higher-order power terms and time-dependent coefficients. Previous efforts toward this goal use largely the numerical integration method with the simple harmonic balance method as a supplement. The numerical integration and the simple harmonic balance method are inefficient for large parametric analysis or with difficulty in improving the solution accuracy. To overcome the weakness of these two methods,we describe formulations of the incremental harmonic balance(IHB) method for periodic forced solutions of such a unique system. Combined with an arc-length continuation technique, the proposed strategy can capture the whole solution branches, both stable and unstable, automatically with any desired accuracy.展开更多
The incremental harmonic balance method was extended to analyze the flutter of systems with multiple structural strong nonlinearities. The strongly nonlinear cubic plunging and pitching stiffness terms were considered...The incremental harmonic balance method was extended to analyze the flutter of systems with multiple structural strong nonlinearities. The strongly nonlinear cubic plunging and pitching stiffness terms were considered in the flutter equations of two-dimensional airfoil. First, the equations were transferred into matrix form, then the vibration process was divided into the persistent incremental processes of vibration moments. And the expression of their solutions could be obtained by using a certain amplitude as control parameter in the harmonic balance process, and then the bifurcation, limit cycle flutter phenomena and the number of harmonic terms were analyzed. Finally, numerical results calculated by the Runge-Kutta method were given to verify the results obtained by the proposed procedure. It has been shown that the incremental harmonic method is effective and precise in the analysis of strongly nonlinear flutter with multiple structural nonlinearities.展开更多
The integration of set-valued ordered rough set models and incremental learning signify a progressive advancement of conventional rough set theory, with the objective of tackling the heterogeneity and ongoing transfor...The integration of set-valued ordered rough set models and incremental learning signify a progressive advancement of conventional rough set theory, with the objective of tackling the heterogeneity and ongoing transformations in information systems. In set-valued ordered decision systems, when changes occur in the attribute value domain, such as adding conditional values, it may result in changes in the preference relation between objects, indirectly leading to changes in approximations. In this paper, we effectively addressed the issue of updating approximations that arose from adding conditional values in set-valued ordered decision systems. Firstly, we classified the research objects into two categories: objects with changes in conditional values and objects without changes, and then conducted theoretical studies on updating approximations for these two categories, presenting approximation update theories for adding conditional values. Subsequently, we presented incremental algorithms corresponding to approximation update theories. We demonstrated the feasibility of the proposed incremental update method with numerical examples and showed that our incremental algorithm outperformed the static algorithm. Ultimately, by comparing experimental results on different datasets, it is evident that the incremental algorithm efficiently reduced processing time. In conclusion, this study offered a promising strategy to address the challenges of set-valued ordered decision systems in dynamic environments.展开更多
A novel incremental nonlinear detection algorithm is presented for Multiple-Input Multiple-Output (MIMO) system. In this algorithm, the received data at multiple receiver antennas are nonlinearly mapped and then sum...A novel incremental nonlinear detection algorithm is presented for Multiple-Input Multiple-Output (MIMO) system. In this algorithm, the received data at multiple receiver antennas are nonlinearly mapped and then summed with weights. The weight coefficients are incrementally computed to avoid direct computation of the inverse of a matrix, which greatly reduce the computational complexity. Simulation and comparison show that the proposed algorithm can obtain better performance of Bit Error Rate (BER) than linear Minimum Mean Square Error (MMSE).展开更多
Warp yarns and weft yarns of plain woven fabric are the principal axes of mate- rial of fabric. They are orthogonal in their original con?guration, but are obliquely crisscross in deformed con?guration in general. I...Warp yarns and weft yarns of plain woven fabric are the principal axes of mate- rial of fabric. They are orthogonal in their original con?guration, but are obliquely crisscross in deformed con?guration in general. In this paper the expressions of incremental components of strain tensor are derived, the non-linear model of woven fabric is linearized physically and its geometric non-linearity survives. The convenience of determining the total deformation is shown by the choice of the coordinate system of the principal axes of the material, with the convergence of the incremental methods illustrated by examples. This incremental model furnishes a basis for numerical simulations of fabric draping and wrinkling based on the micro-mechanical model of fabric.展开更多
This study has focused on developing numerical procedures for the static and dynamic nonlinear analysis of mooring lines. A geometrically nonlinear finite element method using isoparametric cable element with two node...This study has focused on developing numerical procedures for the static and dynamic nonlinear analysis of mooring lines. A geometrically nonlinear finite element method using isoparametric cable element with two nodes is briefly presented on the basis of the total Lagrangian formulation. The static and dynamic equilibrium equations of mooring lines are established. An incremental-iterative method is used to determine the initial static equilibrium state of cable systems under the action of self weights, buoyancy and current. Also the Newmark method is used for dynamic nonlinear analysis of ocean cables. Numerical examples are presented to validate the present numerical method, and examine the effect of various parameters.展开更多
A force-based quadrilateral plate element( 4NQP13) for the analysis of the plate bending problems using large increment method( LIM) was proposed. The LIM, a force-based finite element method( FEM),has been successful...A force-based quadrilateral plate element( 4NQP13) for the analysis of the plate bending problems using large increment method( LIM) was proposed. The LIM, a force-based finite element method( FEM),has been successfully developed for the analysis of truss,beam,frame,and 2D continua problems. In these analyses,LIMcan provide more precise stress results and less computational time consumption compared with displacement-based FEM. The plate element was based on the Mindlin-Reissner plate theory which took into account the transverse shear effects.Numerical examples were presented to study its performance including accuracy and convergence behavior,and the results were compared with the results have been obtained from the displacementbased quadrilateral plate elements and the analytical solutions. The4NQP13 element can analyze the moderately thick plates and the thin plates using LIMand is free from spurious zero energy modes and free from shear locking for thin plate analysis.展开更多
A quantitative analysis of limit cycles and homoclinic orbits, and the bifurcation curve for the Bogdanov-Takens system are discussed. The parameter incremental method for approximate analytical-expressions of these p...A quantitative analysis of limit cycles and homoclinic orbits, and the bifurcation curve for the Bogdanov-Takens system are discussed. The parameter incremental method for approximate analytical-expressions of these problems is given. These analytical-expressions of the limit cycle and homoclinic orbit are shown as the generalized harmonic functions by employing a time transformation. Curves of the parameters and the stability characteristic exponent of the limit cycle versus amplitude are drawn. Some of the limit cycles and homoclinic orbits phase portraits are plotted. The relationship curves of parameters μ and A with amplitude a and the bifurcation diagrams about the parameter are also given. The numerical accuracy of the calculation results is good.展开更多
Rotor system supported by nonlinear bearing such as squeeze film damper(SFD)is widely used in practice owing to its wide range of damping capacity and simplicity in structure.In this paper,an improved and effective In...Rotor system supported by nonlinear bearing such as squeeze film damper(SFD)is widely used in practice owing to its wide range of damping capacity and simplicity in structure.In this paper,an improved and effective Incremental transfer matrix method(ITMM)is first presented by combining ITMM and fast Fourier transform(FFT).Afterwards this method is applied to calculate the dynamic characteristics of a Jeffcott rotor system with SFD.The convergence dificulties incurred caused by strong nonlinearities of SFD has been dealt by adopting a control factor.It is found that for the more general boundary problems where the boundary conditions are not at input and output ends of a chain system,the supplementary equation is necessarily added.Additionally,the Floquet theory is used to analyze the stability and bifurcation type of the obtained periodic solution.The semi-analytical results,including the periodic solutions of the system,the bifurcation points and their types,are in good agreement with the numerical method.Furthermore,the involution mechanism of the quasi-periodic and chaotic motions near the first-order translational mode and the second order bending mode of this system is also clarified by this method with the aid of Floquet theory.展开更多
Aimed at the difficulty in revealing the vibration localization mechanism of mistuned bladed disks by using simple non-linear model,a mechanical model of the bladed disk with random mistuning of hysteretic dry frictio...Aimed at the difficulty in revealing the vibration localization mechanism of mistuned bladed disks by using simple non-linear model,a mechanical model of the bladed disk with random mistuning of hysteretic dry friction damping was established.Then,the incremental harmonic balance method was used to analyze the effects of the parameters of bladed disks,such as the mistuning strength of dry friction force,coupled strength,viscous damping ratio and friction strength,on the forced response of the bladed disks.The results show that the vibrational energy localization phenomenon turns up in the tuned bladed disks if the nonlinear friction damping exists,and the random mistuning of the dry friction force intensifies this kind of vibration localization.展开更多
The aim of this paper is to develop computational models for the ultimate compressive strength analysis of stiffened plate panels with nonuniform thickness.Modeling welding-induced initial deformations and residual st...The aim of this paper is to develop computational models for the ultimate compressive strength analysis of stiffened plate panels with nonuniform thickness.Modeling welding-induced initial deformations and residual stresses was presented with the measured data.Three methods,i.e.,ANSYS finite element method,ALPS/SPINE incremental Galerkin method,and ALPS/ULSAP analytical method,were employed together with existing test database obtained from a full-scale collapse testing of steel-stiffened plate structures.Sensitivity study was conducted with varying the difference in plate thickness to define a representative(equivalent)thickness for plate panels with nonuniform thickness.Guidelines are provided for structural modeling to compute the ultimate compressive strength of plate panels with variable thickness.展开更多
The vibration of a longitudinally moving rectangular plate submersed in an infinite liquid domain is studied analytically with the Rayleigh-Ritz method. The liquid is assumed to be incompressible, inviscid, and irrota...The vibration of a longitudinally moving rectangular plate submersed in an infinite liquid domain is studied analytically with the Rayleigh-Ritz method. The liquid is assumed to be incompressible, inviscid, and irrotational, and the velocity potential is used to describe the fluid velocity in the whole liquid field. The classical thin plate theory is used to derive mechanical energies of the traveling plate. As derivative of transverse displacement with respect to time in the compatibility condition equation exists, an exponential function is introduced to depict the dynamic deformation of the moving plate. It is shown that this exponential function works well with the Rayleigh- Ritz method. A convergence study shows a quick convergence speed for the immersed moving plate. Furthermore, the parametric study is carried out to demonstrate the effect of system parameters including the moving speed, the plate location, the liquid depth, the plate-liquid ratio, and the boundary condition. Results show that the above system parameters have significant influence on the vibration characteristics of the immersed moving plate. To extend the study, the method of added virtual mass incremental (AVMI) factor is used. The results show good agreement with those from the Rayleigh-Ritz method.展开更多
The subharmonic resonance and bifurcations of a clamped-clamped buckled beam under base harmonic excitations are investigated.The nonlinear partial integrodifferential equation of the motion of the buckled beam with b...The subharmonic resonance and bifurcations of a clamped-clamped buckled beam under base harmonic excitations are investigated.The nonlinear partial integrodifferential equation of the motion of the buckled beam with both quadratic and cubic nonlinearities is given by using Hamilton’s principle.A set of second-order nonlinear ordinary differential equations are obtained by spatial discretization with the Galerkin method.A high-dimensional model of the buckled beam is derived,concerning nonlinear coupling.The incremental harmonic balance(IHB)method is used to achieve the periodic solutions of the high-dimensional model of the buckled beam to observe the nonlinear frequency response curve and the nonlinear amplitude response curve,and the Floquet theory is used to analyze the stability of the periodic solutions.Attention is focused on the subharmonic resonance caused by the internal resonance as the excitation frequency near twice of the first natural frequency of the buckled beam with/without the antisymmetric modes being excited.Bifurcations including the saddle-node,Hopf,perioddoubling,and symmetry-breaking bifurcations are observed.Furthermore,quasi-periodic motion is observed by using the fourth-order Runge-Kutta method,which results from the Hopf bifurcation of the response of the buckled beam with the anti-symmetric modes being excited.展开更多
As a force-based finite element method (FEM), large increment method (LIM) has been developed in recent years. It has been shown that LIM provided prominent advantage of parallel computation with high efficiency and l...As a force-based finite element method (FEM), large increment method (LIM) has been developed in recent years. It has been shown that LIM provided prominent advantage of parallel computation with high efficiency and low time consumption for member structural system. To fully utilize its advantage in parallel computation, it is the time to extend LIM to 2D and 3D continua analysis. In this paper, a 2D finite element library with the capability of modeling arbitrary configurations is developed. Some illustrative numerical examples are solved by using the proposed library; the obtained results are compared with those obtained from both traditional displacement-based FEM and analytical solutions, which has clearly shown the advantages of LIM.展开更多
Many displacement-based quadrilateral plate elements based on Mindlin-Reissner plate theory have been proposed to analyze the thin and moderately thick plate problems. However, numerical inaccuracies of some elements ...Many displacement-based quadrilateral plate elements based on Mindlin-Reissner plate theory have been proposed to analyze the thin and moderately thick plate problems. However, numerical inaccuracies of some elements appear since the presence of shear locking and spurious zero energy modes for thin plate problems. To overcome these shortcomings, we employ the large increment method(LIM) for the analyses of the plate bending problems, and propose a force-based 8-node quadrilateral plate(8NQP) element which is based on MindlinReissner plate theory and has no extra spurious zero energy mode. Several benchmark plate bending problems are presented to illustrate the accuracy and convergence of the plate element by comparing with the analytical solutions and displacement-based plate elements. The results show that the 8-node plate element produces fast convergence and accurate stress distributions in both the moderately thick and thin plate bending problems. The plate element is insensitive to mesh distortion and it can avoid the shear locking for thin plate analysis.展开更多
X-ray pulsars offer stable, periodic X-ray pulse sequences that can be used in spacecraft positioning systems. A method using X-ray pulsars to determine the initial orbit of a satellite is presented in this paper. Thi...X-ray pulsars offer stable, periodic X-ray pulse sequences that can be used in spacecraft positioning systems. A method using X-ray pulsars to determine the initial orbit of a satellite is presented in this paper. This method suggests only one detector to be equipped on the satellite and assumes that the detector observes three pulsars in turn. To improve the performance, the use of incremental phase in one observation duration is proposed, and the incremental phase is combined with the time difference of arrival(TDOA). Then, a weighted least squares(WLS) algorithm is formulated to calculate the initial orbit. Numerical simulations are performed to assess the proposed orbit determination method.展开更多
In this paper,a combined optimization of a coupled electricity and gas system is presented.For the electricity network a unit commitment problem with optimization of energy and reserves under a power pool,considering ...In this paper,a combined optimization of a coupled electricity and gas system is presented.For the electricity network a unit commitment problem with optimization of energy and reserves under a power pool,considering all system operational and unit technical constraints is solved.The gas network subproblem is a medium-scale mixed-integer nonconvex and nonlinear programming problem.The coupling constraints between the two networks are nonlinear as well.The resulting mixed-integer nonlinear program is linearized with the extended incremental method and an outer approximation technique.The resulting model is evaluated using the Greek power and gas system comprising fourteen gas-fired units under four different approximation accuracy levels.The results indicate the efficiency of the proposed mixed-integer linear program model and the interplay between computational requirements and accuracy.展开更多
基金the funding support from the National Natural Science Foundation of China (Grant No. 51573156, 51675335)
文摘Based on the principle of piecewise linearization, the incremental forms of microstructure evolution models were integrated into the thermo-mechanical coupled finite element(FE) model to simulate nonlinear microstructure evolution during multi-pass hot deformation. This is an unsteady-state deformation where dynamic recrystallization(DRX), meta-dynamic recrystallization(MDRX), static recrystallization(SRX) and grain growth(GG) take place during hot deformation or deformation interval. The distributions of deformation and microstructure for cylindrical AZ31 sample during single-pass and double-pass hot compressions were quantitatively calculated and compared with the metallographic observation. It is shown that both the deformation and microstructure are non-uniformly distributed due to the presence of friction between the die and the flat end of sample. The average grain size and its standard deviation under the double-pass hot compression are slightly smaller than those under single-pass compression.The simulated average grain sizes agree well with the experiments, which validates that the developed FE model on the basis of incremental forms of microstructure evolution models is reasonable.
文摘By the improvement of Riks’and Crisfield’s arc-length method,the adaptiveparameter incremental method is preasted for predicting the snapping response ofstructures. Its justification is fulfilled. Finally,the effectiveness of this method isdemonstrated by solving the snapping response of spherical caps subjected to centrallydistributed pressures.
基金supported by the National Natural Science Foundation of China (10772202)Doctoral Program Foundation of Ministry of Education of China (20050558032)Guangdong Province Natural Science Foundation (07003680, 05003295)
文摘We have deduced incremental harmonic balance an iteration scheme in the (IHB) method using the harmonic balance plus the Newton-Raphson method. Since the convergence of the iteration is dependent upon the initial values in the iteration, the convergent region is greatly restricted for some cases. In this contribution, in order to enlarge the convergent region of the IHB method, we constructed the zeroth-order deformation equation using the homotopy analysis method, in which the IHB method is employed to solve the deformation equation with an embedding parameter as the active increment. Taking the Duffing and the van der Pol equations as examples, we obtained the highly accurate solutions. Importantly, the presented approach renders a convenient way to control and adjust the convergence.
基金the National Natural Science Foundation of China(Nos.11702215 and11972277)the Natural Science Basic Research Plan in Shaanxi Province of China(Nos.2017JQ5062 and 2018JQ1029)。
文摘Dielectric elastomer(DE) is suitable in soft transducers for broad applications,among which many are subjected to dynamic loadings, either mechanical or electrical or both. The tuning behaviors of these DE devices call for an efficient and reliable method to analyze the dynamic response of DE. This remains to be a challenge since the resultant vibration equation of DE, for example, the vibration of a DE balloon considered here is highly nonlinear with higher-order power terms and time-dependent coefficients. Previous efforts toward this goal use largely the numerical integration method with the simple harmonic balance method as a supplement. The numerical integration and the simple harmonic balance method are inefficient for large parametric analysis or with difficulty in improving the solution accuracy. To overcome the weakness of these two methods,we describe formulations of the incremental harmonic balance(IHB) method for periodic forced solutions of such a unique system. Combined with an arc-length continuation technique, the proposed strategy can capture the whole solution branches, both stable and unstable, automatically with any desired accuracy.
基金Project supported by the Ph. D. Programs Foundation of Ministry of Education of China (No.20050558032) the Natural Science Foundation of Guangdong Province of China (No.05003295) the Foundation of Sun Yat-sen University Advanced Research Center (No.06M8) the Young Teacher Scientific Research Foundation of Sun Sat-sen University (No.1131011)
文摘The incremental harmonic balance method was extended to analyze the flutter of systems with multiple structural strong nonlinearities. The strongly nonlinear cubic plunging and pitching stiffness terms were considered in the flutter equations of two-dimensional airfoil. First, the equations were transferred into matrix form, then the vibration process was divided into the persistent incremental processes of vibration moments. And the expression of their solutions could be obtained by using a certain amplitude as control parameter in the harmonic balance process, and then the bifurcation, limit cycle flutter phenomena and the number of harmonic terms were analyzed. Finally, numerical results calculated by the Runge-Kutta method were given to verify the results obtained by the proposed procedure. It has been shown that the incremental harmonic method is effective and precise in the analysis of strongly nonlinear flutter with multiple structural nonlinearities.
文摘The integration of set-valued ordered rough set models and incremental learning signify a progressive advancement of conventional rough set theory, with the objective of tackling the heterogeneity and ongoing transformations in information systems. In set-valued ordered decision systems, when changes occur in the attribute value domain, such as adding conditional values, it may result in changes in the preference relation between objects, indirectly leading to changes in approximations. In this paper, we effectively addressed the issue of updating approximations that arose from adding conditional values in set-valued ordered decision systems. Firstly, we classified the research objects into two categories: objects with changes in conditional values and objects without changes, and then conducted theoretical studies on updating approximations for these two categories, presenting approximation update theories for adding conditional values. Subsequently, we presented incremental algorithms corresponding to approximation update theories. We demonstrated the feasibility of the proposed incremental update method with numerical examples and showed that our incremental algorithm outperformed the static algorithm. Ultimately, by comparing experimental results on different datasets, it is evident that the incremental algorithm efficiently reduced processing time. In conclusion, this study offered a promising strategy to address the challenges of set-valued ordered decision systems in dynamic environments.
文摘A novel incremental nonlinear detection algorithm is presented for Multiple-Input Multiple-Output (MIMO) system. In this algorithm, the received data at multiple receiver antennas are nonlinearly mapped and then summed with weights. The weight coefficients are incrementally computed to avoid direct computation of the inverse of a matrix, which greatly reduce the computational complexity. Simulation and comparison show that the proposed algorithm can obtain better performance of Bit Error Rate (BER) than linear Minimum Mean Square Error (MMSE).
基金Project supported by the National Natural Science Foundation of China (No. 10272079).
文摘Warp yarns and weft yarns of plain woven fabric are the principal axes of mate- rial of fabric. They are orthogonal in their original con?guration, but are obliquely crisscross in deformed con?guration in general. In this paper the expressions of incremental components of strain tensor are derived, the non-linear model of woven fabric is linearized physically and its geometric non-linearity survives. The convenience of determining the total deformation is shown by the choice of the coordinate system of the principal axes of the material, with the convergence of the incremental methods illustrated by examples. This incremental model furnishes a basis for numerical simulations of fabric draping and wrinkling based on the micro-mechanical model of fabric.
基金supported by the National Natural Science Foundation of China (Grant No.11072052)the National High Technology Research and Development Program of China (863 Program,Grant No.2006AA09A109-3)
文摘This study has focused on developing numerical procedures for the static and dynamic nonlinear analysis of mooring lines. A geometrically nonlinear finite element method using isoparametric cable element with two nodes is briefly presented on the basis of the total Lagrangian formulation. The static and dynamic equilibrium equations of mooring lines are established. An incremental-iterative method is used to determine the initial static equilibrium state of cable systems under the action of self weights, buoyancy and current. Also the Newmark method is used for dynamic nonlinear analysis of ocean cables. Numerical examples are presented to validate the present numerical method, and examine the effect of various parameters.
基金National Natural Science Foundation of China(No.10872128)
文摘A force-based quadrilateral plate element( 4NQP13) for the analysis of the plate bending problems using large increment method( LIM) was proposed. The LIM, a force-based finite element method( FEM),has been successfully developed for the analysis of truss,beam,frame,and 2D continua problems. In these analyses,LIMcan provide more precise stress results and less computational time consumption compared with displacement-based FEM. The plate element was based on the Mindlin-Reissner plate theory which took into account the transverse shear effects.Numerical examples were presented to study its performance including accuracy and convergence behavior,and the results were compared with the results have been obtained from the displacementbased quadrilateral plate elements and the analytical solutions. The4NQP13 element can analyze the moderately thick plates and the thin plates using LIMand is free from spurious zero energy modes and free from shear locking for thin plate analysis.
基金the National Natural Science Foundation of China (No.10672193)
文摘A quantitative analysis of limit cycles and homoclinic orbits, and the bifurcation curve for the Bogdanov-Takens system are discussed. The parameter incremental method for approximate analytical-expressions of these problems is given. These analytical-expressions of the limit cycle and homoclinic orbit are shown as the generalized harmonic functions by employing a time transformation. Curves of the parameters and the stability characteristic exponent of the limit cycle versus amplitude are drawn. Some of the limit cycles and homoclinic orbits phase portraits are plotted. The relationship curves of parameters μ and A with amplitude a and the bifurcation diagrams about the parameter are also given. The numerical accuracy of the calculation results is good.
文摘Rotor system supported by nonlinear bearing such as squeeze film damper(SFD)is widely used in practice owing to its wide range of damping capacity and simplicity in structure.In this paper,an improved and effective Incremental transfer matrix method(ITMM)is first presented by combining ITMM and fast Fourier transform(FFT).Afterwards this method is applied to calculate the dynamic characteristics of a Jeffcott rotor system with SFD.The convergence dificulties incurred caused by strong nonlinearities of SFD has been dealt by adopting a control factor.It is found that for the more general boundary problems where the boundary conditions are not at input and output ends of a chain system,the supplementary equation is necessarily added.Additionally,the Floquet theory is used to analyze the stability and bifurcation type of the obtained periodic solution.The semi-analytical results,including the periodic solutions of the system,the bifurcation points and their types,are in good agreement with the numerical method.Furthermore,the involution mechanism of the quasi-periodic and chaotic motions near the first-order translational mode and the second order bending mode of this system is also clarified by this method with the aid of Floquet theory.
基金Project(2007CB707706) supported by the National Basic Research Program of China
文摘Aimed at the difficulty in revealing the vibration localization mechanism of mistuned bladed disks by using simple non-linear model,a mechanical model of the bladed disk with random mistuning of hysteretic dry friction damping was established.Then,the incremental harmonic balance method was used to analyze the effects of the parameters of bladed disks,such as the mistuning strength of dry friction force,coupled strength,viscous damping ratio and friction strength,on the forced response of the bladed disks.The results show that the vibrational energy localization phenomenon turns up in the tuned bladed disks if the nonlinear friction damping exists,and the random mistuning of the dry friction force intensifies this kind of vibration localization.
文摘The aim of this paper is to develop computational models for the ultimate compressive strength analysis of stiffened plate panels with nonuniform thickness.Modeling welding-induced initial deformations and residual stresses was presented with the measured data.Three methods,i.e.,ANSYS finite element method,ALPS/SPINE incremental Galerkin method,and ALPS/ULSAP analytical method,were employed together with existing test database obtained from a full-scale collapse testing of steel-stiffened plate structures.Sensitivity study was conducted with varying the difference in plate thickness to define a representative(equivalent)thickness for plate panels with nonuniform thickness.Guidelines are provided for structural modeling to compute the ultimate compressive strength of plate panels with variable thickness.
基金Project supported by the National Natural Science Foundation of China(Nos.11302046 and 11672071)the Fundamental Research Funds for the Central Universities(No.N150504003)
文摘The vibration of a longitudinally moving rectangular plate submersed in an infinite liquid domain is studied analytically with the Rayleigh-Ritz method. The liquid is assumed to be incompressible, inviscid, and irrotational, and the velocity potential is used to describe the fluid velocity in the whole liquid field. The classical thin plate theory is used to derive mechanical energies of the traveling plate. As derivative of transverse displacement with respect to time in the compatibility condition equation exists, an exponential function is introduced to depict the dynamic deformation of the moving plate. It is shown that this exponential function works well with the Rayleigh- Ritz method. A convergence study shows a quick convergence speed for the immersed moving plate. Furthermore, the parametric study is carried out to demonstrate the effect of system parameters including the moving speed, the plate location, the liquid depth, the plate-liquid ratio, and the boundary condition. Results show that the above system parameters have significant influence on the vibration characteristics of the immersed moving plate. To extend the study, the method of added virtual mass incremental (AVMI) factor is used. The results show good agreement with those from the Rayleigh-Ritz method.
基金Project supported by the National Natural Science Foundation of China(Nos.11972381 and 11572354)the Fundamental Research Funds for the Central Universities(No.18lgzd08)。
文摘The subharmonic resonance and bifurcations of a clamped-clamped buckled beam under base harmonic excitations are investigated.The nonlinear partial integrodifferential equation of the motion of the buckled beam with both quadratic and cubic nonlinearities is given by using Hamilton’s principle.A set of second-order nonlinear ordinary differential equations are obtained by spatial discretization with the Galerkin method.A high-dimensional model of the buckled beam is derived,concerning nonlinear coupling.The incremental harmonic balance(IHB)method is used to achieve the periodic solutions of the high-dimensional model of the buckled beam to observe the nonlinear frequency response curve and the nonlinear amplitude response curve,and the Floquet theory is used to analyze the stability of the periodic solutions.Attention is focused on the subharmonic resonance caused by the internal resonance as the excitation frequency near twice of the first natural frequency of the buckled beam with/without the antisymmetric modes being excited.Bifurcations including the saddle-node,Hopf,perioddoubling,and symmetry-breaking bifurcations are observed.Furthermore,quasi-periodic motion is observed by using the fourth-order Runge-Kutta method,which results from the Hopf bifurcation of the response of the buckled beam with the anti-symmetric modes being excited.
基金the National Natural Science Foundation of China (No. 10872128)
文摘As a force-based finite element method (FEM), large increment method (LIM) has been developed in recent years. It has been shown that LIM provided prominent advantage of parallel computation with high efficiency and low time consumption for member structural system. To fully utilize its advantage in parallel computation, it is the time to extend LIM to 2D and 3D continua analysis. In this paper, a 2D finite element library with the capability of modeling arbitrary configurations is developed. Some illustrative numerical examples are solved by using the proposed library; the obtained results are compared with those obtained from both traditional displacement-based FEM and analytical solutions, which has clearly shown the advantages of LIM.
基金the National Natural Science Foundation of China(No.10872128)
文摘Many displacement-based quadrilateral plate elements based on Mindlin-Reissner plate theory have been proposed to analyze the thin and moderately thick plate problems. However, numerical inaccuracies of some elements appear since the presence of shear locking and spurious zero energy modes for thin plate problems. To overcome these shortcomings, we employ the large increment method(LIM) for the analyses of the plate bending problems, and propose a force-based 8-node quadrilateral plate(8NQP) element which is based on MindlinReissner plate theory and has no extra spurious zero energy mode. Several benchmark plate bending problems are presented to illustrate the accuracy and convergence of the plate element by comparing with the analytical solutions and displacement-based plate elements. The results show that the 8-node plate element produces fast convergence and accurate stress distributions in both the moderately thick and thin plate bending problems. The plate element is insensitive to mesh distortion and it can avoid the shear locking for thin plate analysis.
基金supported by the National Natural Science Foundation of China(No.61401340)the Natural Science Basic Research Plan in Shaanxi Province of China(No.2016JM6035)+1 种基金the Fundamental Research Funds for the Central Universities,China(No.JB161303)and the Areospace T.T.&C.Innovation Program(No.201515A)
文摘X-ray pulsars offer stable, periodic X-ray pulse sequences that can be used in spacecraft positioning systems. A method using X-ray pulsars to determine the initial orbit of a satellite is presented in this paper. This method suggests only one detector to be equipped on the satellite and assumes that the detector observes three pulsars in turn. To improve the performance, the use of incremental phase in one observation duration is proposed, and the incremental phase is combined with the time difference of arrival(TDOA). Then, a weighted least squares(WLS) algorithm is formulated to calculate the initial orbit. Numerical simulations are performed to assess the proposed orbit determination method.
基金funding through the DFG SFB/Transregio 154, Subprojects A05 and Z01
文摘In this paper,a combined optimization of a coupled electricity and gas system is presented.For the electricity network a unit commitment problem with optimization of energy and reserves under a power pool,considering all system operational and unit technical constraints is solved.The gas network subproblem is a medium-scale mixed-integer nonconvex and nonlinear programming problem.The coupling constraints between the two networks are nonlinear as well.The resulting mixed-integer nonlinear program is linearized with the extended incremental method and an outer approximation technique.The resulting model is evaluated using the Greek power and gas system comprising fourteen gas-fired units under four different approximation accuracy levels.The results indicate the efficiency of the proposed mixed-integer linear program model and the interplay between computational requirements and accuracy.