A theoretical study is presented herein on the pen- etration of a semi-infinite target by a spherical-headed long rod for Yp 〉 S, where Yp is the penetrator strength and S is the static target resistance. For Yp 〉 S...A theoretical study is presented herein on the pen- etration of a semi-infinite target by a spherical-headed long rod for Yp 〉 S, where Yp is the penetrator strength and S is the static target resistance. For Yp 〉 S, depending upon initial impact velocity, there exist three types of penetration, namely, penetration by a rigid long rod, penetration by a deforming non-erosive long rod and penetration by an erosive long rod. If the impact velocity of the penetrator is higher than the hydrodynamic velocity (VH), it will penetrate the target in an erosive mode; if the impact velocity lies between the hydrodynamic velocity (VH) and the rigid body velocity (VR), it will penetrate the target in a deformable mode; if the impact velocity is less than the rigid body velocity (VR), it will penetrate the target in a rigid mode. The critical conditions for the transition among these three penetration modes are proposed. It is demonstrated that the present model predictions correlate well with the experimental observations in terms of depth of penetration (DOP) and the critical transition conditions.展开更多
Compliant bistable mechanisms, devices with two distinct stable equilibrium positions, are used in a variety of applications, such as switches, clasps, and valves. Many kinds of compliant bistable mechanisms were prop...Compliant bistable mechanisms, devices with two distinct stable equilibrium positions, are used in a variety of applications, such as switches, clasps, and valves. Many kinds of compliant bistable mechanisms were proposed and studied during the past decade. Among them, tension-based compliant bistable mechanisms, that incorporate tension pivots as their flexible members, feature in short travel distance and low power consuming. So far, the design of this kind of bistable mechanisms is done using finite element method through trial and error, thus is time-consuming. By treating the tension pivots as fixed-guided segments and their elongation as a spring, we developed a novel three degree-of-freedom (3-DOF) pseudo-rigid-body model (PRBM) for this kind of bistable mechanisms. The principle of virtual work is used to derive the force-deflection relationship of the PRBM model. The comparison between the PRBM results and the experimental results of the force-deflection characteristics shows that the PRBM can predict not only the bistable behavior of the tension-based bistable mechanisms, but also their soft spring-like post-bistable behavior and the spring-like force-deflection characteristics when pulling in the reverse direction from the as-fabricated position, which is called reverse behavior. The 3-DOF PRBM can be used to design and identify tension-based bistable mechanisms. Using the PRBM instead of the trial-and-error method can greatly reduce the development time of this kind of bistable mechanisms.展开更多
Energies of the yrast positive-and negative-parity excited states in140 Xe are reproduced by two different models considering quadrupole-octupole deformations, namely the axial vibrational-rotational model and the tri...Energies of the yrast positive-and negative-parity excited states in140 Xe are reproduced by two different models considering quadrupole-octupole deformations, namely the axial vibrational-rotational model and the triaxial rigid rotor model, and compared with the stable octupole-deformed222 Th. The origin of the energy difference between the opposite parity sequences is considered from two different mechanisms, the vibration in axial deformed energy minima and the rotation considering the effective triaxial deformation. The success of reproducing the data in both the models implies that these two mechanisms are equivalent on some level for the octupole-soft nuclei. By investigating the probability distributions for projection of total angular momentum in the triaxial rigid rotor model, it is found that such an energy difference is associated with the difference of orientation of the rotational axis.展开更多
The paper builds a multi-rigid-body model of human with a 4-rigid-body foot in the 3D CAD software Solidworks, based on human anatomy. By controlling the rotation of the ankle and major joints of human body while walk...The paper builds a multi-rigid-body model of human with a 4-rigid-body foot in the 3D CAD software Solidworks, based on human anatomy. By controlling the rotation of the ankle and major joints of human body while walking, the Kinematic simulation was performed in the dynamics simulation software ADAMS. The paper analyzes the simulate results and points out deficiencies in the current work and the direction of research efforts in future.展开更多
Based on the assumption that the slope bodies are rigid, the dynamic model of the landsiding (forward model) was put forward. According to the dynamic model, the system equations of Kalman filter were constituted. The...Based on the assumption that the slope bodies are rigid, the dynamic model of the landsiding (forward model) was put forward. According to the dynamic model, the system equations of Kalman filter were constituted. The mechanical status of a slope was hence combined with the monitoring data by Kalman filter. The model uncertainties or model errors could also be considered through some fictitious observation equations. Different from existed methods, the presented method can make use for not only the statistic information contained in the data but also the information provided by the mechanical and geological aspect of slopes. At last a numerical example was given out to show the feasibility of the method. [展开更多
In 3D models retrieval, feature description and retrieval of non-rigid model face more complex problems due to the isometry transformation of itself. We introduce the hierarchical combination matching into the feature...In 3D models retrieval, feature description and retrieval of non-rigid model face more complex problems due to the isometry transformation of itself. We introduce the hierarchical combination matching into the feature comparison, and build a map between the divided regions of two models, and then achieve accurate feature matching based on patch-by-patch, which successfully introduces the spatial information into feature matching. Verified by experiment, the 3D model retrieval method proposed in this paper based on hierarchical combination matching can make sure more accurate feature matching, so as to enhance the precision of retrieval.展开更多
The dynamics of a coupled rigid-flexible rocket launcher is reported. The coupled rigid-flexible rocket launcher is divided into two subsystems, one is a system of rigid bodies, the other a flexible launch tube which ...The dynamics of a coupled rigid-flexible rocket launcher is reported. The coupled rigid-flexible rocket launcher is divided into two subsystems, one is a system of rigid bodies, the other a flexible launch tube which can undergo large overall motions spatially. First, the mathematical models for these two subsystems were established respectively. Then the dynamic model for the whole system was obtained by considering the coupling effect between these two subsystems. The approach, which divides a complex system into several simple subsystems first and then obtains the dynamic model for the whole system via combining the existing dynamic models for simple subsystems, can make the modeling procedure efficient and convenient.展开更多
Semi-rigid liquid crystal polymer is a class of liquid crystal polymers different from long rigid rod liquid crystal polymer to which the well-known Onsager and Flory theories are applied. In this paper, three statist...Semi-rigid liquid crystal polymer is a class of liquid crystal polymers different from long rigid rod liquid crystal polymer to which the well-known Onsager and Flory theories are applied. In this paper, three statistical models for the semi-rigid nematic polymer were addressed. They are the elastically jointed rod model, worm-like chain model, and non-homogeneous chain model. The nematic-isotropic transition temperature was examined. The pseudo-second transition temperature is expressed analytically. Comparisons with the experiments were made and the agreements were found.展开更多
Based on Maxwell’s constraint counting theory, rigidity percolation in GexSe1-x glasses occurs when the mean coordination number reaches the value of 2.4. This corresponds to Ge0.20Se0.80 glass. At this composition, ...Based on Maxwell’s constraint counting theory, rigidity percolation in GexSe1-x glasses occurs when the mean coordination number reaches the value of 2.4. This corresponds to Ge0.20Se0.80 glass. At this composition, the number of constraints experienced by an atom equals the number of degrees of freedom in three dimensions. Hence, at this composition, the network changes from a floppy phase to a rigid phase, and rigidity starts to percolate. In this work, we use reverse Monte Carlo (RMC) modeling to model the structure of Ge0.20Se0.80 glass by simulating its experimental total atomic pair distribution function (PDF) obtained via high energy synchrotron radiation. A three-dimensional configuration of 2836 atoms was obtained, from which we extracted the partial atomic pair distribution functions associated with Ge-Ge, Ge-Se and Se-Se real space correlations that are hard to extract experimentally from total scattering methods. Bond angle distributions, coordination numbers, mean coordination numbers and the number of floppy modes were also extracted and discussed. More structural insights about network topology at this composition were illustrated. The results indicate that in Ge0.20Se0.80 glass, Ge atoms break up and cross-link the Se chain structure, and form structural units that are four-fold coordinated (the GeSe4 tetrahedra). These tetrahedra form the basic building block and are connected via shared Se atoms or short Se chains. The extent of the intermediate ranged oscillations in real space (as extracted from the width of the first sharp diffraction peak) was found to be around 19.6 ?. The bonding schemes in this glass are consistent with the so-called “8-N” rule and can be interpreted in terms of a chemically ordered network model.展开更多
A unilateral non-penetration constraint dynamical simulation model withfriction is constructed based on compliant model for mechanical system VP (virtual prototyping)simulation. This model combines computer graphics w...A unilateral non-penetration constraint dynamical simulation model withfriction is constructed based on compliant model for mechanical system VP (virtual prototyping)simulation. This model combines computer graphics with multi-body system dynamics. It avoidshandling multiplicity of solution, such as cases of no solution, multi-solution brought about byfriction during traditional construction of non-penetration constraint based on rigid model. At thesame time, the realism of VE (virtual environment) is improved in process of simulation.Furthermore, the valid condition of rolling and sliding unilateral contact is constituted based onsingular perturbation and linear complementary theory. Finally, the compliant method is verified byan interaction between a multi-legged robot and VE.展开更多
The need to find alternative urban mobility solutions for delivery and transport has led mobility companies to devote enormous resources for research-based solutions to increase vehicle safety. This paper documents a ...The need to find alternative urban mobility solutions for delivery and transport has led mobility companies to devote enormous resources for research-based solutions to increase vehicle safety. This paper documents a virtual approach to investigate the influences of different load states to the vehicle dynamic of light electric vehicle. A model basing on a three-dimensional </span><span style="font-family:""><span style="font-family:Verdana;">multibody system was used, which consists of five bodies. By applying methods of multibody modelling the generalized equations of motion were generated. To </span><span style="font-family:Verdana;">include the behavior within the contact point between road and vehicle a</span><span style="font-family:Verdana;"> simplified tire models was added. The implementation of the equations allowed a first validation of the model via simulations. In a final modeling step the simulation results were interpreted in respect of plausibility. Afterwards, the model was simulated numerically to investigate different load states of the vehicle, by applying constant steering stimuli and variable velocities. In sum, the investigated model approach is useful to identify safety relevant parameters and shows the effects of load states to the vehicle dynamics. Furthermore, it behaves plausibly regarding general vehicle dynamics. These results prove </span><span style="font-family:Verdana;">the general usability of the model for the development controllers and esti</span><span style="font-family:Verdana;">mators in driver assistances systems.展开更多
This paper aims to reveal the depth distribution law of non-limit passive soil pressure on rigid retaining wall that rotates about the top of the wall(rotation around the top(RT) model). Based on Coulomb theory, the d...This paper aims to reveal the depth distribution law of non-limit passive soil pressure on rigid retaining wall that rotates about the top of the wall(rotation around the top(RT) model). Based on Coulomb theory, the disturbance degree theory, as well as the spring-element model, by setting the rotation angle of the wall as the disturbance parameter, we establish both a depth distribution function for sand and a nonlinear depth distribution calculation method for the non-limit passive soil pressure on a rigid retaining wall under the RT model, which is then compared with experiment. The results suggest that under the RT model: the non-limit soil pressure has a nonlinear distribution; the backfill disturbance degree and the lateral soil pressure increase with an increase in the wall rotation angle; and, the points where the resultant lateral soil pressure acts on the retaining wall are less than 2/3 of the height of the wall. The soil pressure predicted by the theoretical calculation put forward in this paper are quite similar to those obtained by the model experiment, which verifies the theoretical value, and the engineering guidance provided by the calculations are of significance.展开更多
The characteristics of a disturbed soil-structure interface were studied based on the variation regularities of the disturbed soil within its mining subsidence area using direct shear tests.The effects of the initial ...The characteristics of a disturbed soil-structure interface were studied based on the variation regularities of the disturbed soil within its mining subsidence area using direct shear tests.The effects of the initial moisture content on the shear strength parameters of the soil-structure interfaces were analyzed.The results indicate that the cohesion of the interface initially increased and then decreased as the initial moisture content increased.In addition,the friction angle of the interface decreased as the initial moisture content increased.A constitutive model of the disturbed soil-structure interface,a rigid-plastic model based on the initial void ratio and saturability(VSRP) model,was established based on the results.In order to validate this model,a finite element analysis of DRS-1 direct shear tests was conducted.The finite element model calculations coincided with the results of the DRS-1 direct shear tests.The proposed model also reflected the nonlinear features of the soil-structure interface.展开更多
For submerged vegetated flow, the velocity profile has two distinctive distributions in the vegetation layer in the lower region and the surface layer in the upper non-vegetated region. Based on a mixing-layer analogy...For submerged vegetated flow, the velocity profile has two distinctive distributions in the vegetation layer in the lower region and the surface layer in the upper non-vegetated region. Based on a mixing-layer analogy, different analytical models have been proposed for the velocity profile in the two layers. This paper evaluates the four analytical models of Klopstra et al., Defina & Bixio, Yang et al. and Nepf against a wide range of independent experimental data available in the literature. To test the applicability and robust of the models, the author used the 19 datasets with various relative depths of submergence, different vegetation densities and bed slopes (1.8 × 10?6 - 4.0 × 10?3). This study shows that none of the models can predict the velocity profiles well for all datasets. The three models except Yang’s model performed reasonably well in certain cases, but Yang’s model failed in most the cases studied. It was also found that the Defina model is almost the same as the Klopstra model, if the same mixing length scale of eddies (λ) is used. Finally, close examination of the mixing length scale of eddies (λ) in the Defina model showed that when λ/h = 1/40(H/h)1/2, this model can predict velocity profiles well for all the datasets used.展开更多
基金supported by the National Natural Science Foundation of China (10872195)
文摘A theoretical study is presented herein on the pen- etration of a semi-infinite target by a spherical-headed long rod for Yp 〉 S, where Yp is the penetrator strength and S is the static target resistance. For Yp 〉 S, depending upon initial impact velocity, there exist three types of penetration, namely, penetration by a rigid long rod, penetration by a deforming non-erosive long rod and penetration by an erosive long rod. If the impact velocity of the penetrator is higher than the hydrodynamic velocity (VH), it will penetrate the target in an erosive mode; if the impact velocity lies between the hydrodynamic velocity (VH) and the rigid body velocity (VR), it will penetrate the target in a deformable mode; if the impact velocity is less than the rigid body velocity (VR), it will penetrate the target in a rigid mode. The critical conditions for the transition among these three penetration modes are proposed. It is demonstrated that the present model predictions correlate well with the experimental observations in terms of depth of penetration (DOP) and the critical transition conditions.
基金supported by National Natural Science Foundation of China(Grant No.50805110)Foundation for Key Program of Ministry of Education of China (Grant No.109145)Fundamental Research Funds for the Central Universities of China (Grant No.JY10000904010)
文摘Compliant bistable mechanisms, devices with two distinct stable equilibrium positions, are used in a variety of applications, such as switches, clasps, and valves. Many kinds of compliant bistable mechanisms were proposed and studied during the past decade. Among them, tension-based compliant bistable mechanisms, that incorporate tension pivots as their flexible members, feature in short travel distance and low power consuming. So far, the design of this kind of bistable mechanisms is done using finite element method through trial and error, thus is time-consuming. By treating the tension pivots as fixed-guided segments and their elongation as a spring, we developed a novel three degree-of-freedom (3-DOF) pseudo-rigid-body model (PRBM) for this kind of bistable mechanisms. The principle of virtual work is used to derive the force-deflection relationship of the PRBM model. The comparison between the PRBM results and the experimental results of the force-deflection characteristics shows that the PRBM can predict not only the bistable behavior of the tension-based bistable mechanisms, but also their soft spring-like post-bistable behavior and the spring-like force-deflection characteristics when pulling in the reverse direction from the as-fabricated position, which is called reverse behavior. The 3-DOF PRBM can be used to design and identify tension-based bistable mechanisms. Using the PRBM instead of the trial-and-error method can greatly reduce the development time of this kind of bistable mechanisms.
基金Supported by the National Natural Science Foundation of China under Grant Nos 11675094 and 11622540the Young Scholars Program of Shandong University under Grant No 2015WHWLJH01
文摘Energies of the yrast positive-and negative-parity excited states in140 Xe are reproduced by two different models considering quadrupole-octupole deformations, namely the axial vibrational-rotational model and the triaxial rigid rotor model, and compared with the stable octupole-deformed222 Th. The origin of the energy difference between the opposite parity sequences is considered from two different mechanisms, the vibration in axial deformed energy minima and the rotation considering the effective triaxial deformation. The success of reproducing the data in both the models implies that these two mechanisms are equivalent on some level for the octupole-soft nuclei. By investigating the probability distributions for projection of total angular momentum in the triaxial rigid rotor model, it is found that such an energy difference is associated with the difference of orientation of the rotational axis.
文摘The paper builds a multi-rigid-body model of human with a 4-rigid-body foot in the 3D CAD software Solidworks, based on human anatomy. By controlling the rotation of the ankle and major joints of human body while walking, the Kinematic simulation was performed in the dynamics simulation software ADAMS. The paper analyzes the simulate results and points out deficiencies in the current work and the direction of research efforts in future.
文摘Based on the assumption that the slope bodies are rigid, the dynamic model of the landsiding (forward model) was put forward. According to the dynamic model, the system equations of Kalman filter were constituted. The mechanical status of a slope was hence combined with the monitoring data by Kalman filter. The model uncertainties or model errors could also be considered through some fictitious observation equations. Different from existed methods, the presented method can make use for not only the statistic information contained in the data but also the information provided by the mechanical and geological aspect of slopes. At last a numerical example was given out to show the feasibility of the method. [
基金Supported by National Nature Science Foundation of China(61379106,61379082,61227802)Shandong Provincial Natural Science Foundation(ZR2013FM036,ZR2015FM011,ZR2015FM022)
文摘In 3D models retrieval, feature description and retrieval of non-rigid model face more complex problems due to the isometry transformation of itself. We introduce the hierarchical combination matching into the feature comparison, and build a map between the divided regions of two models, and then achieve accurate feature matching based on patch-by-patch, which successfully introduces the spatial information into feature matching. Verified by experiment, the 3D model retrieval method proposed in this paper based on hierarchical combination matching can make sure more accurate feature matching, so as to enhance the precision of retrieval.
文摘The dynamics of a coupled rigid-flexible rocket launcher is reported. The coupled rigid-flexible rocket launcher is divided into two subsystems, one is a system of rigid bodies, the other a flexible launch tube which can undergo large overall motions spatially. First, the mathematical models for these two subsystems were established respectively. Then the dynamic model for the whole system was obtained by considering the coupling effect between these two subsystems. The approach, which divides a complex system into several simple subsystems first and then obtains the dynamic model for the whole system via combining the existing dynamic models for simple subsystems, can make the modeling procedure efficient and convenient.
基金The work was supported by the Foundation of State Education Committee of China
文摘Semi-rigid liquid crystal polymer is a class of liquid crystal polymers different from long rigid rod liquid crystal polymer to which the well-known Onsager and Flory theories are applied. In this paper, three statistical models for the semi-rigid nematic polymer were addressed. They are the elastically jointed rod model, worm-like chain model, and non-homogeneous chain model. The nematic-isotropic transition temperature was examined. The pseudo-second transition temperature is expressed analytically. Comparisons with the experiments were made and the agreements were found.
文摘Based on Maxwell’s constraint counting theory, rigidity percolation in GexSe1-x glasses occurs when the mean coordination number reaches the value of 2.4. This corresponds to Ge0.20Se0.80 glass. At this composition, the number of constraints experienced by an atom equals the number of degrees of freedom in three dimensions. Hence, at this composition, the network changes from a floppy phase to a rigid phase, and rigidity starts to percolate. In this work, we use reverse Monte Carlo (RMC) modeling to model the structure of Ge0.20Se0.80 glass by simulating its experimental total atomic pair distribution function (PDF) obtained via high energy synchrotron radiation. A three-dimensional configuration of 2836 atoms was obtained, from which we extracted the partial atomic pair distribution functions associated with Ge-Ge, Ge-Se and Se-Se real space correlations that are hard to extract experimentally from total scattering methods. Bond angle distributions, coordination numbers, mean coordination numbers and the number of floppy modes were also extracted and discussed. More structural insights about network topology at this composition were illustrated. The results indicate that in Ge0.20Se0.80 glass, Ge atoms break up and cross-link the Se chain structure, and form structural units that are four-fold coordinated (the GeSe4 tetrahedra). These tetrahedra form the basic building block and are connected via shared Se atoms or short Se chains. The extent of the intermediate ranged oscillations in real space (as extracted from the width of the first sharp diffraction peak) was found to be around 19.6 ?. The bonding schemes in this glass are consistent with the so-called “8-N” rule and can be interpreted in terms of a chemically ordered network model.
基金This project is supported by National Natural Science Foundation of China (No.60375020, No.50305033, No.50405037)973 Program of China (N0.2002CB312106, No.2004CB719400)863 Program of China (No.2003-AA413310)Special Research Foundation of Ministry of Education for PhD Project in Colleges, China (No.20020335112).
文摘A unilateral non-penetration constraint dynamical simulation model withfriction is constructed based on compliant model for mechanical system VP (virtual prototyping)simulation. This model combines computer graphics with multi-body system dynamics. It avoidshandling multiplicity of solution, such as cases of no solution, multi-solution brought about byfriction during traditional construction of non-penetration constraint based on rigid model. At thesame time, the realism of VE (virtual environment) is improved in process of simulation.Furthermore, the valid condition of rolling and sliding unilateral contact is constituted based onsingular perturbation and linear complementary theory. Finally, the compliant method is verified byan interaction between a multi-legged robot and VE.
文摘The need to find alternative urban mobility solutions for delivery and transport has led mobility companies to devote enormous resources for research-based solutions to increase vehicle safety. This paper documents a virtual approach to investigate the influences of different load states to the vehicle dynamic of light electric vehicle. A model basing on a three-dimensional </span><span style="font-family:""><span style="font-family:Verdana;">multibody system was used, which consists of five bodies. By applying methods of multibody modelling the generalized equations of motion were generated. To </span><span style="font-family:Verdana;">include the behavior within the contact point between road and vehicle a</span><span style="font-family:Verdana;"> simplified tire models was added. The implementation of the equations allowed a first validation of the model via simulations. In a final modeling step the simulation results were interpreted in respect of plausibility. Afterwards, the model was simulated numerically to investigate different load states of the vehicle, by applying constant steering stimuli and variable velocities. In sum, the investigated model approach is useful to identify safety relevant parameters and shows the effects of load states to the vehicle dynamics. Furthermore, it behaves plausibly regarding general vehicle dynamics. These results prove </span><span style="font-family:Verdana;">the general usability of the model for the development controllers and esti</span><span style="font-family:Verdana;">mators in driver assistances systems.
基金financially supported by the National Natural Science Foundation of China (No.51274192)Jiangsu Key Laboratory of Environmental Impact and Structural Safety in Engineering Open Foundation of China (No.JSKL2014K12)Jiangsu Ordinary University Graduate Students Research and Innovation Project of China (No.KYLX-1392)
文摘This paper aims to reveal the depth distribution law of non-limit passive soil pressure on rigid retaining wall that rotates about the top of the wall(rotation around the top(RT) model). Based on Coulomb theory, the disturbance degree theory, as well as the spring-element model, by setting the rotation angle of the wall as the disturbance parameter, we establish both a depth distribution function for sand and a nonlinear depth distribution calculation method for the non-limit passive soil pressure on a rigid retaining wall under the RT model, which is then compared with experiment. The results suggest that under the RT model: the non-limit soil pressure has a nonlinear distribution; the backfill disturbance degree and the lateral soil pressure increase with an increase in the wall rotation angle; and, the points where the resultant lateral soil pressure acts on the retaining wall are less than 2/3 of the height of the wall. The soil pressure predicted by the theoretical calculation put forward in this paper are quite similar to those obtained by the model experiment, which verifies the theoretical value, and the engineering guidance provided by the calculations are of significance.
基金Project(51274192)supported by National Natural Science Foundation of China
文摘The characteristics of a disturbed soil-structure interface were studied based on the variation regularities of the disturbed soil within its mining subsidence area using direct shear tests.The effects of the initial moisture content on the shear strength parameters of the soil-structure interfaces were analyzed.The results indicate that the cohesion of the interface initially increased and then decreased as the initial moisture content increased.In addition,the friction angle of the interface decreased as the initial moisture content increased.A constitutive model of the disturbed soil-structure interface,a rigid-plastic model based on the initial void ratio and saturability(VSRP) model,was established based on the results.In order to validate this model,a finite element analysis of DRS-1 direct shear tests was conducted.The finite element model calculations coincided with the results of the DRS-1 direct shear tests.The proposed model also reflected the nonlinear features of the soil-structure interface.
文摘For submerged vegetated flow, the velocity profile has two distinctive distributions in the vegetation layer in the lower region and the surface layer in the upper non-vegetated region. Based on a mixing-layer analogy, different analytical models have been proposed for the velocity profile in the two layers. This paper evaluates the four analytical models of Klopstra et al., Defina & Bixio, Yang et al. and Nepf against a wide range of independent experimental data available in the literature. To test the applicability and robust of the models, the author used the 19 datasets with various relative depths of submergence, different vegetation densities and bed slopes (1.8 × 10?6 - 4.0 × 10?3). This study shows that none of the models can predict the velocity profiles well for all datasets. The three models except Yang’s model performed reasonably well in certain cases, but Yang’s model failed in most the cases studied. It was also found that the Defina model is almost the same as the Klopstra model, if the same mixing length scale of eddies (λ) is used. Finally, close examination of the mixing length scale of eddies (λ) in the Defina model showed that when λ/h = 1/40(H/h)1/2, this model can predict velocity profiles well for all the datasets used.