首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 140 毫秒
1.
由充液弯管三维振动模型切入,应用动刚度法构建了弯管及直管单元的振动求解方法,进而用于组装求解充液管系的振动,可同时适用于含弯管单元的连续模型或只含直管单元的离散模型;通过算例对比,证明动刚度法比传递矩阵法和有限元法在计算效率和精度上有所提升;与充液L型管道振动实验测得的加速度频响曲线对比,验证了本文对于管道组装的计算方法的有效性,此外还分析了连续模型和离散模型的区别及适用范围。  相似文献   

2.
充液管道系统的模态分析   总被引:25,自引:0,他引:25  
考虑管道与液体之间的泊松耦合与连接耦合,推导了低频情况下的充液直管轴向、横向振动传递矩阵与弯管单元的传递矩阵,引入了弯管弯曲因子,对-“L”形充液管系进行了研究,使用ANSYS有限元分析软件与传递矩阵法计算了两种边界条件下充液、充所管系的模态频率与模态振型,结果表明,气体与管道之间的耦合作用对管系模态的影响可忽略;而对充液管系而言,因液耦合作用不能忽略。  相似文献   

3.
分别基于经典薄板理论和一阶剪切理论研究了沿半径方向变厚度的圆板及圆环板的横向自由振动,将结构离散为若干个等厚度同心圆环单元,在得出圆环单元的精确解后,通过动刚度法组装单元。应用该方法将变厚度圆板退化至等厚度板,与解析解对比验证了计算方法的正确性;用于计算线性或非线性变厚度板,也能与有限元三维解吻合。计算结果表明:基于一阶剪切理论和薄板理论的动刚度法计算等厚度薄板的振动均能取得与解析解完全吻合的数值解;而计算变厚度薄板时则采用基于一阶剪切理论的动刚度法更准确;与有限元法相比,本文采用的动刚度法划分单元少,具有较高的计算效率,尤其在工程中的大型板结构振动方面有较好的应用前景。  相似文献   

4.
管道联接件参数识别的行波法   总被引:2,自引:0,他引:2  
基于结构振动的行波理论,提出了一种管道联接件物理参数识别的新方法.以未充液直管在平面内的弯曲振动为例,建立了入射波和出射波在联接件处的散射关系,并结合波模式幅值系数的识别,得到一组包含待识别参数的非线性方程组.在每一个频率下求解该方程组得到联接件的动刚度,最后采用最小二乘法一次识别得到联接件的物理参数.数值仿真结果表明,行波法在管道联接件参数辨识中,识别精度高,而且采用矩阵条件数小的传感器阵列,可以提高参数辨识的精度.  相似文献   

5.
求解不连续中厚板自由振动的微分容积单元法   总被引:2,自引:0,他引:2  
基于区域叠加原理和微分容积法,发展了一种新型的数值方法——微分容积单元法,用以分析具有不连续几何特征的中厚板的自由振动。根据板的不连续情况将其划分为若干单元,在每个单元内用微分容积法将控制微分方程离散成为一组线性代数方程.在相邻的单元连接处应用位移连续条件和平衡条件,引入边界约束条件后得到一套关于各配点位移的齐次线性代数方程,由此可导出求解系统固有频率的特征方程。本文用子空间迭代法求解特征方程,并以开孔板、混合边界条件板和突变厚度板为例研究了方法的收敛性和计算精度。  相似文献   

6.
本文提出一种用于含分层的双层梁线弹性断裂分析的有限元方法.将上下子梁均模拟为多个子层,采用只有平动位移自由度的新型梁单元,假设单元内的位移沿纵向和横向均线性变化,推导了该单元的单元刚度矩阵.将开裂部分和未开裂部分的子梁进行单元刚度矩阵组装,施加相应的等效结点力,得到整体平衡方程,并结合边界条件进行求解.为验证该方法的有效性和精度,开展非对称双悬臂梁(Asymmetric Double Cantilever Beam, ADCB)和单臂弯曲梁(Single Leg Bending, SLB)试样的断裂分析,利用虚拟裂纹闭合技术(Virtual Crack Closure Technique, VCCT)得到了试样的能量释放率及其分量,并将求得的结果与解析解和二维有限元解进行对比.计算结果表明,相对于传统双层模拟方法,该多层模拟方法能够精确、高效地计算各类梁试样的能量释放率及其分量,并且无需引入界面连续条件.  相似文献   

7.
基于"三心重合"的设计思想和飞行器的姿态需求,建立了含立方项非线性刚度的惯组小系统动力学模型。提出了惯组小系统存在系统动刚度和减振器动刚度的"双层级"概念。利用龙格-库塔法求解渐软非线性系统对正弦扫频激励的响应,得出减振器动刚度存在对激励幅值和激励频率的敏感区域,采用动刚度曲线表征了减振器的非线性软化特性。利用虚弧长延拓法计算惯组小系统的非线性频响函数和传递特性,预示了高量级振动下减振器的动力失稳现象。通过惯组小系统传递特性试验,验证了减振器具有渐软刚度的非线性特性。扫频法计算结果与正弦扫频试验结果的吻合度达到96.5%,检验了构建模型的正确性。所建立的模型工程应用简便,对惯组小系统非线性特性的预示具有较高的精度,可供飞行器姿控系统设计时参考。  相似文献   

8.
弱不连续问题(如含夹杂问题)是固体力学计算中的一类重要问题。高阶有限元方法由于其具有更好的逼近效果,是确保数值解在界面保持较高精度的计算方法之一。但与线性元相比,高阶单元需要更多的计算机存储单元,具有更高的计算复杂性。本文利用两水平算法的思想,将高阶有限元离散系统化归于线性元离散系统的求解,为弱不连续问题高阶有限元离散系统设计了一种新的基于几何与分析信息的代数多重网格(GAMG)法,并应用于圆形求解域含单夹杂问题的高阶有限元离散系统的求解。数值试验结果表明,相比于常用GAMG法,新方法的迭代次数基本不依赖于问题规模、单元阶次以及杨氏模量的间断性,CPU计算时间得到明显改善,具有更好的计算效率和鲁棒性,可大大提高弱不连续问题有限元分析的整体效率。  相似文献   

9.
根据独立连续拓扑变量概念,建立了桁架和平面膜结构拓扑优化的有无复合体模型,从而不引入过滤函数实现拓扑变量在连续型和离散型之间的转换,推导了有无复合体杆单元的面积与膜单元的厚度与重量、单元刚度阵都是“有单元”和“无单元”相应量的线性组合,进而把这一线性关系延拓到许用应力,借助于有无复合体建立了在应力约束下骨架和连续体结构拓扑优化的统一模型,同时提出了求解这一模型的有效算法,获得了令人满意的计算结果。  相似文献   

10.
应用特征线差分法求解耦合瞬变问题时存在难以避免的多波插值问题,会引入较大的插值误差。为了解决此问题,本文提出了一种求解管道耦合水力瞬变模型的Godunov计算格式。首先基于有限体积法对模型进行数值离散,然后采用时空均为二阶精度的三步MUSCL-Hancock方法计算单元界面上的数值通量,同时引入斜率限制器函数来抑制虚假的数值振荡。在计算边界单元时,采用Rankine-Hugoniot条件与边界条件相结合的方法建立边界方程,有效降低了计算的复杂程度。实验与仿真对比表明:本文的计算结果与实验结果吻合较好,激波捕捉准确且无虚假的数值振荡,进而证明了该方法的可行性和有效性。  相似文献   

11.
In this paper, a Hybrid Analytical/Two-Dimensional Finite Element Method (2-D HAFEM) is proposed to analyze wave propagation characteristics of fluid-filled, composite pipes. In the proposed method, a fluid-filled pipe with a constant cross-section is modeled by using a 2-D finite element approximation in the cross-sectional area while an analytical wave solution is assumed in the axial direction. Thus, it makes possible to use a small number of finite elements even for high frequency analyses in a computationally efficient manner. Both solid and fluid elements as well as solid–fluid interface boundary conditions are developed to model the cross-section of the fluid-filled pipe. In addition, an acoustical transfer function (ATF) approach based on the 2-D HAFEM formulation is suggested to analyze a pipe system assembled with multiple pipe sections with different cross-sections. An ATF matrix relating two sets of acoustic wave variables at the ends of each individual pipe section with a constant cross-section is first calculated and the total ATF matrix for the multi-sectional pipe system is then obtained by multiplying all individual ATF matrices. Therefore, the HAFEM-based ATF approach requires significantly low computational resources, in particular, when there are many pipe sections with a same cross-sectional shape since a single 2-D HAFEM model is needed for these pipe sections. For the validation of the proposed method, experimental and full 3-D FE modeling results are compared to the results obtained by using the HAFEM-based ATF procedure.  相似文献   

12.
This paper proposes a new method for investigating the Hopf bifurcation of a curved pipe conveying fluid with nonlinear spring support. The nonlinear equation of motion is derived by forces equilibrium on microelement of the system under consideration. The spatial coordinate of the system is discretized by the differential quadrature method and then the dynamic equation is solved by the Newton-Raphson method. The numerical solutions show that the inner fluid velocity of the Hopf bifurcation point of the curved pipe varies with different values of the parameter,nonlinear spring stiffness. Based on this, the cycle and divergent motions are both found to exist at specific fluid flow velocities with a given value of the nonlinear spring stiffness. The results are useful for further study of the nonlinear dynamic mechanism of the curved fluid conveying pipe.  相似文献   

13.
???????????????????????????о?   总被引:1,自引:0,他引:1  
唐斌 《力学与实践》2009,31(4):32-36
针对连续Bernoulli-Euler和Timoshenko梁单元的动态刚度矩阵,分析了在使用连续梁单元 进行结构动态特性分析中的数值问题. 基于连续梁单元的运动方程,导出了连续 Bernoulli-Euler和Timoshenko梁单元的动态刚度矩阵. 分析了影响动态刚度矩阵中双曲函 数自变量的各个独立变量及其产生的影响,并给出了初估连续梁单元合理长度的方法. 使用 单一连续Bernoulli-Euler和Timoshenko梁单元的动态刚度矩阵分别进行了悬臂梁频响曲线 的数值求解. 研究表明,在合理选择连续梁单元的长度时,大多数工程结构的动态特性分析 中都不会产生数值问题.  相似文献   

14.
For the cases of using the finite curved beam elements and taking the effects of both the shear deformation and rotary inertias into consideration, the literature regarding either free or forced vibration analysis of the curved beams is rare. Thus, this paper tries to determine the dynamic responses of a circular curved Timoshenko beam due to a moving load using the curved beam elements. By taking account of the effect of shear deformation and that of rotary inertias due to bending and torsional vibrations, the stiffness matrix and the mass matrix of the curved beam element were obtained from the force–displacement relations and the kinetic energy equations, respectively. Since all the element property matrices for the curved beam element are derived based on the local polar coordinate system (rather than the local Cartesian one), their coefficients are invariant for any curved beam element with constant radius of curvature and subtended angle and one does not need to transform the property matrices of each curved beam element from the local coordinate system to the global one to achieve the overall property matrices for the entire curved beam structure before they are assembled. The availability of the presented approach has been verified by both the existing analytical solutions for the entire continuum curved beam and the numerical solutions for the entire discretized curved beam composed of the conventional straight beam elements based on either the consistent-mass model or the lumped-mass model. In addition to the typical circular curved beams, a hybrid curved beam composed of one curved-beam segment and two identical straight-beam segments subjected to a moving load was also studied. Influence on the dynamic responses of the curved beams of the slenderness ratio, moving-load speed, shear deformation and rotary inertias was investigated.  相似文献   

15.
In this study, the forced vibration of a curved pipe conveying fluid resting on a nonlinear elastic foundation is considered. The governing equations for the pipe system are formed with the consideration of viscoelastic material, nonlinearity of foundation, external excitation, and extensibility of centre line. Equations governing the in-plane vibration are solved first by the Galerkin method to obtain the static in-plane equilibrium configuration. The out-of-plane vibration is simplified into a constant coefficient gyroscopic system. Subsequently, the method of multiple scales (MMS) is developed to investigate external first and second primary resonances of the out-of-plane vibration in the presence of three-to-one internal resonance between the first two modes. Modulation equations are formed to obtain the steady state solutions. A parametric study is carried out for the first primary resonance. The effects of damping, nonlinear stiffness of the foundation, internal resonance detuning parameter, and the magnitude of the external excitation are investigated through frequency response curves and force response curves. The characteristics of the single mode response and the relationship between single and two mode steady state solutions are revealed for the second primary resonance. The stability analysis is carried out for these plots. Finally, the approximately analytical results are confirmed by the numerical integrations.  相似文献   

16.
It is a new attempt to extend the differential quadrature method (DQM) to stability analysis of the straight and curved centerline pipes conveying fluid. Emphasis is placed on the study of the influences of several parameters on the critical flow velocity. Compared to other methods, this method can more easily deal with the pipe with spring support at its boundaries and asks for much less computing effort while giving aceptable precision in the numerical results. Supported by National Key Project of China (No. PD9521907) and the National Science Foundation of China (No. 19872025).  相似文献   

17.
The nonlinear response characteristics for a dynamic system with a geometric nonlinearity is examined using a multibody dynamics method. The planar system is an initially straight clamped-clamped beam subject to high frequency excitation in the vicinity of its third natural mode. The model includes a pre-applied static axial load, linear bending stiffness and a cubic in-plane stretching force. Constrained flexibility is applied to a multibody method that lumps the beam into N elements for three substructures subjected to the nonlinear partial differential equation of motion and N-1 linear modal constraints. This procedure is verified by d'Alembert's principle and leads to a discrete form of Galerkin's method. A finite difference scheme models the elastic forces. The beam is tuned by the axial force to obtain fourth order internal resonance that demonstrates bimodal and trimodal responses in agreement with low and moderate excitation test results. The continuous Galerkin method is shown to generate results conflicting with the test and multibody method. A new checking function based on Gauss' principle of least constraint is applied to the beam to minimize modal constraint error.  相似文献   

18.
输液曲管平面内振动的波动方法研究   总被引:2,自引:0,他引:2  
采用Flügge曲梁模拟弯曲管道,推导了管内流体的加速度,在总体轴线不可伸长假定的基础上建立了曲管平面内振动的动力学方程;采用波动方法,获得了曲管内振动波的传播和反射矩阵,提出了计算曲管平面内振动固有频率的数值方法。算例分析中,通过计算两端固定半圆形曲管的临界流速并与已有文献结果对比,验证了本文方法的正确性。最后,计算了两端固定半圆形曲管在四种不同流速下的前四阶固有频率,结果表明,管内流速的增大会降低管道的固有频率,当流速增大到某一特定值时,管道的一阶固有频率消失。  相似文献   

19.
赵翔  李思谊  李映辉 《力学学报》2021,53(11):3035-3044
建立了含裂纹损伤的曲梁压电能量俘获系统在强迫振动下的动力学模型. 基于Prescott型压电曲梁力电耦合振动方程的解析解和裂纹截面处的连续性条件, 求解了含裂纹损伤的压电曲梁的格林函数. 根据线性叠加原理, 对含裂纹的力电耦合模型的系统方程解耦, 得到强迫振动下含裂纹损伤的曲梁压电俘能器的输出电压. 在得到模型的强迫振动解析解后, 提出逆方法检测结构中的裂纹损伤, 这一检测方法适用于处于振动状态下的结构. 在数值计算中, 令裂纹深度为零, 通过对比本文的解析解与现有文献中的解析解, 验证了本文解的有效性. 分别分析了含裂纹损伤的压电曲梁的电压响应与裂纹深度、裂纹位置、材料的几何参数以及阻尼之间的关系. 研究结果表明: 裂纹的存在对曲梁式压电俘能器的影响比直梁式更加复杂; 裂纹出现时, 损伤曲梁在健康曲梁的一阶频率值处一定会出现波动并被激励出二阶频率, 此时的二阶频率是开路中健康压电曲梁的一阶频率值; 通过对电压响应的检测可以确定的损伤裂纹的深度和在结构中出现的位置范围; 利用振动问题的解来检测压电曲梁的健康状况是可行且准确的.   相似文献   

20.
两端弹性支承输流管道固有特性研究   总被引:2,自引:1,他引:1  
颜雄  魏莎  毛晓晔  丁虎  陈立群 《力学学报》2022,54(5):1341-1352
输流管道广泛应用于航天航空、石油化工、海洋等重要的工程领域, 其振动特性尤其是系统固有特性一直是国内外学者研究的热点问题. 本文研究了两端弹性支承输流管道横向振动的固有特性, 尤其是在非对称弹性支承下的系统固有特性. 使用哈密顿原理得到了输流管道的控制方程及边界条件, 通过复模态法得到了静态管道的模态函数, 以其作为伽辽金法的势函数和权函数对线性派生系统控制方程进行截断处理. 分析了两端对称支承刚度、两端非对称支承刚度、管道长度以及流体质量比对系统固有频率的影响规律, 重点讨论了管道两端可能形成的非对称支承条件下固有频率的变化规律. 结果表明, 较大的对称支承刚度下管道的第一阶固有频率下降较快; 当管道两端支承刚度变化时, 管道的各阶固有频率在两端支承刚度相等时取得最值; 对于两端非对称支承的管道而言, 两端支承刚度越接近, 第一阶固有频率下降的越快, 而且相应的临界流速越小; 流体的流速越大, 其对两端非对称弹簧支承的管道固有频率的影响更为明显.   相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号