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 PDF

Info

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
Application number
CN202110248281.2A
Other languages
Chinese (zh)
Other versions
CN113051718A (en
Inventor
王世宇
王久霞
夏春花
刘晨
魏振航
王一凡
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin University
Original Assignee
Tianjin University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tianjin University filed Critical Tianjin University
Priority to CN202110248281.2A priority Critical patent/CN113051718B/en
Publication of CN113051718A publication Critical patent/CN113051718A/en
Application granted granted Critical
Publication of CN113051718B publication Critical patent/CN113051718B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force 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

The invention discloses a static analysis method for a radial loaded ring of an extended hypothesis packet topology, which comprises the following steps: firstly, intercepting a micro-segment on a circular ring, and establishing a static model of the circular ring under the action of a single radial load by adopting a Dalnbell principle and a cross-section method; secondly, establishing a distribution function of each static parameter of the ring under the action of a single radial load with an extension hypothesis through the static model; and thirdly, obtaining a function of each static parameter of the grouped topological radial loaded circular ring based on the extension hypothesis by using an overlap method.

Description

一种有延展假设分组拓扑径向受载圆环静力学分析方法A Static Analysis Method for Radially Loaded Circular Ring with Extended Assumption Grouping Topology

技术领域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:

步骤一、在圆环上截取微段,采用达朗贝尔原理,通过截面法,建立单个径向载荷作用下圆环的静力学模型;Step 1, intercepting micro-segments on the ring, adopting D'Alembert's principle, and establishing a static model of the ring under the action of a single radial load through the section method;

步骤二、通过所述静力学模型建立基于有延展假设的单个径向载荷作用下圆环的各个静力学参数的分布函数;Step 2, establishing the distribution function of each static parameter of the ring under the action of a single radial load based on the assumption of extension through the static model;

步骤三、利用叠加法得到基于有延展假设的分组拓扑径向受载圆环的各个静力学参数的函数。Step 3, using the superposition method to obtain the functions of the static parameters of the grouped topology radially loaded ring based on the extension assumption.

进一步的,所述单个径向载荷作用下圆环的静力学模型具体为:Further, the static model of the ring under the action of the single radial load is specifically:

Figure BDA0002964960800000021
Figure BDA0002964960800000021

Figure BDA0002964960800000022
Figure BDA0002964960800000022

Figure BDA0002964960800000023
Figure BDA0002964960800000023

θ为圆环上某质点的角度,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:

Figure BDA0002964960800000024
Figure BDA0002964960800000024

Figure BDA0002964960800000025
Figure BDA0002964960800000025

Figure BDA0002964960800000026
Figure BDA0002964960800000026

Figure BDA0002964960800000031
Figure BDA0002964960800000031

Figure BDA0002964960800000032
Figure BDA0002964960800000032

式中,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:

Figure BDA0002964960800000033
Figure BDA0002964960800000033

Figure BDA0002964960800000034
Figure BDA0002964960800000034

Figure BDA0002964960800000035
Figure BDA0002964960800000035

Figure BDA0002964960800000036
Figure BDA0002964960800000036

Figure BDA0002964960800000037
Figure BDA0002964960800000037

式中,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:

Figure BDA0002964960800000051
Figure BDA0002964960800000051

Figure BDA0002964960800000052
Figure BDA0002964960800000052

Figure BDA0002964960800000053
Figure BDA0002964960800000053

式中,θ为圆环上某质点的角度,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:

Figure BDA0002964960800000061
Figure BDA0002964960800000061

Figure BDA0002964960800000062
Figure BDA0002964960800000062

Figure BDA0002964960800000063
Figure BDA0002964960800000063

Figure BDA0002964960800000064
Figure BDA0002964960800000064

(3)将式(4)-(7)代入式(1)-(3)中,化简可得:(3) Substituting formulas (4)-(7) into formulas (1)-(3) and simplifying can get:

Figure BDA0002964960800000065
Figure BDA0002964960800000065

Figure BDA0002964960800000066
Figure BDA0002964960800000066

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):

Figure BDA0002964960800000067
Figure BDA0002964960800000067

式(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:

Figure BDA0002964960800000068
Figure BDA0002964960800000068

由待定系数法可得a1和b1分别为:According to the undetermined coefficient method, a 1 and b 1 can be obtained as:

Figure BDA0002964960800000069
Figure BDA0002964960800000069

b1=0 (16)b 1 =0 (16)

因此,微分方程的特解为:Therefore, the particular solution of the differential equation is:

Figure BDA0002964960800000071
Figure BDA0002964960800000071

式中,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:

Figure BDA0002964960800000072
Figure BDA0002964960800000072

式中,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:

Figure BDA0002964960800000073
Figure BDA0002964960800000073

Figure BDA0002964960800000074
Figure BDA0002964960800000074

(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:

Figure BDA0002964960800000075
Figure BDA0002964960800000075

Figure BDA0002964960800000076
Figure BDA0002964960800000076

式中,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:

Figure BDA0002964960800000077
Figure BDA0002964960800000077

对于二元高阶微分方程,可利用算子发求解。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:

Figure BDA0002964960800000081
Figure BDA0002964960800000081

根据式(24)消去u可得:Eliminate u according to formula (24):

Figure BDA0002964960800000082
Figure BDA0002964960800000082

由式(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 3v=0 (26)λ v 5 +2λ v 3v =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:

Figure BDA0002964960800000083
Figure BDA0002964960800000083

由待定系数法可得a2和b2分别为:According to the undetermined coefficient method, a 2 and b 2 can be obtained as:

Figure BDA0002964960800000084
Figure BDA0002964960800000084

b2=0 (30)b 2 =0 (30)

因此,微分方程的特解为:Therefore, the particular solution of the differential equation is:

Figure BDA0002964960800000085
Figure BDA0002964960800000085

由特征方程得到的特征值以及微分方程的一个特解,可得径向变形的通解为: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:

Figure BDA0002964960800000086
Figure BDA0002964960800000086

式中,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:

Figure BDA0002964960800000091
Figure BDA0002964960800000091

式中,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

由材料力学知识可知,对于小曲率圆环,任意截面上的转角

Figure BDA0002964960800000092
为:According to the knowledge of material mechanics, for a ring with small curvature, the rotation angle on any section
Figure BDA0002964960800000092
for:

Figure BDA0002964960800000093
Figure BDA0002964960800000093

取圆环θ=0处的截面分析,径向内力为:Taking the section analysis of the ring at θ=0, the radial internal force is:

Figure BDA0002964960800000094
Figure BDA0002964960800000094

圆环在θ=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:

Figure BDA0002964960800000095
Figure BDA0002964960800000095

Figure BDA0002964960800000096
Figure BDA0002964960800000096

(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:

Figure BDA0002964960800000101
Figure BDA0002964960800000101

因此,基于有延展假设的受单个径向载荷时,圆环的径向内力、切向内力、弯矩、径向变形以及切向变形分别为: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:

Figure BDA0002964960800000102
Figure BDA0002964960800000102

Figure BDA0002964960800000103
Figure BDA0002964960800000103

Figure BDA0002964960800000104
Figure BDA0002964960800000104

Figure BDA0002964960800000105
Figure BDA0002964960800000105

Figure BDA0002964960800000111
Figure BDA0002964960800000111

由材料力学中内力与应力的关系及式(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:

Figure BDA0002964960800000112
Figure BDA0002964960800000112

Figure BDA0002964960800000113
Figure BDA0002964960800000113

(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

Figure BDA0002964960800000114
Figure BDA0002964960800000114

由于应力等静力学参数具有矢量性不能直接进行旋转,加减等操作,故而为方便后续操作,需要将分布函数分解成傅里叶级数的形式。傅里叶级数是将满足狄利赫条件的某周期函数展开成三角函数的线性组合。由于在单个径向载荷的作用下,圆环的切向应力和径向应力的周期为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:

Figure BDA0002964960800000121
Figure BDA0002964960800000121

Figure BDA0002964960800000122
Figure BDA0002964960800000122

Figure BDA0002964960800000123
Figure BDA0002964960800000123

Figure BDA0002964960800000124
Figure BDA0002964960800000124

Figure BDA0002964960800000125
Figure BDA0002964960800000125

式中,

Figure BDA0002964960800000126
Figure BDA0002964960800000127
分别为:In the formula,
Figure BDA0002964960800000126
and
Figure BDA0002964960800000127
They are:

Figure BDA0002964960800000128
Figure BDA0002964960800000128

Figure BDA0002964960800000129
Figure BDA0002964960800000129

Figure BDA00029649608000001210
Figure BDA00029649608000001210

Figure BDA00029649608000001211
Figure BDA00029649608000001211

Figure BDA00029649608000001212
Figure BDA00029649608000001212

Figure BDA0002964960800000131
Figure BDA0002964960800000131

Figure BDA0002964960800000132
Figure BDA0002964960800000132

Figure BDA0002964960800000133
Figure BDA0002964960800000133

Figure BDA0002964960800000134
Figure BDA0002964960800000134

Figure BDA0002964960800000135
Figure BDA0002964960800000135

Figure BDA0002964960800000136
Figure BDA0002964960800000136

Figure BDA0002964960800000137
Figure BDA0002964960800000137

Figure BDA0002964960800000138
Figure BDA0002964960800000138

Figure BDA0002964960800000141
Figure BDA0002964960800000141

Figure BDA0002964960800000142
Figure BDA0002964960800000142

由将式(53)-(55)化简并代入式(48)可得:By simplifying formulas (53)-(55) and substituting them into formula (48), we can get:

Figure BDA0002964960800000143
Figure BDA0002964960800000143

由将式(56)-(58)化简并代入式(49)可得:By simplifying equations (56)-(58) and substituting them into equation (49), we can get:

Figure BDA0002964960800000144
Figure BDA0002964960800000144

由将式(59)-(61)化简并代入式(50)可得:By simplifying formulas (59)-(61) and substituting them into formula (50), we can get:

Figure BDA0002964960800000145
Figure BDA0002964960800000145

由将式(62)-(64)化简并代入式(51)可得:By simplifying formulas (62)-(64) and substituting them into formula (51), we can get:

Figure BDA0002964960800000151
Figure BDA0002964960800000151

由将式(65)-(67)化简并代入式(52)可得:By simplifying formulas (65)-(67) and substituting them into formula (52), we can get:

Figure BDA0002964960800000152
Figure BDA0002964960800000152

径向载荷旋转θ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:

Figure BDA0002964960800000153
Figure BDA0002964960800000153

Figure BDA0002964960800000154
Figure BDA0002964960800000154

Figure BDA0002964960800000155
Figure BDA0002964960800000155

Figure BDA0002964960800000156
Figure BDA0002964960800000156

Figure BDA0002964960800000161
Figure BDA0002964960800000161

Figure BDA0002964960800000162
Figure BDA0002964960800000162

式中,θ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个径向载荷的位置角,其中,

Figure BDA0002964960800000163
α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,
Figure BDA0002964960800000163
α 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:

Figure BDA0002964960800000164
Figure BDA0002964960800000164

Figure BDA0002964960800000165
Figure BDA0002964960800000165

Figure BDA0002964960800000166
Figure BDA0002964960800000166

Figure BDA0002964960800000171
Figure BDA0002964960800000171

Figure BDA0002964960800000172
Figure BDA0002964960800000172

Figure BDA0002964960800000173
Figure BDA0002964960800000173

Figure BDA0002964960800000174
Figure BDA0002964960800000174

综上所述,本发明实施例提供了一种分组拓扑径向受载圆环应力叠加的方法。该方法运用达朗贝尔原理在圆环的微段上建立静力学模型,采用截面法得到基于有延展假设的单个径向载荷作用下圆环的各个静力学参数的分布函数,再利用叠加法,得到基于有延展假设的分组拓扑径向受载圆环的各个静力学参数的函数,显著提高了圆环静力学分析的准确性、普适性以及真实性,更好地贴近工程实际。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)

1. A static analysis method for a radial loaded ring of an extended hypothesis packet topology is characterized by comprising the following steps:
firstly, intercepting a micro-segment on a circular ring, and establishing a static model of the circular ring under the action of a single radial load by adopting a Dalnbell principle and a cross-section method;
secondly, establishing a distribution function of each static parameter of the ring under the action of a single radial load with an extension hypothesis through the static model;
the distribution function of each static parameter of the ring under the action of a single radial load based on the extension hypothesis is specifically as follows:
Figure FDA0003889377550000011
Figure FDA0003889377550000012
Figure FDA0003889377550000013
Figure FDA0003889377550000014
Figure FDA0003889377550000015
wherein, theta is the angle of a certain mass point on the circular ring, and R is the circular ringRadius of neutral circle, F ef For radial loading, M bm Is a bending moment, F tfθ As tangential stress, F sfr Radial stress, v radial deformation, u tangential deformation, A cross-sectional area of the ring, and A = bh; h is the radial thickness, b is the axial thickness, E is the Young's modulus, I is the moment of inertia, and μ is the Poisson's ratio;
and thirdly, obtaining a function of each static parameter of the grouped topological radial loaded circular ring based on the extension hypothesis by using an overlap method.
2. The method for statically analyzing the radially loaded ring with the extended hypothesis packet topology according to claim 1, wherein the static model of the ring under the action of the single radial load is specifically:
Figure FDA0003889377550000021
Figure FDA0003889377550000022
Figure FDA0003889377550000023
where θ is the angle of a certain point on the ring, R is the neutral radius of the ring, and F ef For radial loading, F ff For virtual loads, F sf As a radially internal force, F tf As a tangential internal force, M bm Is a bending moment.
3. The method for analyzing statics of the packet topology radial loaded ring with the extended hypothesis according to claim 1, wherein the function of each static parameter of the packet topology radial loaded ring with the extended hypothesis is specifically:
Figure FDA0003889377550000024
Figure FDA0003889377550000025
Figure FDA0003889377550000026
Figure FDA0003889377550000027
Figure FDA0003889377550000028
where θ is the angle of a certain mass point on the ring, R is the neutral radius of the ring, and F ef For radial load, A = bh is the cross-sectional area of the ring, E is Young's modulus, I is the moment of inertia, μ is Poisson's ratio, N 1 Number of groups of radial loads, i 1 =1,2,...N 1 ,N 2 The number of radial loads in each group, N is the total number of radial loads, N = N 1 N 2 ,α i,j Is the position angle of the jth radial load in the ith group, F tfθN For tangential stress of the ring under N radial loads, F sfrN For radial stress of the ring under N radial loads, M bmN Bending moment of the ring under N radial loads, v N For radial deformation of the ring under N radial loads, u N Is the tangential deformation of the ring under N radial loads.
CN202110248281.2A 2021-03-07 2021-03-07 Static analysis method for packet topology radial loaded ring with extension hypothesis Active CN113051718B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (5)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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