有限宽中心圆孔板应力集中系数数值实验
基于ANSYS的有限宽板孔边应力集中分析
4 结论
圆孔附近发生了明显的应力集中现象 ,且孔径越小应力集中越明显 ,应力突然变大的趋势越快 . 参考文献 :
[1 ] 徐芝纶 . 弹性力学简明教程 : 第 2 版 [M] . 北京 : 高等教育出版社 ,1993. [2 ] 王国强 . 实用工程数值模拟技术及其在 ANSYS 上的实践 [M] . 西安 : 西北工业大学出版社 ,1999. [3 ] 刘 波 ,曹晓东 . 平板中心圆孔边应力集中的有限之分析 [J ] . 石油化工设备技术 ,2004 ,25 (5) :20～23. [4 ] 朱晓东 ,覃启东 . 基于 ANSYS 平台含圆孔薄板的应力集中分析 [J ] . 苏州大学学报 ( 工科版) ,2004 ,25 (5) :20～23. [5 ] 张胜明 . 基于有限元软件 ANSYS/ 7. 0 的结构分析 [M] . 北京 : 清华大学出版社 ,2003.
图1 带孔平板
第 1 期 张宁锋 : 基于 ANSYS 的有限宽板孔边应力集中分析 ・35 ・
析中采用八节点实体单元 PLANE82 ,单元属性设置为 Plane stress w/
thk ,弹性模量和泊松比分别为 200 GPa 和 0. 3 ,边界条件为 x = 0 , UX = 0 ; y = 0 , UY = 0 . 在板远端作用有沿 x 轴方向的 100 MPa 的均匀
( 上接第 19 页)
参考文献 :
[1 ] BARTTREPAK,辛晓英 . 虚拟红外防盗报警器 [J ] . 家庭电子 ,2001 , (10) :33～34. [2 ] 魏立君 ,韩华琦 . 1COM S40000 系列 60 种常用集成电路的应用 [M] . 北京 : 人民邮电出版社 ,1995. [3 ] 蔡凡弟 . 门窗监控远距离转发报警器 [J ] . 电子世界 ,1999 , (10) :28～31. [4 ] 赵建华 ,张荷芳 ,李 静 . 门窗监控远距离转发报警器 [J ] . 现代电子技术 ,2001 , (8) :37～39. [5 ] 李建华 . 实用遥控器原理与制作 [M] . 北京 : 人民邮电出版社 ,1996.
平板孔口应力集中的ANSYS有限元分析喻光安
平板孔口应力集中的ANSYS 有限元分析一、开孔的应力集中和应力集中系数容器开孔后使承载截面减小，破坏了原有的应力分布，并产生应力集中，而且接管处容器壳体与接管形成不连续结构而产生边缘应力，这两种因素均使开孔或开孔接管部位的局部应力比壳体的薄膜应力大，这种现象称为开孔的应力集中。
常用应力集中系数t K 来描述接管处的应力集中特性。
未开孔时的名义应力为σ，开孔后按弹性方法计算出最大应力若为max σ，则弹性应力集中系数的定义为σσ/max t =K 。
下面以两向拉伸应力作用下的平板为例，利用ansys 有限元分析得出平板的受力情况，求出t K 的值，并与理论解作分析比较。
二、两向拉伸应力作用下平板的理论分析。
如图所示为无限平板受21σσ≥两向拉伸应力作用，由弹性力学的知识可得A 、B 两点的应力为213σσσ-=A ,12-3σσσ=B比较可得 1211max t -3σσσσσ==K 当σσσ==21时 2-31211max t ===σσσσσK 当σσ=1，σσ212=时 5.20.5-31max t ===σσσσσK三、建立模型。
设有中心带圆孔的长方形平板，板的厚度为0.05m ，圆孔的孔半径r=0.05m,材料的弹性模量E 为2e11，泊松比为0.3，板长度为30m ，宽度为230m ，m N /401=σ，m /202N =σ2σ 平板开小圆孔的应力集中取四分之一薄板，模型如下：对模型进行网格划分并施加荷载，并对圆孔周围的区域进行局部网格划分，划分后的模型。
，Ansys计算后的应力云图如下：由应力云图可知，圆孔处最大应力m N /27.100max =σ 验证公式当m /401N ==σσ，m N /20212==σσ时 50675.24027.1001max t ≈==σσK ，基本符合理论解2.5。
有限大含孔板的应力集中分析
有限大含孔板的应力集中分析作者姓名：李姝指导教师：李英梅副教授、王建祥教授单位名称：理学院专业名称：工程力学东北大学2013年6月Theoretical Analysis of the Stress Concentration in Finite-size Plates with HolesBy Li ShuSupervisor: Associate Professor Li YingmeiProfessor Wang JianxiangNortheastern UniversityJune 2013毕业设计（论文）任务书- i -东北大学毕业设计（论文）摘要有限大含孔板的应力集中分析摘要开孔板被广泛应用于机械制造、土木工程等领域中，由于开孔而产生的应力集中现象会降低结构的承载能力，并导致工程结构在开孔附近区域发生破坏。
目前，对无限大含孔板的理论研究已经比较成熟。
但是，工程结构多为有限尺寸板。
为有效预防开孔有限尺寸板由应力集中导致破坏的现象，本文研究含圆孔或椭圆孔有限板的应力分布等问题，进行了相应的理论分析和数值模拟。
本文采用复变函数法，设定待定复函数，使其预先满足所有的基本方程，运用最小二乘法使其满足物体的边界条件，最终求解应力集中系数。
该方法的收敛性很好，且随着Laurent级数展开项数的增多，本文解的精度会增高。
此外，本文给出了含中心圆孔的细长板在面内剪切载荷作用下，计算应力集中系数的经验公式；含中心椭圆孔（长短轴比为1.2）的细长板在拉伸载荷作用下，计算应力集中系数的经验公式。
本文采用最小二乘法分析含非圆形孔口边沿的应力集中问题，该方法比常用的保角映射法简单，而且也满足工程实际的精度要求。
关键词：含孔有限大板，应力集中系数，最小二乘法- ii -东北大学毕业设计（论文）AbstractTheoretical Analysis of the Stress Concentration in Finite-size Plates with HolesAbstractPlates with holes are widely used in machinery manufacturing, civil engineering and other fields. Stress concentration resulting from holes can reduce the carrying capacity of the structure and lead engineering structures’failure occurring in the vicinity of holes. Currently, there is abundant research about the theoretical analysis of the stress concentration in infinite plates with holes. However, in the application of real engineering cases, most plates are finite. To effectively prevent the destruction of finite-size plates with holes, we study the stress field in finite plates with circular holes or elliptical holes and study the theoretical analysis and numerical simulation analysis.This paper uses the complex expansion method together with least-square method. Set the pending complex function that meets all of the basic equations in advance, use the least-square method so as to satisfy the boundary conditions of the object, and ultimately compute the stress concentration factor. When the number of terms of Laurent-expansion increases, the accuracy of the results is more satisfying. In addition, this paper presents empirical formulas of stress concentration of plates with a central circular hole in the condition of in-plane shear. Simultaneously, this paper also presents empirical formulas of stress concentration of plates with a central elliptical circular hole (the ratio is 1.2) in the condition of simple tension.This paper uses the complex expansion method together with least-square fitting to compute the stress concentration factor in the plates with a central non-circular hole. Compared with conformal mapping method, this method is more simple. This method also meets the accuracy requirements of engineering practice.Keywords: finite plates with holes, stress concentration factor, least-square method- iii -目录任务书 (i)摘要 (ii)Abstract (iii)第1章绪论 (1)1.1 应力集中问题概述 (1)1.2 开口有限板的研究背景和意义 (1)1.3 开口有限板的研究方法 (2)1.3.1 有限元法 (2)1.3.2 实验方法测试 (3)1.3.3 弹性力学分析方法 (3)1.4 论文研究内容 (4)第2章开口有限板应力集中理论分析方法 (5)2.1 应力函数的复变函数表示 (5)2.2 应力的复变函数表示 (6)2.3 合力表达式 (7)2.4 多连通区域里的ϕ和ψ (8)2.5 含孔各向同性有限板应力表达式的确定 (10)2.6 最小二乘法求解待定系数 (12)2.6.1 含圆孔、椭圆孔有限板在面内拉伸载荷作用下 (12)2.6.2 含圆孔、椭圆孔有限板在面内剪切载荷作用下 (12)2.7 本章小结 (13)第3章开口有限板应力集中结果分析 (15)3.1 含圆孔有限板在面内拉伸载荷作用下 (15)3.1.1 级数收敛性分析 (15)3.1.2 计算结果 (16)3.1.3 Heywood经验公式 (17)3.1.4 算例分析 (18)3.1.5 小结 (21)3.2 含圆孔有限板在面内剪切载荷作用下 (21)3.2.1 级数收敛性分析 (21)3.3.2 计算结果 (22)3.2.3 本文所提经验公式 (24)3.2.4 算例分析 (25)3.2.5 小结 (29)3.3 含椭圆孔有限板在面内拉伸载荷作用下 (29)3.3.1 级数收敛性分析 (29)3.3.2 计算结果 (30)3.3.3 本文所提经验公式 (31)3.3.4 算例分析 (33)3.3.5 小结 (35)3.4 含椭圆孔有限板在面内剪切载荷作用下 (35)3.4.1 级数收敛性分析 (36)3.4.2 计算结果 (37)3.4.3 算例分析 (38)3.4.4 小结 (42)第4章总结与展望 (43)4.1 论文总结 (43)4.2 工作展望 (44)参考文献 (45)结束语 (47)附录A：中文译文 (49)附录B：英文原文 (63)第1章 绪论1.1 应力集中问题概述在构件强度设计中所用的基本公式，一般只适用于等截面的情况。
板中孔应力集合ANSYS有限元分析
3.2.4 设置材料属性
执行 Main
Menu→Preprocessor→Material
Models→Structural→Linear→Elastic→Isotropic→输入实常数（在 EX
框中输入 200000，在 PRXY 框中输入 0.3）→OK，如图 3-2 所示。
图 3-2
对全部高中资料试卷电气设备，在安装过程中以及安装结束后进行高中资料试卷调整试验；通电检查所有设备高中资料电试力卷保相护互装作置用调与试相技互术关，系电通，力1根保过据护管生高线产中0不工资仅艺料可高试以中卷解资配决料置吊试技顶卷术层要是配求指置，机不对组规电在范气进高设行中备继资进电料行保试空护卷载高问与中题带资2负料2，荷试而下卷且高总可中体保资配障料置2试时32卷，3各调需类控要管试在路验最习；大题对限到设度位备内。进来在行确管调保路整机敷使组设其高过在中程正资1常料中工试，况卷要下安加与全强过，看度并25工且52作尽22下可护都能1关可地于以缩管正小路常故高工障中作高资；中料对资试于料卷继试连电卷接保破管护坏口进范处行围理整，高核或中对者资定对料值某试，些卷审异弯核常扁与高度校中固对资定图料盒纸试位，卷置编工.写况保复进护杂行层设自防备动腐与处跨装理接置，地高尤线中其弯资要曲料避半试免径卷错标调误高试高等方中，案资要，料求编试技5写、卷术重电保交要气护底设设装。备备置管4高调、动线中试电作敷资高气，设料中课并技3试资件且、术卷料中拒管试试调绝路包验卷试动敷含方技作设线案术，技槽以来术、及避管系免架统不等启必多动要项方高方案中式；资，对料为整试解套卷决启突高动然中过停语程机文中。电高因气中此课资，件料电中试力管卷高壁电中薄气资、设料接备试口进卷不行保严调护等试装问工置题作调，并试合且技理进术利行，用过要管关求线运电敷行力设高保技中护术资装。料置线试做缆卷到敷技准设术确原指灵则导活：。。在对对分于于线调差盒试动处过保，程护当中装不高置同中高电资中压料资回试料路卷试交技卷叉术调时问试，题技应，术采作是用为指金调发属试电隔人机板员一进，变行需压隔要器开在组处事在理前发；掌生同握内一图部线纸故槽资障内料时，、，强设需电备要回制进路造行须厂外同家部时出电切具源断高高习中中题资资电料料源试试，卷卷线试切缆验除敷报从设告而完与采毕相用，关高要技中进术资行资料检料试查，卷和并主检且要测了保处解护理现装。场置设。备高中资料试卷布置情况与有关高中资料试卷电气系统接线等情况，然后根据规范与规程规定，制定设备调试高中资料试卷方案。
理论应力集中系数的有限元求法
图 8 路径上的应力分布
式中 : s———路径长度 ; σds———上述应力沿路径积分值 。
从而该带孔平板在拉伸条件下的理论应力集中
系数 KT 为 :
KT
σ
=
m ax
σ
0
= 2. 8.
173 586
×107 ×106
= 2.
53
·18·
航空兵器 2005年第 3期
图 10 受拉伸的带轴肩轴结构
5 结 论
通过有限元计算以及与传统方法的比较可知 , 利用有限元计算方案求解工程结构的理论应力集中 系数具有不受材料 、结构形状等条件的限制 , 不但 省时 、省力 、节约费用 , 而且通用 、可靠 、精确 , 是 一种非常值得推广的数值分析方法 。
参考文献 :
[ 1 ] 姚卫星. 结构疲劳寿命分析 [M ]. 北京 : 国防工业出版 社 , 2003.
[ 5 ] 王国强. 实用工程数值模拟技术及其在 ANSYS上的实 践 [M ]. 西安 : 西北工业大学 , 1999.
同时该应力集中现象是可设计的 , 在结构的设计初 期 , 为提高产品的疲劳寿命 , 必须尽可能地避免或 降低构件的应力集中 。
随着计算机计算速度的飞速发展 , 以及有限元 计算方法的建立和完善 , 使得通过数值计算的方法 求取理论应力集中系数 ( KT )成为了可能 。有限元法 能够十分精确地逼近工程上的实际结构 , 且计算结 果准确度高 , 使得有限元法在很多复杂结构的强度 、 刚度 、模态等理论分析领域大显身手 。
应力集中系数图 。
尽管工程手册中列举了大量的图表数据 , 但是
由于受实际结构 、受力状态和结构形式的变化 , 手
册中的数据图表并不能全部地反映实际情况 。
实验4孔边应力分布及应力集中系数的测定
实验四 静态多点应变测量——孔边应力分布及应力集中系数的测定一.实验目的1. 掌握静态多点应变测量的方法；2. 学习拟订实验加载方案；3. 学习数据处理及回归分析方法；4. 测定孔边应力分布及应力集中系数。
二.实验设备和器材1. 已贴片的带孔板状拉伸试件；2. DH3818静态电阻应变仪及电脑；3. 万能材料试验机；4. 游标卡尺；三.实验方法和步骤1. 加载方案的制订（a ） 测量试件圆孔处的试件宽度、厚度及圆孔的直径；（b ） 测量每片应变计的贴片位置到试件边缘的距离；（c ） 实验时分四级加载，所加最大载荷不能超过材料的屈服强度，最大载荷P max =A ×σs ×80%/K ，其中A=(b-d)t ，σs =235MPa 为材料的屈服强度，K 为估计的孔边应力集中系数，可以按2.4选取。
初载荷P 0=选用量程×10%。
每级载荷的增量为： 40max P P P -=∆ 以计算结果取整为准。
例：b=60.00mm ，d=20.00mm ，t=6.00mm,σs =235MPa,则A=(b-d)t=(60-20)×6.00=240.00mm 2，P max =A ×σs ×80%/K=240.00×235×0.8/2.4=18.8kN,试验机的选用量程为60kN ，初载荷P 0=选用量程×10%=6kN, 每级载荷的大小为： kN P P P 2.3468.1840max =-=-=∆ 则每级载荷的增量为3kN 。
2. 依次将每片电阻应变计接入DH3818应变仪的AB 桥臂，在公共补偿桥路的BC桥臂中接入1片温度补偿应变计（可在另一已贴片的试件中选择1片）。
3. 按实验三操作DH3818应变仪的方法将应变仪调平衡。
4. 按拟订的加载方案逐级加载，在载荷状态栏中输入载荷值并采集数据，将数据转换为EXCEL 文本形式并保存，然后卸载到零。
有限元作业一带孔平板圆孔应力集中分析。二内六角扳手静力分析。三弹簧质量阻尼系统受谐载荷响应分析
学号：S2*******程序版本：ANSYS 10作业一：带孔平板圆孔应力集中分析问题描述：如右图所示，一个承受单向拉伸的无限大板，在中心位置有一个小圆孔。
材料属性为弹性模量a P E 6101⨯=，泊松比为0，拉伸的均布载荷Pa p 7101⨯=，平板厚度mm t 1=。
ANSYS 10 分析步骤：1. 定义工作文件名：Utility Menu>File>Change Jobname>输入Plate>OK2. 定义工作标题：Utility Menu>File>Change Title>输入The Ansysis of Plate withsmall Circle>OK3. 重新显示：Utility Menu>Plot>Replot4. 设置系统单位制：命令输入窗口，输入命令/UNITS,SI 并回车5. 设置计算类型：ANSYS Main Menu>Preferences>选Structural>OK6. 选择单元类型：ANSYS Main Menu>Preprocessor>Element Type>Add/Edit/Delte>Add>选Solid Quad 4node 42>OK>Options>K3:Plate Strs w/thk>OK>Close7. 定义实常数：ANSYS Main Menu>Preprocessor>Real Constants> Add/Edit/Delte>Add>OK>在THK 输入1 >OK>Close8. 定义材料特性：ANSYS Main Menu>Preprocessor>Material Props> Material Models>双击选Structural>双击Linear>双击Elastic>双击Isotropic>在EX 输入1e6，PRXY 输入0>OK>点击“X”关闭9. 生成平面方板：ANSYS Main Menu>Preprocessor>Modeling>Create>Areas>Rectangle>By2 Corners>输入WP X:0 WP Y:0 Width:10 Height:10 >OK10. 生成圆孔平面：ANSYS Main Menu>Preprocessor>Modeling>Create>Areas>Circle>SolidCircle>输入WP X:5 WP Y:5 Radius:1 >OK11. 布尔运算生成孔：ANSYS Main Menu>Preprocessor>Modeling>Operate>Booleans>Subtract>Areas>选方板>点OK(Multi Entities 窗)>OK(Subtract Areas 窗) 选方板>点NEXT>OK(Multi Entities 窗)>OK(Subtract Areas 窗)12. 网格划分：ANSYS Main Menu>Preprocessor>Meshing>MeshTool>Size Control:Global>set>在NDIV 输入6>OK> MeshTool> Mesh>Pick All>Close(Warning)> Close(MeshTool)13. 施加约束：(1): ANSYS Main Menu>Solution>Define Loads>Apply>Structural>Displacement>OnNodes>点选结构左侧所有节点>OK>Lab2 DOFs:UX,VALUE:0>OK (2):ANSYS Main Menu>Solution>Define Loads>Apply>Structural>Displacement>On Nodes>点选结构左下侧（0,0）节点>OK>Lab2 DOFs:UX,UY,VALUE:0>OK14. 施加均布载荷：ANSYS Main Menu>Solution>Define Loads>Apply>Structural>Pressure>OnLines>点选结构右侧所有节点>OK>VALUE:-1E7> OK>Close15. 分析计算：ANSYS Main Menu>Solution>Solve>Current LS>OK>Yes>Close>关闭文字窗16. 结果显示：ANSYS Main Menu>General Postpro>Plot Results>Deformed Shape>点选Def+undeformed>OK> Plot Results>Contour Plot>Nodal Solu>选Stress 选von Mises stress>Def+undeformed Model>OK17. 退出系统图1 带孔平板变形形状的结果图2带孔平板应力分布的结果作业二：内六角扳手静力分析如右图所示，截面宽度为10mm的内六角扳手，在手柄端部施加扭转力100N，以及垂直向下的力20N，分析在两种载荷的作用下扳手的应力分布。
不同板宽的孔边的应力集中问题
不同板宽的孔边的应力集中问题1 选题目的：对于如图所示的平面圆孔的孔边问题，通过数值实验的方法研究不同板宽的孔边应力集中问题，与弹性力学的解析解进行比较。
给出应力集中系数与相对孔径尺度的关系。
图一 不同板宽的孔边的应力集中问题2 背景：就无限大板宽的孔边应力集中问题，有以下弹性力学的解析解：004020002020040020020200200390)2321(90y )31)(1(2sin 2)31(2cos 2)1(2)31)(1(2cos 2)1(2422242222q R r rR rR q rR rR q rR q rR q rR rR q rR q r r r ===++==+--==+-+=--+-=），（）（分布：轴上有在孔边的θσθσσθττθσθσθθθθθθ3 数值分析我们定义板宽和孔径的相对尺度的特征参数： 0R B=ξ进行研究，具体取值如表：结果如图：图2 30=ξ时的应力分布R B =ξ30 24 20 10 6图3 24=ξ时的应力分布图4 20=ξ时的应力分布图5 10=ξ时的应力分布图6 6=ξ时的应力分布4 应力集中系数如表：5 在matlab 中划出曲线：0R B =ξ30 24 20 10 6 0maxq k σ=3.0143.0323.0783.2134.0806 结论 随着0R B =ξ增大，k 值减小。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
有限宽中心圆孔板应力集中系数数值实验
冯美生，张红珠
辽宁工程技术大学力学与工程科学系，辽宁阜新 （123000）
摘 要：在anays 平台上，采用有限元方法对拉伸有限宽中心圆孔板应力集中问题进行了数值实验，定义了应力集中的特征参数，定量分析特征尺度的变化规律，研究应力集中系数与孔径尺度的关系见图3，并与解析解比较，给出了解析解的适用范围。
关键词: 应力集中，应力集中系数，圆孔，特征尺度，数值实验
1 引言
受力的弹性平面板具有小孔，则孔边的应力将远大于无孔时的应力，也远大于距孔稍远处的应力，这种现象称为孔边应力集中。
应力集中现象是局部现象。
在几倍于孔径以外，应力几乎不受孔的影响，应力的分布情况以及数值都与无孔时相同。
一般来说，集中的程度越高，集中的现象越是局部性的，就是说应力随着与孔的距离增大而越快的趋进于无孔时的应力。
应力集中的程度，首先与孔的形状有关，一般来说，圆孔孔边的集中程度最低。
另外集中系数还与相对孔径尺度有关。
基于ansys 平台，通过数值试验的方法，研究不同板宽，不同孔径时的孔边应力集中问题，并与弹性力学的解析解进行比较，研究应力集中系数与孔径尺度的关系。
2 实例分析
2.1力学模型及假设
如图1所示，平面带孔平板，孔位于板正中，假设板为各向同性完全弹性，板左端固定，右端受均布荷载q 0=10N/mm 作用，长为200mm ，厚为10mm ，泊松比为0.3，E=2.1×1011Pa,板宽和孔径变化，数值实验其应力集中时的特征参数。
定义一个描述板宽与孔径的相对尺度的特征参数，0
B R ε=，定义应力集中系数max 0k q σ=，其中B 为板宽，R 0为孔半径，max σ为孔边最大应力，q 0为均布荷载。
2.2数值实验
在ansys 平台上变化各种ε值，计算相应的k 值，进行相应的数值研究。
整个过程采用
APDL 语言[1]，基于命令流进行参数化处理。
正式试验前，已经用固定板宽和固定孔径的有限元模型在ansys 上进行了严格的精度计算和收敛性效验，网格划分的精度足够高，误差小于1％。
实验一：B=100mm ，R 0分别为50，40，35，30，25，20，15，10，7，6，5，4。
实验二：B=75mm ，R 0分别为35，30，25，20，15，10，7，6，5，4。
实验三：B=50mm ，R 0分别为25，20，15，10，7，6，5，4。
下面给出B=50mm 时，计算的有限元网格图及其ε=3.16时的应力云图，如图2所示。
图2 有限元网格图和应力云图 Fig.2 Meshing picture and stress envelope by FEM
对应不同B 值时，分别计算不同ε值时的k 值，得到三组ε～k 关系曲线，见图3所示。
由图可见，当2<ε≤4.10时，应力集中系数经历一个由无穷大趋向于有限值的过程。
当ε<2时，应力集中系数趋向于无穷大,此时集中系数无意义，因为此时的孔径大于等于板宽，工程上很少见。
ε≥4.10时，即k ≤3.16时，应力集中系数的数值解是稳定的，因为此时的孔径会越变越小，工程上也没这种情况。
进一步由ε～k 图得：三组实验的稳定k 值均满足2.93≤ k ≤3.16，此时的孔径特征尺度在工程上很常见，此时变化的板宽与变化的孔径之间的尺度关系满足4.10≤ε≤9.46，研究这样的ε才最有意义。
3 误差分析
无限大板宽情形下的解析解[2]：max (21)a q b
σ=−，当孔为圆孔时，a=b ，此max 3q σ=，应力集中系数。
数值解与解析解的误差上限为：
3k =max 3100% 5.33%3
k η−=×=， 此时的特征尺度稳定于[4.10,9.46]。
20
40
60
80100
120
图3 不同B ，R 0值时ε～k 关系曲线
Fig.3 ε～k curve with different B and R 0
4 结论
（1）文中提出了应力集中时的特征尺度，通过有限元数值实验验证，对有限宽中心圆孔板的集中系数在某一特征尺度上是稳定的，误差上限满足工程上的需要。
（2）由ε～k 关系曲线看出，文章定义的尺度效应对应力集中系数的影响有一极限值
4.10。
（3）数值实验的过程证明，一定板宽时孔径大小对应力集中系数的影响很小。
参考文献
[1]王国强.实用工程数值模拟技术及其在ansys 上的实践.西安:西北工业大学出版社,2001.
[2]徐芝纶.弹性力学(第二版).北京:人民教育出版社,1982.
Numerical Experiment of Stress Concentration Factors for
Finite-width Center Plate with a Circular Hole
Feng Meshing, Zhang Hongzhu
(Department of mechanics and engineering science of Liaoning Technical
University,Fuxin ,123000,China)
Abstract
Using finite element method on ansys9.0 platform makes numerical experiment with stress concentration problem to finite-width plate with a circular hole in center. Characteristic parameter of stress concentration problem is defined,and variation regular of characteristic scale is quantitative analysed. By comparing numerical solution of concentration factors to analytic solution ,it obtains the application scope of analytic solution.
Keywords: stress concentration; stress concentration factor; circular hole; characteristic scale; numerical experiment
作者简介：冯美生，男，1980年生。
山西晋中人，硕士研究生，从事岩石力学系统稳定性理论、环境岩石力学数值模拟,流变岩石力学方向研究，E-mail ：daokedaoke@ 。