CN113051718B - Static analysis method for packet topology radial loaded ring with extension hypothesis - Google Patents
Static analysis method for packet topology radial loaded ring with extension hypothesis Download PDFInfo
- Publication number
- CN113051718B CN113051718B CN202110248281.2A CN202110248281A CN113051718B CN 113051718 B CN113051718 B CN 113051718B CN 202110248281 A CN202110248281 A CN 202110248281A CN 113051718 B CN113051718 B CN 113051718B
- Authority
- CN
- China
- Prior art keywords
- radial
- ring
- static
- loads
- hypothesis
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000003068 static effect Effects 0.000 title claims abstract description 50
- 238000004458 analytical method Methods 0.000 title claims abstract description 12
- 238000000034 method Methods 0.000 claims abstract description 37
- 238000005315 distribution function Methods 0.000 claims abstract description 12
- 238000005452 bending Methods 0.000 claims description 11
- 230000007935 neutral effect Effects 0.000 claims description 7
- 238000010586 diagram Methods 0.000 description 9
- 239000000463 material Substances 0.000 description 4
- 230000000694 effects Effects 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000005489 elastic deformation Effects 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Complex Calculations (AREA)
Abstract
Description
技术领域technical field
本发明涉及材料力学应力分布与弹性变形领域,尤其涉及一种基于有延展假设的分组拓扑径向受载圆环应力及变形叠加的方法。The invention relates to the field of mechanical stress distribution and elastic deformation of materials, in particular to a method for superimposing stress and deformation of a grouped topology radially loaded circular ring with an extension assumption.
背景技术Background technique
环状结构广泛应用于机械、液压、化工、水利、交通以及航空航天等领域,如齿轮齿圈、汽车轮毂、离心泵叶轮、柱塞泵的配流盘、电机的定转子、陀螺仪以及航空发动机等结构。它们的工作状态可以分为静止、定轴旋转以及多轴旋转。在工程实际中,静态及其动态特性是主要关注的问题,特别是在高负荷、高速度以及高精度要求的情况下。环状结构在承受不平衡的径向力和切向力时,应力分布状态会发生改变,使结构精度下降,振动更加剧烈,甚至缩短使用寿命。因而,为了使结构拥有最佳使用性能,需要首先研究结构在径向力和切向力的作用下的应力分布,使结构获得最佳的强度、刚度和稳定性。The ring structure is widely used in the fields of machinery, hydraulics, chemical industry, water conservancy, transportation and aerospace, such as gear rings, automobile hubs, impellers of centrifugal pumps, valve plates of plunger pumps, stators and rotors of motors, gyroscopes and aeroengines and other structures. Their working states can be divided into stationary, fixed-axis rotation and multi-axis rotation. In engineering practice, static and dynamic characteristics are the main concerns, especially in the case of high load, high speed and high precision requirements. When the annular structure bears unbalanced radial force and tangential force, the stress distribution state will change, which will reduce the structural precision, make the vibration more severe, and even shorten the service life. Therefore, in order to make the structure have the best performance, it is necessary to first study the stress distribution of the structure under the action of radial force and tangential force, so that the structure can obtain the best strength, stiffness and stability.
文献(Sakamoto S,Hirata T,Kobayashi T,et al.Vibration analysisconsidering higher harmonics of electromagnetic forces for rotating electricmachines.J IEEE T Magn,1999,35(3):1662-1665)研究了永磁体产生的径向磁拉力对应力分布的影响。结果表明,径向磁拉力就是致使环形定子产生不平衡应力分布的主要原因,也是导致定子产生振动而辐射噪声的主要原因。Literature (Sakamoto S, Hirata T, Kobayashi T, et al. Vibration analysis considering higher harmonics of electromagnetic forces for rotating electric machines. J IEEE T Magn, 1999, 35 (3): 1662-1665) studied the radial magnetic force produced by permanent magnets. Effect of tension on stress distribution. The results show that the radial magnetic pull is the main reason for the unbalanced stress distribution of the annular stator, and it is also the main reason for the stator to vibrate and radiate noise.
文献(Barber J R,Force and displacement influence functions for thecircular ring.J Journal of Strain Analysis,1978,13(2):77-81)研究了无延展假设条件下超静定圆环内力求解的问题。建立了一种通用方法,使无延展条件下的各种复杂情形都能通过文献所给的分布函数求得。The literature (Barber J R, Force and displacement influence functions for the circular ring. J Journal of Strain Analysis, 1978, 13(2): 77-81) studies the problem of solving the internal force of a statically indeterminate ring under the assumption of no extension. A general method is established, so that various complex situations under the condition of no extension can be obtained through the distribution functions given in the literature.
但值得一提的是,目前,很多学者对圆环进行研究都是基于无延展假设进行的,即假设受载环状结构中性圆的长度保持不变。可实际上,受到外载荷作用的结构必定会产生变形,尤其是在进行静力学研究时,无延展假设不可能满足实际情况。But it is worth mentioning that, at present, many scholars are conducting research on rings based on the assumption of no extension, that is, assuming that the length of the neutral circle of the loaded ring structure remains constant. But in fact, the structure subjected to external loads will inevitably deform, especially in static research, the assumption of no extension cannot meet the actual situation.
发明内容Contents of the invention
本发明的目的是为了克服现有技术中的不足,提供一种有延展假设分组拓扑径向受载圆环静力学分析方法,本发明针对径向受载圆环的静力学分析等问题,基于有延展假设,在圆环的微段上建立静力学模型,采用叠加法来计算基于有延展假设的分组拓扑径向受载圆环的应力分布,使所得结果更加接近于工程实际。The purpose of the present invention is to overcome the deficiencies in the prior art, and to provide a static analysis method for radially loaded circular rings with extended hypothetical grouping topology. The present invention aims at problems such as the static analysis of radially loaded circular rings, based on With the assumption of extension, a static model is established on the micro-segment of the ring, and the superposition method is used to calculate the stress distribution of the radially loaded ring based on the group topology with the assumption of extension, so that the obtained results are closer to the actual engineering.
本发明的目的是通过以下技术方案实现的:The purpose of the present invention is achieved through the following technical solutions:
一种有延展假设分组拓扑径向受载圆环静力学分析方法,包括以下步骤:A static analysis method for radially loaded rings with extended assumptions grouping topology, comprising the following steps:
步骤一、在圆环上截取微段,采用达朗贝尔原理,通过截面法,建立单个径向载荷作用下圆环的静力学模型;
步骤二、通过所述静力学模型建立基于有延展假设的单个径向载荷作用下圆环的各个静力学参数的分布函数;
步骤三、利用叠加法得到基于有延展假设的分组拓扑径向受载圆环的各个静力学参数的函数。
进一步的,所述单个径向载荷作用下圆环的静力学模型具体为:Further, the static model of the ring under the action of the single radial load is specifically:
θ为圆环上某质点的角度,R为圆环中性圆半径,Fef为径向载荷,Fff为虚拟载荷,Fsf为径向内力,Ftf为切向内力,Mbm为弯矩。θ is the angle of a particle on the ring, R is the radius of the neutral circle of the ring, F ef is the radial load, F ff is the virtual load, F sf is the radial internal force, F tf is the tangential internal force, M bm is the bending force moment.
进一步的,所述基于有延展假设的单个径向载荷作用下圆环的各个静力学参数的分布函数具体为:Further, the distribution function of each static parameter of the ring under the action of a single radial load based on the assumption of extension is specifically:
式中,Ftfθ为切向应力,Fsfr为径向应力,v为径向变形,u为切向变形,A(A=bh)为圆环的截面面积,h是径向厚度,b是轴向厚度,E为杨氏模量,I为转动惯量,μ为泊松比。In the formula, F tfθ is the tangential stress, F sfr is the radial stress, v is the radial deformation, u is the tangential deformation, A(A=bh) is the cross-sectional area of the ring, h is the radial thickness, b is Axial thickness, E is Young's modulus, I is moment of inertia, μ is Poisson's ratio.
进一步的,所述基于有延展假设的分组拓扑径向受载圆环的各个静力学参数的函数具体为:Further, the function of each static parameter of the group topology radially loaded ring based on the extension hypothesis is specifically:
式中,N1为径向载荷的组数,N2为每组中径向载荷的个数,N(N=N1 N2)为径向载荷的总个数,αi,j是第i组中第j个径向载荷的位置角,FtfθN为受N个径向载荷时圆环的切向应力,FsfrN为受N个径向载荷时圆环的径向应力,Mbm为受N个径向载荷时圆环的弯矩,vN为受N个径向载荷时圆环的径向变形,uN为受N个径向载荷时圆环的切向变形。In the formula, N 1 is the number of groups of radial loads, N 2 is the number of radial loads in each group, N (N=N 1 N 2 ) is the total number of radial loads, α i,j is the number of radial loads The position angle of the jth radial load in group i, F tfθN is the tangential stress of the ring when it is subjected to N radial loads, F sfrN is the radial stress of the ring when it is subjected to N radial loads, and M bm is Bending moment of the ring when subjected to N radial loads, v N is the radial deformation of the ring when subjected to N radial loads, u N is the tangential deformation of the ring when subjected to N radial loads.
与现有技术相比,本发明的技术方案所带来的有益效果是:Compared with the prior art, the beneficial effects brought by the technical solution of the present invention are:
1.本发明采用达朗贝尔原理基于有延展假设通过截面法在圆环的微段上建立了静力学模型,然后利用算子法求解基于有延展假设的单个径向载荷作用下圆环的各个静力学参数的分布函数,既得到了分布函数,又利用算子法简化了微分方程求解过程;1. The present invention adopts d'Alembert's principle to set up a statics model on the micro-segment of the annulus by the section method based on the assumption of extension, and then utilizes the operator method to solve each of the annulus under the action of a single radial load based on the assumption of extension. The distribution function of the static parameters not only obtains the distribution function, but also simplifies the process of solving the differential equation by using the operator method;
2.本发明利用叠加原理采用叠加法求解基于有延展假设的分组拓扑径向受载圆环的各个静力学参数的函数,为求解复杂情形下受载圆环的各个静力学参数提供了一种新的可供参考的思路、方法;2. The present invention utilizes the principle of superposition and adopts the superposition method to solve the function of each static parameter based on the grouping topology radially loaded ring with extension hypothesis, and provides a method for solving each static parameter of the loaded ring under complex situations New ideas and methods for reference;
3.本发明具有普适、便捷、精确和新颖等特征。根据该方法可研究径向载荷与各个静力学参数之间的关系,比如,径向载荷作用对与圆环切向应力的影响、径向载荷与圆环径向变形的关系等等,亦可求解各种复杂径向载荷作用下各个静力学参数的函数,比如径向应力、切向变形等等,有了这些基本参数以及参数间的关系便可为静力学以及动力学的后续研究都提供便利。3. The present invention has the characteristics of universality, convenience, precision and novelty. According to this method, the relationship between radial load and various static parameters can be studied, for example, the effect of radial load on the tangential stress of the ring, the relationship between radial load and radial deformation of the ring, etc., can also be Solve the functions of various static parameters under the action of various complex radial loads, such as radial stress, tangential deformation, etc. With these basic parameters and the relationship between parameters, it can be used for subsequent research on statics and dynamics. convenient.
附图说明Description of drawings
图1a和图1b为本发明提供的单个径向载荷作用下整环及微段上力的分布示意图;Fig. 1a and Fig. 1b are the distribution schematic diagrams of force on the whole ring and micro-segment under the action of single radial load provided by the present invention;
图2为本发明提供的旋转θk角的基于有延展假设的单个径向载荷作用下整环上力的分布示意图;Fig. 2 is a schematic diagram of the distribution of the force on the whole ring under the action of a single radial load based on the extension assumption of the rotation θ k angle provided by the present invention;
图3a为本发明提供的基于有延展假设的单个径向载荷作用下圆环的切向和径向应力及变形分布的示意图;Figure 3a is a schematic diagram of the tangential and radial stress and deformation distribution of the ring under the action of a single radial load based on the assumption of extension provided by the present invention;
图3b至图3f为本发明提供的基于有延展假设的单个径向载荷作用下圆环的切向和径向应力及变形大小的示意图;Figures 3b to 3f are schematic diagrams of the tangential and radial stress and deformation of the ring under the action of a single radial load based on the assumption of extension provided by the present invention;
图3g为本发明提供的旋转θk角后基于有延展假设的单个径向集中载荷作用下圆环的切向和径向应力及变形分布的示意图;Figure 3g is a schematic diagram of the tangential and radial stress and deformation distribution of the ring under the action of a single radial concentrated load based on the assumption of extension after the rotation θ k angle provided by the present invention;
图4为本发明提供的基于有延展假设的分组拓扑径向受载圆环上力的分布示意图。Fig. 4 is a schematic diagram of the force distribution on the radially loaded ring of the group topology based on the assumption of extension provided by the present invention.
具体实施方式detailed description
以下结合附图和具体实施例对本发明作进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。The present invention will be described in further detail below in conjunction with the accompanying drawings and specific embodiments. It should be understood that the specific embodiments described here are only used to explain the present invention, not to limit the present invention.
本发明实施例提供了一种分组拓扑径向受载圆环应力叠加的方法。该方法运用达朗贝尔原理在圆环的微段上建立静力学模型,采用截面法得到基于有延展假设的单个径向载荷作用下圆环的各个静力学参数的分布函数,再利用叠加法,得到基于有延展假设的分组拓扑径向受载圆环的各个静力学参数的函数,显著提高了圆环静力学分析的准确性、普适性以及真实性,更好地贴近工程实际。本方法也可用于齿轮、齿圈、旋转电机的定转子以及精密仪器中的环形构件等典型周期结构的各个静力学参数的求解。The embodiment of the present invention provides a group topology radially loaded ring stress superposition method. This method uses d'Alembert's principle to establish a static model on the micro-segment of the ring, uses the section method to obtain the distribution function of each static parameter of the ring under the action of a single radial load with the assumption of extension, and then uses the superposition method, The function of each static parameter of the grouped topology radially loaded ring based on the extension assumption is obtained, which significantly improves the accuracy, universality and authenticity of the static analysis of the ring, and is closer to engineering practice. This method can also be used to solve various static parameters of typical periodic structures such as gears, ring gears, stators and rotors of rotating electrical machines, and ring components in precision instruments.
圆环受到分组拓扑径向载荷的作用;该静力学参数的分布函数的叠加方法的基本特征在于:采用叠加法实现圆环的参数求解,具体步骤为:The circular ring is subjected to the radial load of the group topology; the basic feature of the superposition method of the distribution function of the static parameters is: the parameter solution of the circular ring is realized by the superposition method, and the specific steps are as follows:
(1)利用截面法,在圆环的微段上根据力和力矩平衡原理建立单个径向集中力作用下圆环的静力学模型:(1) Using the cross-section method, the static model of the ring under the action of a single radial concentrated force is established on the micro-segment of the ring according to the principle of force and moment balance:
式中,θ为圆环上某质点的角度,R为圆环中性圆半径,Fef为径向集中力,Fff为虚拟支撑力,Fsf为径向内力,Ftf为切向内力,Mbm为弯矩。In the formula, θ is the angle of a particle on the ring, R is the radius of the neutral circle of the ring, F ef is the radial concentrated force, F ff is the virtual support force, F sf is the radial internal force, F tf is the tangential internal force , M bm is the bending moment.
图1a和图1b为单个径向集中力作用下圆环及微段上力的分布图。如图1a所示,圆环中性圆半径是R,圆环径向宽度是h,圆环轴向高度是b。圆环在θ=0处受一个径向集中力Fef的作用,方向向左。在圆环的一周,分布着均匀的虚拟支撑,虚拟支撑产生虚拟力Fff,方向向右。圆环在一个集中力以及均布的虚拟力的作用下维持平衡。为了研究圆环在单个径向集中力作用下的应力分布,在圆环θ(θ∈(0,2π))处截取dθ的微段,采用截面法对微段进行受力分析,如图1b所示,O和O′分别是圆环的几何中心和微段的中点,Fsf、Ftf和Mbm分别为径向内力、切向内力和弯矩。Figure 1a and Figure 1b are the distribution diagrams of the force on the ring and the micro-segment under the action of a single radial concentrated force. As shown in Figure 1a, the radius of the neutral circle of the ring is R, the radial width of the ring is h, and the axial height of the ring is b. The ring is subjected to a radial concentrated force F ef at θ=0, and the direction is leftward. Around the ring, uniform virtual supports are distributed, and the virtual supports generate virtual force F ff , and the direction is rightward. The ring is balanced by a concentrated force and a uniform virtual force. In order to study the stress distribution of the ring under the action of a single radial concentrated force, the micro-segment of dθ is intercepted at the ring θ(θ∈(0,2π)), and the force analysis of the micro-segment is carried out by using the section method, as shown in Figure 1b As shown, O and O' are the geometric center of the ring and the midpoint of the micro-segment, respectively, and F sf , F tf and M bm are the radial internal force, tangential internal force and bending moment, respectively.
(2)由于研究的是圆环的微段,dθ为微量,利用极限的思想,因此,含有微量的三角函数可以化简为:(2) Since the micro-segment of the ring is studied, dθ is a micro-quantity, and the idea of limit is used, therefore, the trigonometric function containing a micro-quantity can be simplified as:
(3)将式(4)-(7)代入式(1)-(3)中,化简可得:(3) Substituting formulas (4)-(7) into formulas (1)-(3) and simplifying can get:
dMbm=RFsfdθ (10)dM bm = RF sf dθ (10)
(4)求径向内力的通解(4) Find the general solution of the radial internal force
由式(8)和(9)可得径向内力Fsf与径向集中力Fef的关系:The relationship between radial internal force F sf and radial concentrated force F ef can be obtained from formulas (8) and (9):
式(11)是径向内力的二阶非齐次微分方程,由式(11)可得特征方程:Equation (11) is the second-order inhomogeneous differential equation of radial internal force, and the characteristic equation can be obtained from Equation (11):
λsf 2+1=0 (12)λ sf 2 +1=0 (12)
式中,λsf是特征方程的特征值,解之得λsf1,2=±i,i为虚数单位。In the formula, λ sf is the eigenvalue of the characteristic equation, and the solution is λ sf1,2 =±i, where i is the imaginary unit.
为了求解微分方程的通解,设该微分方程的一个特解为:In order to solve the general solution of the differential equation, let a special solution of the differential equation be:
Fsf *=θ(a1cosθ+b1sinθ) (13)F sf * =θ(a 1 cosθ+b 1 sinθ) (13)
式中,a1和b1为实数。In the formula, a 1 and b 1 are real numbers.
将式(13)代入式(11)中,化简可得:Substituting formula (13) into formula (11), it can be simplified to get:
由待定系数法可得a1和b1分别为:According to the undetermined coefficient method, a 1 and b 1 can be obtained as:
b1=0 (16)b 1 =0 (16)
因此,微分方程的特解为:Therefore, the particular solution of the differential equation is:
式中,Fsf *为微分方程的特解。In the formula, F sf * is the special solution of the differential equation.
由特征方程得到的特征值以及微分方程的一个特解,可得径向内力的通解为:From the eigenvalue obtained from the characteristic equation and a special solution of the differential equation, the general solution of the radial internal force can be obtained as:
式中,c1和c2为实数。In the formula, c 1 and c 2 are real numbers.
(5)求切向内力及弯矩的通解(5) Find the general solution of tangential internal force and bending moment
由式(8)、(10)和(18)可得:From formulas (8), (10) and (18), we can get:
(6)求径向变形的通解(6) Find the general solution of radial deformation
对于小曲率圆环,由材料力学的知识可知弯矩、切向内力分别与径向变形、切向变形由如下关系:For small curvature rings, it can be known from the knowledge of material mechanics that the bending moment, tangential internal force and radial deformation, tangential deformation have the following relationship:
式中,v为圆环的径向变形,u为圆环的切向变形,E、I、A和μ分别为圆环的弹性模量、转动惯量、截面面积和泊松比。In the formula, v is the radial deformation of the ring, u is the tangential deformation of the ring, E, I, A and μ are the elastic modulus, moment of inertia, cross-sectional area and Poisson's ratio of the ring, respectively.
由式(10)、(11)、(21)和(22)可得二元高阶非齐次微分方程组:From equations (10), (11), (21) and (22), we can obtain binary high-order non-homogeneous differential equations:
对于二元高阶微分方程,可利用算子发求解。For binary high-order differential equations, operators can be used to solve them.
根据算子法的定义,对式(23)进行化简,可得:According to the definition of the operator method, formula (23) is simplified to get:
根据式(24)消去u可得:Eliminate u according to formula (24):
由式(25)可得径向变形的五阶非齐次微分方程特征方程为:From formula (25), the characteristic equation of the fifth-order non-homogeneous differential equation of radial deformation can be obtained as:
λv 5+2λv 3+λv=0 (26)λ v 5 +2λ v 3 +λ v =0 (26)
解之得λv1=0,λv2,3=±i,λv4,5=±i,i为虚数单位。The solution is λ v1 =0, λ v2,3 =±i, λ v4,5 =±i, i is the imaginary unit.
为了求解微分方程的通解,设该微分方程的一个特解为:In order to solve the general solution of the differential equation, let a special solution of the differential equation be:
v*=θ2(a2cosθ+b2 sinθ) (27)v * =θ 2 (a 2 cosθ+b 2 sinθ) (27)
式中,a2和b2为实数,v*为微分方程的特解。。In the formula, a 2 and b 2 are real numbers, and v * is a special solution of the differential equation. .
将式(27)代入式(25)中,化简可得:Substituting formula (27) into formula (25), it can be simplified to get:
由待定系数法可得a2和b2分别为:According to the undetermined coefficient method, a 2 and b 2 can be obtained as:
b2=0 (30)b 2 =0 (30)
因此,微分方程的特解为:Therefore, the particular solution of the differential equation is:
由特征方程得到的特征值以及微分方程的一个特解,可得径向变形的通解为:From the eigenvalue obtained from the characteristic equation and a special solution of the differential equation, the general solution of the radial deformation can be obtained as:
式中,av1~av5均为实数。In the formula, a v1 ~a v5 are all real numbers.
同理根据式(24)消去v,再利用待定系数法可得切向变形的通解为:Similarly, according to formula (24) to eliminate v, and then use the undetermined coefficient method, the general solution of tangential deformation can be obtained as follows:
式中,au1~au5均为实数。In the formula, a u1 ~a u5 are all real numbers.
(7)利用边界条件求各通解的系数(7) Use the boundary conditions to find the coefficients of each general solution
由材料力学知识可知,对于小曲率圆环,任意截面上的转角为:According to the knowledge of material mechanics, for a ring with small curvature, the rotation angle on any section for:
取圆环θ=0处的截面分析,径向内力为:Taking the section analysis of the ring at θ=0, the radial internal force is:
圆环在θ=0处的径向变形为零,即:The radial deformation of the ring at θ=0 is zero, that is:
v(θ=0)=0 (36a)v (θ = 0) = 0 (36a)
圆环的径向变形关于θ=0对称,即:The radial deformation of the ring is symmetrical about θ=0, namely:
v(θ)=v(2π-θ) (36a)v (θ) = v (2π-θ) (36a)
圆环在θ=0和θ=π处的切向变形均为零,即:The tangential deformation of the ring at θ=0 and θ=π is zero, that is:
u(θ=0)=0 (37a)u (θ=0) =0 (37a)
u(θ=π)=0 (37b)u (θ=π) =0 (37b)
圆环的切向变形关于θ=0中心对称,即:The tangential deformation of the ring is symmetrical about the center of θ=0, that is:
u(θ)=-u(2π-θ) (38)u (θ) = -u (2π-θ) (38)
圆环在θ=0和θ=π处的转角均为零,即:The rotation angles of the ring at θ=0 and θ=π are both zero, that is:
(8)求解(8) Solving
将径向变形通解(式(32))、切向变形的通解(式(33))带入式(21)和式(22)中,并分别于式(19)和式(20)对比,得到两组关系式,再联立式(35)~式(39)即可求得各个系数:Bring the general solution of radial deformation (Equation (32)) and the general solution of tangential deformation (Equation (33)) into Equation (21) and Equation (22), and compare them with Equation (19) and Equation (20), respectively, Obtain two sets of relational expressions, and then combine the vertical equations (35) to (39) to obtain each coefficient:
因此,基于有延展假设的受单个径向载荷时,圆环的径向内力、切向内力、弯矩、径向变形以及切向变形分别为:Therefore, when subjected to a single radial load based on the assumption of extension, the radial internal force, tangential internal force, bending moment, radial deformation, and tangential deformation of the ring are respectively:
由材料力学中内力与应力的关系及式(41)-(42)可知,圆环的切向应力和径向应力分别为:From the relationship between internal force and stress in material mechanics and formulas (41)-(42), it can be known that the tangential stress and radial stress of the ring are respectively:
(9)应力的叠加(9) Superposition of stress
以应力为例,如图2所示,圆环在θ=θk处受一个径向载荷Fef的作用,在圆环的一周分布着与载荷相平衡的虚拟力,与图1a相比,径向载荷和虚拟力的大小分别相等,方向分别逆时针旋转θk角。Taking stress as an example, as shown in Figure 2, the ring is subjected to a radial load F ef at θ=θ k , and a virtual force that is balanced with the load is distributed around the ring. Compared with Figure 1a, The magnitude of the radial load and the virtual force are respectively equal, and the direction is rotated counterclockwise by θ k angle respectively.
图3a是在单个径向载荷作用下圆环的切向应力、径向应力、径向变形和切向变形的分布图,其大小如图3b至图3f所示。所取参数均取自下表1-1。图3g是由图3a旋转θk角之后得到的各个静力学参数的分布图。Fig. 3a is a distribution diagram of tangential stress, radial stress, radial deformation and tangential deformation of the ring under a single radial load, and its magnitude is shown in Fig. 3b to Fig. 3f. All parameters are taken from Table 1-1 below. Fig. 3g is a distribution diagram of various static parameters obtained after rotating Fig. 3a by an angle θ k .
表1-1永磁转子基本参数Table 1-1 Basic parameters of permanent magnet rotor
由于应力等静力学参数具有矢量性不能直接进行旋转,加减等操作,故而为方便后续操作,需要将分布函数分解成傅里叶级数的形式。傅里叶级数是将满足狄利赫条件的某周期函数展开成三角函数的线性组合。由于在单个径向载荷的作用下,圆环的切向应力和径向应力的周期为2π,可设式(46)、式(47)、式(43)、式(44)和式(45)的傅里叶级数展开式为:Due to the vector nature of stress and other static parameters, operations such as rotation, addition and subtraction cannot be performed directly. Therefore, for the convenience of subsequent operations, the distribution function needs to be decomposed into the form of Fourier series. The Fourier series is a linear combination that expands a periodic function that satisfies the Dilich condition into a trigonometric function. Since under the action of a single radial load, the period of the tangential stress and radial stress of the ring is 2π, formula (46), formula (47), formula (43), formula (44) and formula (45) can be set )’s Fourier series expansion is:
式中,和分别为:In the formula, and They are:
由将式(53)-(55)化简并代入式(48)可得:By simplifying formulas (53)-(55) and substituting them into formula (48), we can get:
由将式(56)-(58)化简并代入式(49)可得:By simplifying equations (56)-(58) and substituting them into equation (49), we can get:
由将式(59)-(61)化简并代入式(50)可得:By simplifying formulas (59)-(61) and substituting them into formula (50), we can get:
由将式(62)-(64)化简并代入式(51)可得:By simplifying formulas (62)-(64) and substituting them into formula (51), we can get:
由将式(65)-(67)化简并代入式(52)可得:By simplifying formulas (65)-(67) and substituting them into formula (52), we can get:
径向载荷旋转θk角之后,圆环的各个静力学参数的分布也应该旋转相同的角度,即:After the radial load is rotated by θ k angle, the distribution of each static parameter of the ring should also be rotated by the same angle, namely:
式中,θk为第一个集中力旋转到第k个集中力时的旋转角,FtfθN、FsfrN、Mbm、vN、uN分别为第一个集中力旋转θk后得到的切向内力、径向内力、弯矩、径向变形、切向变形。In the formula, θ k is the rotation angle when the first concentrated force rotates to the kth concentrated force, and F tfθN , F sfrN , M bm , v N , u N are obtained after the first concentrated force rotates θ k Tangential internal force, radial internal force, bending moment, radial deformation, tangential deformation.
如图4所示,圆环的一周受到N(N=N1*N2)个径向载荷的作用,这N个径向载荷采用分组拓扑的布置形式,径向载荷被分为N1组,如图Gi1(i1=1,2,...N1)所示,每组N2个,如图Li1,j1(i1=1,2,...N1,j1=1,2,...N2)所示。ψi1描述的是第i1组中第1个径向载荷的位置角,其中,αi1,j1描述的是第i1组中第j1个径向载荷与第i1组中第1个径向载荷间的夹角,其中,αij=(j-1)αi2。利用叠加法,当圆环受到N个径向集中力作用时,设第一个径向集中力作用在θ=0处,圆环的各个静力学参数分别为:As shown in Figure 4, the circle of the ring is affected by N (N=N 1* N 2 ) radial loads, and the N radial loads are arranged in group topology, and the radial loads are divided into N 1 groups , as shown in figure G i1 (i 1 =1,2,...N 1 ), N 2 for each group, as shown in figure L i1,j1 (i 1 =1,2,...N 1 ,j 1 =1,2,...N 2 ). ψ i1 describes the position angle of the first radial load in group i 1 , where, α i1,j1 describes the angle between the j 1st radial load in the i 1st group and the 1st radial load in the i 1st group, where α ij =(j-1)α i2 . Using the superposition method, when the ring is subjected to N radial concentrated forces, the first radial concentrated force is assumed to act at θ=0, and the static parameters of the ring are respectively:
综上所述,本发明实施例提供了一种分组拓扑径向受载圆环应力叠加的方法。该方法运用达朗贝尔原理在圆环的微段上建立静力学模型,采用截面法得到基于有延展假设的单个径向载荷作用下圆环的各个静力学参数的分布函数,再利用叠加法,得到基于有延展假设的分组拓扑径向受载圆环的各个静力学参数的函数,显著提高了圆环静力学分析的准确性、普适性以及真实性,更好地贴近工程实际。To sum up, the embodiment of the present invention provides a method for stress superposition of grouped topology radially loaded circular rings. This method uses d'Alembert's principle to establish a static model on the micro-segment of the ring, uses the section method to obtain the distribution function of each static parameter of the ring under the action of a single radial load with the assumption of extension, and then uses the superposition method, The function of each static parameter of the grouped topology radially loaded ring based on the extension assumption is obtained, which significantly improves the accuracy, universality and authenticity of the static analysis of the ring, and is closer to engineering practice.
本发明实施例对各器件的型号除做特殊说明的以外,其他器件的型号不做限制,只要能完成上述功能的器件均可。In the embodiments of the present invention, unless otherwise specified, the models of the devices are not limited, as long as they can complete the above functions.
本发明并不限于上文描述的实施方式。以上对具体实施方式的描述旨在描述和说明本发明的技术方案,上述的具体实施方式仅仅是示意性的,并不是限制性的。在不脱离本发明宗旨和权利要求所保护的范围情况下,本领域的普通技术人员在本发明的启示下还可做出很多形式的具体变换,这些均属于本发明的保护范围之内。The present invention is not limited to the embodiments described above. The above description of the specific embodiments is intended to describe and illustrate the technical solution of the present invention, and the above specific embodiments are only illustrative and not restrictive. Without departing from the gist of the present invention and the scope of protection of the claims, those skilled in the art can also make many specific changes under the inspiration of the present invention, and these all belong to the protection scope of the present invention.
Claims (3)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110248281.2A CN113051718B (en) | 2021-03-07 | 2021-03-07 | Static analysis method for packet topology radial loaded ring with extension hypothesis |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110248281.2A CN113051718B (en) | 2021-03-07 | 2021-03-07 | Static analysis method for packet topology radial loaded ring with extension hypothesis |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113051718A CN113051718A (en) | 2021-06-29 |
CN113051718B true CN113051718B (en) | 2023-01-03 |
Family
ID=76510169
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110248281.2A Active CN113051718B (en) | 2021-03-07 | 2021-03-07 | Static analysis method for packet topology radial loaded ring with extension hypothesis |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113051718B (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3267339A1 (en) * | 2016-07-06 | 2018-01-10 | Siemens AG Österreich | Method for creating a meta-model for calculating notch stresses |
CN107808048A (en) * | 2017-10-26 | 2018-03-16 | 哈尔滨工程大学 | A kind of vibration analysis method of uniform cross-section girder construction |
CN110569560A (en) * | 2019-08-16 | 2019-12-13 | 天津大学 | A Stress Superposition Method of Mirror Topology Tangentially Loaded Ring |
CN110580383A (en) * | 2019-08-16 | 2019-12-17 | 天津大学 | A Method of Grouping Topological Radially Loaded Circular Stress Superposition |
CN112347591A (en) * | 2020-11-30 | 2021-02-09 | 天津大学 | An internal force analysis and free vibration modeling method of an eccentric rotating annular structure |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102651054B (en) * | 2012-04-10 | 2014-04-30 | 北京航空航天大学 | Probability method of electronic product service life model based on Bayesian theory |
-
2021
- 2021-03-07 CN CN202110248281.2A patent/CN113051718B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3267339A1 (en) * | 2016-07-06 | 2018-01-10 | Siemens AG Österreich | Method for creating a meta-model for calculating notch stresses |
CN107808048A (en) * | 2017-10-26 | 2018-03-16 | 哈尔滨工程大学 | A kind of vibration analysis method of uniform cross-section girder construction |
CN110569560A (en) * | 2019-08-16 | 2019-12-13 | 天津大学 | A Stress Superposition Method of Mirror Topology Tangentially Loaded Ring |
CN110580383A (en) * | 2019-08-16 | 2019-12-17 | 天津大学 | A Method of Grouping Topological Radially Loaded Circular Stress Superposition |
CN112347591A (en) * | 2020-11-30 | 2021-02-09 | 天津大学 | An internal force analysis and free vibration modeling method of an eccentric rotating annular structure |
Non-Patent Citations (1)
Title |
---|
环状周期结构静力学与动力学特性研究;柳金龙;《万方平台在线出版》;20210129;正文第三至四章 * |
Also Published As
Publication number | Publication date |
---|---|
CN113051718A (en) | 2021-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106547957B (en) | A parametric elastic vibration analysis method for rotating annular periodic structures | |
CN110580383B (en) | Grouping topology radial loaded ring stress superposition method | |
Zhu et al. | Investigation of bread-loaf magnet on vibration performance in FSCW PMSM considering force modulation effect | |
CN106528959B (en) | A Simplified Analysis Method for the Natural Frequency and Stability of Rotationally Symmetrical Structures | |
CN110046418B (en) | A Method for Analyzing Vibration Characteristics of Periodic Stator of Permanent Magnet Motor | |
CN109284569B (en) | Combined instability analysis and prediction method for periodic stator of permanent magnet motor | |
Zheng et al. | Rotor balancing for magnetically levitated TMPs integrated with vibration self-sensing of magnetic bearings | |
Liu et al. | Model and experimental verification of a four degrees-of-freedom rotor considering combined eccentricity and electromagnetic effects | |
CN113051718B (en) | Static analysis method for packet topology radial loaded ring with extension hypothesis | |
CN113515843A (en) | Permanent magnet motor rotor strength optimization design method based on multi-dimensional visualization | |
CN102025304A (en) | Method for reducing vibration and noise of permanent-magnet motor | |
CN112953344A (en) | Unbalance vibration compensation control method for rotor of bearingless asynchronous motor | |
CN110569560B (en) | Method for superposing stresses of mirror topology tangential loaded circular ring | |
Li et al. | Analysis on the vibration modes of the electric vehicle motor stator | |
CN105550395A (en) | Turbomachinery single supporting shafting geometric centering installation method | |
Martynenko | Stability analysis of rotor motion in nonlinear systems with passive and active magnetic bearings | |
Zhou et al. | Research on steady-state characteristics of centrifugal pump rotor system with weak nonlinear stiffness | |
Rodriguez et al. | Designing, simulations and experiments of a passive permanent magnet bearing | |
CN114611315A (en) | Calculation and analysis method for unbalanced electromagnetic force of motor | |
Ma et al. | Characteristic analysis and direct measurement for air gap magnetic field of external rotor permanent magnet synchronous motors in electric vehicles | |
CN110569574B (en) | Method for improving out-of-plane vibration stability of rotor by using sine-shaped magnetic poles of permanent magnet motor | |
Xing et al. | Critical speed calculation and modal analysis of rotor for high-speed machine | |
Zhang et al. | Complete control of radial suspension force for bearingless induction motors | |
CN110555254A (en) | Method for reducing out-of-plane vibration instability of permanent magnet motor by grouping magnetic poles | |
CN112861291B (en) | Multistage centrifugal pump annular seal design method based on damping ratio analysis |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |