The surface fracture toughness is an important mechanical parameter for studying the failure behavior of air plasma sprayed(APS)thermal barrier coatings(TBCs).As APS TBCs are typical multilayer porous ceramic material...The surface fracture toughness is an important mechanical parameter for studying the failure behavior of air plasma sprayed(APS)thermal barrier coatings(TBCs).As APS TBCs are typical multilayer porous ceramic materials,the direct applications of the traditional single edge notched beam(SENB)method that ignores those typical structural characters may cause errors.To measure the surface fracture toughness more accurately,the effects of multilayer and porous characters on the fracture toughness of APS TBCs should be considered.In this paper,a modified single edge V-notched beam(MSEVNB)method with typical structural characters is developed.According to the finite element analysis(FEA),the geometry factor of the multilayer structure is recalculated.Owing to the narrower V-notches,a more accurate critical fracture stress is obtained.Based on the Griffith energy balance,the reduction of the crack surface caused by micro-defects is corrected.The MSEVNB method can measure the surface fracture toughness more accurately than the SENB method.展开更多
This paper extends Le van's work to the case of nonlinear problem and the complicated configuration. The wrinkling stress distribution and the pressure effects are also included in our analysis. Pseudo-beam method is...This paper extends Le van's work to the case of nonlinear problem and the complicated configuration. The wrinkling stress distribution and the pressure effects are also included in our analysis. Pseudo-beam method is presented based on the inflatable beam theory to model the inflatable structures as a set of inflatable beam elements with a prestressed state. In this method, the discretized nonlinear equations are given based upon the virtual work principle with a 3-node Timoshenko's beam model. Finite element simulation is performed by using a 3-node BEAM189 element incorporating ANSYS nonlinear program. The pressure effect is equivalent included in our method by modifying beam element cross-section parameters related to pressure. A benchmark example, the bending case of an inflatable cantilever beam, is performed to verify the accuracy of our proposed method. The comparisons reveal that the numerical results obtained with our method are close to open published analytical and membrane finite element results. The method is then used to evaluate the whole buckling and the loadcarrying characteristics of an inflatable support frame subjected to a compression force. The wrinkling stress and region characteristics are also shown in the end. This method gives better convergence characteristics, and requires much less computation time. It is very effective to deal with the whole load-carrying ability analytical problems for large scale inflatable structures with complex configuration.展开更多
The behavior of beams with variable stiffness subjected to the action of variable loadings (impulse or harmonic) is analyzed in this paper using the successive approximation method. This successive approximation metho...The behavior of beams with variable stiffness subjected to the action of variable loadings (impulse or harmonic) is analyzed in this paper using the successive approximation method. This successive approximation method is a technique for numerical integration of partial differential equations involving both the space and time, with well-known initial conditions on time and boundary conditions on the space. This technique, although having been applied to beams with constant stiffness, is new for the case of beams with variable stiffness, and it aims to use a quadratic parabola (in time) to approximate the solutions of the differential equations of dynamics. The spatial part is studied using the successive approximation method of the partial differential equations obtained, in order to transform them into a system of time-dependent ordinary differential equations. Thus, the integration algorithm using this technique is established and applied to examples of beams with variable stiffness, under variable loading, and with the different cases of supports chosen in the literature. We have thus calculated the cases of beams with constant or variable rigidity with articulated or embedded supports, subjected to the action of an instantaneous impulse and harmonic loads distributed over its entire length. In order to justify the robustness of the successive approximation method considered in this work, an example of an articulated beam with constant stiffness subjected to a distributed harmonic load was calculated analytically, and the results obtained compared to those found numerically for various steps (spatial h and temporal τ ¯ ) of calculus, and the difference between the values obtained by the two methods was small. For example for ( h=1/8 , τ ¯ =1/ 64 ), the difference between these values is 17%.展开更多
As an important model in quantum semiconductor devices, the SchrSdinger-Poisson equations have generated widespread interests in both analysis and numerical simulations in recent years. In this paper, we present Gauss...As an important model in quantum semiconductor devices, the SchrSdinger-Poisson equations have generated widespread interests in both analysis and numerical simulations in recent years. In this paper, we present Gaussian beam methods for the numerical simulation of the one-dimensional Schrodinger-Poisson equations. The Gaussian beam methods for high frequency waves outperform the geometrical optics method in that the former are accurate even around caustics. The purposes of the paper are first to develop the Gaussian beam methods, based on our previous methods for the linear SchrSdinger equation, for the Schrodinger-Poisson equations, and then check their validity for this weakly-nonlinear system.展开更多
A novel Eulerian Gaussian beam method was developed in[8]to compute the Schrödinger equation efficiently in the semiclassical regime.In this paper,we introduce an efficient semi-Eulerian implementation of this me...A novel Eulerian Gaussian beam method was developed in[8]to compute the Schrödinger equation efficiently in the semiclassical regime.In this paper,we introduce an efficient semi-Eulerian implementation of this method.The new algorithm inherits the essence of the Eulerian Gaussian beam method where the Hessian is computed through the derivatives of the complexified level set functions instead of solving the dynamic ray tracing equation.The difference lies in that,we solve the ray tracing equations to determine the centers of the beams and then compute quantities of interests only around these centers.This yields effectively a local level set implementation,and the beam summation can be carried out on the initial physical space instead of the phase plane.As a consequence,it reduces the computational cost and also avoids the delicate issue of beam summation around the caustics in the Eulerian Gaussian beam method.Moreover,the semi-Eulerian Gaussian beam method can be easily generalized to higher order Gaussian beam methods,which is the topic of the second part of this paper.Several numerical examples are provided to verify the accuracy and efficiency of both the first order and higher order semi-Eulerian methods.展开更多
A modified alternating direction implicit algorithm is proposed to solve the full-vectorial finite-difference beam propagation method formulation based on H fields. The cross-coupling terms are neglected in the first ...A modified alternating direction implicit algorithm is proposed to solve the full-vectorial finite-difference beam propagation method formulation based on H fields. The cross-coupling terms are neglected in the first sub-step, but evaluated and doubly used in the second sub-step. The order of two sub-steps is reversed for each transverse magnetic field component so that the cross-coupling terms are always expressed in implicit form, thus the calculation is very efficient and stable. Moreover, an improved six-point finite-difference scheme with high accuracy independent of specific structures of waveguide is also constructed to approximate the cross-coupling terms along the transverse directions. The imaginary-distance procedure is used to assess the validity and utility of the present method. The field patterns and the normalized propagation constants of the fundamental mode for a buried rectangular waveguide and a rib waveguide are presented. Solutions are in excellent agreement with the benchmark results from the modal transverse resonance method.展开更多
Chaotic vibrations of flexible non-linear Euler-Bernoulli beams subjected to harmonic load and with various boundary conditions(symmetric and non-symmetric)are studied in this work.Reliability of the obtained result...Chaotic vibrations of flexible non-linear Euler-Bernoulli beams subjected to harmonic load and with various boundary conditions(symmetric and non-symmetric)are studied in this work.Reliability of the obtained results is verified by the finite difference method(FDM)and the finite element method(FEM)with the Bubnov-Galerkin approximation for various boundary conditions and various dynamic regimes(regular and non-regular).The influence of boundary conditions on the Euler-Bernoulli beams dynamics is studied mainly,dynamic behavior vs.control parameters { ωp,q0 } is reported,and scenarios of the system transition into chaos are illustrated.展开更多
The mathematical modeling of a rotating tapered Timoshenko beam with preset and pre-twist angles is constructed. The partial differential equations governing the six degrees, i.e., three displacements in the axial, fl...The mathematical modeling of a rotating tapered Timoshenko beam with preset and pre-twist angles is constructed. The partial differential equations governing the six degrees, i.e., three displacements in the axial, flapwise, and edgewise directions and three cross-sectional angles of torsion, flapwise bending, and edgewise bending, are obtained by the Euler angle descriptions. The power series method is then used to inves- tigate the natural frequencies and the corresponding complex mode functions. It is found that all the natural frequencies are increased by the centrifugal stiffening except the twist frequency, which is slightly decreased. The tapering ratio increases the first transverse, torsional, and axial frequencies, while decreases the second transverse frequency. Because of the pre-twist, all the directions are gyroscopically coupled with the phase differences among the six degrees.展开更多
With the third innovation in science and technology worldwide, China has also experienced thismarvelous progress. Concerning the longwall mining in China, the "masonry beam theory" (MBT) wasfirst proposed in the 1...With the third innovation in science and technology worldwide, China has also experienced thismarvelous progress. Concerning the longwall mining in China, the "masonry beam theory" (MBT) wasfirst proposed in the 1960s, illustrating that the transmission and equilibrium method of overburdenpressure using reserved coal pillar in mined-out areas can be realized. This forms the so-called "121mining method", which lays a solid foundation for development of mining science and technology inChina. The "transfer rock beam theory" (TRBT) proposed in the 1980s gives a further understanding forthe transmission path of stope overburden pressure and pressure distribution in high-stress areas. In thisregard, the advanced 121 mining method was proposed with smaller coal pillar for excavation design,making significant contributions to improvement of the coal recovery rate in that era. In the 21st century,the traditional mining technologies faced great challenges and, under the theoretical developmentspioneered by Profs. Minggao Qian and Zhenqi Song, the "cutting cantilever beam theory" (CCBT) wasproposed in 2008. After that the 110 mining method is formulated subsequently, namely one stope face,after the first mining cycle, needs one advanced gateway excavation, while the other one is automaticallyformed during the last mining cycle without coal pillars left in the mining area. This method can beimplemented using the CCBT by incorporating the key technologies, including the directional presplittingroof cutting, constant resistance and large deformation (CRLD) bolt/anchor supporting systemwith negative Poisson's ratio (NPR) effect material, and remote real-time monitoring technology. TheCCBT and 110 mining method will provide the theoretical and technical basis for the development ofmining industry in China.展开更多
Since the 1960 s, mining science and technology in China has experienced two technical innovations, i.e.the ‘‘Masonry Beam Theory(MBT)" and ‘‘Transfer Rock Beam Theory(TRBT)". Based on those theories, th...Since the 1960 s, mining science and technology in China has experienced two technical innovations, i.e.the ‘‘Masonry Beam Theory(MBT)" and ‘‘Transfer Rock Beam Theory(TRBT)". Based on those theories, the conventional mining method(being called the 121 mining method) was established, consisting of excavating two tunnels with a pillar left for mining a working panel. However, with increasing mining depth,engineering geological disasters in the underground caverns have been frequently encountered. In addition, the use of the coal-pillar mining results in a large amount of coal resources unexploited. In order to address the problems above, the ‘‘Roof Cut Short-Arm Beam Theory(RCSBT), being called the 110 mining method)" was proposed by He Manchao in 2008. The 110 mining method features the mining of one coal seam panel, excavating necessarily only one roadway tunnel and leaving no pillars. Realization of the 110 mining method includes the following steps:(1) directional pre-splitting roof cutting,(2) supporting the roof by using high Constant Resistance Large Deformation bolt/cable(CRLD), and(3) blocking gangue by hydraulic props. This paper presents an overview of the principles, techniques and application of the 110 mining method. Special emphasis is placed on the numerical simulation of the geostress distribution found in the mining panel using the 110 method compared to that of the 121 method. In addition, the stress distribution on the ‘‘short beam" left by the roof cutting when performing the 110 method was also investigated using both numerical simulation and theoretical formulation.展开更多
Large deformation of a cantilever axially functionally graded (AFG) beam subject to a tip load is analytically studied using the homotopy analysis method (HAM). It is assumed that its Young’s modulus varies along the...Large deformation of a cantilever axially functionally graded (AFG) beam subject to a tip load is analytically studied using the homotopy analysis method (HAM). It is assumed that its Young’s modulus varies along the longitudinal direction according to a power law. Taking the solution of the corresponding homogeneous beam as the initial guess and obtaining a convergence region by adjusting an auxiliary parameter, the analytical expressions for large deformation of the AFG beam are provided. Results obtained by the HAM are compared with those obtained by the finite element method and those in the previous works to verify its validity. Good agreement is observed. A detailed parametric study is carried out. The results show that the axial material variation can greatly change the deformed configuration, which provides an approach to control and manage the deformation of beams. By tailoring the axial material distribution, a desired deformed configuration can be obtained for a specific load. The analytical solution presented herein can be a helpful tool for this procedure.展开更多
Based on the method of reverberation ray matrix(MRRM), a reverberation matrix for planar framed structures composed of anisotropic Timoshenko(T) beam members containing completely hinged joints is developed for st...Based on the method of reverberation ray matrix(MRRM), a reverberation matrix for planar framed structures composed of anisotropic Timoshenko(T) beam members containing completely hinged joints is developed for static analysis of such structures.In the MRRM for dynamic analysis, amplitudes of arriving and departing waves for joints are chosen as unknown quantities. However, for the present case of static analysis, displacements and rotational angles at the ends of each beam member are directly considered as unknown quantities. The expressions for stiffness matrices for anisotropic beam members are developed. A corresponding reverberation matrix is derived analytically for exact and unified determination on the displacements and internal forces at both ends of each member and arbitrary cross sectional locations in the structure. Numerical examples are given and compared with the finite element method(FEM) results to validate the present model. The characteristic parameter analysis is performed to demonstrate accuracy of the present model with the T beam theory in contrast with errors in the usual model based on the Euler-Bernoulli(EB) beam theory. The resulting reverberation matrix can be used for exact calculation of anisotropic framed structures as well as for parameter analysis of geometrical and material properties of the framed structures.展开更多
The analytical solution of a viscoelastic continuous beam whose damping characteristics are described in terms of a fractional derivative of arbitrary order was derived by means of the AdoInian decomposition method. T...The analytical solution of a viscoelastic continuous beam whose damping characteristics are described in terms of a fractional derivative of arbitrary order was derived by means of the AdoInian decomposition method. The solution contains arbitrary initial conditions and zero input. For specific analysis, the initial conditions were assumed homogeneous, and the input force was treated as a special process with a particular beam. Two simple cases, step and impulse function responses, were considered respectively. Subsequently, some figures were plotted to show the displacement of the beam under different sets of parameters including different orders of the fractional derivatives.展开更多
The bearing beams and the supporting beams under low velocity impact may be in four different strain stages of deformation depending on the impact intensity and beam structure strength.Based on the different judging c...The bearing beams and the supporting beams under low velocity impact may be in four different strain stages of deformation depending on the impact intensity and beam structure strength.Based on the different judging conditions of deformation stages,the corresponding calculation models are proposed,the calculation formulae for the determination of the impact force and the beam's lateral displacement are obtained.Calculation shows that the beam's total deflection is small when the flexibility of the supporting component is high and the effect of diminishing deflection disappears almost when the stiffness of the supporting component is high.展开更多
The vibration suppression of the finite plate with square steel beams is studied using traveling wave method. The finite plate with square beams is modeled as the coupling systems between the plate flexural motion and...The vibration suppression of the finite plate with square steel beams is studied using traveling wave method. The finite plate with square beams is modeled as the coupling systems between the plate flexural motion and the flexural and torsional motions for the square beams. The vibration response at any position of the coupling structure can be obtained by wave method. Numerical results show that comparing to finite element method (FEM), not only the low frequency but also the medium-high frequency vibration response of the finite plate with square beam can be effectively calculated by wave method. The suppression effect can be increased as the square beam is located at one-third of the length of plate or increasing the height of the beam. The study provides reference for arranged square beams applying to vibration suppression of ship and train structures.展开更多
The characteristics of transverse free vibration of a tapered Timoshenko beam under an axially conservative compression resting on visco-Pasternak foundations are investigated by the interpolating matrix method. The r...The characteristics of transverse free vibration of a tapered Timoshenko beam under an axially conservative compression resting on visco-Pasternak foundations are investigated by the interpolating matrix method. The research is executed in view of a three-parameter foundation which includes the eff ects of the Winkler coeffi cient, Pasternak coeffi cient and damping coeffi cient of the elastic medium. The governing equations of free vibration of a non-prismatic Timoshenko beam under an axially conservative force resting on visco-Pasternak foundations are transformed into ordinary diff erential equations with variable coeffi cients in light of the bending rotation angle and transverse displacement. All the natural frequencies orders together with the corresponding mode shapes of the beam are calculated at the same time, and a good convergence and accuracy of the proposed method is verifi ed through two numerical examples. The infl uences of foundation mechanical characteristics together with rotary inertia and shear deformation on natural frequencies of the beam with diff erent taper ratios are analyzed. A comprehensive parametric numerical study is carried out emphasizing the primary parameters that describe the dynamic property of the beam.展开更多
Based on the multi-rigid body discretization model, namely, finite segment model,a chain multi-rigid-body-hinge-spring system model of a beam was presented, then a nonlinear parametrically exacted vibration equation o...Based on the multi-rigid body discretization model, namely, finite segment model,a chain multi-rigid-body-hinge-spring system model of a beam was presented, then a nonlinear parametrically exacted vibration equation of multi-degrees of freedom system was established using the coordination transformation method, and its resonance fields were derived by the restriction parameter method, that is, the dynamical buckling analysis of the beam. Because the deformation of a beam is not restricted by the discrete model and dynamic equation, the post buckling analysis can be done in above math model. The numerical solutions of a few examples were obtained by direct integrated method, which shows that the mechanical and math model gotten is correct.展开更多
The dynamic characteristics of a beam-cable coupled system are investigated using an improved Chebyshev spectral element method in order to observe the effects of adding cables on the beam. The system is modeled as a ...The dynamic characteristics of a beam-cable coupled system are investigated using an improved Chebyshev spectral element method in order to observe the effects of adding cables on the beam. The system is modeled as a double Timoshenko beam system interconnected by discrete springs. Utilizing Chebyshev series expansion and meshing the system according to the locations of its connections, numerical results of the natural frequencies and mode shapes are obtained using only a few elements, and the results are validated by comparing them with the results of a finite-element method. Then the effects of the cable parameters and layout of connections on the natural frequencies and mode shapes of a fixed-pinned beam are studied. The results show that the modes of a beam-cable coupled system can be classified into two types, beam mode and cable mode, according to the dominant deformation. To avoid undesirable vibrations of the cable, its parameters should be controlled in a reasonable range, or the layout of the connections should be optimized.展开更多
基金Project supported by the National Natural Science Foundation of China(Nos.12172048 and 12027901)the National Science and Technology Major Project of China(Nos.2019-Ⅶ-0007-0147 and 2017-Ⅵ-0020-0093)。
文摘The surface fracture toughness is an important mechanical parameter for studying the failure behavior of air plasma sprayed(APS)thermal barrier coatings(TBCs).As APS TBCs are typical multilayer porous ceramic materials,the direct applications of the traditional single edge notched beam(SENB)method that ignores those typical structural characters may cause errors.To measure the surface fracture toughness more accurately,the effects of multilayer and porous characters on the fracture toughness of APS TBCs should be considered.In this paper,a modified single edge V-notched beam(MSEVNB)method with typical structural characters is developed.According to the finite element analysis(FEA),the geometry factor of the multilayer structure is recalculated.Owing to the narrower V-notches,a more accurate critical fracture stress is obtained.Based on the Griffith energy balance,the reduction of the crack surface caused by micro-defects is corrected.The MSEVNB method can measure the surface fracture toughness more accurately than the SENB method.
基金supported by the Specialized Fund for the Doctoral Program of Higher Education of China (200802131046)China Postdoctoral Science Foundation Funded Major Project (200801290)+1 种基金Development Program of Outstanding Young Teachers in Harbin Institute of Technology (HITQNJS.2008.004)Specialized Fund for Innovation Talents of Science and Technology in Harbin (2008RFQXG057).
文摘This paper extends Le van's work to the case of nonlinear problem and the complicated configuration. The wrinkling stress distribution and the pressure effects are also included in our analysis. Pseudo-beam method is presented based on the inflatable beam theory to model the inflatable structures as a set of inflatable beam elements with a prestressed state. In this method, the discretized nonlinear equations are given based upon the virtual work principle with a 3-node Timoshenko's beam model. Finite element simulation is performed by using a 3-node BEAM189 element incorporating ANSYS nonlinear program. The pressure effect is equivalent included in our method by modifying beam element cross-section parameters related to pressure. A benchmark example, the bending case of an inflatable cantilever beam, is performed to verify the accuracy of our proposed method. The comparisons reveal that the numerical results obtained with our method are close to open published analytical and membrane finite element results. The method is then used to evaluate the whole buckling and the loadcarrying characteristics of an inflatable support frame subjected to a compression force. The wrinkling stress and region characteristics are also shown in the end. This method gives better convergence characteristics, and requires much less computation time. It is very effective to deal with the whole load-carrying ability analytical problems for large scale inflatable structures with complex configuration.
文摘The behavior of beams with variable stiffness subjected to the action of variable loadings (impulse or harmonic) is analyzed in this paper using the successive approximation method. This successive approximation method is a technique for numerical integration of partial differential equations involving both the space and time, with well-known initial conditions on time and boundary conditions on the space. This technique, although having been applied to beams with constant stiffness, is new for the case of beams with variable stiffness, and it aims to use a quadratic parabola (in time) to approximate the solutions of the differential equations of dynamics. The spatial part is studied using the successive approximation method of the partial differential equations obtained, in order to transform them into a system of time-dependent ordinary differential equations. Thus, the integration algorithm using this technique is established and applied to examples of beams with variable stiffness, under variable loading, and with the different cases of supports chosen in the literature. We have thus calculated the cases of beams with constant or variable rigidity with articulated or embedded supports, subjected to the action of an instantaneous impulse and harmonic loads distributed over its entire length. In order to justify the robustness of the successive approximation method considered in this work, an example of an articulated beam with constant stiffness subjected to a distributed harmonic load was calculated analytically, and the results obtained compared to those found numerically for various steps (spatial h and temporal τ ¯ ) of calculus, and the difference between the values obtained by the two methods was small. For example for ( h=1/8 , τ ¯ =1/ 64 ), the difference between these values is 17%.
基金Research Program of China under the grant 2005CB321701supported by a Van Vleck Distinguished Research Prize from University of Wisconsin-Madison
文摘As an important model in quantum semiconductor devices, the SchrSdinger-Poisson equations have generated widespread interests in both analysis and numerical simulations in recent years. In this paper, we present Gaussian beam methods for the numerical simulation of the one-dimensional Schrodinger-Poisson equations. The Gaussian beam methods for high frequency waves outperform the geometrical optics method in that the former are accurate even around caustics. The purposes of the paper are first to develop the Gaussian beam methods, based on our previous methods for the linear SchrSdinger equation, for the Schrodinger-Poisson equations, and then check their validity for this weakly-nonlinear system.
基金This work was partially supported by NSF grant No.DMS-0608720,NSF FRG grant DMS-0757285,NSAF Projects 10676017,NSFC Projects 10971115the National Basic Research Program of China under the grant 2005CB321701.SJin was also supported by a Van Vleck Distinguished Research Prize from University of Wisconsin-Madison.
文摘A novel Eulerian Gaussian beam method was developed in[8]to compute the Schrödinger equation efficiently in the semiclassical regime.In this paper,we introduce an efficient semi-Eulerian implementation of this method.The new algorithm inherits the essence of the Eulerian Gaussian beam method where the Hessian is computed through the derivatives of the complexified level set functions instead of solving the dynamic ray tracing equation.The difference lies in that,we solve the ray tracing equations to determine the centers of the beams and then compute quantities of interests only around these centers.This yields effectively a local level set implementation,and the beam summation can be carried out on the initial physical space instead of the phase plane.As a consequence,it reduces the computational cost and also avoids the delicate issue of beam summation around the caustics in the Eulerian Gaussian beam method.Moreover,the semi-Eulerian Gaussian beam method can be easily generalized to higher order Gaussian beam methods,which is the topic of the second part of this paper.Several numerical examples are provided to verify the accuracy and efficiency of both the first order and higher order semi-Eulerian methods.
文摘A modified alternating direction implicit algorithm is proposed to solve the full-vectorial finite-difference beam propagation method formulation based on H fields. The cross-coupling terms are neglected in the first sub-step, but evaluated and doubly used in the second sub-step. The order of two sub-steps is reversed for each transverse magnetic field component so that the cross-coupling terms are always expressed in implicit form, thus the calculation is very efficient and stable. Moreover, an improved six-point finite-difference scheme with high accuracy independent of specific structures of waveguide is also constructed to approximate the cross-coupling terms along the transverse directions. The imaginary-distance procedure is used to assess the validity and utility of the present method. The field patterns and the normalized propagation constants of the fundamental mode for a buried rectangular waveguide and a rib waveguide are presented. Solutions are in excellent agreement with the benchmark results from the modal transverse resonance method.
文摘Chaotic vibrations of flexible non-linear Euler-Bernoulli beams subjected to harmonic load and with various boundary conditions(symmetric and non-symmetric)are studied in this work.Reliability of the obtained results is verified by the finite difference method(FDM)and the finite element method(FEM)with the Bubnov-Galerkin approximation for various boundary conditions and various dynamic regimes(regular and non-regular).The influence of boundary conditions on the Euler-Bernoulli beams dynamics is studied mainly,dynamic behavior vs.control parameters { ωp,q0 } is reported,and scenarios of the system transition into chaos are illustrated.
基金Project supported by the National Natural Science Foundation of China(Nos.11672007,11402028,11322214,and 11290152)the Beijing Natural Science Foundation(No.3172003)the Key Laboratory of Vibration and Control of Aero-Propulsion System Ministry of Education,Northeastern University(No.VCAME201601)
文摘The mathematical modeling of a rotating tapered Timoshenko beam with preset and pre-twist angles is constructed. The partial differential equations governing the six degrees, i.e., three displacements in the axial, flapwise, and edgewise directions and three cross-sectional angles of torsion, flapwise bending, and edgewise bending, are obtained by the Euler angle descriptions. The power series method is then used to inves- tigate the natural frequencies and the corresponding complex mode functions. It is found that all the natural frequencies are increased by the centrifugal stiffening except the twist frequency, which is slightly decreased. The tapering ratio increases the first transverse, torsional, and axial frequencies, while decreases the second transverse frequency. Because of the pre-twist, all the directions are gyroscopically coupled with the phase differences among the six degrees.
基金supported by the National Natural Science Foundation of China (No. 51404278)the State Key Program of National Natural Science Foundation of China (No. 51134005)
文摘With the third innovation in science and technology worldwide, China has also experienced thismarvelous progress. Concerning the longwall mining in China, the "masonry beam theory" (MBT) wasfirst proposed in the 1960s, illustrating that the transmission and equilibrium method of overburdenpressure using reserved coal pillar in mined-out areas can be realized. This forms the so-called "121mining method", which lays a solid foundation for development of mining science and technology inChina. The "transfer rock beam theory" (TRBT) proposed in the 1980s gives a further understanding forthe transmission path of stope overburden pressure and pressure distribution in high-stress areas. In thisregard, the advanced 121 mining method was proposed with smaller coal pillar for excavation design,making significant contributions to improvement of the coal recovery rate in that era. In the 21st century,the traditional mining technologies faced great challenges and, under the theoretical developmentspioneered by Profs. Minggao Qian and Zhenqi Song, the "cutting cantilever beam theory" (CCBT) wasproposed in 2008. After that the 110 mining method is formulated subsequently, namely one stope face,after the first mining cycle, needs one advanced gateway excavation, while the other one is automaticallyformed during the last mining cycle without coal pillars left in the mining area. This method can beimplemented using the CCBT by incorporating the key technologies, including the directional presplittingroof cutting, constant resistance and large deformation (CRLD) bolt/anchor supporting systemwith negative Poisson's ratio (NPR) effect material, and remote real-time monitoring technology. TheCCBT and 110 mining method will provide the theoretical and technical basis for the development ofmining industry in China.
文摘Since the 1960 s, mining science and technology in China has experienced two technical innovations, i.e.the ‘‘Masonry Beam Theory(MBT)" and ‘‘Transfer Rock Beam Theory(TRBT)". Based on those theories, the conventional mining method(being called the 121 mining method) was established, consisting of excavating two tunnels with a pillar left for mining a working panel. However, with increasing mining depth,engineering geological disasters in the underground caverns have been frequently encountered. In addition, the use of the coal-pillar mining results in a large amount of coal resources unexploited. In order to address the problems above, the ‘‘Roof Cut Short-Arm Beam Theory(RCSBT), being called the 110 mining method)" was proposed by He Manchao in 2008. The 110 mining method features the mining of one coal seam panel, excavating necessarily only one roadway tunnel and leaving no pillars. Realization of the 110 mining method includes the following steps:(1) directional pre-splitting roof cutting,(2) supporting the roof by using high Constant Resistance Large Deformation bolt/cable(CRLD), and(3) blocking gangue by hydraulic props. This paper presents an overview of the principles, techniques and application of the 110 mining method. Special emphasis is placed on the numerical simulation of the geostress distribution found in the mining panel using the 110 method compared to that of the 121 method. In addition, the stress distribution on the ‘‘short beam" left by the roof cutting when performing the 110 method was also investigated using both numerical simulation and theoretical formulation.
基金Project supported by the China Postdoctoral Science Foundation(No.2018M630167)
文摘Large deformation of a cantilever axially functionally graded (AFG) beam subject to a tip load is analytically studied using the homotopy analysis method (HAM). It is assumed that its Young’s modulus varies along the longitudinal direction according to a power law. Taking the solution of the corresponding homogeneous beam as the initial guess and obtaining a convergence region by adjusting an auxiliary parameter, the analytical expressions for large deformation of the AFG beam are provided. Results obtained by the HAM are compared with those obtained by the finite element method and those in the previous works to verify its validity. Good agreement is observed. A detailed parametric study is carried out. The results show that the axial material variation can greatly change the deformed configuration, which provides an approach to control and manage the deformation of beams. By tailoring the axial material distribution, a desired deformed configuration can be obtained for a specific load. The analytical solution presented herein can be a helpful tool for this procedure.
基金Project supported by the Program for New Century Excellent Talents in Universities(NCET)by the Ministry of Education of China(No.NCET-04-0373)
文摘Based on the method of reverberation ray matrix(MRRM), a reverberation matrix for planar framed structures composed of anisotropic Timoshenko(T) beam members containing completely hinged joints is developed for static analysis of such structures.In the MRRM for dynamic analysis, amplitudes of arriving and departing waves for joints are chosen as unknown quantities. However, for the present case of static analysis, displacements and rotational angles at the ends of each beam member are directly considered as unknown quantities. The expressions for stiffness matrices for anisotropic beam members are developed. A corresponding reverberation matrix is derived analytically for exact and unified determination on the displacements and internal forces at both ends of each member and arbitrary cross sectional locations in the structure. Numerical examples are given and compared with the finite element method(FEM) results to validate the present model. The characteristic parameter analysis is performed to demonstrate accuracy of the present model with the T beam theory in contrast with errors in the usual model based on the Euler-Bernoulli(EB) beam theory. The resulting reverberation matrix can be used for exact calculation of anisotropic framed structures as well as for parameter analysis of geometrical and material properties of the framed structures.
基金Project supported by the National Natural Science Foundation of China(Nos.10547124 and 10475055)
文摘The analytical solution of a viscoelastic continuous beam whose damping characteristics are described in terms of a fractional derivative of arbitrary order was derived by means of the AdoInian decomposition method. The solution contains arbitrary initial conditions and zero input. For specific analysis, the initial conditions were assumed homogeneous, and the input force was treated as a special process with a particular beam. Two simple cases, step and impulse function responses, were considered respectively. Subsequently, some figures were plotted to show the displacement of the beam under different sets of parameters including different orders of the fractional derivatives.
基金the National Science Foundation of China(50578158)
文摘The bearing beams and the supporting beams under low velocity impact may be in four different strain stages of deformation depending on the impact intensity and beam structure strength.Based on the different judging conditions of deformation stages,the corresponding calculation models are proposed,the calculation formulae for the determination of the impact force and the beam's lateral displacement are obtained.Calculation shows that the beam's total deflection is small when the flexibility of the supporting component is high and the effect of diminishing deflection disappears almost when the stiffness of the supporting component is high.
基金National Natural Science Foundation of China ( No. 10972065) Natural Science Foundation of Heilongjiang Province of China( No. ZD200905)
文摘The vibration suppression of the finite plate with square steel beams is studied using traveling wave method. The finite plate with square beams is modeled as the coupling systems between the plate flexural motion and the flexural and torsional motions for the square beams. The vibration response at any position of the coupling structure can be obtained by wave method. Numerical results show that comparing to finite element method (FEM), not only the low frequency but also the medium-high frequency vibration response of the finite plate with square beam can be effectively calculated by wave method. The suppression effect can be increased as the square beam is located at one-third of the length of plate or increasing the height of the beam. The study provides reference for arranged square beams applying to vibration suppression of ship and train structures.
基金AHKJT of China under Grant Nos.1708085QE121 and 1808085ME147AHEDU of China under Grant No.TSKJ2017B13
文摘The characteristics of transverse free vibration of a tapered Timoshenko beam under an axially conservative compression resting on visco-Pasternak foundations are investigated by the interpolating matrix method. The research is executed in view of a three-parameter foundation which includes the eff ects of the Winkler coeffi cient, Pasternak coeffi cient and damping coeffi cient of the elastic medium. The governing equations of free vibration of a non-prismatic Timoshenko beam under an axially conservative force resting on visco-Pasternak foundations are transformed into ordinary diff erential equations with variable coeffi cients in light of the bending rotation angle and transverse displacement. All the natural frequencies orders together with the corresponding mode shapes of the beam are calculated at the same time, and a good convergence and accuracy of the proposed method is verifi ed through two numerical examples. The infl uences of foundation mechanical characteristics together with rotary inertia and shear deformation on natural frequencies of the beam with diff erent taper ratios are analyzed. A comprehensive parametric numerical study is carried out emphasizing the primary parameters that describe the dynamic property of the beam.
文摘Based on the multi-rigid body discretization model, namely, finite segment model,a chain multi-rigid-body-hinge-spring system model of a beam was presented, then a nonlinear parametrically exacted vibration equation of multi-degrees of freedom system was established using the coordination transformation method, and its resonance fields were derived by the restriction parameter method, that is, the dynamical buckling analysis of the beam. Because the deformation of a beam is not restricted by the discrete model and dynamic equation, the post buckling analysis can be done in above math model. The numerical solutions of a few examples were obtained by direct integrated method, which shows that the mechanical and math model gotten is correct.
基金supported by the National Basic Research Program of China (Grant 2013CB733004)
文摘The dynamic characteristics of a beam-cable coupled system are investigated using an improved Chebyshev spectral element method in order to observe the effects of adding cables on the beam. The system is modeled as a double Timoshenko beam system interconnected by discrete springs. Utilizing Chebyshev series expansion and meshing the system according to the locations of its connections, numerical results of the natural frequencies and mode shapes are obtained using only a few elements, and the results are validated by comparing them with the results of a finite-element method. Then the effects of the cable parameters and layout of connections on the natural frequencies and mode shapes of a fixed-pinned beam are studied. The results show that the modes of a beam-cable coupled system can be classified into two types, beam mode and cable mode, according to the dominant deformation. To avoid undesirable vibrations of the cable, its parameters should be controlled in a reasonable range, or the layout of the connections should be optimized.