[发明专利]一种修形圆柱齿轮的有限元网格自动生成方法有效

专利信息
申请号: 201410645189.X 申请日: 2014-11-06
公开(公告)号: CN104408241A 公开(公告)日: 2015-03-11
发明(设计)人: 王延忠;唐文;刘鹏 申请(专利权)人: 北京航空航天大学
主分类号: G06F17/50 分类号: G06F17/50;G06T17/30
代理公司: 北京科迪生专利代理有限责任公司 11251 代理人: 杨学明;顾炜
地址: 100191*** 国省代码: 北京;11
权利要求书: 查看更多 说明书: 查看更多
摘要: 发明涉及一种修形圆柱齿轮的有限元网格自动生成方法,用以解决修形圆柱齿轮仿真分析中有限元网格模型建立操作繁杂的问题。该方法基于齿轮修形加工的基本原理,通过计算出修形齿轮的齿廓曲线点以及内部点,形成齿轮有限元模型的网格节点坐标,进而建立起齿轮的有限元网格模型。本发明建立了修形齿轮有限元网格自动生成的方法,可以免去三维建模以及划分网格的繁琐工作,在很短的时间内建立起齿轮有限元模型,大大提高齿轮有限元计算分析的效率和精度。
搜索关键词: 一种 圆柱齿轮 有限元 网格 自动 生成 方法
【主权项】:
一种修形圆柱齿轮的有限元网格自动生成方法,其特征在于实现步骤如下:步骤(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为齿轮分度圆半径,mt为端面模数;步骤(5):根据步骤(4)得到的rpa按公式计算齿廓曲线其中:为齿轮转角,决定齿廓上点的位置,ξ为相邻两齿间夹角的一半rp为齿轮分度圆半径,m为齿轮模数,αt为端面压力角,st为刀具端面上齿距的一半mt为端面模数;步骤(6):计算齿根过渡曲线根据刀具圆角半径顶隙c=c*·m,计算刀具圆角中心坐标<mrow><msub><mi>R</mi><mi>o</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><mi>m</mi><mo>&CenterDot;</mo><msup><msub><mi>h</mi><mi>a</mi></msub><mo>*</mo></msup><mo>+</mo><mi>c</mi><mo>-</mo><msub><mi>r</mi><mi>c</mi></msub><mo>,</mo><mfrac><msub><mi>s</mi><mi>t</mi></msub><mn>2</mn></mfrac><mo>+</mo><mi>m</mi><mo>&CenterDot;</mo><msup><msub><mi>h</mi><mi>a</mi></msub><mo>*</mo></msup><mo>&CenterDot;</mo><mi>tan</mi><msub><mi>&alpha;</mi><mi>t</mi></msub><mo>+</mo><msub><mi>r</mi><mi>c</mi></msub><mo>&CenterDot;</mo><mi>cos</mi><msub><mi>&alpha;</mi><mi>t</mi></msub><mo>)</mo></mrow><mo>,</mo></mrow>在此基础上,计算刀具圆角上任意一点<mrow><msub><mi>R</mi><mi>r</mi></msub><mrow><mo>(</mo><msub><mi>&alpha;</mi><mi>c</mi></msub><mo>)</mo></mrow><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mo>-</mo><msub><mi>R</mi><mi>o</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>r</mi><mi>c</mi></msub><mo>&CenterDot;</mo><mi>sin</mi><msub><mi>&alpha;</mi><mi>c</mi></msub></mtd></mtr><mtr><mtd><mo>-</mo><msub><mi>R</mi><mi>o</mi></msub><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>r</mi><mi>c</mi></msub><mo>&CenterDot;</mo><mi>cos</mi><msub><mi>&alpha;</mi><mi>c</mi></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>步骤(7):根据步骤(2)中的中间参数计算内部分割边界点M1(racosξ,ra sinξ),M2(rfcosξ,rf sinξ),M3(rdcosξ,rd sinξ),步骤(8):根据步骤(5)计算的齿廓曲线步骤(6)计算的齿根过渡曲线步骤(7)计算的内部分割边界点,按步骤(1)输入的网格划分数计算出各点之间的中间点,得到lx=0截面位置的网格点步骤(9):将步骤(8)得到的截面网格点经过旋转运动变换,得到当前lx对应的左侧截面网格点计算公式为:式中:<mrow><msub><mi>M</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>cos</mi><mi>&theta;</mi></mtd><mtd><mo>-</mo><mi>sin</mi><mi>&theta;</mi></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mi>sin</mi><mi>&theta;</mi></mtd><mtd><mi>cos</mi><mi>&theta;</mi></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><msub><mi>l</mi><mi>x</mi></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced></mrow>步骤(10):根据步骤(9)得到的lx截面位置网格点经过对称变换得到lx截面的右侧网格点计算公式为:式中:<mrow><msub><mi>M</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>cos</mi><mrow><mo>(</mo><mn>2</mn><mi>&xi;</mi><mo>+</mo><mn>2</mn><mi>&theta;</mi><mo>)</mo></mrow></mtd><mtd><mi>sin</mi><mrow><mo>(</mo><mn>2</mn><mi>&xi;</mi><mo>+</mo><mn>2</mn><mi>&theta;</mi><mo>)</mo></mrow></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mi>sin</mi><mrow><mo>(</mo><mn>2</mn><mi>&xi;</mi><mo>+</mo><mn>2</mn><mi>&theta;</mi><mo>)</mo></mrow></mtd><mtd><mo>-</mo><mi>cos</mi><mrow><mo>(</mo><mn>2</mn><mi>&xi;</mi><mo>+</mo><mn>2</mn><mi>&theta;</mi><mo>)</mo></mrow></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced></mrow>与步骤(9)中得到的整合,得到lx截面对应的网格点步骤(11):重复步骤(4)~步骤(10)直到所有lx对应的截面点均计算完成,各截面网格点坐标计算完成后,将相邻8个节点定义为网格单元,按照有限元元件前处理文件格式存储节点和单元信息,利用本文件,即可快速建立起齿轮的有限元网格模型。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于北京航空航天大学,未经北京航空航天大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/patent/201410645189.X/,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top