含边端半圆孔有限宽平板的理论应力集中系数与孔径板宽的变化关系汇总

1、有限元作业报告Ji含边端半圆孔有限宽平板的理论应力集中系数与孔径/板宽的变化关系2013年11月I目录目录.第1章1.21.31.4第2章2.1目的和意义应力集中理论应力集中系数 平面问题的应力集中 利用有限元分析. 软件介绍MSC.PATRAN2.2第3章3.13.2第4章MSC.NASTRAN基本模型和数据基本模型.基本数据.计算分析.学习体会和建议11第1章目的和意义第1章目的和意义1.1 应力集中应力集中是应力在固体局部区域内显著增高的现象。多出现于尖角、孔洞、缺 口、沟槽以及有刚性约束处及其邻域。应力集中会引起脆性材料断裂；使物体产生 疲劳裂纹。在应力集中区域，应力的最大值（峰值应力

2、）与物体的几何形状和加载 方式等因素有关。局部增高的应力值随与峰值应力点的间距的增加而迅速衰减。由 于峰值应力往往超过屈服极限而造成应力的重新分配，所以，实际的峰值应力常低 于按弹性力学计算出的理论峰值应力。1.2 理论应力集中系数工程上，用应力集中系数来表示应力增高的程度。它是峰值应力b max与不考虑应力集中时适当选取的基准应力 bn的比值，恒大于1且与载荷大小无关。峰值应力 可根据弹性力学的理论或有限元计算得到，也可根据光弹性实验或其他实验应力分 析方法测定。其中基准应力是人为规定的，取值方式不唯一，在后面的分析中给出 两种常用表达方式。本文将用a来表示理论应力集中系数，其表达式为：cm

3、ax1.3 平面问题的应力集中由弹性理论可知，对于平面问题的应力集中问题，应力分布只取决于构建的形 状以及边界上的载荷条件，而与材料的弹性常数无关，于是在实验应力分析中，利 用力学性质不同的材料制作模型，以替代真实构件进行应力分析，对平面问题是精 确的。1.4利用有限元分析早起理论应力集中系数通过复杂的计算，可以得到精确解。随着实验手段的发 展，如电测法、光弹性法、散斑干涉法、云纹法等实验手段（见实验应力分析）均 可测出物体的应力集中。随着科技的进步，计算机和有限元法以及边界元法的迅速 发展，为寻找应力集中的数值解开辟了新途径。本文利用P ATRANi行有限元方法的数值计算便是一种快速有效的方

4、法。1第2章软件介绍第2章软件介绍2.1 MSC .P ATRANPat ran 20113旳IStfwartoaiww Iw-and険 比叭 叭 tk*!询通少'pwndef Th?隔use伽ETF CcupufBOon Hih« UnhMSwmflEFar «帕削i 冰祈ifaddpnei他ucp X Rnf luuJirAboteno »iMr« 2011 WSCjE治讪All rtgHlD pwHZ” tECHi C图 2-1 MSC .P ATRAN 软件MSC.PATRA最早由美国宇航局（NASA倡导开发的,是工业领域最著名的并行框架

5、式有限元前后处理及分析系统，其开放式、多功能的体系结构可将工程设计、工程分析、结果评估、用户化身和交互图形界面集于一身，构成一个完整CAE集成环境。MSC.PATRA提供了功能全面,方便灵活的可满足各种分析精度要求的复杂有限元的建模能力。其综合全面、先进的网格划分技术，为用户根据不同的几何模型 提供了多种不同的生成和定义有限元模型工具，包括：多种网格划分器、有限元模 型的编辑处理、单元设定、任意梁截面建模、边界和载荷定义、及交互式计算结果 后处理。MSC. PATRA提供了多种计算分析结果可视化工具，帮助工程师灵活、快速地理解结构在载荷作用下复杂的行为，如结构受力、变形、温度场、疲劳寿命、流

6、体流动等。分析的结果同时可与其它有限元程序联合使用。2.2 MSC.NASTRAN作为世界CAE工业标准及最流行的大型通用结构有限元分析软件,MSC.NASTRA的分析 功能覆盖了绝大多数工程应用领域，并为用户提供了方便的模 块化功能选项,MSC.NASTRAN的主要功能模块有：基本分析模块（含静力、 模态、 屈曲、热应力、流固耦合及数据库管理等）。动力学分析模块、热传导模块、非线性分析模块、设计灵敏度分析及优化模块、超单元分析模块、气动弹性分析 模块、DMAP用户开发工具模块及高级对称分析模块。除模块化外,MSC.NASTRAN还按解题规模分成10,000节点到无限节点，用户引进时可根据自身

7、的经费状况和功 能需求灵活地选择不同的模块和不同的解题规模，以最小的经济投入取得最大效益。MSC.NASTRAN MSC勺相关产品拥有统一的数据库管理，一旦用户需要可方便地 进行模块或解题规模扩充，不必有任何其它的担心。MSC.NASTRA以每年一个小版本，每两年一个大版本的速度更新，用户可不断 获得当今CAE发展的最新技术用于其产品设计。 目前MSC.NASTRA的最新版本是 1999年发布的V70.5版。新版本中无论在设计优化、 P单元、热传导、 非线性 还是在数值算法、 性能、文档手册等方面均有大幅度的改进或突出的新增功能。第4章计算分析第3章基本模型和数据3.1 基本模型本报告建立的模

8、型为含边端半圆孔有限宽平板，并且是两边都有的对称情况。 其形状如下：h MxJb - ck<fliiK.viflwpoil - dtfiuh_g sup - EntityX兮牛心15图3-1含边端半圆孔有限宽平板模型3.2 基本数据基本数据如下：长度为 200mm宽度为100mm厚度为3mm选用钢，定义弹性模量E为200 GPa泊松比为0.3，1mm,3mm,5mm10mm, 20mm，30mm, 40mm，48mm平板的几何特性：平板的材料特性：半圆孔的半径分别选：在两边施加等大的载荷10N/mm平板被拉伸。建模的过程中，会随着孔径的变化，在孔附近撒种的数量进行变化，目的是为 了达到一

9、定的精确度，提高结果的准确性。第4章计算分析这里我们将利用模型的1/4进从图中可以看出我们要分析的模型具有对称性， 行分析，从而可以得出整体的结果。事实上如果不采用这种方法，在进行分析时得 不到结果。利用对称性，可以得到两个位移约束条件。如图，对于左面来说，没有X方向的位移；对于上面来说，没有丫方向的位移。 在施加约束约束时也要固定 Z方向的位移（否则可以发生Z方向的平动，得不出结 果来）。有趣的是这个1/4模型与研究中心圆孔有限宽平板的模型是一样的，不过由于 位移约束条件的差异，结果自然是有区别的。CD i C3 1 E3-眉 rOl.db - dtFiulL*i«4»r

10、l - d«*HtLgroup - EntityTOTWimm图4-1 1/4模型的示意图建模的过程首先是几何作图，这个模型的建立有多种方法。先画出矩形平面，在上面挖去1/4的圆，即扇形便可得到模型。也可先画出整体模型，再用Break的方法逐渐得到1/4的部分。也可建立体模型，这样就不用在后面定义厚度了。下面划分网格，网格的值设为 2,在圆弧周围撒种可以根据不同情况确定不同 的数目，这样可提高计算的精度。施加位移、应力约束条件，定义材料，赋予属性。模型建立就完成了。利用相 似的方法完成其他几种情况的建模，计算并分析结果。得到的图如下："! rilT.db-d*<iiA

11、.=twT»T«E机初'Htran 2011 1 >-S 2fl:Cv u&FTfige DeTaiJt.Al DtetiL 5ubiase Stress T&reor. .wn Mrses AJtZlVi nnliHfOrm 闊茴A l SIMif: SubC=p-/i. 0$伪&0的曲怜.17&1$191«1|-131BgggE5 92000尋那中COD.cacl>C4«+囱.4.2+DU«a.a-MO- 科斗MC打*McL!M0-| 谒期L daLCLFriw 阳沾费怨呂嗣2Mini 6

12、&+K10mnd :朋 曲fadm 佃hort H郎16皿)3国嗣1和图 4-2 R=1mmpBtranS&il仔片叶M和:制ISFtrnge 2eTaiJt ji生仙:二jbcrStr碑 T謂JS 心凶 El卜 TThnni ：S绑MdT 71 +COdW-字 r-2-rIIJ 匚 Al Slyhu ”宙民 CttplyOttfwriH. wster*!SiSiflprflSlAI *(£j«0AdgfaulL*rng6j MaKStSittODiNd USMin37t4M| <ffind 146daraJl.DarofffflhortHgM6W)3

13、信嗣甥图 4-3 R=3mm图 4-4 R=5mmi nHdh - drtiJMil/ijiZp - Er4<rPgtTMi诲FMnge 5C" 5FAULT.A2;5TaftcSubcase. Slrsse-Iersor wn M&es 曲二1口Hfortu 3F&UU.血"Sttli右 SkibCA開.口儈跑鴛riKi闻 TrSnslAhuPfllss+wo6护M叶I3 91-hKIC5 £e-Mi卜I2 E 牡 ML”1 9盼閒4DM Mm1 fis网仙Min5 6HM2<ffind I drauH.EisfcirrahortMa

14、x 1 77-003 ifflMd 2图 4-5 R=10mm了 nlF.db driai?i_Titwi»T dei/ijroup = EfPEtranSOll 1了円叶 19M：已 违FtTige Cef, A :S1a&cZjjtxiaM Stress Terser, .wn MZlCd 茴汕JtjhcCi：pl£Kriitfir5 Tranril?:Ml1 n 牡 rm：-：1 0沪M： gSMWD a9wmcS52tML|7创曲 日込D0«66 6>*MC =JS?*MC- 旷肛MO詁 止的M片 2 了甘 Moy aroMcioLirhgp

15、DCS!国Nd IISdefFMhk 1.03-hMin6 2B-X£Nd I da 31 _ M<h他 h WhH#2 0dM3 国Nd?图 4-6 R=20mmFTnga CelaJt. Aa L$tBfne Subeese Stress Tecsor von Mrs AtZl CrifciM 閔aifl/l St#llG 虫山;门忡 Tr.'isUU-JAilirtBBltaaHBa 1曲墜】1谢Ml10 -H;H;L »»»K»» SI *AH.畧冊 -Pffiria图 4-7 R=30mm(fcraJi.Dsr

16、wdhort E'dti driai?i_-<wi»T Eiviih*图 4-8 R=40mmW Hh5 _-_5 2-JI 76 ItUoIPtranSC-ll 17flDrl3 2O:0e £1Ftifige CeTaUt.h1r£ 丸予CufdM rS-'fiiJr Al SlightTratl riiVilEeiMM'5 0d*MI3 3Er+M J2创+MP牡W'斷h on血*_FrTigp JI tMi 国Nd I儁 l7?-M2®Nd I图 4-9 R=48mm冃 rO2.db: 31=1 II 回 I

17、 £2LEGENDResultl: Stress Tensor, VONM5,0+001-50+001-3,GD+DD1-2,70+001-1.80+001a.DD+DDO-&51I+DII0I JtHIW 出 55*(W13 上 D+DII 込 25M(l&1(ktW1Coord 0.2上图是对整个左面的应力分布图，可以看出在靠近应力集中地地方，应力变化很剧烈，远离应力集中的地方，应力变化不大则比较平滑。应力集中的地方是左面 上圆孔的边界处。下面计算理论应力集中系数，给出两种表达式：F匚 d是板厚冲max* = 其中 cnc“ - maxot =pcn其中FB(B-2Rdd旳21?用的关系« RF面列出表:圆孔半 径 R/mm1351020304048(T max5.928.239.289.6810.311.81984.12R/B0.020.060.10.20.40.60.80.96(T n3.333.333.333.333.333.333.333.33(T 'n3.403.553.704.175.558.3316.6783.33a1.782.472.782.903.093.545.725.231a1.742.322.512