CN109918712A - 一种基于p型有限元法和围线积分法求解复合型应力强度因子的计算方法 - Google Patents
一种基于p型有限元法和围线积分法求解复合型应力强度因子的计算方法 Download PDFInfo
- Publication number
- CN109918712A CN109918712A CN201910063707.XA CN201910063707A CN109918712A CN 109918712 A CN109918712 A CN 109918712A CN 201910063707 A CN201910063707 A CN 201910063707A CN 109918712 A CN109918712 A CN 109918712A
- Authority
- CN
- China
- Prior art keywords
- finite element
- contour
- field
- stress
- model
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000004364 calculation method Methods 0.000 title claims abstract description 12
- 239000002131 composite material Substances 0.000 title claims abstract description 11
- 238000006073 displacement reaction Methods 0.000 claims abstract description 32
- 238000005336 cracking Methods 0.000 claims abstract description 15
- 239000011159 matrix material Substances 0.000 claims description 18
- 239000000463 material Substances 0.000 claims description 13
- 208000037656 Respiratory Sounds Diseases 0.000 claims description 9
- 230000017105 transposition Effects 0.000 claims description 5
- 230000007704 transition Effects 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 5
- 238000002203 pretreatment Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000005477 standard model Effects 0.000 description 1
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明涉及一种基于p型有限元法和围线积分法求解复合型应力强度因子的计算方法,属于断裂力学领域,特别是对含裂纹构件断裂进行计算。步骤1构建含裂纹构件有限元模型;步骤2将含裂纹构件有限元模型采用p型有限元法计算得到含裂纹构件的位移场、应力场和应变场;步骤3采用围线积分法选择包含裂纹尖端的一条围线,根据步骤2得到的该围线上的位移值和应力值导出裂纹处的应力强度因子;步骤4判断所得的应力强度因子是否满足精度要求,如未满足精度要求,提高插值多项式的阶次,返回步骤2。本发明将p型有限元法应用于断裂力学领域,并结合围线积分法,得到一种求解复合型应力强度因子的自适应计算方法。
Description
技术领域
本发明涉及一种基于p型有限元法和围线积分法求解复合型应力强度因子的计算方法,属于断裂力学领域,特别是对含裂纹构件断裂进行计算。
背景技术
应力强度因子是判断含裂纹构件断裂及计算裂纹扩展速度的重要参数,在裂纹体分析中占据着重要地位。目前已有众多理论和数值解法求解应力强度因子,如解析法、混合法、权函数法、线弹簧法、有限元法、边界元法等。其中有限元法由于其强大的建模能力,可充分利用计算机的计算能力,能够适用于各种复杂的几何情况,能够在各种工程问题中获得较高的精度,成为求解应力强度因子较为有效的一种方法。
有限元解按其结构可分为三类:h型;p型;hp型有限元法,即分别通过加密网格、提高插值多项式阶次或两者同时进行来提高有限元解的精度。相比传统的有限元法,p型有限元方法具有网格划分少、收敛速度快、计算精度高、前处理少等优点。
围线积分法基于Betti功互等定理,只需知道远场积分路径上的位移和应力应变,就能计算出混合型裂纹的应力强度因子,导出应力强度因子时具有超收敛性。
使用有限元法对断裂力学问题进行研究的大多是传统的h型有限元法或者基于传统有限元的发展(如扩展有限元法等),因而存在如下问题:
(1)传统的h型有限元缺乏有效的误差估计模式,计算精度的控制较为依赖研究人员的经验。为获得较好的计算精度,需要多次重新划分网格,来判断有限元解是否收敛;
(2)传统的h型有限元基于低阶的插值函数,在模拟裂纹尖端的高梯度应力应变场时存在先天不足,为获得较高精度,需要大量的加密网格,前处理及计算成本较大;
(3)由于裂尖应力场的奇异性,基于传统有限元导出应力强度因子时,由于裂尖应力应变场的误差较大,导出应力强度因子时精度较低;或者要获得较高精度应力强度因子时需要增大前处理和计算成本。
本发明基于国家自然科学基金(资助编号:51769011),提出了将p型有限元法应用于断裂力学领域,并结合围线积分法来求解复合型应力强度因子。
发明内容
针对上述现有技术存在的问题及不足,本发明提供一种基于p型有限元法和围线积分法求解复合型应力强度因子的计算方法。本发明将p型有限元法应用于断裂力学领域,并结合围线积分法,得到一种求解复合型应力强度因子的自适应计算方法。本发明
一种基于p型有限元法和围线积分法求解复合型应力强度因子的计算方法,按以下步骤进行:
步骤1、构建含裂纹构件有限元模型,该步骤包括:
根据裂纹构件的实际情况,采集构件的几何参数、材料参数和荷载参数。
根据几何参数建立几何模型,进行几何模型的建立。
进行网格划分:依据几何模型建立有限元网格。由于裂纹尖端的奇异性,裂纹尖端的网格应适当加密。裂纹外围应力应变场的梯度较小,网格可以较为稀疏。
进行材料设定:根据采集的材料参数建立有限元的材料模型,将对应的材料模型赋予相应的单元。
进行载荷与边界条件的施加:根据采集的荷载参数(位移参数和荷载参数),在相应边界上施加相应的边界条件;
步骤2、将步骤1构建含裂纹构件有限元模型采用p型有限元法计算得到含裂纹构件的位移场、应力场和应变场,
具体步骤包括:
步骤2.1、求解得到结构位移列阵a
根据方程:Ka=F (1),其中K=∑eGTKeG,结构整体刚度矩阵,Ke为单元刚度矩阵;F=∑eGTfe,结构结点载荷列阵;G为单位节点自由度和结构节点自由度的转换矩阵,GT为G的转置,上标T表示矩阵的转置;e表示单个单元;a为结构位移列阵;
再根据公式单元刚度矩阵Ke=∫ΩBTDBdΩ (2),
单元等效节点载荷列阵
单元内部节点力
外部节点力
上述公式(2)至(5)中,Ω表示在单元内部;B=LNI,L为微分算子,NI为插值函数矩阵或形函数矩阵,I为单元的标号,Γt为单元的外部边界;D为应力矩阵,b为体力,为荷载边界条件,分别由步骤1中的材料属性以及荷载和位移边界条件确定;
将公式(2)至(5)以及NI为p型有限元法的插值多项式代入到(1)中,通过求解线性方程组(1)得到结构位移列阵a;
NI采用p型有限元法的插值多项式,典型的二维p型有限元形函数基于勒让德正交多项式,以二维四边形单元为例(附图1所示),在坐标系(η,ξ)下,点p1、p2、p3和p4分为单元的四个顶点,Γ1、Γ2、Γ3和Γ4分为单元的四条边,形函数的构造形式如下所示:
p≥1,点模式基函数与常规的拉格朗日型基函数时一致的,以四边形四节点单元为例按如下展开:
p≥2,边模式基函数:
边Γ1(ξ=-1,-1≤η≤1)上的形函数可按如下公式给出(中1为边号,i为插值多项式的阶次):
式中:
这里,Pn(t)是阶数为n≥0的勒让德多项式
类似的,可以如下定义与边Γk(2≤k≤4)相关联的形状函数:
p≥4时内部模式基函数:
本发明采用的p型有限元方法,其插值多项式的阶次p=1,…,∞,提升插值多项式时低阶的刚度矩阵可继续沿用,只需计算高阶部分,避免了刚度矩阵低阶部分的重复计算,具有良好的承袭性,节约了前处理的成本。
步骤2.2、根据公式aI=Ga求解得到单元位移列阵aI;
步骤2.3、根据位移场公式u=NIaI,应变场公式ε=Lu,应力场公式σ=Dε=DBaI,求解得到位移场u,应变场ε,应力场σ。
步骤3、采用围线积分法选择包含裂纹尖端的一条围线,根据步骤2得到的该围线上的位移值和应力值导出裂纹处的应力强度因子;
本发明采用围线积分法导出应力强度因子。考虑如附图2所示的弹性裂纹场Ω,有局部笛卡尔坐标系(x1,x2,x3)和局部极坐标系(r,θ,z),z轴与x3重合,用u(x1,x2,x3)和u(r,θ,z)分别表示局部笛卡尔和局部极坐标系下导出域Ωs中裂纹尖端的应变场。其中KⅠ、KⅡ、KⅢ分别为Ⅰ型、Ⅱ型、Ⅲ型应力强度因子(即分别为张开型、滑开型和撕开型裂纹应力强度因子),位移场u下计算得到围线上牵引力向量用T(u)表示:
在围线上(6)
其中公式(6)中,为位移场u中计算得到的应力张量,nj为围线上外方向法向量分量。
在围线积分法中,导出域Ωs退化到域Ωs的一个平面片上,如附图3所示。应力强度因子KⅠ、KⅡ和KⅢ便可通过围线和上的线积分得到:
式(7)至(9)中:和为导出函数,是裂尖附近弹性解渐近展开的表达式在相应模式下负的特征值,由以下公式给出:
其中:G为剪切强度模量,κ为克罗斯夫常数,平面应力状态下,κ=(3-ν)/(1+ν),平面应变状态下,κ=3-4ν,ν为泊松比;
和分别为虚位移场和下计算得到的牵引力向量,同公式(6);p3和p4为裂纹面上预先给定的牵引力荷载,由步骤1结构的模型参数确定;这里u为步骤2中计算得到的位移场。
选择一条包含裂纹尖端的围线,由上述公式导出应力强度因子KⅠ、KⅡ和KⅢ。围线Γ2的选择是很灵活的,不必接近裂纹前端。由于裂纹前端最内层单元中数值解的精度较低,此时导出的应力强度因子误差较大,因而选择围线Γ2时,一般包含裂尖附近的最内层单元。
步骤4、判断步骤3所得的应力强度因子是否满足精度要求,如未满足精度要求,提高插值多项式的阶次,返回步骤2:
对标准模型,将数值计算结构与理论解或实验结果进行对比,判断数值计算结果是否满足精度要求。
对非标准模型,没有理论解或者实验结果进行对比的情况下,依次提高插值多项式阶次,再由围线积分法导出应力强度因子。观察后得到的结果是否与前面计算的结果偏差是否在可接受范围内。
精度不满足要求,或所得结果直接偏差不在可接受范围内。再次返回步骤二,重新计算位移场和应力应变场,导出应力强度因子,直至数值计算结果满足精度要求
所述步骤3中的围线为一条远离裂尖且包含裂尖的围线,且由于裂尖外层单元受奇异性的影响,也应当包含裂纹外层的一层单元。
上述同一个标号代表同一个参数值,上标T表示相应该矩阵的转置。
本发明的有益效果是:
(1)本方法能减少计算成本,提高计算效率、收敛速率和精度。
(2)本方法计算方法能直接应用在实际工程问题中,对含裂纹构件断裂进行计算,且计算精度高。
附图说明
图1是本发明标准四边形母单元的示意图;
图2是本发明裂纹尖端邻域及全局坐标、局部坐标系的示意图;
图3为裂纹尖端退化后的导出域及局部笛卡尔坐标系、极坐标系的示意图;
图4是本发明实施例1边缘裂纹构件的结构示意图;
图5为实施例1边缘裂纹网格示意图;
图6为实施例1采用p型有限元法计算有限元模型得到的位移场云图;
图7为实施例1采用p型有限元法计算有限元模型得到的应力场云图;
具体实施方式
下面结合附图和具体实施方式,对本发明作进一步说明。
实施例1
如图4所示一边缘裂纹板,其中a=5为裂缝宽度,板宽b=10,板长h=20,板顶端受竖向单位应力σ=1.0作用,底部左端固定,底部受竖直固定约束。板的材料参数:杨氏模量E=107Pa,泊松比ν=0.3,结构处于平面应变状态。
步骤1、构建含裂纹构件有限元模型:
首先根据结构的几何参数构建几何模型、定义裂纹;
根据结构的几何特征划分网格,模型的网格划分如附图5所示,外围采用稀疏的网格划分,裂纹尖端处于板中央,采用了双层网格加密,最内层网格尺寸是其外层单元的0.15倍。
裂纹尖端奇异点的存在,应力梯度较大,这里按0.15的几何因子对裂纹尖端进行了两个网格加密,即内层网格尺寸为外层网格尺寸的0.15倍。本发明中,除了裂纹尖端网格尺寸比较密集之外,其他地方网格网格非常稀疏,p型有限元方法提升阶次的过程中,提高了稀疏网格处的计算精度。
根据采集的材料参数建立有限元的材料模型,将相应的材料模型赋予对应的单元;
施加结构的位移边界条件以及荷载边界条件:在底部施加竖向的固定约束,左端点施加横向固定约束(位移边界);顶部施加竖直向上的单位应力(σ=1.0)(荷载边界)。
步骤2、将步骤1构建含裂纹构件有限元模型采用p型有限元法计算得到含裂纹构件的位移场、应力场和应变场,设定有限元计算的插值多项式阶次;这里直接计算了p=1~8时的应力场和位移场。
通过本步骤得到的位移场云图如附图6所示,应力场云图如附图7所示(p≥6时,有限元解的能量范数误差较低(下表中给出),云图变化很小,仅给出了p=1~6时的位移云图和应力云图)。
步骤3、采用围线积分法选择包含裂纹尖端的一条围线,根据步骤2得到的该围线上的位移值和应力值导出裂纹处的应力强度因子。
选择以裂纹尖端为圆心,以r为半径的一个圆作为围线,利用前面提高的应力强度因子导出公式计算应力强度因子。这里依次导出了p=1~8时应力强度因子的值。
r的选择没有特殊的要求,由于裂纹尖端奇异性的存在,当围线处于最内层单元中时,导出的应力强度因子振荡较大,因而围线以包含最内层单元为佳,继续增大围线的半径r,应力强度因子的误差变化较小。
本模型具有经验解,由经验公式计算此时应力强度因子的经验解。
上式中:为经验公式。a/b≤0.6时,函数F:
当a/b=0.5时,KI=11.2018。
基于p型有限元求解该构件的位移场、应力应变场通过围线积分法导出应力强度因子,下表给出了p=1,2,…,8时有限元解的能力范数误差和对应的情形下导出的应力强度因子,及相对误差
基于p型有限元法,随插值多项式阶次p的增加,有限元解的能力范数误差逐渐收敛。基于有限元解,采用围线积分法导出应力强度因子,应力强度因子的解也随着应力场精度提高而提高。
本次计算的模型,结构简单,网格划分较少,所需的计算成本较低,因而直接计算了p=1~8时的所有情形。
当计算的模型结构复杂,网格划分较多时,可依次提高插值多少式的阶次,并依次导出计算的应力强度因子,继续计算至应力强度因子的变化较小时,可认为是比较理想的结果。
依次提高插值多项式时,低阶插值多项式的部分不用重复计算,只需计算高阶插值多项式的部分,重新组成高阶的刚度矩阵。
基于上表可看出,本发明计算应力强度因子时,具有自适应性,前处理少,计算收敛快,采用较低的计算成本,便可获得较高精度的数值解。
实例2
模型如实例1所示,改变裂纹尺寸a。
实例1中,当插值多项式阶次p=5时,应力强度因子的精度低于1%。
在实例2中,在相似网格条件下,选择插值多项式p=6,确保应力强度因子精度满足要求。
建立模型、划分网格、材料设定、施加边界条件后,直接计算p=6时的位移场、应力场、应变场,采用围线积分法应力强度因子。改变裂纹长度,继续计算当a/b=0.1,0.2,…,0.6时的应力强度因子,具体数据如下表:
表中KΙ为使用本方法计算的结果,KⅠ解析为经验公式计算的结果。
由上表可看出改变裂纹长度,网格条件略有不同,对比经验解,此时应力强度因子的精度较高,且具有良好的数值稳定性。
以上结合附图对本发明的具体实施方式作了详细说明,但是本发明并不限于上述实施方式,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下作出各种变化。
Claims (3)
1.一种基于p型有限元法和围线积分法求解复合型应力强度因子的计算方法,其特征在于:按以下步骤进行:
步骤1、构建含裂纹构件有限元模型,该步骤包括:含裂纹构件有限元模型中进行几何模型的建立;进行网格划分;进行材料设定;进行载荷与边界条件的施加;
步骤2、将步骤1构建含裂纹构件有限元模型采用p型有限元法计算得到含裂纹构件的位移场、应力场和应变场,
步骤3、采用围线积分法选择包含裂纹尖端的一条围线,根据步骤2得到的该围线上的位移值和应力值导出裂纹处的应力强度因子;
步骤4、判断步骤3所得的应力强度因子是否满足精度要求,如未满足精度要求,提高插值多项式的阶次,返回步骤2。
2.根据权利要求1所述的基于p型有限元法和围线积分法求解复合型应力强度因子的计算方法,其特征在于:所述步骤2具体步骤包括:
步骤2.1、求解得到结构位移列阵a
根据方程:Ka=F (1),其中K=∑eGTKeG,结构整体刚度矩阵,Ke为单元刚度矩阵;F=∑eGTfe,结构结点载荷列阵;G为单位节点自由度和结构节点自由度的转换矩阵,GT为G的转置,上标T表示矩阵的转置;e表示单个单元;a为结构位移列阵;
再根据公式单元刚度矩阵Ke=∫ΩBTDBdΩ (2),
单元等效节点载荷列阵
单元内部节点力
外部节点力
上述公式(2)至(5)中,Ω表示在单元内部;B=LNI,L为微分算子,NI为插值函数矩阵或形函数矩阵,I为单元的标号,Γt为单元的外部边界;D为应力矩阵,b为体力,为荷载边界条件,分别由步骤1中的材料属性以及荷载和位移边界条件确定;
将公式(2)至(5)以及NI为p型有限元法的插值多项式代入到(1)中,通过求解线性方程组(1)得到结构位移列阵a;
步骤2.2、根据公式aI=Ga求解得到单元位移列阵aI;
步骤2.3、根据位移场公式u=NIaI,应变场公式ε=Lu,应力场公式σ=Dε=DBaI,求解得到位移场u,应变场ε和应力场σ。
3.根据权利要求1所述的基于p型有限元法和围线积分法求解复合型应力强度因子的计算方法,其特征在于:所述步骤3中的围线为一条远离裂尖且包含裂尖的围线,且由于裂尖外层单元受奇异性的影响,也应当包含裂纹外层的一层单元。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910063707.XA CN109918712B (zh) | 2019-01-23 | 2019-01-23 | 一种基于p型有限元法和围线积分法求解复合型应力强度因子的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910063707.XA CN109918712B (zh) | 2019-01-23 | 2019-01-23 | 一种基于p型有限元法和围线积分法求解复合型应力强度因子的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109918712A true CN109918712A (zh) | 2019-06-21 |
CN109918712B CN109918712B (zh) | 2022-07-05 |
Family
ID=66960534
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910063707.XA Active CN109918712B (zh) | 2019-01-23 | 2019-01-23 | 一种基于p型有限元法和围线积分法求解复合型应力强度因子的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109918712B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110866350A (zh) * | 2019-09-16 | 2020-03-06 | 合肥通用机械研究院有限公司 | 一种裂纹尖端单元部分增强的扩展有限元方法 |
CN111916160A (zh) * | 2020-07-13 | 2020-11-10 | 烟台南山学院 | 一种计算材料裂纹尖端应力场系数的方法 |
CN112182489A (zh) * | 2020-10-14 | 2021-01-05 | 杭州电子科技大学 | 一种基于偏微分方程求解的二维高阶网格生成方法 |
CN112966421A (zh) * | 2021-03-16 | 2021-06-15 | 昆明理工大学 | 使用p型有限元法计算薄板结构屈曲载荷因子和相应屈曲形状的计算方法 |
CN113742960A (zh) * | 2021-08-10 | 2021-12-03 | 昆明理工大学 | 一种三维平面曲边裂纹断裂参数高阶有限元数值模拟方法 |
CN114169209A (zh) * | 2021-12-24 | 2022-03-11 | 哈尔滨工业大学 | 一种提取非均匀材料混合型弥散裂纹尖端参量的数值方法 |
CN114266776A (zh) * | 2022-03-03 | 2022-04-01 | 南京航空航天大学 | 一种应用复合裂纹位移场函数的数字图像相关方法 |
CN116629079A (zh) * | 2023-07-21 | 2023-08-22 | 北京大学 | 混合有限元空间构造及求解线弹性力学问题的方法及装置 |
CN117057166A (zh) * | 2023-10-11 | 2023-11-14 | 合肥通用机械研究院有限公司 | 应力集中部位裂纹自由表面处应力强度因子的计算方法 |
CN117150822A (zh) * | 2023-10-30 | 2023-12-01 | 中南大学 | 界面裂纹的热力耦合应力强度因子计算方法及系统 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070272018A1 (en) * | 2006-05-24 | 2007-11-29 | Honeywell International Inc. | Determination of remaining useful life of gas turbine blade |
CN101788425A (zh) * | 2010-02-09 | 2010-07-28 | 浙江工业大学 | 一种结构件复合型裂纹前缘应力强度因子分离和分布的确定方法 |
US20110191074A1 (en) * | 2010-02-03 | 2011-08-04 | Kellogg Brown & Root Llc | Systems and Methods for Performing Stress Intensity Factor Calculations Using Non-Singluar Finite Elements |
CN102332046A (zh) * | 2011-09-30 | 2012-01-25 | 北京工业大学 | 一种齿轮裂纹扩展模拟的小波扩展有限元仿真分析方法 |
CN103020426A (zh) * | 2012-11-23 | 2013-04-03 | 北京航空航天大学 | 一种矩形板中心斜裂纹疲劳扩展寿命预测的简化方法 |
US20130090902A1 (en) * | 2010-06-28 | 2013-04-11 | Yao Yao | Method and System for Modeling Fractures in Ductile Rock |
CN103955604A (zh) * | 2014-04-11 | 2014-07-30 | 南京航空航天大学 | 一种含裂纹金属梯度材料剩余强度预测方法 |
CN104239588A (zh) * | 2013-06-18 | 2014-12-24 | 深圳市网蓝实业有限公司 | 一种基于超单元构造方法的结构功能分析方法 |
CN105808884A (zh) * | 2016-03-30 | 2016-07-27 | 北京航空航天大学 | 一种基于分形理论的有界不确定性平面裂纹应力强度因子上下界的预测方法 |
CN108349218A (zh) * | 2015-11-05 | 2018-07-31 | 康宁股份有限公司 | 具有确定模量对比的层压玻璃制品及其形成方法 |
-
2019
- 2019-01-23 CN CN201910063707.XA patent/CN109918712B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070272018A1 (en) * | 2006-05-24 | 2007-11-29 | Honeywell International Inc. | Determination of remaining useful life of gas turbine blade |
US20110191074A1 (en) * | 2010-02-03 | 2011-08-04 | Kellogg Brown & Root Llc | Systems and Methods for Performing Stress Intensity Factor Calculations Using Non-Singluar Finite Elements |
CN101788425A (zh) * | 2010-02-09 | 2010-07-28 | 浙江工业大学 | 一种结构件复合型裂纹前缘应力强度因子分离和分布的确定方法 |
US20130090902A1 (en) * | 2010-06-28 | 2013-04-11 | Yao Yao | Method and System for Modeling Fractures in Ductile Rock |
CN102332046A (zh) * | 2011-09-30 | 2012-01-25 | 北京工业大学 | 一种齿轮裂纹扩展模拟的小波扩展有限元仿真分析方法 |
CN103020426A (zh) * | 2012-11-23 | 2013-04-03 | 北京航空航天大学 | 一种矩形板中心斜裂纹疲劳扩展寿命预测的简化方法 |
CN104239588A (zh) * | 2013-06-18 | 2014-12-24 | 深圳市网蓝实业有限公司 | 一种基于超单元构造方法的结构功能分析方法 |
CN103955604A (zh) * | 2014-04-11 | 2014-07-30 | 南京航空航天大学 | 一种含裂纹金属梯度材料剩余强度预测方法 |
CN108349218A (zh) * | 2015-11-05 | 2018-07-31 | 康宁股份有限公司 | 具有确定模量对比的层压玻璃制品及其形成方法 |
CN105808884A (zh) * | 2016-03-30 | 2016-07-27 | 北京航空航天大学 | 一种基于分形理论的有界不确定性平面裂纹应力强度因子上下界的预测方法 |
Non-Patent Citations (4)
Title |
---|
TRACEY D M: "Finite elements for determination of crack tip elastic stress intensity factors", 《ENGINEERING FRACTURE MECHANICS》 * |
杨晓翔 等: "求解混合型裂纹应力强度因子的围线积分法", 《计算结构力学及其应用》 * |
王周宏 等: "建立p型有限元刚度矩阵的快速方法及分析", 《华中理工大学学报》 * |
陆洋春: "基于p型有限元法研究断裂力学问题", 《中国水运》 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110866350A (zh) * | 2019-09-16 | 2020-03-06 | 合肥通用机械研究院有限公司 | 一种裂纹尖端单元部分增强的扩展有限元方法 |
CN111916160A (zh) * | 2020-07-13 | 2020-11-10 | 烟台南山学院 | 一种计算材料裂纹尖端应力场系数的方法 |
CN112182489B (zh) * | 2020-10-14 | 2023-11-24 | 杭州电子科技大学 | 一种基于偏微分方程求解的二维高阶网格生成方法 |
CN112182489A (zh) * | 2020-10-14 | 2021-01-05 | 杭州电子科技大学 | 一种基于偏微分方程求解的二维高阶网格生成方法 |
CN112966421A (zh) * | 2021-03-16 | 2021-06-15 | 昆明理工大学 | 使用p型有限元法计算薄板结构屈曲载荷因子和相应屈曲形状的计算方法 |
CN112966421B (zh) * | 2021-03-16 | 2023-12-26 | 昆明理工大学 | 使用p型有限元法计算薄板结构屈曲载荷因子和相应屈曲形状的计算方法 |
CN113742960A (zh) * | 2021-08-10 | 2021-12-03 | 昆明理工大学 | 一种三维平面曲边裂纹断裂参数高阶有限元数值模拟方法 |
CN114169209A (zh) * | 2021-12-24 | 2022-03-11 | 哈尔滨工业大学 | 一种提取非均匀材料混合型弥散裂纹尖端参量的数值方法 |
CN114266776A (zh) * | 2022-03-03 | 2022-04-01 | 南京航空航天大学 | 一种应用复合裂纹位移场函数的数字图像相关方法 |
CN114266776B (zh) * | 2022-03-03 | 2022-06-10 | 南京航空航天大学 | 一种应用复合裂纹位移场函数的数字图像相关方法 |
CN116629079A (zh) * | 2023-07-21 | 2023-08-22 | 北京大学 | 混合有限元空间构造及求解线弹性力学问题的方法及装置 |
CN116629079B (zh) * | 2023-07-21 | 2024-01-23 | 北京大学 | 混合有限元空间构造及求解线弹性力学问题的方法及装置 |
CN117057166B (zh) * | 2023-10-11 | 2023-12-26 | 合肥通用机械研究院有限公司 | 应力集中部位裂纹自由表面处应力强度因子的计算方法 |
CN117057166A (zh) * | 2023-10-11 | 2023-11-14 | 合肥通用机械研究院有限公司 | 应力集中部位裂纹自由表面处应力强度因子的计算方法 |
CN117150822A (zh) * | 2023-10-30 | 2023-12-01 | 中南大学 | 界面裂纹的热力耦合应力强度因子计算方法及系统 |
CN117150822B (zh) * | 2023-10-30 | 2024-02-06 | 中南大学 | 界面裂纹的热力耦合应力强度因子计算方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN109918712B (zh) | 2022-07-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109918712A (zh) | 一种基于p型有限元法和围线积分法求解复合型应力强度因子的计算方法 | |
Lee et al. | On error estimation and adaptive refinement for element free Galerkin method: Part I: stress recovery and a posteriori error estimation | |
CN103838852B (zh) | 一种快速查找多块结构化网格对接关系的方法 | |
CN106646645B (zh) | 一种重力正演加速方法 | |
GAFFNEY, JR et al. | Euler calculations for wings using Cartesian grids | |
Cosentino et al. | Numerical optimization design of advanced transonic wing configurations | |
CN104408773B (zh) | 一种结构化网格非匹配界面插值的方法 | |
Kenway et al. | Aerodynamic shape optimization of the CRM configuration including buffet-onset conditions | |
CN104809692A (zh) | 一种等高线生成dem的内插方法 | |
CN112966398A (zh) | 一种基于应力分布的Voronoi多孔梯度结构生成方法 | |
CN110321571A (zh) | 一种蜂窝板壳结构的力学参数数值提取方法 | |
Murman et al. | A vortex wake capturing method for potential flow calculations | |
Wang et al. | Mixed element type unstructured grid generation and its application to viscous flow simulation | |
CN113012286B (zh) | 一种基于密集点云数据构建道路dem的方法 | |
CN108804791A (zh) | 一种适用于埋入式进气道布局的飞行器参数化方法 | |
CN109325257B (zh) | 一种快速生成飞行载荷分析模型的建模方法 | |
Ward et al. | Hybrid prismatic/tetrahedral grid generation for complex 3-D geometries | |
CN113722965A (zh) | 一种基于积分-广义有限差分数值离散算子的断裂模拟方法 | |
CN109087391A (zh) | 一种三维场景下的军标标绘方法 | |
CN115392032A (zh) | 一种gis-mpm无缝集成的动态三维地质模型构建方法 | |
CN104598682B (zh) | 一种等效结点荷载的通用精确积分计算方法 | |
CN111191395A (zh) | 一种嵌套模型建模方法及设备 | |
Paiva et al. | Approximating implicit curves on triangulations with affine arithmetic | |
CN107766673A (zh) | 一种半径可控的参数化三维前缘钝化设计方法 | |
RUMSEY | A computational analysis of flow separation over five different airfoil geometries at high angles-of-attack |
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 |