弹塑性力学有限元上机报告(abaqus版).docxVIP

有限元上机作业系所：同济大学土木工程学院 专业： 建筑与土木工程 . 学号：####### . 姓名：## . 指导教师：### .二零一五年一月二十八日第1题（原题号5）通过Abaqus有限元软件模拟具有小圆孔的平板均匀拉伸产生的应力集中问题，验证弹性力学中的结论以及对比材料只考虑弹性和考虑弹塑性的区别。如图为平面模型，小孔半径远小于平面板的边长，两边受均匀拉应力q.1 问题描述一个承受拉力的平板，在其中心位置有一个小圆孔，平板的边长为100mm，中心小圆孔半径r=5mm，平板厚度1mm, 拉伸荷载P=100Mpa。如图1.1所示.要求分析圆孔应力集中处的Mises应力和小孔处的应力集中系数.图1.1 平板示意图运用弹性力学方法求解由于小圆孔的存在，必然对板内应力的分布发生影响：一是集中性，空口附近的应力远大于较远处的应力，且最大应力和最小应力均发生在孔边上；二是局部性，由于开孔引起的应力扰动，主要发生在距孔边1.5倍孔口尺寸的范围内，在离空边的较远处，这种影响显著地减少。由弹性力学知识知，在距圆孔圆心处的径向正应力、环向正应力、切应力分别为：沿着孔边，当时，小圆孔处的应力为： (度)030456090-q0q2q3q沿着y轴，，环向正应力为： (度)r2r3r4r3q1.22q1.07q1.04q可见孔边最大应力比无孔时提高了3倍，应力集中系数k=3，如图1.2所示。图1.2 1.2 Abaqus有限元分析材料特性：弹性模量E=2.1e11Pa，fy=418e6Pa，泊松比拉伸荷载：=1e8Pa平板尺寸：长度=100mm，宽度=100mm，厚度：1mm，圆孔直径d=5mm量纲：(N、m、Pa)建模要点：1）通过简单的分析可以知道该问题属于平面应力问题；基于结构荷载的对称性，可以之取模型的1/4/进行分析。利用Abaqus软件生成薄板应力图，根据板的应力图可以很直观地看到板的应力变化情况，各处的Mises应力值，可求得应力集中因子。1.2.1 材料只考虑弹性时基于结构和荷载的对称性，只取结构的1/4进行分析，如图1.3所示：图1.3带孔平板的分析模型创建部件Modeling Space选择2D Planar，Type为Ddformable，Base Feature选择shell；圆弧半径为5mm，矩形边长为50mm。如图1.4所示：图1.4 创建part 划分网格在主菜单中选择Seed→Edge，将边界上的单元数设置为8（如图1.5所示）；在主菜单中选择Mesh→Controls，将Techniques设为Structured；在主菜单中选择Mesh→Element Type，弹出Element Type对话框，考虑到对于应力集中问题，使用二次单元可以提高应力结果的精度，将Geometric Order（几何阶次）设为Quadratic（二次单元），并取消对Reduced integration（减缩积分）的选择（如图1.7所示），最后生成网格（如图1.6所示）。图1.5边界单元数设置图1.6生成网格图1.7设置单元类型创建材料和截面属性创建材料特性：弹性模量E= 2.1*1011Pa,泊松比0.3；创建截面属性：Category（截面类型）为Solid（实体），Type（类型）为Homogeneous（匀质的）,Plane stress/strain thickness为1；将创建好的截面赋给部件。装配部件整个分析模型是一个装配件，前面在Part中创建的部件在Assembly功能模块中装配起来。设置分析步ABAQUS中默认有一个初始分析步（initial），可以在其中施加边界条件。再创建一个名为Load的分析步，用来施加荷载。Procedure type: General,选择Static,General。定义边界条件和荷载在初始分析步（initial）中创建两个名为BC-1和BC-2的边界条件，如图1.8所示；定义平板左边上的对称边界条件：根据对称性可知，左侧边界U1=UR2=UR3=0，选中左侧边界，定义边界条件，选中XSYMM，如图1.9所示；定义平板底边上的边界条件：平板底边上的边界条件为U2=UR1=UR3=0，选中底部边界，定义边界条件，选中YSYMM，如图1.10所示；荷载类型选择Pressure,然后选中平板右侧的边界线，输入荷载，Magnitude后面输入-1e8,如图1.11所示边界条件和荷载定义以后的模型如图1.12。图1.8创建边界图1.9定义平板左边上的对称边界条件图1.10定义平板左边上的对称边界条件图1.11施加荷载图1.12位移边界条件和荷载提交作业后处理计算得到平板的Mises应