ANSYS上机算例孔板应力集中问题
- 格式:pptx
- 大小:322.42 KB
- 文档页数:16
【静⼒分析】AnsysWorkBench验证应⼒集中系数应⼒集中是在零件上普遍存在的现象,有各种各样的应⼒集中被计算出来,这次使⽤板中央有孔的零件来计算应⼒集中系数。
新建静⼒分析(Static Structural)导⼊模型后,在DM中打开模型,右键Import1 > Generate,⽣成模型。
使⽤Manual(⼿动⽅式)⽣成中⾯。
点击Tools > Mid-Surface > 选择模型的两个⾯ > Apply(系统默认使⽤⼿动⽅式);右键Midsurf1 > Generate,⽣成中⾯。
使⽤⾃动⽅式⽣成中⾯。
⽣成中⾯可以使问题简化,求解更加快速。
点击Tools > Mid-Surface > Selection Method :Aotumatic;在Min和Max Threshold处分别填1和3,代表最⼩壁厚为1mm,最⼤壁厚为3mm;Find Face Pairs Now:Yes,⽴即查找⾯对,此时在Face Pairs处出现数字1,表⽰已经⾃动找到⼀个⾯对。
右键Midsurf2 > Generate,⽣成中⾯。
退出DM,进⼊静⼒分析。
划分⽹格。
设置⽹格尺⼨右键Mesh > Insert > Sizing > 选中中⾯ > Apply;Element Size:0.5mm。
设置映射⽹格右键Mesh > Insert > Face Meshing > 选中中⾯ > Apply右键Mesh > Generate,⽣成⽹格。
设置边界条件设置固定约束点击A5 > Supports > Fixed Support > 选中中⾯的⼀个短边 > Apply设置⼒点击A5 > Loads > Force > Magnitude :100N;Direction:选择⼀个长边,将板拉伸。
平板孔口应力集中的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。
带孔平板的应力分布及应力集中系数的计算一、问题重述计算带孔平板的应力分布及应力集中系数。
二、模型的建立与计算在ANSYS中建立模型,材料的设置属性如下分析类型为结构(structural),材料为线弹性(Linear Elastic),各向同性(Isotropic)。
弹性模量、泊松比的设定均按照题目要求设定,以N、cm为标准单位,实常数设置中设板厚为1。
采用solid 4 node 42板单元,Element Behavior设置为Plane strs w/thk。
建立模型时先建立完整模型,分别用单元尺度为5cm左右的粗网格和单元尺度为2cm左右的细网格计算。
然后取四分之一模型计算比较精度,为了使粗细网格单元数与完整模型接近,四分之一模型分别用单元尺度为2.5cm左右的粗网格和单元尺度为1cm左右的细网格计算。
(1) 完整模型的计算①粗网格单元网格的划分及约束荷载的施加如图(单元尺度为5cm)约束施加时在模型左侧边界所有节点上只施加x方向的约束,即令UX=0,在左下角节点上施加x、y两个方向的约束,即UX =0、UY=0。
荷载施加在右侧边界上,大小为100。
对模型进行分析求解得到:节点应力云图(最大值222.112)单元应力云图(最大值256.408)可看出在孔周围有应力集中现象,其余地方应力分布较为均匀,孔上部出现最大应力。
②细网格单元网格的划分及约束荷载的施加如图(单元尺度为2cm)约束及荷载的施加方法如前,对模型进行分析求解得到:节点应力云图(最大值272.484)单元应力云图(最大值285.695)(2) 取1/4模型的计算①粗网格单元网格的划分及约束荷载的施加如图(单元尺度为2.5cm)=0,在下约束施加时在模型左侧边界所有节点上只施加x方向的约束,即UX侧边界所有节点上只施加y方向的约束,即U=0。
荷载施加在右侧边界上,大小Y为100。
对模型进行分析求解得到:节点应力云图(最大值251.333)单元应力云图(最大值268.888)②细网格单元网格的划分及约束荷载的施加如图(单元尺度为1cm)约束及荷载的施加方法如前,对模型进行分析求解得到:节点应力云图(最大值290.478)单元应力云图(最大值297.137)(3) 计算结果比较下面按照弹性力学理论求解带孔平板的应力集中系数。
ANSYS静力学分析APDL建模实例-应力集中————————————————————————————————作者:————————————————————————————————日期:计算分析模型如图所示, 习题文件名: scf材料参数:E=205GPa, v = 0.3力载:4500N注意单位的一致性:使用N, mm, MPa单位制建模教程在ANSYS工作文件夹内新建“stress concentration factor”目录,以存放模型文件。
注意定期保存文件,注意不可误操作,一旦误操作,不可撤销。
1.1 进入ANSYS开始→程序→ANSYS 14.5→Mechanical APDL Product Launcher14.5→然后在弹出的启动界面输入相应的working directory及文件名scf如通过Mechanical APDL 14.5进入,则进入预设的working directoryworking directory必须设置在电脑最后一个分区(因为教学用电脑只有最后一个分区不受系统保护)至此ANSYS静力学分析模块启动,ANSYS在“stress concentration factor”目录下自动创建了.log、.err等必要的文件。
2.2设置计算类型ANSYS Main Menu: Preferences →select Structural → OK2.3选择单元类型ANSYS Main Menu: Preprocessor →Element Type→Add/Edit/Delete →Add →select Solid Quad 4 node 182 →OK (back to Element Types window)→ Options… →select K3: Plane Strs w/thk →OK→Close (the Element Type window)2.4定义实常数ANSYS Main Menu: Preprocessor →Real Constants →Add/Edit/Delete →Add →OK →THK 1.2 →OK2.5定义材料参数ANSYS Main Menu: Prepr ocessor →Material Props →Material Models →Structural→Linear →Elastic→Isotropic→input EX:205e3, PRXY:0.3→ OK2.6生成几何模型✓生成特征点(8个)ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints →In Active CS→依次输入四个点的坐标:input:1(0,0,0) ,2(75,0,0) ,3(75,4.5,0) ,4(120,4.5,0) ,5(120,19.5,0),6(75, 19.5,0) ,7(75, 24,0) ,8(0, 24,0)→Apply/OK(开始点Apply,最后一个点OK)Tips:如何用ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints →On Working Plane →又该如何操作才能生成同样的点??✓直线(8条)ANSYS Main Menu: Preprocessor →Modeling →Create →Lines →Lines →Straight Lines→跳出对话框,用鼠标(左键)依次选择点1、2生成直线1,依次类推生成直线2-8。
学号: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,分析在两种载荷的作用下扳手的应力分布。
孔的应力集中问题的数值模拟(陶钧21028039)问题描述:工程结构中,常根据需要设置一些孔口,孔口附近的的应力将远大于无孔时的应力,也远大于距离孔口较远处的应力,这种现象称为孔口应力集中。
孔口的尺寸远小于弹性体尺寸。
此问题属于平面应力问题。
针对本问题,先分别建立了具有圆孔、方孔、椭圆孔的板模型,分别施加了相同的载荷,用ansys软件对其进行了模拟。
模拟过程:一、圆孔1.建立有限元模型1.1工作文件名:jisuanji;修改标题:yuankongyingli1.2单元类型:strctural solid ,quad 8node 82;材料属性:ex为2e11,prxy为0.31.3创建有限元模型,长为100,宽50,孔直径为10。
1.4生成有限元网格设定网格大小为0.5,自由划分2.加载和求解2.1施加均布载荷1000帕3.结果后处理及查看结果显示,在r=5,90=ϕ时,应力最大,最大值为3143帕,约为施加载荷的3倍。
与实际结果相符。
二、方孔1.建立有限元模型1.1工作文件名:jisuanji2;修改标题:圆孔应力1.2单元类型:strctural solid ,quad 8node 82;材料属性:ex为2e11,prxy为0.31.3创建有限元模型,长为100,宽50,正方形孔,孔边长为5。
1.4生成有限元网格设定网格大小为0.5,自由划分2.加载和求解2.1施加均布载荷1000帕3.结果后处理及查看结果显示,最大应力在方孔四角处,应力最大值为2696帕,对比圆孔处的应力集中,此应力应该大于上例中的3143。
此计算结果与实际不相符。
因此我们重新划分网格,进行计算。
设定网格大小为0.2,自由划分。
设置相同的载荷,经计算得到下列结果:结果表明,最大应力为3961帕,高于3143帕,结果与实际相符。
我们同样可以预测,如果对网格进行进一步细化,会得到更加精确的结果。
在一定的计算机能力下,分区划分网格就显得尤其重要。
不同板宽的孔边应力集中问题摘要:应用ANSYS数值模拟的方法(二维和三维)研究了含圆孔有限宽度薄板孔边应力集中问题,分析表明:平板圆孔应力集中系数的收敛性与网格划分的密度有关;应力集中系数与宽径比及长宽比有关;三维状态的内部的应力集中比二维强烈。
关键词:平面圆孔;应力集中;ANSYS;三维有限元1.引言设受力弹性体具有小孔,则孔边应力将远大于无孔时的应力,也远大于距孔稍远处的应力。
这种现象称为孔边应力集中。
孔边应力集中是局部现象,不是由于截面减小了一些而应力有所增大,而是由于开孔后发生的应力扰动所引起的。
圆孔孔边的应力可以用较简单的数学工具进行分析。
图1 平板圆[]孔如图1所示的具有小圆孔的平板,对于无限大板宽的孔边应力集中问题,有以下弹性力学解析解:在孔边的y轴上有分布:然而,实际工程上所涉及的主要是有限板宽的孔边应力集中问题,以上解析解能否适用及适用条件还值得研究。
本文就图1所示有限板宽的孔边应力集中问题,通过ANSYS软件计算其应力分布情况,采用二维模型,讨论在选取合适的网格情况下,不同的长宽比的应力集中系数变化规律及其与宽径比的关系;然后采用三维模型计算分析,与二维模型计算结果进行比较。
2.计算模型由于图1所示矩形薄板几何荷载的对称性,可选用1/4薄板作为有限元模型,坐标原点位于圆孔中心,圆孔半径R=5cm为定值,取不同的宽度和长度进行比较。
分析中采用八节点实体单元PLANE82,单元属性设置为Plane stress w/thk,弹性模量和泊松比分别为200GPa和0.3,边界条件为x=0,UX=0;y=0,UY=0。
在板远端作用有沿x轴方向的q0=1MPa的均匀分布拉力。
为了便于分析比较,定义宽径比,应力集中系数,长宽比,网格划分密度( =1时为初始网格密度,如图2所示;当 =2时,表示网格密度为初始的网格密度的2倍)。
划分的模型如图2所示。
图2 平板圆孔网格模型(网格密度 =1,长宽比 =5,宽径比 =6)3.数值模拟在同样的材料以及同样的荷载作用下,应力集中系数不仅与宽径比有关,还与网格密度以及长宽比有关。