In this paper,the parametric equations with multipliers of nonholonomic nonconservative sys- tems in the event space are established,their properties are studied,and their explicit formulation is obtained. And then th...In this paper,the parametric equations with multipliers of nonholonomic nonconservative sys- tems in the event space are established,their properties are studied,and their explicit formulation is obtained. And then the field method for integrating these equations is given.Finally,an example illustrating the appli- cation of the integration method is given.展开更多
An analytical technique, namely the homotopy analysis method （HAM）, is used to solve problems of nonlinear oscillations with parametric excitation. Unlike perturbation methods, HAM is not dependent on any small phys...An analytical technique, namely the homotopy analysis method （HAM）, is used to solve problems of nonlinear oscillations with parametric excitation. Unlike perturbation methods, HAM is not dependent on any small physical parameters at all, and thus valid for both weakly and strongly nonlinear problems. In addition, HAM is different from all other analytic techniques in providing a simple way to adjust and control convergence region of the series solution by means of an auxiliary parameter h. In the present paper, a periodic analytic approximations for nonlinear oscillations with parametric excitation are obtained by using HAM, and the results are validated by numerical simulations.展开更多
The objective of the paper is to develop a new algorithm for numerical solution of dynamic elastic-plastic strain hardening/softening problems. The gradient dependent model is adopted in the numerical model to overcom...The objective of the paper is to develop a new algorithm for numerical solution of dynamic elastic-plastic strain hardening/softening problems. The gradient dependent model is adopted in the numerical model to overcome the result mesh-sensitivity problem in the dynamic strain softening or strain localization analysis. The equations for the dynamic elastic-plastic problems are derived in terms of the parametric variational principle, which is valid for associated, non-associated and strain softening plastic constitutive models in the finite element analysis. The precise integration method, which has been widely used for discretization in time domain of the linear problems, is introduced for the solution of dynamic nonlinear equations. The new algorithm proposed is based on the combination of the parametric quadratic programming method and the precise integration method and has all the advantages in both of the algorithms. Results of numerical examples demonstrate not only the validity, but also the advantages of the algorithm proposed for the numerical solution of nonlinear dynamic problems.展开更多
Applying the parametric derivation method, Peierls energy and Peierls stress are calculated with a non-sinusoidal force law in the lattice theory, while the results obtained by the power-series expansion according to ...Applying the parametric derivation method, Peierls energy and Peierls stress are calculated with a non-sinusoidal force law in the lattice theory, while the results obtained by the power-series expansion according to sinusoidal law can be deduced as a limiting case of non- sinusoidal law. The simplified expressions of Peierls energy and Peierls stress are obtained for the limit of wide and narrow. Peierls energy and Peierls stress decrease monotonically with the factor of modification of force law. Present results can be used expediently for prediction of the correct order of magnitude of Peierls stress for materials.展开更多
Traditional procedures to treat vibrations of gyroscopic continua involve direct application of perturbation methods to a system with both a strong gyroscopic term and other weakly coupled terms.In this study,a gyrosc...Traditional procedures to treat vibrations of gyroscopic continua involve direct application of perturbation methods to a system with both a strong gyroscopic term and other weakly coupled terms.In this study,a gyroscopic modes decoupling method is used to obtain an equivalent system with decoupled gyroscopic modes having only weak couplings.Taking the axially moving string as an example,the instability boundaries in the vicinity of parametric resonances are detected using both the traditional coupled gyroscopic system and our system with decoupled gyroscopic modes,and the results are compared to show the advantages and disadvantages of each method.展开更多
The paper studies the parametric stochastic roll motion in the random waves.The differential equation of the ship parametric roll under random wave is established with considering the nonlinear damping and ship speed....The paper studies the parametric stochastic roll motion in the random waves.The differential equation of the ship parametric roll under random wave is established with considering the nonlinear damping and ship speed.Random sea surface is treated as a narrow-band stochastic process,and the stochastic parametric excitation is studied based on the effective wave theory.The nonlinear restored arm function obtained from the numerical simulation is expressed as the approximate analytic function.By using the stochastic averaging method,the differential equation of motion is transformed into Ito’s stochastic differential equation.The steady-state probability density function of roll motion is obtained,and the results are validated with the numerical simulation and model test.展开更多
The optimal condition and its geometrical characters of the least square adjustment were proposed. Then the relation between the transformed surface and least squares was discussed. Based on the above, a non iterative...The optimal condition and its geometrical characters of the least square adjustment were proposed. Then the relation between the transformed surface and least squares was discussed. Based on the above, a non iterative method, called the fitting method of pseudo polynomial, was derived in detail. The final least squares solution can be determined with sufficient accuracy in a single step and is not attained by moving the initial point in the view of iteration. The accuracy of the solution relys wholly on the frequency of Taylor’s series. The example verifies the correctness and validness of the method. [展开更多
This article presents the Parametric Iteration Method (PIM) for finding optimal control and its corresponding trajectory of linear systems. Without any discretization or transformation, PIM provides a sequence of func...This article presents the Parametric Iteration Method (PIM) for finding optimal control and its corresponding trajectory of linear systems. Without any discretization or transformation, PIM provides a sequence of functions which converges to the exact solution of problem. Our emphasis will be on an auxiliary parameter which directly affects on the rate of convergence. Comparison of PIM and the Variational Iteration Method (VIM) is given to show the preference of PIM over VIM. Numerical results are given for several test examples to demonstrate the applicability and efficiency of the method.展开更多
Since the research of flare slamming prediction is seldom when parametric rolling happens, we present an efficient approximation method for flare slamming analysis of large container ships in parametric rolling condit...Since the research of flare slamming prediction is seldom when parametric rolling happens, we present an efficient approximation method for flare slamming analysis of large container ships in parametric rolling conditions. We adopt a 6-DOF weakly nonlinear time domain model to predict the ship motions of parametric rolling conditions. Unlike previous flare slamming analysis, our proposed method takes roll motion into account to calculate the impact angle and relative vertical velocity between ship sections on the bow flare and wave surface. We use the Wagner model to analyze the slamming impact forces and the slamming occurrence probability. Through numerical simulations, we investigate the maximum flare slamming pressures of a container ship for different speeds and wave conditions. To further clarify the mechanism of flare slamming phenomena in parametric rolling conditions, we also conduct real-time simulations to determine the relationship between slamming pressure and 3-DOF motions, namely roll, pitch, and heave.展开更多
An axisymmetrical unit cell model was used to represent a bimodal Al alloy that was composed of both nano-grained （NG） and coarse-grained （CG） aluminum. Effects of microstructural and materials parameters on tensi...An axisymmetrical unit cell model was used to represent a bimodal Al alloy that was composed of both nano-grained （NG） and coarse-grained （CG） aluminum. Effects of microstructural and materials parameters on tensile properties of bimodal AI alloy were investigated by finite element method （FEM）. The parameters analyzed included aspect ratios of CG Al and the unit cell, volume fraction of CG Al （VFCG）, and yield strength and strain hardening exponent of CG Al. Aspect ratios of CG Al and the unit cell have no significant influence on tensile stress-strain response of the bimodal Al alloy. This phenomenon derives from the similarity in elastic modulus and coefficient of thermal expansion between CG AI and NG Al. Conversely, tensile properties of bimodal Al alloy are extremely sensitive to VFCG, yield strength and strain hardening exponent of CG Al. Specifically, as VFCG increases, both yield strength and ultimate tensile strength （UTS） of the bimodal Al alloy decreases, while uniform strain of bimodal AI alloy increases. In addition, an increase in yield strength of CG Al results in an increase in both yield stress and UTS of bimodal AI alloy and a decrease in uniform strain of bimodal Al alloy. The lower capability in lowering the increase of stress concentration in NG Al due to a higher yield strength of CG Al causes the lower uniform strain of the bimodal AI alloy. When strain hardening exponent of CG Al increases, 0.2% yield stress, UTS, and uniform strain of the bimodal Al alloy increases. This can be attributed to the increased work-hardening ability of CG Al with a higher strain hardening exponent.展开更多
The Voronoi cell finite element method （VCFEM） is adopted to overcome the limitations of the classic displacement based finite element method in the numerical simulation of heterogeneous materials. The parametric va...The Voronoi cell finite element method （VCFEM） is adopted to overcome the limitations of the classic displacement based finite element method in the numerical simulation of heterogeneous materials. The parametric variational principle and quadratic programming method are developed for elastic-plastic Voronoi finite element analysis of two-dimensional problems. Finite element formulations are derived and a standard quadratic programming model is deduced from the elastic-plastic equations. Influence of microscopic heterogeneities on the overall mechanical response of heterogeneous materials is studied in detail. The overall properties of heterogeneous materials depend mostly on the size, shape and distribution of the material phases of the microstructure. Numerical examples are presented to demonstrate the validity and effectiveness of the method developed.展开更多
Optimization analysis and computational fluid dynamics （CFDs） have been applied simultaneously, in which a parametric model plays an important role in finding the optimal solution. However, it is difficult to create...Optimization analysis and computational fluid dynamics （CFDs） have been applied simultaneously, in which a parametric model plays an important role in finding the optimal solution. However, it is difficult to create a parametric model for a complex shape with irregular curves, such as a submarine hull form. In this study, the cubic Bezier curve and curve-plane intersection method are used to generate a solid model of a parametric submarine hull form taking three input parameters into account： nose radius, tail radius, and length-height hull ratio （L/H）. Application program interface （API） scripting is also used to write code in the ANSYS DesignModeler. The results show that the submarine shape can be generated with some variation of the input parameters. An example is given that shows how the proposed method can be applied successfully to a hull resistance optimization case. The parametric design of the middle submarine type was chosen to be modified. First, the original submarine model was analyzed, in advance, using CFD. Then, using the response surface graph, some candidate optimal designs with a minimum hull resistance coefficient were obtained. Further, the optimization method in goal-driven optimization （GDO） was implemented to find the submarine hull form with the minimum hull resistance coefficient （Ct）. The minimum C, was obtained. The calculated difference in （7, values between the initial submarine and the optimum submarine is around 0.26%, with the C, of the initial submarine and the optimum submarine being 0.001 508 26 and 0.001 504 29, respectively. The results show that the optimum submarine hull form shows a higher nose radius （rn） and higher L/H than those of the initial submarine shape, while the radius of the tail （r1） is smaller than that of the initial shape.展开更多
This paper mainly proposes a new C-XSC （C- for eXtended Scientific Computing） software for the symmetric single step method and relaxation method for computing an enclosure for the solution set and compares the meth...This paper mainly proposes a new C-XSC （C- for eXtended Scientific Computing） software for the symmetric single step method and relaxation method for computing an enclosure for the solution set and compares the methods with others＇ and then makes some modifications and finally, examples illustrating the applicability of the proposed methods are given.展开更多
This paper aims to evaluate the stochastic response of steel columns subjected to blast loads using the modified single degree of freedom(MSDOF)method,which assessed towards the conventional single degree of freedom(S...This paper aims to evaluate the stochastic response of steel columns subjected to blast loads using the modified single degree of freedom(MSDOF)method,which assessed towards the conventional single degree of freedom(SDOF)and the experimentally validated Finite Element(FE)methods(LSDYNA).For this purpose,special atten-tion is given to calculating the response of H-shaped steel columns under blast.The damage amount is determined based on the support rotation criterion,which is expressed as a function of their maximum lateral mid-span dis-placement.To account for uncertainties in input parameters and obtain the failure probability,the Monte Carlo simulation(MCS)method is employed,complemented by the Latin Hypercube Sampling(LHS)method to reduce the number of simulations.A parametric analysis is hence performed to examine the effect of several input pa-rameters(including both deterministic and probabilistic parameters)on the probability of column damage as a function of support rotation.First,the MSDOF method confirms its higher accuracy in estimating the probability of column damage due to blast,compared to the conventional SDOF.The collected results also show that un-certainties of several input parameters have significant effects on the column behavior.In particular,geometric parameters(including cross-sectional characteristics,boundary conditions and column length)have major effect on the corresponding column response,in the same way of input blast load parameters and material properties.展开更多
基金The Project is supported by the National Natural Science Foundation of China
文摘In this paper,the parametric equations with multipliers of nonholonomic nonconservative sys- tems in the event space are established,their properties are studied,and their explicit formulation is obtained. And then the field method for integrating these equations is given.Finally,an example illustrating the appli- cation of the integration method is given.
文摘An analytical technique, namely the homotopy analysis method （HAM）, is used to solve problems of nonlinear oscillations with parametric excitation. Unlike perturbation methods, HAM is not dependent on any small physical parameters at all, and thus valid for both weakly and strongly nonlinear problems. In addition, HAM is different from all other analytic techniques in providing a simple way to adjust and control convergence region of the series solution by means of an auxiliary parameter h. In the present paper, a periodic analytic approximations for nonlinear oscillations with parametric excitation are obtained by using HAM, and the results are validated by numerical simulations.
文摘The objective of the paper is to develop a new algorithm for numerical solution of dynamic elastic-plastic strain hardening/softening problems. The gradient dependent model is adopted in the numerical model to overcome the result mesh-sensitivity problem in the dynamic strain softening or strain localization analysis. The equations for the dynamic elastic-plastic problems are derived in terms of the parametric variational principle, which is valid for associated, non-associated and strain softening plastic constitutive models in the finite element analysis. The precise integration method, which has been widely used for discretization in time domain of the linear problems, is introduced for the solution of dynamic nonlinear equations. The new algorithm proposed is based on the combination of the parametric quadratic programming method and the precise integration method and has all the advantages in both of the algorithms. Results of numerical examples demonstrate not only the validity, but also the advantages of the algorithm proposed for the numerical solution of nonlinear dynamic problems.
基金Project supported by the National Natural Science Foundation of China (No.10774196)the Science Foundation Project of CQ CSTC (No.2006BB4156)Chongqing University Postgraduates'Science and Innovation Fund (No.2007A1A0030240).
文摘Applying the parametric derivation method, Peierls energy and Peierls stress are calculated with a non-sinusoidal force law in the lattice theory, while the results obtained by the power-series expansion according to sinusoidal law can be deduced as a limiting case of non- sinusoidal law. The simplified expressions of Peierls energy and Peierls stress are obtained for the limit of wide and narrow. Peierls energy and Peierls stress decrease monotonically with the factor of modification of force law. Present results can be used expediently for prediction of the correct order of magnitude of Peierls stress for materials.
基金the National Natural Science Foundation of China(Grants 11772009,11672007)the Beijing Municipal Natural Science Foundation(Grant 3172003).
文摘Traditional procedures to treat vibrations of gyroscopic continua involve direct application of perturbation methods to a system with both a strong gyroscopic term and other weakly coupled terms.In this study,a gyroscopic modes decoupling method is used to obtain an equivalent system with decoupled gyroscopic modes having only weak couplings.Taking the axially moving string as an example,the instability boundaries in the vicinity of parametric resonances are detected using both the traditional coupled gyroscopic system and our system with decoupled gyroscopic modes,and the results are compared to show the advantages and disadvantages of each method.
基金the State Administration of Science,Technology and Industry for National Defense of China(Grant No.B2420132001).
文摘The paper studies the parametric stochastic roll motion in the random waves.The differential equation of the ship parametric roll under random wave is established with considering the nonlinear damping and ship speed.Random sea surface is treated as a narrow-band stochastic process,and the stochastic parametric excitation is studied based on the effective wave theory.The nonlinear restored arm function obtained from the numerical simulation is expressed as the approximate analytic function.By using the stochastic averaging method,the differential equation of motion is transformed into Ito’s stochastic differential equation.The steady-state probability density function of roll motion is obtained,and the results are validated with the numerical simulation and model test.
文摘The optimal condition and its geometrical characters of the least square adjustment were proposed. Then the relation between the transformed surface and least squares was discussed. Based on the above, a non iterative method, called the fitting method of pseudo polynomial, was derived in detail. The final least squares solution can be determined with sufficient accuracy in a single step and is not attained by moving the initial point in the view of iteration. The accuracy of the solution relys wholly on the frequency of Taylor’s series. The example verifies the correctness and validness of the method. [
文摘This article presents the Parametric Iteration Method (PIM) for finding optimal control and its corresponding trajectory of linear systems. Without any discretization or transformation, PIM provides a sequence of functions which converges to the exact solution of problem. Our emphasis will be on an auxiliary parameter which directly affects on the rate of convergence. Comparison of PIM and the Variational Iteration Method (VIM) is given to show the preference of PIM over VIM. Numerical results are given for several test examples to demonstrate the applicability and efficiency of the method.
基金supported by the ChinaMinistry of Education Key Research Project "KSHIP-II Project"(Knowledge-based Ship Design Hyper-Integrated Platform)No.GKZY010004the National Key Basic Research Program of China No.2014CB046804
文摘Since the research of flare slamming prediction is seldom when parametric rolling happens, we present an efficient approximation method for flare slamming analysis of large container ships in parametric rolling conditions. We adopt a 6-DOF weakly nonlinear time domain model to predict the ship motions of parametric rolling conditions. Unlike previous flare slamming analysis, our proposed method takes roll motion into account to calculate the impact angle and relative vertical velocity between ship sections on the bow flare and wave surface. We use the Wagner model to analyze the slamming impact forces and the slamming occurrence probability. Through numerical simulations, we investigate the maximum flare slamming pressures of a container ship for different speeds and wave conditions. To further clarify the mechanism of flare slamming phenomena in parametric rolling conditions, we also conduct real-time simulations to determine the relationship between slamming pressure and 3-DOF motions, namely roll, pitch, and heave.
基金supported by the Office of Naval Re-search, contract N00014-03-C-0163, monitored by Rod Pe-terson.
文摘An axisymmetrical unit cell model was used to represent a bimodal Al alloy that was composed of both nano-grained （NG） and coarse-grained （CG） aluminum. Effects of microstructural and materials parameters on tensile properties of bimodal AI alloy were investigated by finite element method （FEM）. The parameters analyzed included aspect ratios of CG Al and the unit cell, volume fraction of CG Al （VFCG）, and yield strength and strain hardening exponent of CG Al. Aspect ratios of CG Al and the unit cell have no significant influence on tensile stress-strain response of the bimodal Al alloy. This phenomenon derives from the similarity in elastic modulus and coefficient of thermal expansion between CG AI and NG Al. Conversely, tensile properties of bimodal Al alloy are extremely sensitive to VFCG, yield strength and strain hardening exponent of CG Al. Specifically, as VFCG increases, both yield strength and ultimate tensile strength （UTS） of the bimodal Al alloy decreases, while uniform strain of bimodal AI alloy increases. In addition, an increase in yield strength of CG Al results in an increase in both yield stress and UTS of bimodal AI alloy and a decrease in uniform strain of bimodal Al alloy. The lower capability in lowering the increase of stress concentration in NG Al due to a higher yield strength of CG Al causes the lower uniform strain of the bimodal AI alloy. When strain hardening exponent of CG Al increases, 0.2% yield stress, UTS, and uniform strain of the bimodal Al alloy increases. This can be attributed to the increased work-hardening ability of CG Al with a higher strain hardening exponent.
基金Project supported by the National Natural Science Foundation of China(Nos.10225212, 10421002 and 10332010)the NCET Program provided by the Ministry of Education and the National Key Basic Research Special Foundation of China (No.2005CB321704)
文摘The Voronoi cell finite element method （VCFEM） is adopted to overcome the limitations of the classic displacement based finite element method in the numerical simulation of heterogeneous materials. The parametric variational principle and quadratic programming method are developed for elastic-plastic Voronoi finite element analysis of two-dimensional problems. Finite element formulations are derived and a standard quadratic programming model is deduced from the elastic-plastic equations. Influence of microscopic heterogeneities on the overall mechanical response of heterogeneous materials is studied in detail. The overall properties of heterogeneous materials depend mostly on the size, shape and distribution of the material phases of the microstructure. Numerical examples are presented to demonstrate the validity and effectiveness of the method developed.
基金Supported by the Ministry of Research,Technology,and Higher Education Republic of Indonesia,through the Budget Implementation List(DIPA)of Diponegoro University,Grant No.DIPA-023.04.02.189185/2014,December 05,2013
文摘Optimization analysis and computational fluid dynamics （CFDs） have been applied simultaneously, in which a parametric model plays an important role in finding the optimal solution. However, it is difficult to create a parametric model for a complex shape with irregular curves, such as a submarine hull form. In this study, the cubic Bezier curve and curve-plane intersection method are used to generate a solid model of a parametric submarine hull form taking three input parameters into account： nose radius, tail radius, and length-height hull ratio （L/H）. Application program interface （API） scripting is also used to write code in the ANSYS DesignModeler. The results show that the submarine shape can be generated with some variation of the input parameters. An example is given that shows how the proposed method can be applied successfully to a hull resistance optimization case. The parametric design of the middle submarine type was chosen to be modified. First, the original submarine model was analyzed, in advance, using CFD. Then, using the response surface graph, some candidate optimal designs with a minimum hull resistance coefficient were obtained. Further, the optimization method in goal-driven optimization （GDO） was implemented to find the submarine hull form with the minimum hull resistance coefficient （Ct）. The minimum C, was obtained. The calculated difference in （7, values between the initial submarine and the optimum submarine is around 0.26%, with the C, of the initial submarine and the optimum submarine being 0.001 508 26 and 0.001 504 29, respectively. The results show that the optimum submarine hull form shows a higher nose radius （rn） and higher L/H than those of the initial submarine shape, while the radius of the tail （r1） is smaller than that of the initial shape.
文摘This paper mainly proposes a new C-XSC （C- for eXtended Scientific Computing） software for the symmetric single step method and relaxation method for computing an enclosure for the solution set and compares the methods with others＇ and then makes some modifications and finally, examples illustrating the applicability of the proposed methods are given.
文摘This paper aims to evaluate the stochastic response of steel columns subjected to blast loads using the modified single degree of freedom(MSDOF)method,which assessed towards the conventional single degree of freedom(SDOF)and the experimentally validated Finite Element(FE)methods(LSDYNA).For this purpose,special atten-tion is given to calculating the response of H-shaped steel columns under blast.The damage amount is determined based on the support rotation criterion,which is expressed as a function of their maximum lateral mid-span dis-placement.To account for uncertainties in input parameters and obtain the failure probability,the Monte Carlo simulation(MCS)method is employed,complemented by the Latin Hypercube Sampling(LHS)method to reduce the number of simulations.A parametric analysis is hence performed to examine the effect of several input pa-rameters(including both deterministic and probabilistic parameters)on the probability of column damage as a function of support rotation.First,the MSDOF method confirms its higher accuracy in estimating the probability of column damage due to blast,compared to the conventional SDOF.The collected results also show that un-certainties of several input parameters have significant effects on the column behavior.In particular,geometric parameters(including cross-sectional characteristics,boundary conditions and column length)have major effect on the corresponding column response,in the same way of input blast load parameters and material properties.
