CN112464400B - 基于库伦摩擦和黏性摩擦的径向驻波型超声波电机转矩转速特性的计算方法 - Google Patents
基于库伦摩擦和黏性摩擦的径向驻波型超声波电机转矩转速特性的计算方法 Download PDFInfo
- Publication number
- CN112464400B CN112464400B CN202011313598.1A CN202011313598A CN112464400B CN 112464400 B CN112464400 B CN 112464400B CN 202011313598 A CN202011313598 A CN 202011313598A CN 112464400 B CN112464400 B CN 112464400B
- Authority
- CN
- China
- Prior art keywords
- rotor
- standing wave
- elastic blade
- load
- torque
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 38
- 238000000034 method Methods 0.000 claims abstract description 36
- 238000003860 storage Methods 0.000 claims description 12
- 238000005452 bending Methods 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 10
- 238000009434 installation Methods 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 8
- 238000012790 confirmation Methods 0.000 claims description 3
- 230000007547 defect Effects 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 229920001967 Metal rubber Polymers 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000000919 ceramic Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000013433 optimization analysis Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02N—ELECTRIC MACHINES NOT OTHERWISE PROVIDED FOR
- H02N2/00—Electric machines in general using piezoelectric effect, electrostriction or magnetostriction
- H02N2/10—Electric machines in general using piezoelectric effect, electrostriction or magnetostriction producing rotary motion, e.g. rotary motors
-
- 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
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- General Electrical Machinery Utilizing Piezoelectricity, Electrostriction Or Magnetostriction (AREA)
Abstract
本发明公开了基于库伦摩擦和黏性摩擦的径向驻波型超声波电机转矩转速特性的计算方法,包括:根据径向驻波型超声波电机的基本参数和预设的转子转速N的值,计算与转子转速N相对应的电机输出转矩T load;判断T load的大小,若T load>0,重置转子转速N的值并计算T load,若T load≤0,结束计算;根据计算过程中所有转子转速N及其相对应的电机输出转矩T load,确定径向驻波型超声波电机的转矩转速特性;径向驻波型超声波电机的基本参数包括弹性叶片与转子间的库伦摩擦系数μ d和黏性摩擦系数f v。本发明能够弥补采用库伦摩擦定律导致计算结果误差大的不足,更准确地计算出径向驻波型超声波电机的转矩转速特性。
Description
技术领域
本发明涉及基于库伦摩擦和黏性摩擦的径向驻波型超声波电机转矩转速特性的计算方法,属于超声波电机技术领域。
背景技术
超声波电机利用压电陶瓷的逆压电效应,将高频交流电转换为金属弹性体的超声振动,然后通过摩擦传动将超声振动转化为转子或动子的旋转或直线运动。与传统电磁电机相比,超声波电机具有结构简单、断电自锁、动态响应快、控制精度高、电磁兼容性好等优点。
超声波电机主要有行波型和驻波型两种。行波超声波电机与驻波超声波电机的接触机理不同。在行波超声波电机中,定转子处于连续接触状态,定转子之间的接触区域随定子波峰的变化而变化。相反,在驻波超声波电机中,定子和转子处于间歇接触状态。在对驻波型超声波电机的研究中,以采用板式弯曲振动的驻波型超声波电机较为常见,且运动形式以直线运动为主,而对采用径向伸缩振动的驻波型超声波电机的研究较少。
目前对径向驻波型超声波电机转矩转速特性的计算,常基于库伦摩擦定律。库伦摩擦定律虽然简单,但无法准确描述摩擦界面的诸多特性,导致计算结果往往存在较大的误差。
发明内容
本发明的目的在于克服现有技术中的不足,提供基于库伦摩擦和黏性摩擦的径向驻波型超声波电机转矩转速特性的计算方法,将黏性摩擦与库伦摩擦定律相结合,以弥补单纯采用库伦摩擦定律导致计算结果误差大的不足,能够更准确地计算出径向驻波型超声波电机的转矩转速特性。为达到上述目的,本发明是采用下述技术方案实现的:
第一方面,本发明提供了径向驻波型超声波电机转矩转速特性的计算方法,包括:
根据径向驻波型超声波电机的基本参数和预设的转子转速N的值,计算与转子转速N相对应的电机输出转矩Tload;
判断与转子转速N相对应的电机输出转矩Tload的大小,若与转子转速N相对应的电机输出转矩Tload>0,重置转子转速N的值并计算与转子转速N相对应的电机输出转矩Tload,若与转子转速N相对应的电机输出转矩Tload≤0,结束计算;
根据计算过程中所有转子转速N以及与转子转速N相对应的电机输出转矩Tload,确定径向驻波型超声波电机的转矩转速特性;
所述径向驻波型超声波电机的基本参数包括弹性叶片与转子间的库伦摩擦系数μd和黏性摩擦系数fv。
结合第一方面,进一步地,所述径向驻波型超声波电机转矩转速特性的计算方法为基于库伦摩擦和黏性摩擦的径向驻波型超声波电机转矩转速特性的计算方法。
结合第一方面,进一步地,与转子转速N相对应的电机输出转矩Tload通过以下公式计算得到:
式(1)中,n为弹性叶片的数量,r为转子的内半径,T为驱动电源的周期,μd为弹性叶片与转子间的库伦摩擦系数,Fp为使弹性叶片产生形变的压力,t1为等速时间点,α为弹性叶片形变后与定子径向的夹角,fv为黏性摩擦系数,Vm为弹性叶片的切向速度幅值,Vr为转子的切向线速度。
结合第一方面,进一步地,等速时间点t1通过以下公式计算得到:
式(2)中,T为驱动电源的周期,Vm为弹性叶片的切向速度幅值,Vr为转子的切向线速度,通过以下公式计算得到:
式(3)中,r为转子的内半径,N为转子转速的值。
结合第一方面,进一步地,弹性叶片的切向速度幅值Vm为弹性叶片的运动参数,通过以下公式计算得到:
式(4)中,A为定子一阶径向振动的振幅,T为驱动电源的周期,β为弹性叶片的安装倾斜角,α为弹性叶片形变后与定子径向的夹角,为弹性叶片产生形变的弯曲角。
结合第一方面,进一步地,弹性叶片形变后与定子径向的夹角α、弹性叶片产生形变的弯曲角和使弹性叶片产生形变的压力Fp为弹性叶片的形变参数;
弹性叶片形变后与定子径向的夹角α,通过以下公式计算得到:
式(5)中,r为转子的内半径,l为弹性叶片的长度,r1为定子的外半径;
弹性叶片产生形变的弯曲角通过以下公式计算得到:
式(6)中,β为弹性叶片的安装倾斜角;
使弹性叶片产生形变的压力Fp,通过以下公式计算得到:
式(7)中,E为弹性叶片的弹性模量,b为弹性叶片的宽度,h为弹性叶片的厚度。
结合第一方面,进一步地,重置转子转速N的值并计算与转子转速N相对应的电机输出转矩Tload包括:
将转子转速N的值加上步进常数Q,得到重置后的转子转速N的值;
根据径向驻波型超声波电机的基本参数和重置后的转子转速N的值,计算与重置后的转子转速N相对应的电机输出转矩Tload。
结合第一方面,优选地,步进常数Q的取值范围为0<Q≤1。
结合第一方面,进一步地,驱动电源的周期T、定子一阶径向振动的振幅A、弹性叶片的长度l、弹性叶片的安装倾斜角β、弹性叶片的数量n、弹性叶片的弹性模量E、弹性叶片的宽度b、弹性叶片的厚度h、定子的外半径r1和转子的内半径r为所述径向驻波型超声波电机的基本参数。
结合第一方面,进一步地,根据计算过程中所有转子转速N以及与转子转速N相对应的电机输出转矩Tload,得到反映径向驻波型超声波电机的转矩转速特性的特性曲线。
第二方面,本发明提供了基于库伦摩擦和黏性摩擦的径向驻波型超声波电机转矩转速特性的计算装置,包括:
计算模块:用于根据径向驻波型超声波电机的基本参数和预设的转子转速N的值,计算与转子转速N相对应的电机输出转矩Tload;
判断模块:用于判断与转子转速N相对应的电机输出转矩Tload的大小,若与转子转速N相对应的电机输出转矩Tload>0,重置转子转速N的值并重新计算与转子转速N相对应的电机输出转矩Tload,若与转子转速N相对应的电机输出转矩Tload≤0,结束计算;
转矩转速特性确认模块:根据计算过程中所有转子转速N以及与转子转速N相对应的电机输出转矩Tload,确定径向驻波型超声波电机的转矩转速特性。
第三方面,本发明提供了基于库伦摩擦和黏性摩擦的径向驻波型超声波电机转矩转速特性的计算装置,所述装置包括处理器及存储介质;
所述存储介质用于存储指令;
所述处理器用于根据所述指令进行操作以执行第一方面任一项所述方法的步骤。
第四方面,本发明提供了一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现第一方面任一项所述方法的步骤。
与现有技术相比,本发明实施例所提供的基于库伦摩擦和黏性摩擦的径向驻波型超声波电机转矩转速特性的计算方法所达到的有益效果包括:
本发明提供的根据径向驻波型超声波电机的基本参数和预设的转子转速N的值,计算与转子转速N相对应的电机输出转矩Tload,径向驻波型超声波电机的基本参数包括弹性叶片与转子间的库伦摩擦系数和黏性摩擦系数,基于库伦摩擦系数和黏性摩擦系数计算径向驻波型超声波电机的输出转矩,计算效率高,计算准确度高,修正了现有的计算方法中在摩擦力计算方面的不足,在电机的设计分析阶段更准确地计算径向驻波型超声波电机的转矩转速特性;
本发明的径向驻波型超声波电机的转矩转速特性的计算方法,不仅可用于径向驻波型超声波电机的结构设计和优化分析,计算不同尺寸的径向驻波型超声波电机的转矩转速特性,还可以用于对径向驻波型超声波电机的接触机理的分析。
附图说明
图1是本发明基于库伦摩擦和黏性摩擦的径向驻波型超声波电机转矩转速特性的计算方法的流程图;
图2是本发明基于库伦摩擦和黏性摩擦的径向驻波型超声波电机转矩转速特性曲线的示意图。
具体实施方式
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
实施例一:
如图1所示,基于库伦摩擦和黏性摩擦的径向驻波型超声波电机转矩转速特性的计算方法,包括如下步骤:
步骤1:获取径向驻波型超声波电机的基本参数。
径向驻波型超声波电机的基本参数包括:驱动电源的周期T、定子一阶径向振动的振幅A、弹性叶片的长度l、弹性叶片的安装倾斜角β、弹性叶片的数量n、弹性叶片的弹性模量E、弹性叶片的宽度b、弹性叶片的厚度h、弹性叶片与转子间的库伦摩擦系数μd、黏性摩擦系数fv、定子的外半径r1、转子的内半径r。
步骤2:计算弹性叶片的形变参数和运动参数。
弹性叶片的形变参数包括:弹性叶片形变后与定子径向的夹角α、弹性叶片产生形变的弯曲角使弹性叶片产生形变的压力Fp;弹性叶片的运动参数包括:弹性叶片的切向速度幅值Vm。
具体地,弹性叶片形变后与定子径向的夹角α,通过以下公式计算得到:
式(1)中,r为转子的内半径,l为弹性叶片的长度,r1为定子的外半径。
弹性叶片产生形变的弯曲角通过以下公式计算得到:
式(2)中,β为弹性叶片的安装倾斜角。
使弹性叶片产生形变的压力Fp,通过以下公式计算得到:
式(3)中,E为弹性叶片的弹性模量,b为弹性叶片的宽度,h为弹性叶片的厚度。
弹性叶片的切向速度幅值Vm,通过以下公式计算得到:
式(4)中,A为定子一阶径向振动的振幅,T为驱动电源的周期,β为弹性叶片的安装倾斜角,α为弹性叶片形变后与定子径向的夹角,为弹性叶片产生形变的弯曲角。
步骤3:设转子转速N的值为0。
步骤4:计算转子的切向线速度Vr,通过以下公式计算得到:
式(5)中,r为转子的内半径,N为转子转速的值。
步骤5:计算等速时间点t1,通过以下公式计算得到:
式(6)中,T为驱动电源的周期,Vm为弹性叶片的切向速度幅值,Vr为转子的切向线速度。
步骤6:计算与转子转速N相对应的电机输出转矩Tload,通过以下公式计算得到:
式(7)中,n为弹性叶片的数量,r为转子的内半径,T为驱动电源的周期,μd为弹性叶片与转子间的库伦摩擦系数,Fp为使弹性叶片产生形变的压力,t1为等速时间点,α为弹性叶片形变后与定子径向的夹角,fv为黏性摩擦系数,Vm为弹性叶片的切向速度幅值,Vr为转子的切向线速度。
步骤7:判断与转子转速N相对应的电机输出转矩Tload的大小,若与转子转速N相对应的电机输出转矩Tload>0,重置转子转速N的值并计算与转子转速N相对应的电机输出转矩Tload,若与转子转速N相对应的电机输出转矩Tload≤0,结束计算。
重置转子转速N的值并计算与转子转速N相对应的电机输出转矩Tload包括:
将转子转速N的值加上步进常数Q,得到重置后的转子转速N的值;
根据径向驻波型超声波电机的基本参数和重置后的转子转速N的值,计算与重置后的转子转速N相对应的电机输出转矩Tload。
具体地,步进常数Q的取值范围为0<Q≤1。
步骤8:根据计算过程中所有转子转速N以及与转子转速N相对应的电机输出转矩Tload,确定径向驻波型超声波电机的转矩转速特性。
实施例二:
本实施例是基于实施例一提供的基于库伦摩擦和黏性摩擦的径向驻波型超声波电机转矩转速特性的计算方法的具体应用场景,以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
以一台外直径为32mm的径向驻波型超声波电机为例,结合图1对各步骤进行详细的描述。
步骤1:获取径向驻波型超声波电机的基本参数。
该径向驻波型超声波电机的基本参数为:驱动电源的周期T=1.3514×10-5s、定子一阶径向振动的振幅A=0.6μm、弹性叶片的长度l=5mm、弹性叶片的安装倾斜角β=27°、弹性叶片的数量n=16、弹性叶片的弹性模量E=2.11×1011N/m2、弹性叶片的宽度b=5mm、弹性叶片的厚度h=0.1mm、弹性叶片与转子间的库伦摩擦系数μd=0.32、黏性摩擦系数fv=2.0、定子的外半径r1=11mm、转子的内半径r=15mm。
步骤2:计算弹性叶片的形变参数和运动参数。
计算得到弹性叶片的形变参数:弹性叶片形变后与定子径向的夹角α=30.6834°、弹性叶片产生形变的弯曲角使弹性叶片产生形变的压力Fp=3.2447N;
计算得到弹性叶片的运动参数:弹性叶片的切向速度幅值Vm=0.1022m/s。
步骤3:设转子转速N的值为0r/min。
步骤4:计算得到转子的切向线速度Vr=0.0016×Nm/s。
步骤5:计算等速时间点t1:当转子的转速N=0r/min时,等速时间点t1=3.3784×10-6s。
步骤6:计算与转子转速N相对应的电机输出转矩Tload:当转子转速N=0r/min时,与转子转速N相对应的电机输出转矩Tload=0.3951Nm。
步骤7:判断与转子转速N相对应的电机输出转矩Tload的大小,在本实施例中电机输出转矩Tload>0,重置转子转速N的值并计算与转子转速N相对应的电机输出转矩Tload,直到与转子转速N相对应的电机输出转矩Tload≤0,结束计算。
重置转子转速N的值包括:
将转子转速N的值加上步进常数Q,在本实施例中,Q=1.0r/min,得到重置后的转子转速N的值为1.0r/min;
根据径向驻波型超声波电机的基本参数和重置后的转子转速N的值,计算与重置后的转子转速N相对应的电机输出转矩Tload=0.3873Nm。
具体地如下表所示,重置转子转速N的值后,若Tload>0,则继续重置转子转速N的值,转子转速N的值依次取0、1.0、2.0、3.0……直到Tload≤0,停止计算。
表1转子转速N与其相对应的电机输出转矩Tload
步骤8:根据计算过程中所有转子转速N以及与转子转速N相对应的电机输出转矩Tload,确定径向驻波型超声波电机的转矩转速特性。
如图2所示,通过在MATLAB编写的基于库伦摩擦和黏性摩擦的径向驻波型超声波电机转矩转速特性的计算方法的程序中手动输入径向驻波型超声波电机的基本参数,可以画出经过循环计算所得出的基于库伦摩擦+黏性摩擦的径向驻波型超声波电机转矩转速特性曲线。
实施例三:
本发明实施例提供基于库伦摩擦和黏性摩擦的径向驻波型超声波电机转矩转速特性的计算装置,包括如下模块:
计算模块:用于根据径向驻波型超声波电机的基本参数和预设的转子转速N的值,计算与转子转速N相对应的电机输出转矩Tload;
判断模块:用于判断与转子转速N相对应的电机输出转矩Tload的大小,若与转子转速N相对应的电机输出转矩Tload>0,重置转子转速N的值并重新计算与转子转速N相对应的电机输出转矩Tload,若与转子转速N相对应的电机输出转矩Tload≤0,结束计算;
转矩转速特性确认模块:用于根据计算过程中所有转子转速N以及与转子转速N相对应的电机输出转矩Tload,确定径向驻波型超声波电机的转矩转速特性。
实施例四:
本发明实施例还提供基于库伦摩擦和黏性摩擦的径向驻波型超声波电机转矩转速特性的计算装置,包括处理器及存储介质;
所述存储介质用于存储指令;
所述处理器用于根据所述指令进行操作以执行实施例一所述方法的步骤。
实施例五:
本发明实施例还提供一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现实施例一所述方法的步骤。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (9)
1.径向驻波型超声波电机转矩转速特性的计算方法,其特征在于,包括:
确定径向驻波型超声波电机的基本参数,基于径向驻波型超声波电机的基本参数和预设的转子转速N的值,计算与转子转速N相对应的电机输出转矩Tload;其中,与转子转速N相对应的电机输出转矩Tload通过以下公式计算得到:
式(1)中,n为弹性叶片的数量,r为转子的内半径,T为驱动电源的周期,μd为弹性叶片与转子间的库伦摩擦系数,Fp为使弹性叶片产生形变的压力,t1为等速时间点,α为弹性叶片形变后与定子径向的夹角,fv为黏性摩擦系数,Vm为弹性叶片的切向速度幅值,Vr为转子的切向线速度;
判断与转子转速N相对应的电机输出转矩Tload的大小,若与转子转速N相对应的电机输出转矩Tload>0,重置转子转速N的值并计算与转子转速N相对应的电机输出转矩Tload,若与转子转速N相对应的电机输出转矩Tload≤0,结束计算;
根据计算过程中所有转子转速N及其相对应的电机输出转矩Tload,确定径向驻波型超声波电机的转矩转速特性;
所述径向驻波型超声波电机的基本参数包括弹性叶片与转子间的库伦摩擦系数μd和黏性摩擦系数fv。
2.根据权利要求1所述的径向驻波型超声波电机转矩转速特性的计算方法,其特征在于,等速时间点t1通过以下公式计算得到:
式(2)中,T为驱动电源的周期,Vm为弹性叶片的切向速度幅值,Vr为转子的切向线速度,通过以下公式计算得到:
式(3)中,r为转子的内半径,N为转子转速的值。
3.根据权利要求2所述的径向驻波型超声波电机转矩转速特性的计算方法,其特征在于,弹性叶片的切向速度幅值Vm为弹性叶片的运动参数,通过以下公式计算得到:
式(4)中,A为定子一阶径向振动的振幅,T为驱动电源的周期,β为弹性叶片的安装倾斜角,α为弹性叶片形变后与定子径向的夹角,为弹性叶片产生形变的弯曲角。
4.根据权利要求3所述的径向驻波型超声波电机转矩转速特性的计算方法,其特征在于,弹性叶片形变后与定子径向的夹角α、弹性叶片产生形变的弯曲角和使弹性叶片产生形变的压力Fp为弹性叶片的形变参数;
弹性叶片形变后与定子径向的夹角α,通过以下公式计算得到:
式(5)中,r为转子的内半径,l为弹性叶片的长度,r1为定子的外半径;
弹性叶片产生形变的弯曲角通过以下公式计算得到:
式(6)中,β为弹性叶片的安装倾斜角;
使弹性叶片产生形变的压力Fp,通过以下公式计算得到:
式(7)中,E为弹性叶片的弹性模量,b为弹性叶片的宽度,h为弹性叶片的厚度。
5.根据权利要求1所述的径向驻波型超声波电机转矩转速特性的计算方法,其特征在于,重置转子转速N的值并计算与转子转速N相对应的电机输出转矩Tload包括:
将转子转速N的值加上步进常数Q,得到重置后的转子转速N的值;
根据径向驻波型超声波电机的基本参数和重置后的转子转速N的值,计算与重置后的转子转速N相对应的电机输出转矩Tload。
6.根据权利要求1所述的径向驻波型超声波电机转矩转速特性的计算方法,其特征在于,驱动电源的周期T、定子一阶径向振动的振幅A、弹性叶片的长度l、弹性叶片的安装倾斜角β、弹性叶片的数量n、弹性叶片的弹性模量E、弹性叶片的宽度b、弹性叶片的厚度h、定子的外半径r1和转子的内半径r为所述径向驻波型超声波电机的基本参数。
7.径向驻波型超声波电机转矩转速特性的计算装置,其特征在于,包括:
计算模块:用于根据径向驻波型超声波电机的基本参数和转子转速N的值,计算与转子转速N相对应的电机输出转矩Tload;其中,与转子转速N相对应的电机输出转矩Tload通过以下公式计算得到:
式(1)中,n为弹性叶片的数量,r为转子的内半径,T为驱动电源的周期,μd为弹性叶片与转子间的库伦摩擦系数,Fp为使弹性叶片产生形变的压力,t1为等速时间点,α为弹性叶片形变后与定子径向的夹角,fv为黏性摩擦系数,Vm为弹性叶片的切向速度幅值,Vr为转子的切向线速度;
判断模块:用于判断与转子转速N相对应的电机输出转矩Tload的大小,若与转子转速N相对应的电机输出转矩Tload>0,重置转子转速N的值并重新计算与转子转速N相对应的电机输出转矩Tload,若与转子转速N相对应的电机输出转矩Tload≤0,结束计算;
转矩转速特性确认模块:用于根据计算过程中所有转子转速N以及与转子转速N相对应的电机输出转矩Tload,确定径向驻波型超声波电机的转矩转速特性。
8.径向驻波型超声波电机转矩转速特性的计算装置,其特征在于,包括处理器及存储介质;
所述存储介质用于存储指令;
所述处理器用于根据所述指令进行操作以执行根据权利要求1~6任一项所述方法的步骤。
9.计算机可读存储介质,其上存储有计算机程序,其特征在于,该程序被处理器执行时实现权利要求1~6任一项所述方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011313598.1A CN112464400B (zh) | 2020-11-20 | 2020-11-20 | 基于库伦摩擦和黏性摩擦的径向驻波型超声波电机转矩转速特性的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011313598.1A CN112464400B (zh) | 2020-11-20 | 2020-11-20 | 基于库伦摩擦和黏性摩擦的径向驻波型超声波电机转矩转速特性的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112464400A CN112464400A (zh) | 2021-03-09 |
CN112464400B true CN112464400B (zh) | 2024-02-13 |
Family
ID=74799486
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011313598.1A Active CN112464400B (zh) | 2020-11-20 | 2020-11-20 | 基于库伦摩擦和黏性摩擦的径向驻波型超声波电机转矩转速特性的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112464400B (zh) |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07333084A (ja) * | 1994-06-03 | 1995-12-22 | Yaskawa Electric Corp | 制御対象の特性を表す関数の定数の同定・補償方法 |
JPH0815058A (ja) * | 1994-06-30 | 1996-01-19 | Yaskawa Electric Corp | 電動機駆動系の負荷定数測定方法 |
CN102255554A (zh) * | 2011-06-27 | 2011-11-23 | 上海师范大学 | 一种旋转行波超声波电机的摩擦补偿的速度控制方法 |
CN106602960A (zh) * | 2015-10-20 | 2017-04-26 | 上海三菱电梯有限公司 | 电机转动惯量辨识方法和装置及电机控制装置 |
CN106612087A (zh) * | 2015-10-20 | 2017-05-03 | 上海三菱电梯有限公司 | 电机转矩系数辨识方法和装置及电机控制装置和判断装置 |
CN106672717A (zh) * | 2015-11-11 | 2017-05-17 | 上海三菱电梯有限公司 | 电梯载荷转矩电流校正方法及装置和电梯控制方法 |
CN108027651A (zh) * | 2015-09-11 | 2018-05-11 | 莱恩斯特里姆技术有限公司 | 自动估计机械系统中的惯性、库伦摩擦和粘性摩擦的方法 |
CN110208691A (zh) * | 2018-02-28 | 2019-09-06 | 株式会社安川电机 | 一种电机负荷参数检测方法及设备 |
CN111177941A (zh) * | 2020-01-03 | 2020-05-19 | 成都卡诺普自动化控制技术有限公司 | 一种机器人摩擦力辨识方法 |
CN111684706A (zh) * | 2018-11-12 | 2020-09-18 | 深圳配天智能技术研究院有限公司 | 伺服系统阻力特性获取方法、控制方法和系统及存储装置 |
CN111817633A (zh) * | 2020-06-05 | 2020-10-23 | 浙江工业大学 | 一种永磁同步电机的机械参数辨识方法 |
-
2020
- 2020-11-20 CN CN202011313598.1A patent/CN112464400B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07333084A (ja) * | 1994-06-03 | 1995-12-22 | Yaskawa Electric Corp | 制御対象の特性を表す関数の定数の同定・補償方法 |
JPH0815058A (ja) * | 1994-06-30 | 1996-01-19 | Yaskawa Electric Corp | 電動機駆動系の負荷定数測定方法 |
CN102255554A (zh) * | 2011-06-27 | 2011-11-23 | 上海师范大学 | 一种旋转行波超声波电机的摩擦补偿的速度控制方法 |
CN108027651A (zh) * | 2015-09-11 | 2018-05-11 | 莱恩斯特里姆技术有限公司 | 自动估计机械系统中的惯性、库伦摩擦和粘性摩擦的方法 |
CN106602960A (zh) * | 2015-10-20 | 2017-04-26 | 上海三菱电梯有限公司 | 电机转动惯量辨识方法和装置及电机控制装置 |
CN106612087A (zh) * | 2015-10-20 | 2017-05-03 | 上海三菱电梯有限公司 | 电机转矩系数辨识方法和装置及电机控制装置和判断装置 |
CN106672717A (zh) * | 2015-11-11 | 2017-05-17 | 上海三菱电梯有限公司 | 电梯载荷转矩电流校正方法及装置和电梯控制方法 |
CN110208691A (zh) * | 2018-02-28 | 2019-09-06 | 株式会社安川电机 | 一种电机负荷参数检测方法及设备 |
CN111684706A (zh) * | 2018-11-12 | 2020-09-18 | 深圳配天智能技术研究院有限公司 | 伺服系统阻力特性获取方法、控制方法和系统及存储装置 |
CN111177941A (zh) * | 2020-01-03 | 2020-05-19 | 成都卡诺普自动化控制技术有限公司 | 一种机器人摩擦力辨识方法 |
CN111817633A (zh) * | 2020-06-05 | 2020-10-23 | 浙江工业大学 | 一种永磁同步电机的机械参数辨识方法 |
Non-Patent Citations (3)
Title |
---|
A single-modal linear ultrasonic motor based on multi vibration modes of PZT ceramics;Danhong Lu等;《Ultrasonics》;第107卷;1-9 * |
Control design for ultrasonic motors with dynamic friction interface;Carlos Canudas-de-Wit;《IFAC Proceedings Volumes》;第32卷(第2期);970-974 * |
数控机床变频主传动系统的空载能量参数特性研究;胡韶华等;《计算机集成制造系统》;第28卷(第2期);1-6 * |
Also Published As
Publication number | Publication date |
---|---|
CN112464400A (zh) | 2021-03-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2018205490A1 (zh) | 控制机器人转动关节运动的方法和装置及机器人 | |
CN102904512A (zh) | 一种步进电机的控制系统及控制方法 | |
CN107449533B (zh) | 一种谐波减速器动态扭转刚度分析装置与方法 | |
CN101080623A (zh) | 用于识别回转轴的受损轴承的方法和设备 | |
CN102255554B (zh) | 一种旋转行波超声波电机的摩擦补偿的速度控制方法 | |
CN107947629A (zh) | 双压电叠堆驱动型精密旋转驱动装置 | |
CN112464400B (zh) | 基于库伦摩擦和黏性摩擦的径向驻波型超声波电机转矩转速特性的计算方法 | |
CN107873122B (zh) | 电动机控制系统 | |
CN112464458B (zh) | 考虑摩擦系数修正的径向驻波型超声波电机转矩转速特性的计算方法 | |
CN107345802A (zh) | 一种行星减速器间隙动态特性的分析装置与方法 | |
CN102906991A (zh) | 无传感器扭转模态阻尼系统和方法 | |
CN105406761B (zh) | 基于输入电压非恒等幅值的转速控制系统及方法 | |
Hu et al. | Research on prediction method of transmission accuracy of harmonic drive | |
CN114374338A (zh) | 基于定子齿结构和动态摩擦力的环形行波超声波电机转矩转速特性的计算方法 | |
CN114374337A (zh) | 考虑定子齿结构的环形行波型超声波电机转矩转速特性的计算方法 | |
CN110726511A (zh) | 一种带旋转机构的惯导系统的质心计算方法 | |
CN114374340A (zh) | 基于动态摩擦力的环形行波超声波电机转矩转速特性的计算方法 | |
CN112395711B (zh) | 一种内啮合齿轮副六自由度动力学模型建模方法 | |
CN114861338A (zh) | 风力发电机组参数的确定方法、装置及电子设备 | |
CN112100849B (zh) | 改进的径向驻波型超声波电机转矩转速特性计算方法 | |
JP5126044B2 (ja) | 車両の制御装置 | |
CN114374339A (zh) | 环形行波型超声波电机转矩转速特性的计算方法及系统 | |
CN111188793B (zh) | 一种离心压气机叶轮分流叶片周向角设计方法及叶轮 | |
CN106787940A (zh) | 一种改进的超声波电机反步自适应伺服控制方法 | |
CN108282107B (zh) | 一种径向驻波型超声波电机转矩转速特性的计算方法 |
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 |