CN106971036A - 减小复杂几何形状孔附近应力集中的方法 - Google Patents
减小复杂几何形状孔附近应力集中的方法 Download PDFInfo
- Publication number
- CN106971036A CN106971036A CN201710177560.8A CN201710177560A CN106971036A CN 106971036 A CN106971036 A CN 106971036A CN 201710177560 A CN201710177560 A CN 201710177560A CN 106971036 A CN106971036 A CN 106971036A
- Authority
- CN
- China
- Prior art keywords
- hole
- stress concentration
- geometric shapes
- ring
- complex geometric
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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
-
- 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/06—Power analysis or power optimisation
Abstract
本发明公开了一种减小复杂几何形状孔附近应力集中的方法,对于复杂几何形状孔,在孔周焊接衬入杨氏模量沿法向连续递增变化的功能梯度材料加强环,此加强环的形状与孔的形状一致。此加强环可以有效缓解孔附近的应力集中。
Description
技术领域
本发明具体涉及一种减小复杂几何形状孔附近应力集中的方法。
背景技术
机械工程领域,为了减轻重量、节约成本或是结构设计需求,各种几何形状的孔洞结构普遍存在。然而,由于孔的存在破坏了结构的几何连续型,在孔附近会产生显著的应力集中,这极大地削弱了材料结构的强度。
功能梯度材料(FGM)由于微观组织结构的连续变化,在缓解应力集中方面具有显著优点。最近,Sburlati等人为了缓解圆孔附近的应力集中,提出在孔周嵌入功能梯度加强环并建立了相应模型,分析论证了该方法的可行性。
但事实上,工程中除了圆形孔洞,椭圆形、正方形、三角形等其它一般几何形状的孔洞结构也被广泛应用,并且通常而言,这些几何形状孔附近的应力集中要比圆孔附近的应力集中更加显著。因此,急需一种能够减小复杂几何形状孔附近应力集中的解决方案。
发明内容
本发明的目的在于克服现有技术中的不足,提供了一种减小复杂几何形状孔附近应力集中的方法,通过在孔周焊接衬入杨氏模量沿法向连续递增变化的功能梯度材料加强环,可以有效缓解孔附近的应力集中。
为解决上述技术问题,本发明提供了一种减小复杂几何形状孔附近应力集中的方法,其特征是,对于复杂几何形状孔,在孔周焊接衬入杨氏模量沿法向连续递增变化的功能梯度材料加强环,此加强环的形状与孔的形状一致。
进一步的,几何形状孔包括圆形孔、椭圆形孔、三角形孔、正方形孔和矩形孔。
进一步的,在映射平面ζ中,加强环厚度t的范围取为0.05rN≤t≤0.2rN,其中rN为ζ平面几何孔半径。
进一步的,在映射平面ζ中,杨氏模量递变规律公式为:
其中,Ep表示板内的杨氏模量,r为ζ平面中半径。
进一步的,梯度指数值n的取值范围为0<n≤5。
进一步的,应力集中的计算过程包括以下步骤:
步骤S1,在z平面中,利用分层均匀化方法将杨氏模量沿法向连续变化的功能梯度材料加强环分解为N层厚度相等的薄环;
步骤S2,利用保角变换函数从z平面映射到ζ平面,原z平面的N层均质薄环转变为ζ平面的N层同心圆环;根据映射关系,建立ζ平面环及板内应力场表达式及ζ平面中的每个圆环内复势函数表达式;
步骤S3,计算在ζ平面中的每个圆环内复势函数,基于ζ平面环及板内应力场的表达式求得环和板内的应力场。
与现有技术相比,本发明所达到的有益效果是:1)通过在孔周焊接衬入杨氏模量沿法向连续递增变化的功能梯度材料加强环,可以有效缓解孔附近的应力集中。2)利用保角变换技术及分层均匀化方法,简化了含功能梯度材料加强环的任意几何形状孔附近应力集中计算过程,方法计算量小、解的精度高、便于工程应用。
附图说明
图1是含功能梯度材料加强环的任意几何形状孔受远场均布载荷作用;
图2是从z平面到ζ平面映射变换;
图3是单向拉伸作用下,圆形孔周环向应力分布图;
图4是单向拉伸作用下,椭圆形孔周环向应力分布图;
图5是单向拉伸作用下,三角形孔周环向应力分布图;
图6是单向拉伸作用下,正方形孔周环向应力分布图;
图7是单向拉伸作用下,矩形孔周环向应力分布图;
图8是纯剪切作用下,圆形孔周环向应力分布图;
图9是纯剪切作用下,椭圆形孔周环向应力分布图;
图10是纯剪切作用下,三角形孔周环向应力分布图;
图11是纯剪切作用下,正方形孔周环向应力分布图;
图12是纯剪切作用下,矩形孔周环向应力分布图。
具体实施方式
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
本发明的一种减小复杂几何形状孔附近应力集中的方法为,对于复杂几何形状孔,在孔周焊接衬入杨氏模量沿法向连续递增变化的功能梯度材料加强环,此加强环的形状与孔的形状一致。
此复杂几何形状孔包括圆形孔、椭圆形孔、三角形孔、正方形孔和矩形孔。但并不限于此5种几何形状。功能梯度材料加强环(或简称加强环)是由功能梯度材料制成的加强环,文中所述的加强环或功能梯度材料加强环均是相同含义。本发明通过在孔周衬入杨氏模量沿法向连续递增变化的功能梯度材料加强环,可以减小孔附近应力集中。
为了更好的分析功能梯度材料加强环对任意几何形状孔附近应力集中的影响,对功能梯度材料加强环对任意几何形状孔附近应力集中影响进行建模,下面对建模过程的具体过程进行详细描述,如下:
如图1所示,含任意几何形状孔的无限大均质材料板,在孔周衬有功能梯度材料加强环,板和环所占的域分别用Ωp和Ωg表示,环的内外边界视为封闭光滑曲线Γ0和ΓN。功能梯度材料加强环内的杨氏模量沿着孔的法线方向连续变化,分析板在远场均布载荷作用下孔周应力场。利用分层均匀化方法将材料参数连续变化的功能梯度材料环分解为N层厚度相等的薄环,每个薄环分别用Ω(1),Ω(2),…Ω(j),…Ω(N)表示。当薄环层数N取得足够大时,任意一个薄环Ω(j)内的材料参数可近似视为常数。
对于各向同性均质材料,在直角坐标系(x,y)中,二维问题应力分量的复变函数表达式为:
其中,σx、σy、σxy为应力分量,和ψ(z)为复势函数,i表示虚数单位。
平面问题的场方程和边界条件复变函数表达式为:
其中,Xn、Yn、u、v分别表示边界上力和位移分量,此外对于平面应力问题κ=(3-v)/(1+v),G=E/(2(1+v));平面应变问题,须将式中的E换为E/(1-v2),v换为v/(1-v),E和v分别表示相应计算域内的弹性模量和泊松比。
为了方便求解,从z平面到ζ平面单位圆映射变换,引入如下保角变换函数:
其中,R>0,mn为实常数,-1<mn·n<1。表1列出了几种常见几何形状孔的映射函数系数mn的值。同时由于参数R只影响孔的尺寸,对于本发明所研究的无限大板而言,它不会影响应力分布。因此,文中取R=1。值得说明的是,在数学中各种几何形状孔的映射函数已被广泛研究,所以几乎对于任意几何形状的孔都可以通过查找文献手册来确定上式中的R和mn(参见文献Complex variable methods in elasticity)。限于篇幅,此处不再一一列举。
表1不同几何形状孔的映射函数系数
利用上述保角变换函数,在z平面的光滑边界Γ0和ΓN将被映射到ζ平面的同心圆环内外边界L0和LN,其中圆环的内外半径分别为r0,rN,如图2所示。与此同时,原z平面的N层均质薄环Ω(j)也将转变为ζ平面的N层同心圆环S(j)。任意几何孔形状都能转换为ζ平面的同心圆环。
为了便于表示,引入和ψ(ζ)=ψ[ω(ζ)],这样在ζ平面式(1)-(4)可以重新表示为:
σθ+σr=4Re[Φ(ζ)], (6)
其中,Ψ(ζ)=ψ′(ζ)/ω′(ζ)。
在ζ平面中的每个圆环S(j)内,复势函数可以表示为,
其中,j=1,2,…N,为未知系数。
在板域Sp中的复势函数可以表示为:
其中,B1,B2为取决于远场外载荷的已知常数
在j=1环的内周(即孔边),由于孔边自由,所以有
Xn=Yn=0. (14)根据上式,式(8)可以重新写为:
其中,ζ=r0σ=r0eiθ。
由于式(15)中为沿圆周L0的连续函数,所以它可展开为傅里叶级数:
其中:
将式(10)和(11)代入式(15),并且比较方程两边σ±k同次幂的系数得
其中k=1,2,…M。k的取值在以上公式中理论上是到无穷大,但是在现有计算中,需要后续联立方程组求解应力,因此要取有限值M。此M是设定值,其取值与计算精度有关,在实际模拟计算应力时,通常取到20即可。
另一方面,ζ平面中任一层j和j+1环间的界面上,有应力和位移连续条件:
u(j)=u(j+1),v(j)=v(j+1). (20)
基于式(8)和(9),式(19)和(20)可以重新表示为
其中,ζ=rjσ=rjeiθ,而连续函数同样可以展开为傅里叶级数:
其中:
与前述过程类似,将式(10)-(13)代入式(21)和(22),可以获得4N×M个线性方程。将这些方程与式(17)和(18)联立构成一组含有4N×M+2×M个方程的线性方程组,此方程组中同时含有4N×M+2×M个未知系数 和所以通过求解该方程组可以完全确定所有未知系数,进而可以获得每一层环内的所有复势函数。最后,利用式(6)和(7)可以求得环和板内的应力场。
从以上应力场分析过程可得,利用保角变换技术及分层均匀化方法,简化了含功能梯度材料加强环的任意几何形状孔附近应力集中计算过程,方法计算量小、解的精度高、便于工程应用。
实施例
对于不同几何形状孔周的FGM加强环,其法向变化的杨氏模量都可统一的转换为映射平面ζ中的径向变化情况,本实施例中,取杨氏模量为以下递变规律
其中,Ep表示板内的杨氏模量。梯度指数值n的取值范围为0<n≤5。此外,假设FGM加强环内的泊松比保持不变,取作νg=νp=0.3。为了讨论杨氏模量的影响,分别取四个不同的梯度指数值n=0,0.5,2和5。其中,n=0为孔周不含FGM加强环的一个特殊情况。
另外,对于加强环的厚度,从理论上来说,加强环的厚度越大,则其减小应力集中的效果越好,但是在实际工程应用中,加强环的厚度有一定的工程要求(如孔大小的要求),至少加强环不能填满整个孔。经试验,映射平面ζ中加强环厚度t的范围取为0.05rN≤t≤0.2rN。
映射到直角坐标中,加强环的厚度视孔的形状而定。对于圆孔,若孔半径为r,则加强环厚度t的范围取为0.05r≤t≤0.2r。
对于椭圆孔,若椭圆孔长短轴分别为a和b,则加强环厚度t的范围取为0.05b≤t≤0.2b。
对于三角形孔,若三角形孔边长为a,则加强环厚度t的范围取为0.05a≤t≤0.2a。
对于正方形孔,若正方形孔边长为a,则加强环厚度t的范围取为0.05a≤t≤0.2a。
对于矩形孔,若矩形孔长短边分别为a和b,则加强环厚度t的范围取为0.05b≤t≤0.2b。
在本实施例中,取映射平面ζ中加强环的厚度t=rN/6,即半径比r0/rN=5/6。由于单向拉伸及纯剪切这两种载荷是工程领域最基本、最典型的外力作用类型。本实施例中分别试验在这两种作用力下各几何形状孔周环向应力分布情况。
单向拉伸作用下孔周环向应力,利用公式(6)和(7)求得的应力场包含r和θ两个变量,对于孔周,只要再令r=r0,即可获得孔周应力随θ的变化。此时,单向拉伸作用下各种形状孔周的环向应力分布如图3至图7所示。图3是单向拉伸作用下,圆形孔周环向应力分布图;图4是单向拉伸作用下,椭圆形孔周环向应力分布图;图5是单向拉伸作用下,三角形孔周环向应力分布图;图6是单向拉伸作用下,正方形孔周环向应力分布图;图7是单向拉伸作用下,矩形孔周环向应力分布图。由图可以发现,首先对于不同形状的孔,应力分布特点明显不同。椭圆孔和矩形孔周的应力集中分别高于圆孔和正方形孔,而其中应力集中最为显著的是三角形孔。其次由于载荷的对称性,这些形状孔周的应力分布都关于x轴对称。另外,最为重要的是梯度指数对于所有形状孔周的应力分布都具有很大影响,当指数值增加时,应力集中显著降低。
远场纯剪切载荷作用下,孔周环向应力与ζ平面的应力存在角度对应关系,如图2中的角度λ和θ,通过查找ζ平面中的θ角度,即可确定孔周对应角度λ的应力。图8-12反映了在远场纯剪切载荷作用下不同几何形状孔周的环向应力分布图。图8是纯剪切作用下,圆形孔周环向应力分布图;图9是纯剪切作用下,椭圆形孔周环向应力分布图;图10是纯剪切作用下,三角形孔周环向应力分布图;图11是纯剪切作用下,正方形孔周环向应力分布图;图12是纯剪切作用下,矩形孔周环向应力分布图。与单向拉伸情况不同,此时的应力分布不再关于x轴对称,但是对于圆孔和正方形孔,应力分布转为关于λ=135°和315°方向对称。同时也发现,纯剪切作用下梯度指数对于不同形状孔周的应力分布仍有较大影响,应力集中随着指数值的增加明显降低。
表2给出了图3-12中不同几何形状孔周的最大应力值及所在位置。根据表中数据可以清晰看出,n=5时的最大应力值明显小于n=0(即不含FGM加强环的均质开孔板的结果)。此结论表明通过合理选择加强环内法向杨氏模量的递变规律,可以有效缓解各种几何形状孔附近的应力集中。
表2各几何形状孔在不同梯度指数中最大环向应力值及所在位置
结果表明,当杨氏模量的梯度指数值增加时,各种形状孔周的应力集中都显著降低。由此可以得出,通过在孔周衬入FGM加强环并合理选择加强环内材料参数的递变规律,可以有效缓解任意几何形状孔附近的应力集中。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变型,这些改进和变型也应视为本发明的保护范围。
Claims (6)
1.减小复杂几何形状孔附近应力集中的方法,其特征是,对于复杂几何形状孔,在孔周焊接衬入杨氏模量沿法向连续递增变化的功能梯度材料加强环,此加强环的形状与孔的形状一致。
2.根据权利要求1所述的减小复杂几何形状孔附近应力集中的方法,其特征是,复杂几何形状孔包括圆形孔、椭圆形孔、三角形孔、正方形孔和矩形孔。
3.根据权利要求1所述的减小复杂几何形状孔附近应力集中的方法,其特征是,在映射平面ζ中,加强环厚度t的范围取为0.05rN≤t≤0.2rN,其中rN为ζ平面几何孔半径。
4.根据权利要求1所述的减小复杂几何形状孔附近应力集中的方法,其特征是,在映射平面ζ中,杨氏模量递变规律公式为:
其中,Ep表示板内的杨氏模量,r为ζ平面中半径。
5.根据权利要求4所述的减小复杂几何形状孔附近应力集中的方法,其特征是,梯度指数值n的取值范围为0<n≤5。
6.根据权利要求1所述的减小复杂几何形状孔附近应力集中的方法,其特征是,应力集中的计算过程包括以下步骤:
步骤S1,在z平面中,利用分层均匀化方法将杨氏模量沿法向连续变化的功能梯度材料加强环分解为N层厚度相等的薄环;
步骤S2,利用保角变换函数从z平面映射到ζ平面,原z平面的N层均质薄环转变为ζ平面的N层同心圆环;根据映射关系,建立ζ平面环及板内应力场表达式及ζ平面中的每个圆环内复势函数表达式;
步骤S3,计算在ζ平面中的每个圆环内复势函数,基于ζ平面环及板内应力场的表达式求得环和板内的应力场。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710177560.8A CN106971036A (zh) | 2017-03-23 | 2017-03-23 | 减小复杂几何形状孔附近应力集中的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710177560.8A CN106971036A (zh) | 2017-03-23 | 2017-03-23 | 减小复杂几何形状孔附近应力集中的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106971036A true CN106971036A (zh) | 2017-07-21 |
Family
ID=59329828
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710177560.8A Pending CN106971036A (zh) | 2017-03-23 | 2017-03-23 | 减小复杂几何形状孔附近应力集中的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106971036A (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030029910A1 (en) * | 2001-08-07 | 2003-02-13 | Goretta Kenneth C. | Joining of advanced materials by plastic deformation |
CN102561532A (zh) * | 2011-12-16 | 2012-07-11 | 南京航空航天大学 | 一种功能梯度泡沫混凝土保温材料及其制备方法 |
CN104028861A (zh) * | 2014-05-29 | 2014-09-10 | 哈尔滨工业大学 | 一种金属-陶瓷功能梯度材料电火花加工逐级自适应模糊控制方法 |
-
2017
- 2017-03-23 CN CN201710177560.8A patent/CN106971036A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030029910A1 (en) * | 2001-08-07 | 2003-02-13 | Goretta Kenneth C. | Joining of advanced materials by plastic deformation |
CN102561532A (zh) * | 2011-12-16 | 2012-07-11 | 南京航空航天大学 | 一种功能梯度泡沫混凝土保温材料及其制备方法 |
CN104028861A (zh) * | 2014-05-29 | 2014-09-10 | 哈尔滨工业大学 | 一种金属-陶瓷功能梯度材料电火花加工逐级自适应模糊控制方法 |
Non-Patent Citations (1)
Title |
---|
QUANQUAN YANG: "Reduction of the stress concentration around an elliptic hole by using a functionally graded layer", 《ACTA MECHANICA》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ferreira et al. | Analysis of laminated shells by a sinusoidal shear deformation theory and radial basis functions collocation, accounting for through-the-thickness deformations | |
Yeh | Vibration analysis of sandwich rectangular plates with magnetorheological elastomer damping treatment | |
Carrera et al. | Finite element analysis of structures through unified formulation | |
Sahoo et al. | A new trigonometric zigzag theory for buckling and free vibration analysis of laminated composite and sandwich plates | |
Versino et al. | C0 triangular elements based on the Refined Zigzag Theory for multilayer composite and sandwich plates | |
Sahoo et al. | A new inverse hyperbolic zigzag theory for the static analysis of laminated composite and sandwich plates | |
Katariya et al. | Numerical evaluation of transient deflection and frequency responses of sandwich shell structure using higher order theory and different mechanical loadings | |
Reddy et al. | Bending analysis of laminated composite plates using finite element method | |
Ren et al. | The finite element analysis of homogeneous and laminated composite plates using a simple higher order theory | |
Osman et al. | Free vibration analysis of laminated composite beams using finite element method | |
Ton That et al. | Nonlinear bending analysis of functionally graded plates using SQ4T elements based on twice interpolation strategy | |
Onyeka | Direct analysis of critical lateral load in a thick rectangular plate using refined plate theory | |
Mamutov et al. | Electrohydraulic forming of low volume and prototype parts: Process design and practical examples | |
Nguyen et al. | A novel finite element formulation for static bending analysis of functionally graded porous sandwich plates | |
Shi et al. | Free vibration analysis of the unified functionally graded shallow shell with general boundary conditions | |
Pavan et al. | Static, free vibrational and buckling analysis of laminated composite beams using isogeometric collocation method | |
Mantari | General recommendations to develop 4-unknowns quasi-3D HSDTs to study FGMs | |
Sinha et al. | On-demand contactless programming of nonlinear elastic moduli in hard magnetic soft beam based broadband active lattice materials | |
Dawe et al. | Buckling of composite plates and plate structures using the spline finite strip method | |
Wagner et al. | An adaptive strategy for the multi-scale analysis of plate and shell structures with elasto-plastic material behaviour | |
Zanichelli et al. | A novel implementation of the LDEM in the ansys LS-DYNA finite element code | |
Fu et al. | Multiple-step incremental air-bending forming of high-strength sheet metal based on simulation analysis | |
Rango et al. | Static and dynamic analysis of thick laminated plates using enriched macroelements | |
CN106971036A (zh) | 减小复杂几何形状孔附近应力集中的方法 | |
Nguyen-Van et al. | Large deflection analysis of plates and cylindrical shells by an efficient four-node flat element with mesh distortions |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170721 |