CN104408241B - 一种修形圆柱齿轮的有限元网格自动生成方法 - Google Patents
一种修形圆柱齿轮的有限元网格自动生成方法 Download PDFInfo
- Publication number
- CN104408241B CN104408241B CN201410645189.XA CN201410645189A CN104408241B CN 104408241 B CN104408241 B CN 104408241B CN 201410645189 A CN201410645189 A CN 201410645189A CN 104408241 B CN104408241 B CN 104408241B
- Authority
- CN
- China
- Prior art keywords
- gear
- point
- tooth
- grid
- finite element
- 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.)
- Expired - Fee Related
Links
Landscapes
- Gears, Cams (AREA)
Abstract
本发明涉及一种修形圆柱齿轮的有限元网格自动生成方法,用以解决修形圆柱齿轮仿真分析中有限元网格模型建立操作繁杂的问题。该方法基于齿轮修形加工的基本原理,通过计算出修形齿轮的齿廓曲线点以及内部点,形成齿轮有限元模型的网格节点坐标,进而建立起齿轮的有限元网格模型。本发明建立了修形齿轮有限元网格自动生成的方法,可以免去三维建模以及划分网格的繁琐工作,在很短的时间内建立起齿轮有限元模型,大大提高齿轮有限元计算分析的效率和精度。
Description
技术领域
本发明涉及齿轮仿真的技术领域,具体涉及一种修形圆柱齿轮的有限元网格自动生成方法,适用于齿轮仿真分析中有限元网格自动生成。
背景技术
目前国内外对于齿轮有限元仿真分析的前处理,一般是先通过三维建模软件建立齿轮的三维几何模型,然后将齿轮模型导入有限元分析软件中进行网格划分,最后根据仿真工况条件定义载荷及边界条件。为了保证有限元仿真分析的计算精度并顾及计算效率,需要建立网格分布合理的齿轮网格模型,在有限元软件中进行齿轮的划分网格操作费时费力,需要进行线、面、体网格划分定义多个步骤,迫切需要一种能够根据给定基本参数自动生成修形齿轮网格模型的方法。
发明内容
本发明要解决的技术问题为:通过给定齿轮的基本参数、修形参数、网格划分参数,便可以自动生成齿轮的三维有限元网格模型,免去了在三维建模软件中的建模及在有限元软件中网格划分的复杂操作流程,大大提高了仿真分析前处理的效率。
本发明采用的技术方案是:一种修形圆柱齿轮的有限元网格自动生成方法,其特征在于实现步骤如下:
步骤(1):已知齿轮基本参数、齿廓修形参数、齿向修形参数、网格参数四组参数,基本参数包括:齿数z、模数m、压力角α、螺旋角β、齿宽B、变位系数x、齿轮安装中心距a、法向侧隙jbn;齿廓修形参数包括:齿顶修形量Δa、齿顶修形高度la、齿顶修形指数ba,齿根修形量Δf、齿根修形高度lf、齿根修形指数bf;齿向修形参数包括:左侧修形量Δl、左侧修形长度ll、左侧修形指数bl,右侧修形量Δr、右侧修形长度lr、右侧修形指数br;网格参数包括:齿廓上部网格数n1、齿廓中部网格数n2、齿廓下部网格数n3、齿根圆角网格数n4、齿根底部网格数n5、齿厚网格数n6、齿底网格数n7、齿向左侧网格数n8、齿向中部网格数n9、齿向右侧网格数n10(网格参数的划分区间位置如附图2所示);
步骤(2):根据步骤(1)中的基本参数计算中间参数,齿顶圆半径,ra=rp+ha,齿根圆半径rf=rp-hf,齿顶高ha=(ha *+x-Δy)·m,齿根高hf=(ha *+c*-x)·m,其中的Δy=x1+x2-y,a’为实际中心距,a为理论中心距齿高系数ha *=1,顶隙系数c*=0.25;
步骤(3):根据齿向网格数将齿轮沿齿向方向分割为n个截面n=n8+n9+n10+1,分别求取n个截面上的网格点,每个截面的位置由lx表示,齿向修形量随着lx的变化而改变,计算截面对应的齿向修形量Δ:当时,当时,Δ=0;当时,
步骤(4):由步骤(3)中的齿向修形量Δ计算其中齿轮分度圆半径端面模数
步骤(5):根据步骤(4)得到的rpa计算齿廓曲线计算公式为其中:
——齿轮转角,决定齿廓上点的位置
ξ——相邻两齿间夹角的一半
与齿轮转角相关的刀具点计算方程:
式中
rp——齿轮分度圆半径,
m——齿轮模数
αt——端面压力角,
st——刀具端面上齿距的一半,
mt——端面模数,
齿廓上距齿轮圆心半径为r的任一点对应的及刀具上点的坐标(x,y):
式中:
r——齿廓上任一点距齿轮中心的距离
rpa——砂轮中心与齿轮圆心的距离,控制齿向修形量的大小
xn——齿轮法向变位系数
Δ为齿廓上距齿轮圆心半径为r的任一点对应的齿廓修形量,当r∈[ra,ra-la]时,当r∈[ra-la,rf+lf]时,Δ=0;当r∈[rf+lf,rf]时,
步骤(6):计算齿根过渡曲线,根据基本参数中的顶隙和压力角计算刀具圆角半径顶隙c=c*·m,计算刀具圆角中心坐标根据刀具圆角中心坐标求得刀具圆角上任意一点通过坐标变换将刀具圆角坐标转变为齿根过渡曲线坐标:经过步骤(5)(6)后可以得到齿廓曲线各点(如附图3a所示);
步骤(7):根据步骤(2)中的中间参数计算内部分割边界点M1(racosξ,rasinξ),M2(rfcosξ,rfsinξ),M3(rdcosξ,rdsinξ), (各分割点位置如图2所示,计算完成后得到各边界点,如图3b所示);
步骤(8):根据步骤(5)计算的齿廓曲线步骤(6)计算的齿根过渡曲线步骤(7)计算的内部分割边界点,按步骤(1)输入的网格划分数计算出截面分割边界点(如附图3c所示),将分割边界点之间进行等分,得到包含内部网格点的lx=0截面位置的网格节点(如附图3d所示);
步骤(9):将步骤(8)得到的截面网格节点经过旋转运动变换,得到当前lx对应的左侧截面网格节点计算公式为:
式中:
θ——螺旋运动角,每个截面对应一个θ,
p——螺旋参数,由齿轮螺旋角β决定,
步骤(10):根据步骤(9)得到的lx截面位置网格点经过对称变换得到lx截面的右侧网格节点计算公式为:
将与步骤(9)中得到的整合,得到lx截面对应的截面点(如附图3e所示);
步骤(11):重复步骤(4)~(10)直到所有lx对应的截面点均计算完成(如附图4);各截面网格点坐标计算完成后,通过相邻8个节点之间连线,构成由多个六面体网格单元组成的齿轮有限元模型(附图5、6),按照有限元元件前处理文件格式存储节点和单元信息,利用该电子文件,可以快速建立起齿轮的有限元网格模型。
本发明的原理:基于齿轮修形加工的基本原理,计算出修形齿轮的齿廓曲线点以及内部点,形成齿轮有限元模型的网格节点坐标,从而建立齿轮的有限元网格模型。
本发明与现有技术相比的有益效果是:目前修形齿轮建模操作复杂,有限元分析前处理操作复杂、专业技术要求高、效率低,本发明建立了修形齿轮有限元网格自动生成的方法,可以免去三维建模以及划分网格的繁琐工作,在很短的时间内建立起齿轮有限元模型,大大提高齿轮有限元计算分析的效率和精度。
附图说明
图1为本发明的方法流程图;
图2为齿轮截面划分示意图;
图3为齿轮界面点生成过程示意图;
图4为单齿节点计算结果示意图;
图5为单齿有限元网格模型示意图;
图6为完整齿轮有限元网格模型示意图。
具体实施方式
下面结合附图和具体实施例进一步说明本发明。
本发明是针对圆柱齿轮有限元建模而设计的快速建模方法,以模数为5的齿轮为例,建模流程图如图所示。本实施实例以表参数中的小轮为对象,建模的实现步骤如下:
步骤(1):给定基本参数
步骤(2):根据基本参数和安装参数,结合齿高系数、顶隙系数,计算齿轮的分度圆、齿顶圆、齿根圆半径。
ha *=1
c*=0.25
Δy=x1+x2-y=0.0281638
ha=(ha *+x-Δy)·m=6.2342
hf=(ha *+c*-x)·m=4.875
ra=rp+ha=58.1957
rf=rp-hf=47.0865
步骤(3):计算齿向修形量,对于中间截面,以lx=0为例,Δ=0;
对于左侧截面,以lx=-0.018m为例,
对于右侧截面,以lx=0.018m为例,
步骤(4):根据齿向修形量计算不同截面刀具中心与齿轮中心的距离:
Δ=0截面,
Δ=0.01截面,
步骤(5):根据齿廓上距齿轮圆心半径为r的任一点对应的计算齿廓点坐标, 以r=51.04248为例:
步骤(6):计算齿根过渡曲线,根据基本参数中的顶隙和压力角计算刀具圆角半径顶隙c=c*·m=1.25,计算刀具圆角中心坐标:
通过坐标变换将刀具圆角坐标转变为齿根过渡曲线坐标,计算公式为其中:
刀具圆角坐标
步骤(7):根据齿轮参数计算分割边界点坐标:
rd=ra-2(ha+hf)=35.977343
M1(racosξ,rasinξ)=(57.311582,10.105578)
M2(rfcosξ,rfsinξ)=(46.310974,8.165874)
M3(rdcosξ,rdsinξ)=(35.430766,6.247400)
上式中xm,ym的计算方法:当αc,取最大值时,利用步骤(7)中过渡曲线公式求得(xm,ym)=(47.025396,0.499455);
步骤(8):根据步骤(5)计算的齿廓曲线步骤(6)计算的齿根过渡曲线步骤(7)计算的内部分割边界点,按步骤(1)输入的网格划分数计算出各点之间的中间点,得到lx=0截面位置的网格点
步骤(9):将步骤(8)得到的截面网格点经过旋转运动变换,得到当前lx对应的左侧截面网格点计算公式为:
式中:
θ——螺旋运动角,每个截面对应一个θ,
p——螺旋参数,由齿轮螺旋角β决定,
步骤(10):根据步骤(9)得到的lx截面位置网格点经过对称变换得到lx截面的右侧网格点计算公式为:
将与步骤(9)中得到的整合变得到lx截面对应的截面点
步骤(11):根据齿向各截面的rpa,重复步骤(4)~(10),计算出其他截面点,获得一个齿的所有节点,将相邻八个节点连接成为单元,建立齿轮的有限元网格模型。
综上,通过以上流程,可以建立建立齿轮的有限元网格模型,将本流程编写成计算程序,可实现齿轮有限元网格模型的快速自动建立,极大简化齿轮有限元分析前处理所需要的建模及网格划分工作。本发明适用于圆柱齿轮有限元分析前处理过程中的有限元网格划分,为齿轮啮合性能有限元分析提供了一种快速的建模方法。
本发明未详细阐述部分属于本领域公知技术。凡采用等同变换或者等效替换而形成的技术方案,均落在本发明权利保护范围之内。
Claims (3)
1.一种修形圆柱齿轮的有限元网格自动生成方法,其特征在于实现步骤如下:
步骤(1):已知齿轮基本参数、齿廓修形参数、齿向修形参数、网格参数四组参数,齿轮基本参数包括:齿数z、模数m、压力角α、螺旋角β、齿宽B、变位系数x、齿轮安装中心距a、法向侧隙jbn;齿廓修形参数包括:齿顶修形量Δa、齿顶修形高度la、齿顶修形指数ba,齿根修形量Δf、齿根修形高度lf、齿根修形指数bf;齿向修形参数包括:左侧修形量Δl、左侧修形长度ll、左侧修形指数bl,右侧修形量Δr、右侧修形长度lr、右侧修形指数br;网格参数包括:齿廓上部网格数n1、齿廓中部网格数n2、齿廓下部网格数n3、齿根圆角网格数n4、齿根底部网格数n5、齿厚网格数n6、齿底网格数n7、齿向左侧网格数n8、齿向中部网格数n9、齿向右侧网格数n10;
步骤(2):根据步骤(1)中的齿轮基本参数计算中间参数,中间参数包括齿顶圆半径ra、齿根圆半径rf、齿顶高ha、齿根高hf;
步骤(3):根据齿向网格数将齿轮沿齿向方向分割为n个截面n=n8+n9+n10+1,分别求取每个截面的网格点,每个截面的位置由lx表示,
步骤(4):由步骤(3)中的齿向修形量计算其中rp为齿轮分度圆半径,rpa为刀具轴线与齿轮轴线的距离,x为变位系数,Δ为左侧修形量Δl或右侧修形量Δr,mt为端面模数;
步骤(5):根据步骤(4)得到的rpa按公式计算齿廓曲线其中:
为齿轮转角,决定齿廓上点的位置,ξ为相邻两齿间夹角的一半
m为齿轮模数,αt为端面压力角,st为刀具端面上齿距的一半z为齿轮基本参数中的齿数;
步骤(6):计算齿根过渡曲线根据刀具圆角半径顶隙c=c*·m,计算刀具圆角中心坐标在此基础上,计算刀具圆角上任意一点齿高系数ha *=1,顶隙系数c*=0.25,αc刀具圆角坐标,R(x,y)为刀具圆角中心坐标;
步骤(7):根据步骤(2)中的中间参数计算内部分割边界点M1(racosξ,rasinξ),M2(rfcosξ,rfsinξ),M3(rdcosξ,rd sinξ), 其中,rd=ra-2(ha+hf),及分别为相应内部分割边界点的相位角,M1,M2,M3,C,E,D,I,H,G为分割边界点及其坐标;
步骤(8):根据步骤(5)计算的齿廓曲线步骤(6)计算的齿根过渡曲线步骤(7)计算的内部分割边界点,按步骤(1)输入的网格划分数计算出各点之间的中间点,得到lx=0截面位置的网格点
步骤(9):将步骤(8)得到的lx=0截面位置网格点经过旋转运动变换,得到当前lx对应的左侧截面位置网格点计算公式为:
式中:
θ为螺旋运动角;
步骤(10):根据步骤(9)得到的当前lx对应的左侧截面位置网格点经过对称变换得到lx截面的右侧网格点计算公式为:
式中:
将与步骤(9)中得到的整合,得到lx截面对应的网格点
步骤(11):重复步骤(4)~步骤(10)直到所有lx对应的截面点均计算完成,各截面网格点坐标计算完成后,将相邻8个节点定义为网格单元,按照有限元元件前处理文件格式存储节点和单元信息,利用本文件,即可快速建立起齿轮的有限元网格模型。
2.根据权利要求1所述的一种修形圆柱齿轮的有限元网格自动生成方法,其特征在于:所述步骤(5)中齿廓曲线上点的位置由齿轮转角表示,其计算方法为:
其中x,y的计算方法如下:
r为所求齿廓点距齿轮中心的距离,rpa为刀具轴线与齿轮轴线的距离,控制齿向修形量的大小,xn为齿轮法向变位系数。
3.根据权利要求1所述的一种修形圆柱齿轮的有限元网格自动生成方法,其特征在于:所述步骤(5)中的齿向修形量△随着lx的变化而改变,计算方法如下:当时,当时,△=0;当时,b为齿向修形参数中的左侧修形指数bl或右侧修形指数br。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410645189.XA CN104408241B (zh) | 2014-11-06 | 2014-11-06 | 一种修形圆柱齿轮的有限元网格自动生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410645189.XA CN104408241B (zh) | 2014-11-06 | 2014-11-06 | 一种修形圆柱齿轮的有限元网格自动生成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104408241A CN104408241A (zh) | 2015-03-11 |
CN104408241B true CN104408241B (zh) | 2017-05-10 |
Family
ID=52645872
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410645189.XA Expired - Fee Related CN104408241B (zh) | 2014-11-06 | 2014-11-06 | 一种修形圆柱齿轮的有限元网格自动生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104408241B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104992027B (zh) * | 2015-07-17 | 2018-05-01 | 中国船舶重工集团公司第七�三研究所 | 人字齿轮有限元网格自动化建模的方法 |
CN106704543B (zh) * | 2016-12-05 | 2019-04-02 | 贵州群建精密机械有限公司 | 一种通过对发动机正时齿轮修形来降低发动机噪音的方法 |
CN107133417B (zh) * | 2017-05-25 | 2020-04-10 | 扬州大学 | 圆弧圆柱蜗杆的齿形检测方法 |
CN108488353B (zh) * | 2018-03-06 | 2021-06-04 | 上汽通用汽车有限公司 | 用于直齿锥齿轮的修形方法 |
CN110750922A (zh) * | 2019-09-11 | 2020-02-04 | 中南大学 | 含腹板结构的直齿轮有限元接触模型快速建模方法 |
CN111161416B (zh) * | 2019-12-11 | 2023-08-29 | 北京互时科技股份有限公司 | 根据模型形状信息精确调整模型显示优先级的方法和系统 |
CN114818181B (zh) * | 2022-04-21 | 2023-08-25 | 哈尔滨工程大学 | 基于齿廓线直齿圆柱齿轮有限元网格自动生成方法和设备 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102567596A (zh) * | 2012-03-05 | 2012-07-11 | 西北工业大学 | 一种面齿轮轮齿参数化建模方法 |
CN104318017A (zh) * | 2014-10-22 | 2015-01-28 | 江苏理工学院 | 一种非对称圆柱直齿轮副的建模方法 |
-
2014
- 2014-11-06 CN CN201410645189.XA patent/CN104408241B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102567596A (zh) * | 2012-03-05 | 2012-07-11 | 西北工业大学 | 一种面齿轮轮齿参数化建模方法 |
CN104318017A (zh) * | 2014-10-22 | 2015-01-28 | 江苏理工学院 | 一种非对称圆柱直齿轮副的建模方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104408241A (zh) | 2015-03-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104408241B (zh) | 一种修形圆柱齿轮的有限元网格自动生成方法 | |
CN102567596B (zh) | 一种面齿轮轮齿参数化建模方法 | |
CN103034751B (zh) | 一种便于实现参数化的渐开线直齿锥齿轮变曲率精确修鼓方法 | |
CN103197605B (zh) | 基于标准模板和表达式驱动的齿轮滚削cnc自动编程方法 | |
CN110008594B (zh) | 一种行星传动人字齿轮有限元网格自动建模与装配方法 | |
CN103678818A (zh) | 一种双圆弧修正摆线齿轮的精确建模方法 | |
CN103942396A (zh) | 一种含齿向误差的斜齿轮精确建模方法 | |
CN104217063A (zh) | 用于钢轨轧制孔型的仿真设计方法 | |
CN109027186B (zh) | 转向器变速比齿条的离散化包络设计方法 | |
CN102489785A (zh) | 一种基于离散数据螺旋锥齿轮的配对建模与加工方法 | |
CN105181327A (zh) | 一种裂纹轮齿啮合刚度计算的方法 | |
CN104317251A (zh) | 基于Obrechkoff算法的三次NURBS曲线实时插补方法 | |
CN104318017A (zh) | 一种非对称圆柱直齿轮副的建模方法 | |
CN106295015B (zh) | 一种渐开线直齿圆柱齿轮副的齿廓修形方法及与其配套的专用参数化cad系统 | |
CN102198543B (zh) | 齿轮建模方法和齿轮建模装置 | |
CN103065020A (zh) | 一种渐开线齿轮精指形铣刀cad 制图应用方法 | |
CN102581384B (zh) | 一种基于等切削面积的插齿加工方法 | |
CN103927409A (zh) | 基于有限元分析的混凝土烟囱应力分析系统 | |
CN102262696A (zh) | 一种用于直齿或斜齿的圆柱齿轮的建模方法 | |
CN103942397A (zh) | 基于幂函数的修形齿轮数字化建模方法 | |
CN105156634B (zh) | 渐开线直齿齿轮修形的优化方法 | |
CN105404738A (zh) | 一种用于指导齿轮应力检测的齿根应力分析方法 | |
CN115186493A (zh) | 一种小模数渐开线圆柱齿轮精确建模方法 | |
CN111005865B (zh) | 一种外啮合直齿轮泵困油面积精确测量方法 | |
CN110009746B (zh) | 反应堆燃料组件的带有边界层的六面体网格自动生成方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170510 Termination date: 20201106 |