A forward recursive formulation based on corotational frame is proposed for flexible planar beams with large displacement.The traditional recursive formulation has been successfully used for flexible mutibody dynamics...A forward recursive formulation based on corotational frame is proposed for flexible planar beams with large displacement.The traditional recursive formulation has been successfully used for flexible mutibody dynamics to improve the computational efficiency based on floating frame,in which the assumption of small strain and deflection is adopted.The proposed recursive formulation could be used for large displacement problems based on the corotational frame.It means that the recursive scheme is used not only for adjacent bodies but also for adjacent beam elements.The nodal relative rotation coordinates of the planar beam are used to obtain equations with minimal generalized coordinates in present formulation.The proposed formulation is different from absolute nodal coordinate formulation and the geometrically exact beam formulation in which the absolute coordinates are used.The recursive scheme and minimal set of dynamic equations lead to a high computational efficiency in numerical integration.Numerical examples are carried out to demonstrate the accuracy and validity of this formulation.For all of the examples,the results of the present formulation are in good agreement with results obtained using commercial software and the published results.Moreover,it is shown that the present formulation is more efficient than the formulation in ANSYS based on GEBF.展开更多
A Bernoulli-Euler beam mechanism for static analysis of large displacement,large rotation but small strain planar tapered beam structures is proposed using the Updated Lagrangian formulation and the moving coordinate ...A Bernoulli-Euler beam mechanism for static analysis of large displacement,large rotation but small strain planar tapered beam structures is proposed using the Updated Lagrangian formulation and the moving coordinate method.The object beam is the tapered one whose profile is assumed to be varying linearly.From the governing differential equation of lateral deflection including second-order effects by beam-column theory,the geometric nonlinear tangent elemental stiffness matrix is derived.The nonlinear effect of the bending distortions on the axial action is considered to manifest itself as an axial change in length.The aforementioned stiffness matrix is amended,by developing the auxiliary stiffness of bowing effect.The moving coordinate method is employed for obtaining the large displacement total equilibrium equations,and the hinged-hinged moving coordinate system is constructed at the last updated configuration.The multiple load steps Newton-Raphson iteration is adopted for the solution of the nonlinear equations.The validity and efficiency of the proposed method are shown by solving various typical numerical examples.展开更多
A study to estimate land surface movement caused by large surface excavations in sedimentary strata is presented.In stratified or jointed strata the stress relief driven movement adjacent to large excavations can be s...A study to estimate land surface movement caused by large surface excavations in sedimentary strata is presented.In stratified or jointed strata the stress relief driven movement adjacent to large excavations can be significantly larger than expected.High lateral stresses measured in Australia and other places around the world indicate that the ratio of horizontal to vertical stress has been particularly high at shallow depths.The in situ strata is in compression and during excavation,stress is relieved towards the opening causing strata movement.Large excavations such as,open cut mines or highway cuttings,can initiate an extensive horizontal slide of surface layers towards the excavation.These ground movements can be damaging to surface structures such as water storage dams and large buildings.Based on stress measurements at shallow depths in Australian coal mines the study presented here calculates and models the extent of potential ground movement along the bedding surface adjacent to large excavations and provides a new prediction tool of land movement at the excavation boundary that can benefit the geotechnical practitioners in the mining industry.展开更多
The variational functional of the Hellinger-Reissner variational principle for the large displacement problem of a thin shallow shell with an arbitrary shape is first established. Then the functional of the modified p...The variational functional of the Hellinger-Reissner variational principle for the large displacement problem of a thin shallow shell with an arbitrary shape is first established. Then the functional of the modified principle suitable for the finite element method is derived. In the functional only two independent variables, the deflection w and the stress function F are included. The displacement expressions in the middle surface on the boundary of the shell is also derived by means of the previous two variables.展开更多
In this paper, a nonlinear mathematical model for analyzing dynamical response to the large deformation of piles with initial displacements is firstly established with the arc-coordinate, and it is a set of nonlinear ...In this paper, a nonlinear mathematical model for analyzing dynamical response to the large deformation of piles with initial displacements is firstly established with the arc-coordinate, and it is a set of nonlinear integral-differential equa- tions, in which, the Winkeler model is used to simulate the resistance of the soil to the pile. Secondly, a set of new auxiliary functions are introduced. The differential-integral equations are transformed into a set of nonlinear differential equations, and the differential quadrature method (DQM) and the finite difference method (FDM) are applied to discretize the set of nonlinear equations in the spatial and time domains, respectively. Then, the Newton-Raphson method is used to solve the set of discretization algebraic equations at each time step. Finally, numerical examples are presented, and the dynamical re- sponses to the deformation of piles, including configuration, bending moment and shear force, are graphically illuminated. In calculation, two types of initial displacements and dynamical loads are applied, and the effects of parameters on the dynamical responses of piles are analyzed in detail.展开更多
The fundamental equations for analysis of a straight slender rod which deform at large deflections of order of several times of the dimension of the cross section of the rod are derived by using convecled coordinate s...The fundamental equations for analysis of a straight slender rod which deform at large deflections of order of several times of the dimension of the cross section of the rod are derived by using convecled coordinate system. In accordance with the practice of oil drilling, an effective method of three dimensional static analysis of bottom-hole assemblies is simply described. Errors in reference [8 ] are pointed out.展开更多
In order to further deepen the understanding of seismic wave propagation characteristics induced by the large volume airgun source,experimental data from multiple fixed excitation points in Fujian Province were used t...In order to further deepen the understanding of seismic wave propagation characteristics induced by the large volume airgun source,experimental data from multiple fixed excitation points in Fujian Province were used to obtain the equivalent single excitation high signal-to-noise ratio velocity and displacement records through linear stacking and simulation techniques.Then the peak displacements of different epicentral distances were used to calculate the equivalent magnitude of the airgun source excitation at different fixed excitation points so as to establish the attenuation relationship between equivalent magnitude,epicenter distance and velocity peak.Our results show that:①Within 270 km of epicentral distance,for the large-volume airgun’s single shot,the peak velocity range is about 700-4 nm/s,and the peak displacement range is about 200.0-0.2 nm;②The equivalent magnitude of the P-wave from the airgun source with a total capacity of 8,000 in 3 is 0.181-0.760,and the equivalent magnitude of the S-wave is 0.294-0.832.By contrast,the equivalent magnitude of the P-wave from the airgun source with a total capacity of 12,000 in 3 is 0.533-0.896,and the equivalent magnitude of the S-wave is 0.611-0.946.The S-wave energy is greater than the P-wave energy,and the excitation efficiency varies greatly with different excitation environment;③The peak velocity increases with the equivalent magnitude,and decreases with the epicentral distance.The vertical component of the P-wave peak velocity is the largest among those three components,while the S-wave has the smallest vertical component and similar horizontal components.Hence,our research can provide an important basis for the quantitative judgment of the seismic wave propagation distance using the airgun and the design of the observation system in deep exploration or monitoring with airgun.展开更多
A force-based quadrilateral plate element( 4NQP13) for the analysis of the plate bending problems using large increment method( LIM) was proposed. The LIM, a force-based finite element method( FEM),has been successful...A force-based quadrilateral plate element( 4NQP13) for the analysis of the plate bending problems using large increment method( LIM) was proposed. The LIM, a force-based finite element method( FEM),has been successfully developed for the analysis of truss,beam,frame,and 2D continua problems. In these analyses,LIMcan provide more precise stress results and less computational time consumption compared with displacement-based FEM. The plate element was based on the Mindlin-Reissner plate theory which took into account the transverse shear effects.Numerical examples were presented to study its performance including accuracy and convergence behavior,and the results were compared with the results have been obtained from the displacementbased quadrilateral plate elements and the analytical solutions. The4NQP13 element can analyze the moderately thick plates and the thin plates using LIMand is free from spurious zero energy modes and free from shear locking for thin plate analysis.展开更多
基金Projects(11772188,11132007,11202126)supported by the National Natural Science Foundation of ChinaProject(11ZR1417000)supported by the Natural Science Foundation of Shanghai,China
文摘A forward recursive formulation based on corotational frame is proposed for flexible planar beams with large displacement.The traditional recursive formulation has been successfully used for flexible mutibody dynamics to improve the computational efficiency based on floating frame,in which the assumption of small strain and deflection is adopted.The proposed recursive formulation could be used for large displacement problems based on the corotational frame.It means that the recursive scheme is used not only for adjacent bodies but also for adjacent beam elements.The nodal relative rotation coordinates of the planar beam are used to obtain equations with minimal generalized coordinates in present formulation.The proposed formulation is different from absolute nodal coordinate formulation and the geometrically exact beam formulation in which the absolute coordinates are used.The recursive scheme and minimal set of dynamic equations lead to a high computational efficiency in numerical integration.Numerical examples are carried out to demonstrate the accuracy and validity of this formulation.For all of the examples,the results of the present formulation are in good agreement with results obtained using commercial software and the published results.Moreover,it is shown that the present formulation is more efficient than the formulation in ANSYS based on GEBF.
基金National Key Technology R & D Program,China (No.2006BAJ12B03-2)
文摘A Bernoulli-Euler beam mechanism for static analysis of large displacement,large rotation but small strain planar tapered beam structures is proposed using the Updated Lagrangian formulation and the moving coordinate method.The object beam is the tapered one whose profile is assumed to be varying linearly.From the governing differential equation of lateral deflection including second-order effects by beam-column theory,the geometric nonlinear tangent elemental stiffness matrix is derived.The nonlinear effect of the bending distortions on the axial action is considered to manifest itself as an axial change in length.The aforementioned stiffness matrix is amended,by developing the auxiliary stiffness of bowing effect.The moving coordinate method is employed for obtaining the large displacement total equilibrium equations,and the hinged-hinged moving coordinate system is constructed at the last updated configuration.The multiple load steps Newton-Raphson iteration is adopted for the solution of the nonlinear equations.The validity and efficiency of the proposed method are shown by solving various typical numerical examples.
文摘A study to estimate land surface movement caused by large surface excavations in sedimentary strata is presented.In stratified or jointed strata the stress relief driven movement adjacent to large excavations can be significantly larger than expected.High lateral stresses measured in Australia and other places around the world indicate that the ratio of horizontal to vertical stress has been particularly high at shallow depths.The in situ strata is in compression and during excavation,stress is relieved towards the opening causing strata movement.Large excavations such as,open cut mines or highway cuttings,can initiate an extensive horizontal slide of surface layers towards the excavation.These ground movements can be damaging to surface structures such as water storage dams and large buildings.Based on stress measurements at shallow depths in Australian coal mines the study presented here calculates and models the extent of potential ground movement along the bedding surface adjacent to large excavations and provides a new prediction tool of land movement at the excavation boundary that can benefit the geotechnical practitioners in the mining industry.
文摘The variational functional of the Hellinger-Reissner variational principle for the large displacement problem of a thin shallow shell with an arbitrary shape is first established. Then the functional of the modified principle suitable for the finite element method is derived. In the functional only two independent variables, the deflection w and the stress function F are included. The displacement expressions in the middle surface on the boundary of the shell is also derived by means of the previous two variables.
基金supported by the National Natural Science Foundation of China (Grant No.50278051)Shanghai Pujiang Program(Grant No.07pj14073)
文摘In this paper, a nonlinear mathematical model for analyzing dynamical response to the large deformation of piles with initial displacements is firstly established with the arc-coordinate, and it is a set of nonlinear integral-differential equa- tions, in which, the Winkeler model is used to simulate the resistance of the soil to the pile. Secondly, a set of new auxiliary functions are introduced. The differential-integral equations are transformed into a set of nonlinear differential equations, and the differential quadrature method (DQM) and the finite difference method (FDM) are applied to discretize the set of nonlinear equations in the spatial and time domains, respectively. Then, the Newton-Raphson method is used to solve the set of discretization algebraic equations at each time step. Finally, numerical examples are presented, and the dynamical re- sponses to the deformation of piles, including configuration, bending moment and shear force, are graphically illuminated. In calculation, two types of initial displacements and dynamical loads are applied, and the effects of parameters on the dynamical responses of piles are analyzed in detail.
文摘The fundamental equations for analysis of a straight slender rod which deform at large deflections of order of several times of the dimension of the cross section of the rod are derived by using convecled coordinate system. In accordance with the practice of oil drilling, an effective method of three dimensional static analysis of bottom-hole assemblies is simply described. Errors in reference [8 ] are pointed out.
基金the National Natural Science Foundation of China(41790461)
文摘In order to further deepen the understanding of seismic wave propagation characteristics induced by the large volume airgun source,experimental data from multiple fixed excitation points in Fujian Province were used to obtain the equivalent single excitation high signal-to-noise ratio velocity and displacement records through linear stacking and simulation techniques.Then the peak displacements of different epicentral distances were used to calculate the equivalent magnitude of the airgun source excitation at different fixed excitation points so as to establish the attenuation relationship between equivalent magnitude,epicenter distance and velocity peak.Our results show that:①Within 270 km of epicentral distance,for the large-volume airgun’s single shot,the peak velocity range is about 700-4 nm/s,and the peak displacement range is about 200.0-0.2 nm;②The equivalent magnitude of the P-wave from the airgun source with a total capacity of 8,000 in 3 is 0.181-0.760,and the equivalent magnitude of the S-wave is 0.294-0.832.By contrast,the equivalent magnitude of the P-wave from the airgun source with a total capacity of 12,000 in 3 is 0.533-0.896,and the equivalent magnitude of the S-wave is 0.611-0.946.The S-wave energy is greater than the P-wave energy,and the excitation efficiency varies greatly with different excitation environment;③The peak velocity increases with the equivalent magnitude,and decreases with the epicentral distance.The vertical component of the P-wave peak velocity is the largest among those three components,while the S-wave has the smallest vertical component and similar horizontal components.Hence,our research can provide an important basis for the quantitative judgment of the seismic wave propagation distance using the airgun and the design of the observation system in deep exploration or monitoring with airgun.
基金National Natural Science Foundation of China(No.10872128)
文摘A force-based quadrilateral plate element( 4NQP13) for the analysis of the plate bending problems using large increment method( LIM) was proposed. The LIM, a force-based finite element method( FEM),has been successfully developed for the analysis of truss,beam,frame,and 2D continua problems. In these analyses,LIMcan provide more precise stress results and less computational time consumption compared with displacement-based FEM. The plate element was based on the Mindlin-Reissner plate theory which took into account the transverse shear effects.Numerical examples were presented to study its performance including accuracy and convergence behavior,and the results were compared with the results have been obtained from the displacementbased quadrilateral plate elements and the analytical solutions. The4NQP13 element can analyze the moderately thick plates and the thin plates using LIMand is free from spurious zero energy modes and free from shear locking for thin plate analysis.