CN102332046B - 一种齿轮裂纹扩展模拟的小波扩展有限元仿真分析方法 - Google Patents
一种齿轮裂纹扩展模拟的小波扩展有限元仿真分析方法 Download PDFInfo
- Publication number
- CN102332046B CN102332046B CN 201110297653 CN201110297653A CN102332046B CN 102332046 B CN102332046 B CN 102332046B CN 201110297653 CN201110297653 CN 201110297653 CN 201110297653 A CN201110297653 A CN 201110297653A CN 102332046 B CN102332046 B CN 102332046B
- Authority
- CN
- China
- Prior art keywords
- tooth root
- unit
- crackle
- crack
- node
- 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
Images
Abstract
本发明涉及针对裂纹扩展的有限元分析法,是一种采用了小波扩展这一新的有限分析单元的数值分析方法。该方法基于计算机辅助制图软件建立含有任意裂纹故障的齿轮啮合模型;将模型导入到有限元分析软件ABAQUS中,结合ABAQUS网格划分和数据输出功能获得有限元网格的几何数据;然后根据所得数据应用数学辅助计算软件编制程序计算小波扩展单元的单元刚度矩阵,进而根据所分析结构的网格排列集成结构的整体刚度矩阵;引入边界约束条件和载荷后,求解有限元方程,得到裂纹扩展的数值解。本发明不仅可以跟踪裂纹的生长状况,还能解决应力高度集中带来的困难,具有较高的计算精度和较高的计算效率,为机械设备故障诊断的研究带来方便。
Description
技术领域
本发明涉及针对裂纹扩展的有限元分析法,是一种采用了小波扩展这一新的有限分析单元的数值分析方法。
背景技术
由于齿轮裂纹尖端区域位移与应力都含有 奇异性(其中r为裂纹尖端附近某一点的极坐标矢径),在裂纹中奇异点附近的解梯度大,还会发生突变,因此给传统有限元方法求数值解造成困难。
传统有限元方法采用连续函数作为形函数,要求在单元内部形状函数连续且材料性能不能跳跃,对于处理像裂纹这样的不连续问题时,需要将裂纹面设置为单元的边,裂尖点设置为单元的结点,而且要在裂尖点附近不连续体的奇异场内进行高密度网格划分以及在模拟裂纹扩展时需要不断的进行网格的重新划分,这使得有限元程序设计相当复杂,且效率极低。在有限元框架内发展起来的扩展有限元法(XFEM),以解决不连续问题为着眼点,在处理裂纹问题是虽然不需要对结构内存在的几何或者物理界面进行网格剖分,但是对于工程中裂纹等奇异性问题的处理还有困难。将小波分析理论引入有限元分析的小波有限元方法因其具有的多尺度、多分辨的特性,可以解决因裂纹尖端区域奇异性带来的困难,但却不能反映裂纹的生长情况。
发明内容
本发明为了方便有效地解决齿根裂纹尖端应力应变奇异性给传统有限元带来的困难,同时得到裂纹的实时生长情况,提出了一种采取小波扩展这一 新型有限分析单元的数值分析方法。本发明将两种有限元方法的优势结合,借鉴小波形函数更加准确的逼近特性和扩展有限元无需重新划分网格的特点,形成一种新的有限元数值方法。本发明不仅可以跟踪裂纹的生长状况,还能解决应力高度集中带来的困难,具有较高的计算精度和较高的计算效率。此外采用此方法模拟齿根裂纹扩展能够为机械设备故障诊断的研究带来方便。
为实现上述目的,本发明的技术方案如下:
一种齿轮裂纹扩展模拟的小波扩展有限元仿真分析方法,该方法基于有限元软件ABAQUS和计算机辅助制图以及数学辅助计算软件构建的平台,对齿轮裂纹扩展进行分析,其特征在于包括以下步骤:
1)建立含齿根裂纹的齿根裂纹模型:
应用计算机辅助制图以及数学辅助计算软件绘制给定齿轮模数和齿数的齿轮模型,将所建立的齿轮模型导入到有限元软件ABAQUS中自动划分网格生成含有H个单元的齿轮模型,H为正整数,并假设在齿根区域内L个单元内存在裂纹,L为自然数,且L<H;
2)获得齿根裂纹模型内各个单元的节点坐标:
根据步骤1)中建立的齿根裂纹模型,应用有限元软件ABAQUS命令流输出功能输出齿根裂纹模型中H个单元的各个节点坐标;
3)根据所建立的齿根裂纹模型计算出H个单元的小波单元刚度矩阵:
将步骤2)齿根裂纹模型中所含有的H个单元的各个节点坐标数据导出,应用平面问题的势能泛函数得到小波单元刚度矩阵;
小波单元刚度矩阵计算式:
分块子矩阵为:
Ke,3=(Ke,2)T
其中a、b、c、d分别为平面单元求解域Ωe={x,y|x∈[a,b],y∈[c,d]}内x与y坐标取值的上下区间;x与y坐标值为步骤2)中所输出H个单元的各个节点坐标值;ξ,η为自然坐标,ξ=(x-a)/(b-a),η=(y-c)/(d-c),则标准求解域Ωs={ξ,η|ξ,η∈[0,1]};E为弹性模量;μ为泊松比;t为厚度;以上各分块子矩阵中Φ1和Φ2分别为m阶j尺度下的一维区间B样条小波尺度函数向量(详细请参看何正嘉 陈雪峰 著《小波有限元理论及其工程应用》科学出版社 第36、171页);
其中ξ,η为自然坐标,ξ=(x-a)/(b-a),η=(y-c)/(d-c),ξ,η取 值范围为齿根裂纹模型中含有齿根裂纹的L个单元内部区域的自然坐标值;s代表含有齿根裂纹的L个单元内各个节点编号;Ns(ξ,η)代表含有齿根裂纹的L个单元内部区域对应的小波形函数,Ns(ξ,η)=Φ(Re)-1,其中 Φ1和Φ2与步骤3)中相同,分别为m阶j尺度下的一维区间B样条小波尺度函数向量;Re为单位矩阵;us、vs表示含有齿根裂纹的L个单元内各个节点的常规自由度;as1和as2为含有齿根裂纹的L个单元内每一个节点的附加自由度;Hs(ξ,η)为反映含有齿根裂纹的L个单元的位移不连续性函数,计算公式如下:
5)应用扩展有限元改进齿根裂纹模型中齿根区域内含有齿根裂纹的L个单元的小波单元位移场函数之后,计算齿根区域内含有齿根裂纹的L个单元的小波扩展单元刚度矩阵(详细请参看:卡坦 著 韩来彬 译《MATLAB有限元分析与应用》清华大学出版社 第十三章),计算公式如下:
t为厚度;Be为应变矩阵(参考文献 杨珺 平面裂纹扩展分析的扩展有限 元法[D].南京航空航天大学航空宇航学院,2007.第31页);D为弹性矩阵, E为弹性模量;μ为泊松比;J为雅可比矩阵(详
细请参看冷纪桐 赵军 著《有限元技术基础》化学工业出版社 第84页);Ωs为标准求解域Ωs={ξ,η|ξ,η∈[0,1]};
6)计算齿根裂纹模型的结构刚度矩阵;
根据步骤3)、4)、5)对齿根裂纹模型中H个小波单元刚度矩阵的计算以及齿根区域内含有齿根裂纹的L个单元的小波扩展单元刚度矩阵的计算,按照齿轮裂纹模型中H个单元的各个节点次序集成结构刚度矩阵K(详细请参看胡于进王璋奇著《有限元分析及应用》清华大学出版社第59页)。
7)建立有限元方程求解节点常规自由度和附加自由度的数值解:
基于步骤6)中齿根裂纹模型结构刚度矩阵K的计算,建立以下有限元求解方程:
式中u、a分别为齿根裂纹模型中各个节点的常规自由度与附加自由度向量;Pu、Pa分别为与模型中常规自由度向量和附加自由度向量相对应的载荷向量,解方程可求出各个节点常规自由度和附加自由度的数值解;
8)计算齿根裂纹模型中裂纹的张开位移和裂纹开裂角度:
基于步骤7)求解有限元方程得到的齿根裂纹模型中附加自由度的数值解,可求齿根裂纹模型中裂纹张开位移,计算式如下:
w代表裂纹张开位移值; 为裂纹上任一点的法向量;s代表含有齿根裂纹的L个单元内各个节点编号;ns为齿根区域内含有齿根裂纹的L个单元包含的所有改进节点的数量;Ns(ξ,η)代表含有齿根裂纹的L个单元内部区域对应的小波形函数,计算裂纹张开位移时Ns(ξ,η)自变量取值为齿根裂纹面上,即当 时对应的数值;as为含有齿根裂纹的L个单元内每一个节点的附加自由度向量,每个该向量中含有步骤4)中所述as1和as2两个附加自由度值;
裂纹的裂纹开裂角度θ1可由下式计算:
式中σ1为单元的最大主应力,计算公式如下:
以上两式中σx、σy、τxy分别表示单元应力在横向、纵向以及切向的数值分量(详细请参看冷纪桐 赵军 著《有限元技术基础》化学工业出版社 第53页),根据步骤7)中求解有限元方程得到的齿根裂纹模型中各个节点常规自由度数值解计算;
以上所述为应用小波扩展有限元方法模拟齿根裂纹扩展的仿真方法,本发明的流程图如图1所示。
本发明的有益效果是:利用计算机辅助制图比如Pro/E在齿轮建模方面的优势,建立精确的齿轮啮合模型,避免了有限元软件在建模,特别是其在有 复杂裂纹故障的非正常齿轮建模时的缺憾;再结合ABAQUS强大的公式计算和数据输出功能,获得所划分网格的坐标数据;应用数学辅助计算软件如MATLAB编制程序计算出小波扩展这种单元的结构刚度矩阵,最后引入边界条件和载荷求解有限元方程,得到齿根裂纹扩展的模拟结果,还可根据需要改变差值函数的阶次获得不同的求解精度,这为齿轮系统故障的机理研究提供了准确、可靠的理论基础。
附图说明
图1是本发明的流程图;
图2是本实施例工作流程示意图;
图3是本实施例齿轮副二维实体模型示意图;
图4是本实施例含齿根裂纹轮齿二维实体模型示意图;
图5是本实施例在ABAQUS中划分网格后的齿轮有限元模型示意图;
图6是本实施例要改进的单元1与单元10的示意图;
图7是本实施例所中单元10对应的数值解;
图8是本实施例裂纹张开位移随裂纹长度变化的曲线;
图9是是本实施例模拟的裂纹扩展效果图;
具体实施方式
下面具体结合附图与实例对本发明作进一步的说明。
如图2所示,是本实施例工作流程示意图。主要由3大部分组成,详细的步骤如下。
一、利用Pro/E绘图软件绘制模数为2,齿数为19的齿轮与另一大齿轮啮合模型如图3所示。再根据需求可在轮齿上添加任意形式的裂纹故障,可 以改变裂纹的位置,大小,形态等特征,仿真出更符合实际情况的复杂裂纹故障。图4所示为假设一个在齿根处有图示裂纹故障的轮齿。
二、利用ABAQUS划分有限元网格并输出网格节点坐标:
1)将用Pro/E绘制的齿轮模型导入ABAQUS并划分生成含有22个有限单元的齿轮模型;假设裂纹故障存在于单元1与单元10中,为方便计算,实例中仅用到了含有裂纹故障的一个齿而且采用了较为粗糙的网格划分,利用有限元软件的命令显示各个单元编号和节点编号,如图5所示为划分网格后显示单元和节点编号的齿轮有限元模型;
2)在进行步骤1)的同时,ABAQUS会根据菜单操作的步骤自动生成对应的命令流文件,此命令流文件包含了菜单操作时的所有信息。本发明中,利用ABAQUS的计算功能和结果数据写入文本命令,即可得到整个结构中各个节点的几何坐标数据。
三、应用MATLAB计算所有22个单元的单元刚度矩阵以及齿根裂纹模型的结构刚度矩阵并建立有限元方程:
2)对含有齿根裂纹的两个小波单元位移场即单元1与单元10进行改进,如图6所示为需要改进的单元1与单元10的示意图,图中斜线表示裂纹,方形符号表示经过改进的裂纹贯穿单元节点,方程η=0.8ξ+0.1为裂纹在自然坐标下的表达式,(0,0.3)、(1,0.9)为裂纹与单元10的边界交点坐标,改进后 的单元位移场表达式为:
其中(ξ,η)代表单元1与单元10内部区域的自然坐标值;s代表单元1与单元10内各个节点编号;Hs(ξ,η)反映齿根裂纹的位移不连续性函数;改进的主要内容是在每一个原有节点常规自由度的基础上加入了一个附加自由度;含有齿根裂纹的单元1与单元10内每一个节点的附加自由度为as1和as2;
3)然后将扩展有限元对单元位移场的改进引入到单元刚度矩阵的积分方程 中得到相应的两个小波扩展单元刚度矩阵,t为厚度;Be为改进单元位移场后所得到的应变矩阵(参考文献 杨珺 平面裂纹扩展分析的扩展有限元法[D].南京航空航天大学航空宇航学院,2007.第31页),弹性矩阵 E为弹性模量;μ为泊松比;J为雅可比矩阵;Ωs={ξ,η|ξ,η∈[0,1]}为标准求解域。
4)计算完成所有22个单元的刚度矩阵之后,根据结构的单元和节点排列次序构造出整体的结构刚度矩阵K(详细请参看胡于进 王璋奇 著《有限元分析及应用》清华大学出版社 第59页)。
5)引入边界约束条件和载荷,单元1和单元10含有初始裂纹,限制齿轮两侧及底边节点的自由度,并在节点4施加载荷求解。基于本实施例步骤4)中齿根裂纹模型结构刚度矩阵K的计算,建立以下有限元求解方程:
式中u、a分别为齿根裂纹模型中各个节点的常规自由度与附加自由度向量:u=[u1 v1…ur vr]T a=[a11 a12…ar1 ar2]T,模型节点总数为r=33;u1 v1、ur vr分别代表模型中第一个节点常规自由度值和最后一个节点常规自由度值;a11 a12、ar1 ar2分别代表模型中第一个节点的附加自由度值和最后一个节点的附加自由度值;Pu、Pa为与模型中常规自由度向量和附加自由度向量相对应的载荷向量,利用MATLAB解以上方程可得到齿根裂纹模型中所有节点常规自由度和附加自由度的数值解;
6)根据5)中所求节点附加自由度数值解计算单元1和单元10中裂纹的张开位移,以单元10为例, 为自然坐标下裂纹上的任意法向量,节点29、21、3、25的小波形函数值分别为N1=0.75、N2=0.05、N3=0.85、N4=0.25,各节点的附加自由度值分别为(-0.1714,0.0375)、(-0.2013,-0.1112)、(-0.0566,0.0030)、(-0.0656,-0.0826),将以上数据代入公式 中即可得到如图7所示扩展位移数值。
由此例可知,附加自由度是与常规自由度类似的节点位移,在每一个常规自由度上增加一个附加自由度,从而可以反映裂纹的张开位移,这就是对含有齿根裂纹的单元1和单元10的小波单元位移场改进之后所产生的效果;图7为根据齿根单元10中常规自由度和附加自由度求得的裂纹张开位移,图8为裂纹张开位移随裂纹长度的变化曲线;
7)根据5)中所求节点常规自由度数值解计算单元22中的单元应力在横向、纵向及切向上的数值分量σx、σy、τxy(详细请参看冷纪桐赵军著《有限元技术基础》化学工业出版社 第53页);然后根据下列公式计算裂纹在单元22中的扩展角:
式中σ1为单元的最大主应力,计算公式如下:
得到裂纹的扩展角度为θ1=32.46度;裂纹在单元22中的扩展方向以及裂纹扩展的整个过程如图9所示。
Claims (1)
1.一种齿轮裂纹扩展模拟的小波扩展有限元仿真分析方法,该方法基于有限元软件ABAQUS和计算机辅助制图以及数学辅助计算软件构建的平台,对齿轮裂纹扩展进行分析,其特征在于包括以下步骤:
1)建立含齿根裂纹的齿根裂纹模型:
应用计算机辅助制图以及数学辅助计算软件绘制给定齿轮模数和齿数的齿轮模型,将所建立的齿轮模型导入到有限元软件ABAQUS中自动划分网格生成含有H个单元的齿轮模型,H为正整数,并假设在齿根区域内L个单元内存在裂纹,L为自然数,且L<H;
2)获得齿根裂纹模型内各个单元的节点坐标:
根据步骤1)中建立的齿根裂纹模型,应用有限元软件ABAQUS命令流输出功能输出齿根裂纹模型中H个单元的各个节点坐标;
3)根据所建立的齿根裂纹模型计算出H个单元的小波单元刚度矩阵:
将步骤2)齿根裂纹模型中所含有的H个单元的各个节点坐标数据导出,应用平面问题的势能泛函数得到小波单元刚度矩阵;
小波单元刚度矩阵计算式:
分块子矩阵为:
Ke,3=(Ke,2)T
其中a、b、c、d分别为平面单元求解域Ωe={x,y|x∈[a,b],y∈[c,d]}内x与y坐标取值的上下区间;x与y坐标值为步骤2)中所输出H个单元的各个节点坐标值;ξ,η为自然坐标,ξ=(x-a)/(b-a),η=(y-c)/(d-c),则标准求解域Ωs={ξ,η|ξ,η∈[0,1]};E为弹性模量;μ为泊松比;t为厚度;以上各分块子矩阵中Φ1和Φ2分别为m阶j尺度下的一维区间B样条小波尺度函数向量;
其中ξ,η为自然坐标,ξ=(x-a)/(b-a),η=(y-c)/(d-c),ξ,η取值范围为齿根裂纹模型中含有齿根裂纹的L个单元内部区域的自然坐标值;s代表含有齿根裂纹的L个单元内各个节点编号;Ns(ξ,η)代表含有齿根裂纹的L个单元内部区域对应的小波形函数,Ns(ξ,η)=Φ(Re)-1,其中 Φ1和Φ2与步骤3)中相同,分别为m阶j尺度下的一维区间B样条小波尺度函数向量;Re为单位矩阵;us、vs表示含有齿根裂纹的L个单元内各个节点的常规自由度; as1和as2为含有齿根裂纹的L个单元内每一个节点的附加自由度;Hs(ξ,η)为反映含有齿根裂纹的L个单元的位移不连续性函数,计算公式如下:
5)应用扩展有限元改进齿根裂纹模型中齿根区域内含有齿根裂纹的L个单元的小波单元位移场函数之后,计算齿根区域内含有齿根裂纹的L个单元的小波扩展单元刚度矩阵,计算公式如下:
6)计算齿根裂纹模型的结构刚度矩阵;
根据步骤3)、4)、5)对齿根裂纹模型中H个小波单元刚度矩阵的计算以及齿根区域内含有齿根裂纹的L个单元的小波扩展单元刚度矩阵的计算,按照齿轮裂纹模型中H个单元的各个节点次序集成结构刚度矩阵K;
7)建立有限元方程求解节点常规自由度和附加自由度的数值解:
基于步骤6)中齿根裂纹模型结构刚度矩阵K的计算,建立以下有限元求解方程:
式中u、a分别为齿根裂纹模型中各个节点的常规自由度与附加自由度向量;Pu、Pa分别为与模型中常规自由度向量和附加自由度向量相对应的载荷向量,解方程可求出各个节点常规自由度和附加自由度的数值解;
8)计算齿根裂纹模型中裂纹的张开位移和裂纹开裂角度:
基于步骤7)求解有限元方程得到的齿根裂纹模型中附加自由度的数值解,可求齿根裂纹模型中裂纹张开位移,计算式如下:
w代表裂纹张开位移值; 为裂纹上任一点的法向量;s代表含有齿根裂纹的L个单元内各个节点编号;ns为齿根区域内含有齿根裂纹的L个单元包含的所有改进节点的数量;Ns(ξ,η)代表含有齿根裂纹的L个单元内部区域对应的小波形函数,计算裂纹张开位移时Ns(ξ,η)自变量取值为齿根裂纹面上,即当 时对应的数值;as为含有齿根裂纹的L个单元内每一个节点的附加自由度向量,
裂纹的裂纹开裂角度θ1可由下式计算:
式中σ1为单元的最大主应力,计算公式如下:
以上两式中σx、σy、τxy分别表示单元应力在横向、纵向以及切向的数值分量,根据步骤7)中求解有限元方程得到的齿根裂纹模型中各个节点常规自由度数值解计算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110297653 CN102332046B (zh) | 2011-09-30 | 2011-09-30 | 一种齿轮裂纹扩展模拟的小波扩展有限元仿真分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110297653 CN102332046B (zh) | 2011-09-30 | 2011-09-30 | 一种齿轮裂纹扩展模拟的小波扩展有限元仿真分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102332046A CN102332046A (zh) | 2012-01-25 |
CN102332046B true CN102332046B (zh) | 2012-12-19 |
Family
ID=45483822
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110297653 Expired - Fee Related CN102332046B (zh) | 2011-09-30 | 2011-09-30 | 一种齿轮裂纹扩展模拟的小波扩展有限元仿真分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102332046B (zh) |
Families Citing this family (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103412982B (zh) * | 2013-07-20 | 2016-04-06 | 北京工业大学 | 多次降温消除弹性变形损失的预应力建立数值模拟方法 |
CN104574472B (zh) * | 2014-12-31 | 2017-10-31 | 北京大学 | 基于嵌入网格的固体碎裂模拟和动画方法 |
CN104915489A (zh) * | 2015-06-01 | 2015-09-16 | 西南交通大学 | 一种分析裂纹萌生的新型二次插值有限元法 |
US10671768B2 (en) * | 2016-07-01 | 2020-06-02 | The Boeing Company | Finite element modeling and analysis of crack propagation in multiple planes of a structure |
CN106407620B (zh) * | 2016-11-17 | 2019-08-06 | 中国人民解放军济南军区72465部队 | 一种基于abaqus的工程结构响应面随机有限元分析处理方法 |
CN107420523A (zh) * | 2017-09-14 | 2017-12-01 | 东北大学 | 一种具有齿面裂纹缺陷的斜齿轮副啮合刚度计算方法 |
CN107766682B (zh) * | 2017-12-29 | 2020-09-29 | 大连理工大学 | 一种梁结构的静力学分析方法和系统 |
CN108388695A (zh) * | 2018-01-26 | 2018-08-10 | 昌宇应力技术(上海)有限公司 | 根据齿根残余应力分布定量估算齿轮局部疲劳强度的方法 |
CN108389630A (zh) * | 2018-02-08 | 2018-08-10 | 西南交通大学 | 集成纳米hap分布影响的牙釉质宏观弹性模量计算方法 |
CN109165404B (zh) * | 2018-07-16 | 2021-04-23 | 四川大学 | 一种扩展光滑无网格伽辽金法 |
CN109271655B (zh) * | 2018-07-23 | 2021-05-11 | 南京航空航天大学 | 一种基于非对称有限元算法的材料尺度效应分析方法 |
CN108763841B (zh) * | 2018-07-24 | 2022-04-05 | 北京航空航天大学青岛研究院 | 一种基于对偶边界元和应变能优化分析的弹性断裂仿真方法 |
CN109509224B (zh) * | 2018-11-09 | 2022-04-19 | 西南科技大学 | 一种平面材料小孔应力集中极坐标化展示方法 |
CN109918712B (zh) * | 2019-01-23 | 2022-07-05 | 昆明理工大学 | 一种基于p型有限元法和围线积分法求解复合型应力强度因子的计算方法 |
CN110288569B (zh) * | 2019-05-27 | 2023-05-23 | 中国飞机强度研究所 | 一种基于水平集方法的裂纹近场变形分析区域确定方法 |
CN110400362B (zh) * | 2019-07-12 | 2023-02-24 | 重庆大学 | 一种基于图像的abaqus二维裂纹建模方法、系统及计算机可读存储介质 |
CN110688717B (zh) * | 2019-09-25 | 2022-08-16 | 西北工业大学 | 齿轮齿条结构的可靠性分析方法及装置 |
CN113155651B (zh) * | 2021-03-30 | 2022-05-10 | 东北大学 | 基于疲劳试验与仿真的直齿轮齿根裂纹扩展规律分析方法 |
CN113373299B (zh) * | 2021-06-07 | 2023-03-10 | 广东工业大学 | 一种抑制裂纹或改变裂纹走向的方法及其装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101625710B (zh) * | 2009-08-03 | 2011-01-05 | 北京工业大学 | 含裂纹故障的变工况单级重载齿轮箱振动响应仿真分析方法 |
-
2011
- 2011-09-30 CN CN 201110297653 patent/CN102332046B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN102332046A (zh) | 2012-01-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102332046B (zh) | 一种齿轮裂纹扩展模拟的小波扩展有限元仿真分析方法 | |
Bacciocchi et al. | Vibration analysis of variable thickness plates and shells by the generalized differential quadrature method | |
Olshevskiy et al. | Three-dimensional solid brick element using slopes in the absolute nodal coordinate formulation | |
Bediz | A spectral-Tchebychev solution technique for determining vibrational behavior of thick plates having arbitrary geometry | |
CN106844818A (zh) | 基于粗糙表面的直齿轮三维接触刚度计算方法 | |
CN102096736B (zh) | 一种基于渐近变分法的复合材料层合板仿真及优化方法 | |
Castonguay et al. | Application of high-order energy stable flux reconstruction schemes to the Euler equations | |
CN103838913B (zh) | 曲线箱梁弯桥的有限单元法 | |
CN104317985B (zh) | 一种基于界带有限元和拉格朗日坐标的流体仿真方法 | |
Zhang et al. | Extract of stress intensity factors on honeycomb elements by the numerical manifold method | |
CN105955928A (zh) | 基于cfd预报船舶阻力的计算方法 | |
Pi et al. | First order sensitivity analysis of flexible multibody systems using absolute nodal coordinate formulation | |
CN104281730A (zh) | 一种大转动变形的板壳结构动响应的有限元分析方法 | |
CN109992824A (zh) | 一种任意形状实心弹性板横向振动半解析分析方法 | |
CN106354954A (zh) | 一种基于叠层基函数的三维力学模态仿真模拟方法 | |
Wang et al. | Position-dependent stability prediction for multi-axis milling of the thin-walled component with a curved surface | |
Lavrentiev et al. | Hardware/Software Solution for Low Power Evaluation of Tsunami Danger | |
Archibald et al. | Multiwavelet discontinuous Galerkin-accelerated Exact Linear Part (ELP) method for the shallow-water equations on the cubed sphere | |
CN103177168A (zh) | 一种齿轮齿面摩擦因数的反求方法 | |
Xu et al. | A new method and application of full 3d numerical simulation for hydraulic fracturing horizontal fracture | |
CN106909699A (zh) | 基于Galerkin条形传递函数的薄板振动特性分析方法 | |
Carrington et al. | Three-dimensional ALE-FEM method for fluid flow in domains with moving boundaries part 1: algorithm description | |
Chinthala et al. | A novel strong-form random differential quadrature method to compute the stress intensity factor in fracture mechanics | |
Chen et al. | New direction of computational fluid dynamics and its applications in industry | |
Guasch et al. | A stabilized Arbitrary Lagrangian Eulerian Finite Element Method for the mixed wave equation with application to diphthong production |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20121219 Termination date: 20150930 |
|
EXPY | Termination of patent right or utility model |