This paper proposes a new step-by-step Chebyshev space-time spectral method to analyze the force vibration of functionally graded material structures.Although traditional space-time spectral methods can reduce the acc...This paper proposes a new step-by-step Chebyshev space-time spectral method to analyze the force vibration of functionally graded material structures.Although traditional space-time spectral methods can reduce the accuracy mismatch between tem-poral low-order finite difference and spatial high-order discre tization,the ir time collocation points must increase dramatically to solve highly oscillatory solutions of structural vibration,which results in a surge in computing time and a decrease in accuracy.To address this problem,we introduced the step-by-step idea in the space-time spectral method.The Chebyshev polynomials and Lagrange's equation were applied to derive discrete spatial goverming equations,and a matrix projection method was used to map the calculation results of prev ious steps as the initial conditions of the subsequent steps.A series of numerical experiments were carried out.The results of the proposed method were compared with those obtained by traditional space-time spectral methods,which showed that higher accuracy could be achieved in a shorter computation time than the latter in highly oscillatory cases.展开更多
Reducing the peak actuating force(PAF)and parasitic displacement is of high significance for improving the performance of compliant parallel mechanisms(CPMs).In this study,a 2-DOF 4-4R compliant parallel pointing mech...Reducing the peak actuating force(PAF)and parasitic displacement is of high significance for improving the performance of compliant parallel mechanisms(CPMs).In this study,a 2-DOF 4-4R compliant parallel pointing mechanism(4-4R CPPM)was used as the object,and the actuating force of the mechanism was optimized through redundant actuation.This was aimed at minimizing the PAF and parasitic displacement.First,a kinetostatic model of the redundantly actuated 4-4R CPPM was established to reveal the relationship between the input forces/displacements and the output displacements of the mobile platform.Subsequently,based on the established kinetostatic model,methods for optimizing the actuating force distribution with the aim of minimizing the PAF and parasitic displacement were introduced successively.Second,a simulated example of a mobile platform’s spatial pointing trajectory validated the accuracy of the kinetostatic model.The results show a less than 0.9%relative error between the analytical and finite element(FE)results,and the high consistency indicates the accuracy of the kinetostatic model.Then,the effectiveness of the method in minimizing the PAF and parasitic displacement was validated using two simulated examples.The results indicate that compared with the non-redundant actuation case,the PAF of the mechanism could be reduced by up to 50%,and the parasitic displacement was reduced by approximately three-four orders of magnitude by means of redundant actuation combined with the optimal distribution of the actuating force.As expected,with the reduction in parasitic displacement,the FE-results of the output angular displacements(θ_(x) andθ_(z))of the mobile platform were closer to the target oscillation trajectory.This further verified that the reduction in parasitic displacement is indeed effective in improving the motion accuracy of the mechanism.The advantage of this proposed method is that it reduces the PAF and parasitic displacement from the perspective of the actuating force control strategy,without the requirement of structural changes to the original mechanism.展开更多
The influence of supercooled melt forced lamina flow on microsegregation was investigated. The concentration distribution at solid-liquid boundary of binary alloy Ni-Cu was simulated using phase field model coupled wi...The influence of supercooled melt forced lamina flow on microsegregation was investigated. The concentration distribution at solid-liquid boundary of binary alloy Ni-Cu was simulated using phase field model coupled with flow field. The microsegregation, concentration maximum value, boundary thickness of concentration near upstream dendrite and normal to flow dendrite, and downstream dendrite were studied quantitatively in the case of forced lamia flow. The simulation results show that solute field and flow field interact complexly. Compared with melt without flow, in front of upstream dendrite tip, the concentration boundary thickness is the lowest and the concentration maximum value is the smallest for melt with flow. However, in front of downstream dendrite tip, the results are just the opposite. The zone of poor Cu in upstream dendrite where is the most severely microsegregation and shrinkage cavity is wider and the concentration is lower for melt with flow than that without flow.展开更多
A mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic gr...A mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic growth was simulated from undercooled nickel melt under the forced flow. The simulation results show that the asymmetry behavior of the dendritic growth is caused by the forced flow. When the flow velocity is less than the critical value, the asymmetry of dendrite is little influenced by the forced flow. Once the flow velocity reaches or exceeds the critical value, the controlling factor of dendrite growth gradually changes from thermal diffusion to convection. With the increase of the flow velocity, the deflection angle towards upstream direction of the primary dendrite stem becomes larger. The effect of the dendrite growth on the flow field of the melt is apparent. With the increase of the dendrite size, the vortex is present in the downstream regions, and the vortex region is gradually enlarged. Dendrite tips appear to remelt. In addition, the adaptive finite element method can reduce CPU running time by one order of magnitude compared with uniform grid method, and the speed-up ratio is proportional to the size of computational domain.展开更多
Using the concept of the base forces, a new finite element method （base force element method, BFEM） based on the complementary energy principle is presented for accurate modeling of structures with large displacemen...Using the concept of the base forces, a new finite element method （base force element method, BFEM） based on the complementary energy principle is presented for accurate modeling of structures with large displacements and large rotations. First, the complementary energy of an element is described by taking the base forces as state variables, and is then separated into deformation and rotation parts for the case of large deformation. Second, the control equations of the BFEM based on the complementary energy principle are derived using the Lagrange multiplier method. Nonlinear procedure of the BFEM is then developed. Finally, several examples are analyzed to illustrate the reliability and accuracy of the BFEM.展开更多
Atomic force microscopy（AFM） is a commonly used technique for graphene thickness measurement.However, due to surface roughness caused by graphene itself and variation introduced in AFM measurement, graphene thicknes...Atomic force microscopy（AFM） is a commonly used technique for graphene thickness measurement.However, due to surface roughness caused by graphene itself and variation introduced in AFM measurement, graphene thickness is difficult to be accurately determined by AFM. In this paper, a histogram method was used for reliable measurements of graphene thickness using AFM. The influences of various measurement parameters in AFM analysis were investigated. The experimental results indicate that significant deviation can be introduced using various order of flatten and improperly selected measurement parameters including amplitude setpoint and drive amplitude. At amplitude setpoint of 100 mV and drive amplitude of 100 m V, thickness of 1 layer（1L）, 2 layers（2L） and 4 layers（4L） graphene were measured.The height differences for 1L, 2L and 4L were 1.51 ± 0.16 nm, 1.92 ± 0.13 nm and 2.73 ± 0.10 nm, respectively. By comparing these values, thickness of single layer graphene can be accurately determined to be0.41 ± 0.09 nm.展开更多
Cantilever casting concrete arch bridge using form traveller has a broad application prospect.However,it is difficult to obtain reasonable initial cable force in construction stage.In this study,stress balance and inf...Cantilever casting concrete arch bridge using form traveller has a broad application prospect.However,it is difficult to obtain reasonable initial cable force in construction stage.In this study,stress balance and influence matrix methods were developed to determine the initial cable force of cantilever casting concrete arch bridge.The stress balance equation and influence matrix of arch rib critical section were established,and the buckle cable force range was determined by the allowable stress of arch rib critical section.Then a group of buckle cable forces were selected and substituted into the stress balance equation,and the reasonable initial buckle cable force was determined through iteration.Based on the principle of force balance,the initial anchor cable force was determined.In an engineering application example,it is shown that the stress balance and influence matrix methods for the determination of initial cable force are feasible and reliable.The initial cable forces of arch rib segments only need to be adjusted once in the corresponding construction process,which improves the working efficiency and reduces the construction risk.It is found that the methods have great advantages for determining initial cable force in cantilever casting construction process of concrete arch bridge.展开更多
This paper presents a combined application of the finite element method （FEM） and the differential quadrature method （DQM） to vibration and buckling problems of rectangular plates. The proposed scheme combines the...This paper presents a combined application of the finite element method （FEM） and the differential quadrature method （DQM） to vibration and buckling problems of rectangular plates. The proposed scheme combines the geometry flexibility of the FEM and the high accuracy and efficiency of the DQM. The accuracy of the present method is demonstrated by comparing the obtained results with those available in the literature. It is shown that highly accurate results can be obtained by using a small number of finite elements and DQM sample points. The proposed method is suitable for the problems considered due to its simplicity and potential for further development.展开更多
Each joint of hydraulic drive quadruped robot is driven by the hydraulic drive unit(HDU),and the contacting between the robot foot end and the ground is complex and variable,which increases the difficulty of force con...Each joint of hydraulic drive quadruped robot is driven by the hydraulic drive unit(HDU),and the contacting between the robot foot end and the ground is complex and variable,which increases the difficulty of force control inevitably.In the recent years,although many scholars researched some control methods such as disturbance rejection control,parameter self-adaptive control,impedance control and so on,to improve the force control performance of HDU,the robustness of the force control still needs improving.Therefore,how to simulate the complex and variable load characteristics of the environment structure and how to ensure HDU having excellent force control performance with the complex and variable load characteristics are key issues to be solved in this paper.The force control system mathematic model of HDU is established by the mechanism modeling method,and the theoretical models of a novel force control compensation method and a load characteristics simulation method under different environment structures are derived,considering the dynamic characteristics of the load stiffness and the load damping under different environment structures.Then,simulation effects of the variable load stiffness and load damping under the step and sinusoidal load force are analyzed experimentally on the HDU force control performance test platform,which provides the foundation for the force control compensation experiment research.In addition,the optimized PID control parameters are designed to make the HDU have better force control performance with suitable load stiffness and load damping,under which the force control compensation method is introduced,and the robustness of the force control system with several constant load characteristics and the variable load characteristics respectively are comparatively analyzed by experiment.The research results indicate that if the load characteristics are known,the force control compensation method presented in this paper has positive compensation effects on the load characteristics variation,i.e.,this method decreases the effects of the load characteristics variation on the force control performance and enhances the force control system robustness with the constant PID parameters,thereby,the online PID parameters tuning control method which is complex needs not be adopted.All the above research provides theoretical and experimental foundation for the force control method of the quadruped robot joints with high robustness.展开更多
The wheel-rail force measurement is of great importance to the condition monitoring and safety evaluation of railway vehicles. In this paper, an improved indirect method for wheel-rail force measurement is proposed to...The wheel-rail force measurement is of great importance to the condition monitoring and safety evaluation of railway vehicles. In this paper, an improved indirect method for wheel-rail force measurement is proposed to evaluate the running safety of railway vehicles. In this method, the equilibrium equations of a suspended wheelset are derived and the wheel-rail forces are then be obtained from measured suspension and inertia forces. This indirect method avoids structural modifications to the wheelset and is applicable to the long-term operation of railway vehicles. As the wheel-rail lateral forces at two sides of the wheelset are difficult to separate, a new derailment criterion by combined use of wheelset derailment coefficient and wheel unloading ratio is proposed. To illustrate its effectiveness, the indirect method is applied to safety evaluation of rail- way vehicles in different scenarios, such as the cross wind safety of a high-speed train and the safety of a metro vehicle with hunting motions. Then, the feasibility of using this method to identify wheel-rail forces for low-floor light rail vehicles with resilient wheels is discussed. The values identified by this method is compared with that by Simpack simulation for the same low-floor vehicle, which shows a good coincidence between them in the time domain of the wheelset lateral force and the wheel-rail vertical force. In addition, use of the method to determine the high-frequency wheel-rail interaction forces reveals that it is possible to identify the high-frequency wheel-rail forces through the accelerations on the axle box.展开更多
Combining Dual Reciprocity Method （DRM） with Hybrid Boundary Node Method （HBNM）, the Dual Reciprocity Hybrid Boundary Node Method （DRHBNM） is developed for three-dimensional linear elasticity problems with body ...Combining Dual Reciprocity Method （DRM） with Hybrid Boundary Node Method （HBNM）, the Dual Reciprocity Hybrid Boundary Node Method （DRHBNM） is developed for three-dimensional linear elasticity problems with body force. This method can be used to solve the elasticity problems with body force without domain integral, which is inevitable by HBNM. To demonstrate the versatility and the fast convergence of this method, some numerical examples of 3-D elasticity problems with body forces are examined. The computational results show that the present method is effective and can be widely applied in solving practical engineering problems.展开更多
Rolling force and rolling moment are prime process parameter of external spline cold rolling.However,the precise theoretical formulae of rolling force and rolling moment are still very fewer,and the determination of t...Rolling force and rolling moment are prime process parameter of external spline cold rolling.However,the precise theoretical formulae of rolling force and rolling moment are still very fewer,and the determination of them depends on experience.In the present study,the mathematical models of rolling force and rolling moment are established based on stress field theory of slip-line.And the isotropic hardening is used to improve the yield criterion.Based on MATLAB program language environment,calculation program is developed according to mathematical models established.The rolling force and rolling moment could be predicted quickly via the calculation program,and then the reliability of the models is validated by FEM.Within the range of module of spline m=0.5-1.5 mm,pressure angle of reference circleα=30.0°-45.0°,and number of spline teeth Z=19-54,the rolling force and rolling moment in rolling process(finishing rolling is excluded)are researched by means of virtualizing orthogonal experiment design.The results of the present study indicate that:the influences of module and number of spline teeth on the maximum rolling force and rolling moment in the process are remarkable;in the case of pressure angle of reference circle is little,module of spline is great,and number of spline teeth is little,the peak value of rolling force in rolling process may appear in the midst of the process;the peak value of rolling moment in rolling process appears in the midst of the process,and then oscillator weaken to a stable value.The results of the present study may provide guidelines for the determination of power of the motor and the design of hydraulic system of special machine,and provide basis for the farther researches on the precise forming process of external spline cold rolling.展开更多
The character of the deformation zone on pair cross rolls is different from that of the regular 4-high mill. Considering comprehensive influences of normal stress and shear stress in rolling direction, width direction...The character of the deformation zone on pair cross rolls is different from that of the regular 4-high mill. Considering comprehensive influences of normal stress and shear stress in rolling direction, width direction and thickness direction, the rolling force calculation model of the PC （pair crossed） hot strip mill is built. Taking the entry stress and the exit stress of strip as the boundary conditions, the longitudinal and transverse distribution of rolling force is worked out by the differential method. Then the total rolling force is calculated, and calculated results are verified by experimental data.展开更多
Numerical simulation based on a new regularized phase field model was presented to simulate the dendritic shape of a non-isothermal alloy with strong anisotropy in a forced flow. The simulation results show that a cry...Numerical simulation based on a new regularized phase field model was presented to simulate the dendritic shape of a non-isothermal alloy with strong anisotropy in a forced flow. The simulation results show that a crystal nucleus grows into a symmetric dendrite in a free flow and into an asymmetry dendrite in a forced flow. As the forced flow velocity is increased, both of the promoting effect on the upstream arm and the inhibiting effects on the downstream and perpendicular arms are intensified, and the perpendicular arm tilts to the upstream direction. With increasing the anisotropy value to 0.14, all of the dendrite arms tip velocities are gradually stabilized and finally reach their relative saturation values. In addition, the effects of an undercooling parameter and a forced compound flow on the faceted dendrite growth were also investigated.展开更多
In order to evaluate the seismic stability of reinforced soil walls against bearing capacity failure,the seismic safety factor of reinforced soil walls was determined by using pseudo-dynamic method,and calculated by c...In order to evaluate the seismic stability of reinforced soil walls against bearing capacity failure,the seismic safety factor of reinforced soil walls was determined by using pseudo-dynamic method,and calculated by considering different parameters,such as horizontal and vertical seismic acceleration coefficients,ratio of reinforcement length to wall height,back fill friction angle,foundation soil friction angle,soil reinforcement interface friction angle and surcharge.The parametric study shows that the seismic safety factor increases by 24-fold when the foundation soil friction angle varies from 25°to 45°,and increases by 2-fold when the soil reinforcement interface friction angle varies from 0 to 30°.That is to say,the bigger values the foundation soil and/or soil reinforcement interface friction angles have,the safer the reinforced soil walls become in the seismic design.The results were also compared with those obtained from pseudo-static method.It is found that there is a higher value of the safety factor by the present work.展开更多
This paper presents the elastic and plastic deformation of the steel helmet with coldextrusion moulding. The plastic streamline of the plastic mould-making process for ellipse thinplate is described. The distribution ...This paper presents the elastic and plastic deformation of the steel helmet with coldextrusion moulding. The plastic streamline of the plastic mould-making process for ellipse thinplate is described. The distribution of slip-line is established based on the plastic streamline. Theextrusion force of plastic moulding of the steel helmet is calculated by using of slip-line method.Furthermore, an applied example is given.展开更多
In loosely coupled or large-scale problems with high dominance ratios,slow fission source convergence can take extremely long time,reducing Monte Carlo(MC)criticality calculation efficiency.Although various accelerati...In loosely coupled or large-scale problems with high dominance ratios,slow fission source convergence can take extremely long time,reducing Monte Carlo(MC)criticality calculation efficiency.Although various acceleration methods have been developed,some methods cannot reduce convergence times,whereas others have been limited to specific problem geometries.In this study,a new fission source convergence acceleration(FSCA)method,the forced propagation(FP)method,has been proposed,which forces the fission source to propagate and accelerate fission source convergence.Additionally,some stabilization techniques have been designed to render the method more practical.The resulting stabilized method was then successfully implemented in the MC transport code,and its feasibility and effectiveness were tested using the modified OECD/NEA,one-dimensional slab benchmark,and the Hoogenboom full-core problem.The comparison results showed that the FP method was able to achieve efficient FSCA.展开更多
A new analytical method is proposed to analyze the force acting on a rectangular oscillating buoy due to linear waves.In the method a new analytical expression for the diffraction velocity potential is obtained first ...A new analytical method is proposed to analyze the force acting on a rectangular oscillating buoy due to linear waves.In the method a new analytical expression for the diffraction velocity potential is obtained first by use of theeigenfunction expansion method and then the wave excitation force is calculated by use of the known incident wavepotential and the diffraction potential. Compared with the classical analytical method, it can be seen that the presentmethod is simpler for a two-dimensional problem due to the comparable effort needed for the computation ofdiffraction potential and for that of radiated potential. To verify the correctness of the method, a classical example inthe reference is recomputed and the obtained results are in good accordance with those by use of other methods,which shows that the present method is correct.展开更多
Nonlinear dynamic equation is a common engineering model.There is not precise analytical solution for most of nonlinear differential equations.These nonlinear differential equations should be solved by using approxima...Nonlinear dynamic equation is a common engineering model.There is not precise analytical solution for most of nonlinear differential equations.These nonlinear differential equations should be solved by using approximate methods.Classical perturbation methods such as LP method,KBM method,multi-scale method and the averaging method on weakly nonlinear vibration system is effective,while the strongly nonlinear system is difficult to apply.Approximate solutions of primary resonance for forced Duffing equation is investigated by means of homotopy analysis method （HAM）.Different from other approximate computational method,the HAM is totally independent of small physical parameters,and thus is suitable for most nonlinear problems.The HAM provides a great freedom to choose base functions of solution series,so that a nonlinear problem may be approximated more effectively.The HAM provides us a simple way to adjust and control the convergence region of the series solution by means of an auxiliary parameter and the auxiliary function.Therefore,HAM not only may solve the weakly non-linear problems but also may be suitable for the strong non-linear problem.Through the approximate solution of forced Duffing equation with cubic non-linearity,the HAM and fourth order Runge-Kutta method of numerical solution were compared,the results show that the HAM not only can solve the steady state solution,but also can calculate the unsteady state solution,and has the good computational accuracy.展开更多
Ship collision on bridge is a dynamic process featured by high nonlinearity and instantaneity. Calculating ship-bridge collision force typically involves either the use of design-specification-stipulated equivalent st...Ship collision on bridge is a dynamic process featured by high nonlinearity and instantaneity. Calculating ship-bridge collision force typically involves either the use of design-specification-stipulated equivalent static load, or the use of finite element method (FEM) which is more time-consuming and requires supercomputing resources. In this paper, we proposed an alternative approach that combines FEM with artificial neural network (ANN). The radial basis function neural network (RBFNN) employed for calculating the impact force in consideration of ship-bridge collision mechanics. With ship velocity and mass as the input vectors and ship collision force as the output vector, the neural networks for different network parameters are trained by the learning samples obtained from finite element simulation results. The error analyses of the learning and testing samples show that the proposed RBFNN is accurate enough to calculate ship-bridge collision force. The input-output relationship obtained by the RBFNN is essentially consistent with the typical empirical formulae. Finally, a special toolbox is developed for calculation efficiency in application using MATLAB software.展开更多
基金supported by the Advance Research Project of Civil Aerospace Technology(Grant No.D020304)National Nat-ural Science Foundation of China(Grant Nos.52205257 and U22B2083).
文摘This paper proposes a new step-by-step Chebyshev space-time spectral method to analyze the force vibration of functionally graded material structures.Although traditional space-time spectral methods can reduce the accuracy mismatch between tem-poral low-order finite difference and spatial high-order discre tization,the ir time collocation points must increase dramatically to solve highly oscillatory solutions of structural vibration,which results in a surge in computing time and a decrease in accuracy.To address this problem,we introduced the step-by-step idea in the space-time spectral method.The Chebyshev polynomials and Lagrange's equation were applied to derive discrete spatial goverming equations,and a matrix projection method was used to map the calculation results of prev ious steps as the initial conditions of the subsequent steps.A series of numerical experiments were carried out.The results of the proposed method were compared with those obtained by traditional space-time spectral methods,which showed that higher accuracy could be achieved in a shorter computation time than the latter in highly oscillatory cases.
基金Supported by Key Project of Hubei Provincial Department of Education Research Program(Grant No.D20211401).
文摘Reducing the peak actuating force(PAF)and parasitic displacement is of high significance for improving the performance of compliant parallel mechanisms(CPMs).In this study,a 2-DOF 4-4R compliant parallel pointing mechanism(4-4R CPPM)was used as the object,and the actuating force of the mechanism was optimized through redundant actuation.This was aimed at minimizing the PAF and parasitic displacement.First,a kinetostatic model of the redundantly actuated 4-4R CPPM was established to reveal the relationship between the input forces/displacements and the output displacements of the mobile platform.Subsequently,based on the established kinetostatic model,methods for optimizing the actuating force distribution with the aim of minimizing the PAF and parasitic displacement were introduced successively.Second,a simulated example of a mobile platform’s spatial pointing trajectory validated the accuracy of the kinetostatic model.The results show a less than 0.9%relative error between the analytical and finite element(FE)results,and the high consistency indicates the accuracy of the kinetostatic model.Then,the effectiveness of the method in minimizing the PAF and parasitic displacement was validated using two simulated examples.The results indicate that compared with the non-redundant actuation case,the PAF of the mechanism could be reduced by up to 50%,and the parasitic displacement was reduced by approximately three-four orders of magnitude by means of redundant actuation combined with the optimal distribution of the actuating force.As expected,with the reduction in parasitic displacement,the FE-results of the output angular displacements(θ_(x) andθ_(z))of the mobile platform were closer to the target oscillation trajectory.This further verified that the reduction in parasitic displacement is indeed effective in improving the motion accuracy of the mechanism.The advantage of this proposed method is that it reduces the PAF and parasitic displacement from the perspective of the actuating force control strategy,without the requirement of structural changes to the original mechanism.
基金Project (10964004) supported by the National Natural Science Foundation of ChinaProject (20070731001) supported by Research Fund for the Doctoral Program of ChinaProject (096RJZA104) supported by the Natural Science Foundation of Gansu Province,China
文摘The influence of supercooled melt forced lamina flow on microsegregation was investigated. The concentration distribution at solid-liquid boundary of binary alloy Ni-Cu was simulated using phase field model coupled with flow field. The microsegregation, concentration maximum value, boundary thickness of concentration near upstream dendrite and normal to flow dendrite, and downstream dendrite were studied quantitatively in the case of forced lamia flow. The simulation results show that solute field and flow field interact complexly. Compared with melt without flow, in front of upstream dendrite tip, the concentration boundary thickness is the lowest and the concentration maximum value is the smallest for melt with flow. However, in front of downstream dendrite tip, the results are just the opposite. The zone of poor Cu in upstream dendrite where is the most severely microsegregation and shrinkage cavity is wider and the concentration is lower for melt with flow than that without flow.
基金Projects(51161011,11364024)supported by the National Natural Science Foundation of ChinaProject(1204GKCA065)supported by the Key Technology R&D Program of Gansu Province,China+1 种基金Project(201210)supported by the Fundamental Research Funds for the Universities of Gansu Province,ChinaProject(J201304)supported by the Funds for Distinguished Young Scientists of Lanzhou University of Technology,China
文摘A mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic growth was simulated from undercooled nickel melt under the forced flow. The simulation results show that the asymmetry behavior of the dendritic growth is caused by the forced flow. When the flow velocity is less than the critical value, the asymmetry of dendrite is little influenced by the forced flow. Once the flow velocity reaches or exceeds the critical value, the controlling factor of dendrite growth gradually changes from thermal diffusion to convection. With the increase of the flow velocity, the deflection angle towards upstream direction of the primary dendrite stem becomes larger. The effect of the dendrite growth on the flow field of the melt is apparent. With the increase of the dendrite size, the vortex is present in the downstream regions, and the vortex region is gradually enlarged. Dendrite tips appear to remelt. In addition, the adaptive finite element method can reduce CPU running time by one order of magnitude compared with uniform grid method, and the speed-up ratio is proportional to the size of computational domain.
基金supported by the China Postdoctoral Science Foundation Funded Project (20080430038) the Funding Project for Academic Human Resources Development in Institutions of Higher Learning Under the Jurisdiction of Beijing Municipality (05004999200602)
文摘Using the concept of the base forces, a new finite element method （base force element method, BFEM） based on the complementary energy principle is presented for accurate modeling of structures with large displacements and large rotations. First, the complementary energy of an element is described by taking the base forces as state variables, and is then separated into deformation and rotation parts for the case of large deformation. Second, the control equations of the BFEM based on the complementary energy principle are derived using the Lagrange multiplier method. Nonlinear procedure of the BFEM is then developed. Finally, several examples are analyzed to illustrate the reliability and accuracy of the BFEM.
基金supported by the Program of National Key Technology R&D Program of the Ministry of Science and Technology of China(2011BAK15B04)
文摘Atomic force microscopy（AFM） is a commonly used technique for graphene thickness measurement.However, due to surface roughness caused by graphene itself and variation introduced in AFM measurement, graphene thickness is difficult to be accurately determined by AFM. In this paper, a histogram method was used for reliable measurements of graphene thickness using AFM. The influences of various measurement parameters in AFM analysis were investigated. The experimental results indicate that significant deviation can be introduced using various order of flatten and improperly selected measurement parameters including amplitude setpoint and drive amplitude. At amplitude setpoint of 100 mV and drive amplitude of 100 m V, thickness of 1 layer（1L）, 2 layers（2L） and 4 layers（4L） graphene were measured.The height differences for 1L, 2L and 4L were 1.51 ± 0.16 nm, 1.92 ± 0.13 nm and 2.73 ± 0.10 nm, respectively. By comparing these values, thickness of single layer graphene can be accurately determined to be0.41 ± 0.09 nm.
基金Projects(51478049,51778068)supported by the National Natural Science Foundation of ChinaProject(14JJ2075,2019JJ40301)supported by the Hunan Natural Science Foundation of China+1 种基金Project(17A010)supported by the Scientific Research Fund of Hunan Provincial Education Department of ChinaProject(2017GK4034)supported by the Major Technological Achievements Transformation Program of Hunan Strategic Emerging Industries of China
文摘Cantilever casting concrete arch bridge using form traveller has a broad application prospect.However,it is difficult to obtain reasonable initial cable force in construction stage.In this study,stress balance and influence matrix methods were developed to determine the initial cable force of cantilever casting concrete arch bridge.The stress balance equation and influence matrix of arch rib critical section were established,and the buckle cable force range was determined by the allowable stress of arch rib critical section.Then a group of buckle cable forces were selected and substituted into the stress balance equation,and the reasonable initial buckle cable force was determined through iteration.Based on the principle of force balance,the initial anchor cable force was determined.In an engineering application example,it is shown that the stress balance and influence matrix methods for the determination of initial cable force are feasible and reliable.The initial cable forces of arch rib segments only need to be adjusted once in the corresponding construction process,which improves the working efficiency and reduces the construction risk.It is found that the methods have great advantages for determining initial cable force in cantilever casting construction process of concrete arch bridge.
文摘This paper presents a combined application of the finite element method （FEM） and the differential quadrature method （DQM） to vibration and buckling problems of rectangular plates. The proposed scheme combines the geometry flexibility of the FEM and the high accuracy and efficiency of the DQM. The accuracy of the present method is demonstrated by comparing the obtained results with those available in the literature. It is shown that highly accurate results can be obtained by using a small number of finite elements and DQM sample points. The proposed method is suitable for the problems considered due to its simplicity and potential for further development.
基金Supported by National Key Basic Research Program of China(973 Program,Grant No.2014CB046405)State Key Laboratory of Fluid Power and Mechatronic Systems(Zhejiang University)Open Fund Project(Grant No.GZKF-201502)Hebei Military and Civilian Industry Development Funds Projects of China(Grant No.2015B060)
文摘Each joint of hydraulic drive quadruped robot is driven by the hydraulic drive unit(HDU),and the contacting between the robot foot end and the ground is complex and variable,which increases the difficulty of force control inevitably.In the recent years,although many scholars researched some control methods such as disturbance rejection control,parameter self-adaptive control,impedance control and so on,to improve the force control performance of HDU,the robustness of the force control still needs improving.Therefore,how to simulate the complex and variable load characteristics of the environment structure and how to ensure HDU having excellent force control performance with the complex and variable load characteristics are key issues to be solved in this paper.The force control system mathematic model of HDU is established by the mechanism modeling method,and the theoretical models of a novel force control compensation method and a load characteristics simulation method under different environment structures are derived,considering the dynamic characteristics of the load stiffness and the load damping under different environment structures.Then,simulation effects of the variable load stiffness and load damping under the step and sinusoidal load force are analyzed experimentally on the HDU force control performance test platform,which provides the foundation for the force control compensation experiment research.In addition,the optimized PID control parameters are designed to make the HDU have better force control performance with suitable load stiffness and load damping,under which the force control compensation method is introduced,and the robustness of the force control system with several constant load characteristics and the variable load characteristics respectively are comparatively analyzed by experiment.The research results indicate that if the load characteristics are known,the force control compensation method presented in this paper has positive compensation effects on the load characteristics variation,i.e.,this method decreases the effects of the load characteristics variation on the force control performance and enhances the force control system robustness with the constant PID parameters,thereby,the online PID parameters tuning control method which is complex needs not be adopted.All the above research provides theoretical and experimental foundation for the force control method of the quadruped robot joints with high robustness.
基金supported by the National Natural Science Foundation of China (Grant No. U1334206 and No. 51475388)Science & Technology Development Project of China Railway Corporation (Grant No. J012-C)
文摘The wheel-rail force measurement is of great importance to the condition monitoring and safety evaluation of railway vehicles. In this paper, an improved indirect method for wheel-rail force measurement is proposed to evaluate the running safety of railway vehicles. In this method, the equilibrium equations of a suspended wheelset are derived and the wheel-rail forces are then be obtained from measured suspension and inertia forces. This indirect method avoids structural modifications to the wheelset and is applicable to the long-term operation of railway vehicles. As the wheel-rail lateral forces at two sides of the wheelset are difficult to separate, a new derailment criterion by combined use of wheelset derailment coefficient and wheel unloading ratio is proposed. To illustrate its effectiveness, the indirect method is applied to safety evaluation of rail- way vehicles in different scenarios, such as the cross wind safety of a high-speed train and the safety of a metro vehicle with hunting motions. Then, the feasibility of using this method to identify wheel-rail forces for low-floor light rail vehicles with resilient wheels is discussed. The values identified by this method is compared with that by Simpack simulation for the same low-floor vehicle, which shows a good coincidence between them in the time domain of the wheelset lateral force and the wheel-rail vertical force. In addition, use of the method to determine the high-frequency wheel-rail interaction forces reveals that it is possible to identify the high-frequency wheel-rail forces through the accelerations on the axle box.
文摘Combining Dual Reciprocity Method （DRM） with Hybrid Boundary Node Method （HBNM）, the Dual Reciprocity Hybrid Boundary Node Method （DRHBNM） is developed for three-dimensional linear elasticity problems with body force. This method can be used to solve the elasticity problems with body force without domain integral, which is inevitable by HBNM. To demonstrate the versatility and the fast convergence of this method, some numerical examples of 3-D elasticity problems with body forces are examined. The computational results show that the present method is effective and can be widely applied in solving practical engineering problems.
基金supported by National Natural Science Foundation of China(Grant No.50675145)Shanxi Provincial Key Project of Science and Technology of China(Grant No.2006031147)+1 种基金Shanxi Provincial Innovation Project for Graduate Students of China(Grant No.20061027)Shanxi Provincial Key Project for Studied-abroad Returnee of China
文摘Rolling force and rolling moment are prime process parameter of external spline cold rolling.However,the precise theoretical formulae of rolling force and rolling moment are still very fewer,and the determination of them depends on experience.In the present study,the mathematical models of rolling force and rolling moment are established based on stress field theory of slip-line.And the isotropic hardening is used to improve the yield criterion.Based on MATLAB program language environment,calculation program is developed according to mathematical models established.The rolling force and rolling moment could be predicted quickly via the calculation program,and then the reliability of the models is validated by FEM.Within the range of module of spline m=0.5-1.5 mm,pressure angle of reference circleα=30.0°-45.0°,and number of spline teeth Z=19-54,the rolling force and rolling moment in rolling process(finishing rolling is excluded)are researched by means of virtualizing orthogonal experiment design.The results of the present study indicate that:the influences of module and number of spline teeth on the maximum rolling force and rolling moment in the process are remarkable;in the case of pressure angle of reference circle is little,module of spline is great,and number of spline teeth is little,the peak value of rolling force in rolling process may appear in the midst of the process;the peak value of rolling moment in rolling process appears in the midst of the process,and then oscillator weaken to a stable value.The results of the present study may provide guidelines for the determination of power of the motor and the design of hydraulic system of special machine,and provide basis for the farther researches on the precise forming process of external spline cold rolling.
基金Item Sponsored by ("863"Program) of China National High Technology Research and Development Program(2009AA04Z143)
文摘The character of the deformation zone on pair cross rolls is different from that of the regular 4-high mill. Considering comprehensive influences of normal stress and shear stress in rolling direction, width direction and thickness direction, the rolling force calculation model of the PC （pair crossed） hot strip mill is built. Taking the entry stress and the exit stress of strip as the boundary conditions, the longitudinal and transverse distribution of rolling force is worked out by the differential method. Then the total rolling force is calculated, and calculated results are verified by experimental data.
基金Project(11102164)supported by the National Natural Science Foundation of ChinaProject(G9KY101502)supported by NPU Foundation for Fundamental Research,China
文摘Numerical simulation based on a new regularized phase field model was presented to simulate the dendritic shape of a non-isothermal alloy with strong anisotropy in a forced flow. The simulation results show that a crystal nucleus grows into a symmetric dendrite in a free flow and into an asymmetry dendrite in a forced flow. As the forced flow velocity is increased, both of the promoting effect on the upstream arm and the inhibiting effects on the downstream and perpendicular arms are intensified, and the perpendicular arm tilts to the upstream direction. With increasing the anisotropy value to 0.14, all of the dendrite arms tip velocities are gradually stabilized and finally reach their relative saturation values. In addition, the effects of an undercooling parameter and a forced compound flow on the faceted dendrite growth were also investigated.
文摘In order to evaluate the seismic stability of reinforced soil walls against bearing capacity failure,the seismic safety factor of reinforced soil walls was determined by using pseudo-dynamic method,and calculated by considering different parameters,such as horizontal and vertical seismic acceleration coefficients,ratio of reinforcement length to wall height,back fill friction angle,foundation soil friction angle,soil reinforcement interface friction angle and surcharge.The parametric study shows that the seismic safety factor increases by 24-fold when the foundation soil friction angle varies from 25°to 45°,and increases by 2-fold when the soil reinforcement interface friction angle varies from 0 to 30°.That is to say,the bigger values the foundation soil and/or soil reinforcement interface friction angles have,the safer the reinforced soil walls become in the seismic design.The results were also compared with those obtained from pseudo-static method.It is found that there is a higher value of the safety factor by the present work.
文摘This paper presents the elastic and plastic deformation of the steel helmet with coldextrusion moulding. The plastic streamline of the plastic mould-making process for ellipse thinplate is described. The distribution of slip-line is established based on the plastic streamline. Theextrusion force of plastic moulding of the steel helmet is calculated by using of slip-line method.Furthermore, an applied example is given.
基金supported by the National Natural Science Foundation of China(Nos.11775126,11545013,11605101)the Young Elite Scientists Sponsorship Program by CAST(No.2016QNRC001)+1 种基金Science Challenge Project by MIIT of China(No.TZ2018001)Tsinghua University,Initiative Scientific Research Program。
文摘In loosely coupled or large-scale problems with high dominance ratios,slow fission source convergence can take extremely long time,reducing Monte Carlo(MC)criticality calculation efficiency.Although various acceleration methods have been developed,some methods cannot reduce convergence times,whereas others have been limited to specific problem geometries.In this study,a new fission source convergence acceleration(FSCA)method,the forced propagation(FP)method,has been proposed,which forces the fission source to propagate and accelerate fission source convergence.Additionally,some stabilization techniques have been designed to render the method more practical.The resulting stabilized method was then successfully implemented in the MC transport code,and its feasibility and effectiveness were tested using the modified OECD/NEA,one-dimensional slab benchmark,and the Hoogenboom full-core problem.The comparison results showed that the FP method was able to achieve efficient FSCA.
基金This work Was supported by the High Tech Research and Development(863)Program of China under Grant No.2003AA5 16010the Chinese Academy of Science Pilot Project of the National Knowledge Innovation Program under Grant No.KGCX2-SW-305Chinese National Science Fund for Distinguished Young Scholars under Grant No.50125924.
文摘A new analytical method is proposed to analyze the force acting on a rectangular oscillating buoy due to linear waves.In the method a new analytical expression for the diffraction velocity potential is obtained first by use of theeigenfunction expansion method and then the wave excitation force is calculated by use of the known incident wavepotential and the diffraction potential. Compared with the classical analytical method, it can be seen that the presentmethod is simpler for a two-dimensional problem due to the comparable effort needed for the computation ofdiffraction potential and for that of radiated potential. To verify the correctness of the method, a classical example inthe reference is recomputed and the obtained results are in good accordance with those by use of other methods,which shows that the present method is correct.
基金supported by Fundamental Research Funds for the Central Universities of China (Grant No. N090405009)
文摘Nonlinear dynamic equation is a common engineering model.There is not precise analytical solution for most of nonlinear differential equations.These nonlinear differential equations should be solved by using approximate methods.Classical perturbation methods such as LP method,KBM method,multi-scale method and the averaging method on weakly nonlinear vibration system is effective,while the strongly nonlinear system is difficult to apply.Approximate solutions of primary resonance for forced Duffing equation is investigated by means of homotopy analysis method （HAM）.Different from other approximate computational method,the HAM is totally independent of small physical parameters,and thus is suitable for most nonlinear problems.The HAM provides a great freedom to choose base functions of solution series,so that a nonlinear problem may be approximated more effectively.The HAM provides us a simple way to adjust and control the convergence region of the series solution by means of an auxiliary parameter and the auxiliary function.Therefore,HAM not only may solve the weakly non-linear problems but also may be suitable for the strong non-linear problem.Through the approximate solution of forced Duffing equation with cubic non-linearity,the HAM and fourth order Runge-Kutta method of numerical solution were compared,the results show that the HAM not only can solve the steady state solution,but also can calculate the unsteady state solution,and has the good computational accuracy.
基金the National Natural Science Foundation of China (No. 50778131)the National key Technology R&D Pro-gram, Ministry of Science and Technology (No. 2006BAG04B01), China
文摘Ship collision on bridge is a dynamic process featured by high nonlinearity and instantaneity. Calculating ship-bridge collision force typically involves either the use of design-specification-stipulated equivalent static load, or the use of finite element method (FEM) which is more time-consuming and requires supercomputing resources. In this paper, we proposed an alternative approach that combines FEM with artificial neural network (ANN). The radial basis function neural network (RBFNN) employed for calculating the impact force in consideration of ship-bridge collision mechanics. With ship velocity and mass as the input vectors and ship collision force as the output vector, the neural networks for different network parameters are trained by the learning samples obtained from finite element simulation results. The error analyses of the learning and testing samples show that the proposed RBFNN is accurate enough to calculate ship-bridge collision force. The input-output relationship obtained by the RBFNN is essentially consistent with the typical empirical formulae. Finally, a special toolbox is developed for calculation efficiency in application using MATLAB software.
