板结构有限元分析实例详解

2、分析类型静力分析3、问题描述板壳问题4、ANSYS单位m N Pa5、单元SHELL63：4节点（每个节点6个自由度）6、材料弹性模量和泊松比7、实常数厚度8、建模取1/4模型41进入ANSYS程序→ANSYSED 9.0→ANSYS Product Launcher→change the working directory into yours →input Initial jobname:board→Run2 设置计算类型ANSYS Main Menu: Preferences →select Structural →OK3 定义材料参数ANSYS Main Menu: Preprocessor →Material Props→Material Models →Structural →Linear →Elastic→Isotropic →input EX:2e11, PRXY:0.3 →OK54 选择单元类型ANSYS Main Menu:Preprocessor →Element Type→Add/Edit/Delete→Add →select Shell Elastic 4node 63→OK (back to Element Types window) →Close (the Element Type window)65、定义实常数ANSYS Main Menu: Preprocessor →Real Constants…→Add…→select Type 1→OK→input Shell thickness at node I :0.01→OK →Close (the Real Constants Window)76、创建矩形ANSYS Main Menu:Preprocessor →Modeling →Create →Areas→Rectangle→By Dimensions →依次按下图输入：→OK87、划分单元ANSYS Main Menu: Preprocessor →Meshing →Mesh Tool →OK(Size Controls) Global: Set →input NDIV:20→Mesh Tool →Mesh :select Areas→Shape:Quad→Free→Mesh →Pick All →Close( the Mesh Tool window)98、施加固定边界ANSYS Main Menu: Solution→Define Loads →Apply→Structural →Displacement→On Lines→拾取上边和右边→OK→select Lab2:ALL DOF →OK109、施加对称约束边界ANSYS Main Menu: Solution→Define Loads →Apply→Structural →Displacement →Symmetry B.C →On Lines→拾取下边和左边→OK1111 分析计算ANSYS Main Menu: Solution →Solve →Current LS→OK(to closethe solve Current Load Step window) →OK12、挠度ANSYS Main Menu: General Postproc→Plot Results→Contour Plot→Nodal Solu→select: DOF solution, Z-Component of displacement→OK13按右侧的等侧视图141617板壳单元的压力载荷面1 (I-J-K-L)，（底部，+Z方向）面2 (I-J-K-L)，（顶部，-Z方向）面3 (J-I)，面 4 (K-J)，面 5 (L-K)，面 6 (I-L)19显示壳单元的厚度ANSYS命令菜单栏：PlotCtrls>Style >Size and Style→弹出对话框中（Display of element项设置为On)控制壳单元的结果输出ANSYS Main Menu: General Postproc→Options for OutpTop layer: 顶面Middle layer:中面Bottom layer：底面2021作业三要求1、属于力学的那类问题？2、单位制；单元类型；单元描述；实常数；材料参数3、划分网格的方案，施加载荷和边界条件4、计算结果的体现：MISIS 应力和位移，需要知道最大值以及位置。

板结构有限元分析实例详解板结构是一种常见的结构形式，广泛应用于建筑、航空航天、机械、电子等领域。

板结构的特点是结构主要由板和边界构件组成，受到外加载荷作用时，产生弯曲和剪切变形。

为了评估板结构的强度和稳定性，可以使用有限元分析方法进行分析。

本文将以一座大跨度板结构为例，详解板结构有限元分析的步骤及其相关实例。

首先，我们需要对板结构进行几何建模。

通常情况下，板结构可以简化为二维平面问题。

我们可以使用专业的有限元分析软件，如ANSYS、ABAQUS等，进行几何建模。

在建模过程中，需要确定结构的几何形状、边界条件、加载方式等参数。

以一块长方形板作为例子，我们可以在软件中创建一个二维平面，并定义板的几何尺寸和材料属性。

接下来，我们需要对板结构进行网格划分。

有限元分析方法将结构划分为许多小的单元，然后对每个单元进行分析计算。

在板结构分析中，常用的单元类型包括矩形单元、三角形单元、四边形单元等。

我们可以根据实际需要选择适当的单元类型和网格密度，并利用软件自动生成板结构的网格。

然后，我们需要为板结构定义边界条件。

边界条件包括支撑条件和加载条件两个方面。

支撑条件描述了板结构受力的边界，通常包括固定支撑、滑动支撑、自由支撑等情况。

加载条件描述了外力或外载荷施加在板结构上的方式和大小。

在我们的例子中，假设板结构的四个边界均为固定支撑，我们可以在软件中设置相应的边界条件。

之后，我们需要为板结构定义材料属性。

板结构的材料属性包括弹性模量、泊松比、密度等参数。

这些参数描述了板结构在受力时的材料性能和特性。

我们需要根据实际的材料情况，为板结构指定合适的材料属性，并在软件中进行设置。

最后，我们可以对板结构进行有限元分析计算。

在软件中，我们可以选择合适的求解器和分析方法，进行结构的静力分析、动力分析、稳定性分析等。

通过有限元分析，我们可以得到板结构在受力状态下的变形、应力分布、应变分布等结果。

总之，通过板结构的有限元分析，我们可以对结构的强度、稳定性、振动等性能进行评估和优化。

一、平面3节点三角形单元分析的算例如图所示为一矩形薄平板，在右端部受集中力F=10000N 作用，材料常数为：弹性模量1E 二1 107 Pa 、泊松比，板的厚度为t 二0.1m ，试按平面应力问题计算各个节点位3移及支座反力。

閹4-^20右瑞部受集中力作用的平面问题(高深梁)解:(1) 结构的离散化与编号形单元。

载荷F 按静力等效原则向节点 1节点2移置等效。

约束的支反力列阵： R =「0 0 0 0 R X 3 &3 R x4 Ry 4T其中(R X 3, R y3)和(R x4,R y4)分别为节点3和节点4的两个方向的支反力。

(2)各个单元的描述当两个单元取图示中的局部编码(i,j,m)时，其单元刚度矩阵完全相同，即a=X j y m —X m y j , b i=y j —y m , C i=X j —X m a j =X m y i -xy m , b j =y m -y i , c^x^X i a m =X j y m —X m y j , b m =y i —y j , c^x^X j对该结构进行离散，单元编号及节点编号如图4-20(b)所示，即有二个3节点三角 节点位移列阵:q - g v u 2 v 2U 3 V 3 U 4 V 4 r1-^b r b s 十 ---- C r C s2 1 - J亠 -------- b r C s 2^C r b s十丁农I1 - J GG 2brbs"b r C sk iik ⑴,（2）k jjkmma)制赳描述 b)冇限元分折模型节点外载列阵:k ijk jmk mj1 00132I23-12T1I3023434323220042一』—112427433212d413L 333T(3)建立整体刚度方程按单元的位移自由度所对应的位置进行组装可以得到整体刚度矩阵,成k二k⑴k⑵具体写出单元刚度矩阵的各个子块在总刚度矩阵中的对应位置如下该组装过程可以写代入整体刚度方程Kq =P中，有(5)支反力的计算将所求得的节点位移式代入总刚度方程中，可求得支反力如下9Et/ 2 4 、 “(-U i ——V i + — V 2) = —2 F 32 3 3 9Et 2 1 4(U v U 2)»0.07F 32 3 3 3 9Et 2 (-U 2 V 2) =2F32 3 9Et 2 1 (u 2 v 2) = 1.07F7 亍474 3 n — 3-1—3斗1 3 q21-4一 一 □0 33 3 34 r74 "q—0 □-13 3 7亍31341-4 033 33 I2 4 "7'4"-1 7 0亍T7亍r 141 32— —□=-4 33J3 3 i hi BW-Wa-ivKaH4q7 4 0 □ -1——-—33 3 331 2413■^= -^^a--433373 J9Z732(4) 边界条件的处理及刚度方程求解该问题的位移边界条件为 U 3 = 0,V 3 二 0, U 4 二 0,V4 -0 将其代入上式中，划去已知节点位移对应的第5行至第8行(列)，有13由上式可求出节点位移如下[U i w u 2 v 2]Et131.88 -8.99-1.50 -8.42 TR X 4 R y432 3 3、MATLAB —平面3节点三角形单元分析的算例(Triangle2D3Node)解：(1)结构的离散化与编号将结构离散为二个3节点三角形单元，单元编号及节点编号如图4-20(b)所示。

基于ANSYS的木板床的有限元分析姚美良（昆明理工大学建筑工程学院，云南昆明650500）摘要：木板床是我们生活中的必需品。

本文主要利用ANSYS来分析不同重量的人躺在床上的时候，木板床的应力应变；并通过不同的网格划分进行同一目标的计算。

对比两种不同的情况，发现网格的优劣，对于计算结果的影响较大；找到木板床整个结构的薄弱点，并且通过对比不同种材料属性的木头的床体应变，为木床选材提供了部分依据。

关键词：木板床；ANSYS；有限元；对比分析中图分类号：TP319 文献标识码：A 文章编号：Finite element analysis of wood bed based on ANSYSYAO Meiliang（Faculty of Civil Engineering, Kunming University of Science and Technology, Kunming, Yunnan, 650500, China）Abstract：Board bed is a necessity in our life. This paper mainly uses ANSYS to analyze the stress strain of the wood bed when people lying in bed with different weight, and calculate the same target by different grid. The advantages and disadvantages of grid influence on the calculation results, find plank bed the whole structure of the weak point and through the comparison of different kinds of material attributes of wood bed body strain, as the material of the wooden bed provides part of the basis of the comparison of two different cases.Key words：Wood bed ；ANSYS；Finite element；Comparative analysis0引言木板床在我们的生活中非常常见，而且是我们的生活必须品，我们每天都会接触到。

1.工程概况：在一个展览馆设计中，对一个4.8m跨的折板楼梯按上述思路进行分析设计，板厚t=100mm，取到了折板投影跨度的1/48，简化算法至少要取到160mm以上，仅此一项砼量节省了近40%。

配筋量更是不到简化算法配筋的一半。

所用软件为midas/gen（V.7.3.0）。

一、楼梯剖面：二、荷载标准值（kN/㎡）：2四、梯板设计：4.1 梯板折角处M=5.8 kN.m：计算面积As：205 mm2；采用方案： 8@150；实配面积：335 mm2。

4.2平板跨中处M=1.9 kN.m计算面积As：69 mm2；采用方案： 8@200；实配面积：251 mm2。

4.3斜板跨中处M=4.3 kN.m计算面积As：151 mm2；采用方案： 8@200；实配面积：251 mm2。

4.4斜板支座处M=8.6 kN.m计算面积As：310 mm2；采用方案： 8@150；实配面积：335 mm2。

由上述分析我们能够得出如下结论：1、斜板计算跨中弯矩时可以按等投影跨度的PB模型计算。

M PBmax=29.2kN.m；M XBmax=31.1kN.m(PB弯矩理论解：Mmax=q*L^2/8=(1.2*8+1.4*3.5)*42/8=29kN.m，即数值解与理论解吻合较好，表明有限元计算时的剖分尺度合理，分析结果可信)2、折板（含上、下折板）不适合简化为等投影长度的PB计算：1）对于折板计算，弯折处均有负弯矩存在，若简化为PB便不能充分体现,造成局部的不安全；2）对于板底正弯矩，简化为PB时，在某些情况下计算弯矩是其实际弯矩的5倍有余，比如：M CT2/4=5.4 kN.m; M PB=29.2kN.m，则：M PB /M CT2/4=29.2/5.4=5.43）板厚取值不宜按投影跨度的1/30取值，应根据折板的实际形态有区别地对待。

要不然不仅造成不必要的浪费，而且毫无意义地增加结构自重。

根据弯折位置及弯折角度的不同，有时可取到投影跨度的1/50（详见附件工程案例1）。

板的屈曲分析为证明ansys 程序计算的有效性，特分析了两个个有代表性的算例。

算例1四边简支的矩形板，边长a =600mm ，b =600mm ，t =10mm ，E =206000 N/m m 2 ν=0.3，两对边受到均匀分布的压力。

图1.1 薄板有限元模型约束图根据简支模型和及均布压力屈曲模式建立板模型，下面是对模型的说明。

(1) 单元：采用三维壳体单元shell63。

(2) 单元划分情况：沿长度a 方向为100份，沿宽度b 方向为100份。

(3) 约束情况：四边为简支。

(4) 加载情况：两对边受单位均布压力XY Z图1.2加载图最终利用ansys 计算出的临界屈曲应力为196.67Mpa 。

根据板的稳定理论可知：在a/b ≤ 2的情况下F =π2D b 2 b a +a b2D =Et 312 1−ν2根据计算得到σ =206.66Mpa Ansys 计算结果的相对误差为4.8%图1.3 应力云图算例2XYZSEP 3 2012MNMXXY ZSEP 4 2012三边简支的矩形板，边长a =600mm ，b =600mm ，t =10mm ，E =206000 N/m m 2 ν=0.25，两对边受到单位均匀分布的压力。

图1.4单元剖分与约束图根据简支模型和及均布压力屈曲模式建立板模型，下面是对模型的说明。

(1) 单元：采用三维壳体单元shell63。

(2) 单元划分情况：沿长度a 方向为100份，沿宽度b 方向为100份。

(3) 约束情况：三边为简支。

(4) 加载情况：两对边受单位均布压力图1.5加载图XYZXYZSEP 3 2012最终利用ansys 计算出的临界屈曲应力为72.42Mpa 。

根据板的稳定理论可知：F =k π2D b 2D =Et 312 1−ν2当a/b =1时，k=1.44 根据计算得到σ = 74.4 Mpa Ansys 计算结果的相对误差为3.9%图1.6 应力云图加筋板算例1四边简支的矩形板，边长a =3000mm ，b =1500mm ，t =10mm ，E =206000 N/m m 2ν=0.3，两对边受到均匀分布的压力，中间加一道纵向扁钢。

1问题描述某周边简支非均匀的矩形（或圆形）板在均布载荷作用下挠度过大。

结合实际，提出集中改进设计方案，并进行对比分析。

2.问题分析不均匀板有两种主要的情况，结构不均匀和材料不均匀，结构不均匀是指板的厚度不是常量，材料不均匀体现在板的弹性模量和泊松比是变化的。

另外，有的板可以是以上两种情况的混合情形。

不均匀板与均匀板的有限元问题有哪些差别呢？下面从均匀板问题推导出非均匀板有限元问题的解决方法。

2.1应力应变先以结构不均匀板为例来讨论。

假设一矩形板长为2，宽为2，厚度沿x ，y 不均匀，由一函数()h ,h x y =描述，但仍然符合薄板假设。

对于均匀板，显然h 是一个常数。

设挠度为()=x,y ωω，则板内应变向量可以表示为{}2222211==z 12x x y y xy xy x z y x y ρεεεωεγγ⎧⎫⎧⎫∂⎪⎪⎪⎪∂⎪⎪⎪⎪⎧⎫⎪⎪⎪⎪∂⎪⎪⎪⎪⎪⎪=-⎨⎬⎨⎬⎨⎬∂⎪⎪⎪⎪⎪⎪⎩⎭⎪⎪⎪⎪∂⎪⎪⎪⎪∂∂⎪⎪⎪⎪⎩⎭⎩⎭应力应变关系为{}1p z D σρ⎧⎫⎡⎤=⎨⎬⎣⎦⎩⎭弯矩扭矩矩阵{}{}()()h ,2h ,2x y x y M zdz σ-=⎰这里就体现出不均匀板和均匀板的区别了。

积分完毕后，可以得到{}[]1M D ρ⎧⎫=⎨⎬⎩⎭其中薄板的弯曲系数矩阵[]()()()321,1012101/2Eh x y D μμμμ⎡⎤⎢⎥=⎢⎥-⎢⎥-⎣⎦是关于薄板总体坐标的函数，所以对各个分单元都是不同的。

各单元的弯曲系数矩阵可以采用单元中心处的代替。

那么就可以得出一系列的弯曲系数矩阵[]D ei 。

如果单元划分得足够细，是可以代替真实解的。

2.2单元分析可以将板分为边长为0.25的矩形小单元，每一个单元都是一样的。

对于任何一个单元的节点，都有3项独立的位移{}i i i xi i yi i w w w y w x δθθ⎧⎫⎪⎪⎪⎪⎧⎫⎪⎪⎛⎫∂⎪⎪⎪⎪==⎨⎬⎨⎬⎪∂⎝⎭⎪⎪⎪⎪⎩⎭⎪⎪∂⎛⎫⎪⎪- ⎪∂⎪⎪⎝⎭⎩⎭位移模式()22312345672233389101112,w x y x y x xy y x x y xy y x y xy αααααααααααα=+++++++++++形状函数矩阵是一个112⨯的行向量()[],kl mn N x y N N N N =⎡⎤⎣⎦其中222222222222222211128111111i i i i i i i i i i i i i x x y y x x y y x y N a b a b a b x x y y y y x x y y x x y x a b b a b a ⎛⎫⎡⎛⎫⎛⎫=++++--⎡⎤ ⎪⎪⎪⎣⎦⎢⎝⎭⎝⎭⎣⎝⎭⎤⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫++--++-⎥ ⎪⎪ ⎪ ⎪ ⎪⎪⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭⎥⎦(),,,i k l m n =单元刚度矩阵[][][][]1212ee TS k B D B dxdy ⨯=⎰很明显，积分式中包含了弹性系数矩阵，而不同单元的弹性系数矩阵是不同的，所以，即便单元划分相同，得到的单元刚度矩阵也不同。

板中圆孔的应力集中问题：如图所示为一个承受单向拉伸的无限大板，在其中心位置有一个小圆孔。

材料属性为弹性模量E= Pa泊松比为0.3,拉伸载荷q=1000Pa平板厚度t=0.1.201、定义工作名和工作标题（1）定义工作文件名：在弹出的Cha nge Job name对话框中输入Plate。

选择New log and error files复选框，单击OK按钮。

（2）定义工作标题：在弹出的的Change Title对话框中输入The analysis of plate stress with small circle单击OK按钮。

（3）重新显示：执行replot命令。

2、定义单元类型和材料属性（1）选择单元类型：在弹出的Element Type中，单击Add按钮，弹出所示对话框，选择Structural Solid 和Quad 8node 82选项，单击OK ,然后 单击close 。

(2) 设置材料属性：在弹出的 defi ne material models behavior 窗口中，双 击 structural/linear/elastic/isotropic 选项，弹出 linear isotropic material properties for material number 1 对话框，EX 和 PRXY 分别输入 2e11 和0.3,单击OK,执行exit 命令Lintar Isotropic f^lattrial Proptrtits for f^lattrial Numbtr 1Add Temperature Delete Temperature□K（3）保存数据：单击SAVE_DB 按钮。

3、创建几何模型PRXYCancel Linejr Isotropic Properties for Material Nuinbcr 1T1Help（1）生成一个矩形面：执行相应操作弹出create recta ngle by dime nsio ns对话框，输入数据，单击OK,显示一个矩形。

板结构有限元分析实例详解1：带孔平板结构静力分析本节介绍带孔平板结构静力分析问题，同时介绍布尔操作的基本用法。

8.3.1 问题描述与分析有孔的矩形平板，左侧边缘固定，长400mm，宽200 mm，厚度为10 mm，圆孔在板的正中心，半径为40 mm，左侧全约束，右侧边缘均布应力1MPa，如图8.7所示。

求板的变形、位移及应力变化情况。

（材料的材料属性为：弹性模量为300000 MPa，剪切模量为0.31。

）图8.7 带孔的矩形平板由于小孔处边缘不规则，本文采用PLANE82高阶平面单元进行分析。

8.3.2 求解过程8.3.2.1 定义工作目录及文件名启动ANSYS Mechanical APDL Product Launcher窗口，如图8.8所示。

在License下拉选框中选择ANSYS Multiphysics产品，在Working Directory输入栏中输入工作目录：C:\ANSYS12.0 Structural Finite Elements Analysis and Practice\Chapter8\8-1，在Job Name一栏中输入工作文件名：Chapter8-1。

以上参数设置完毕后，单击Run按钮运行ANSYS。

图8.8 ANSYS设置窗口菜单可以先在目标文件位置建立工作目录，然后单击Browse按钮选择工作目录；也可以通过单击Browse按钮选择工作文件名。

8.3.2.2 定义单元类型和材料属性选择Main Menu>Preferences命令，出现Preferences for GUI Filtering对话框，如图8.9所示，在Individual discipline(s) to show in the GUI中勾选Structural，过滤掉ANSYS GUI菜单中与结构分析无关的选项，单击OK按钮关闭该对话框。

图8.9 Preferences for GUI Filtering对话框选择Main Menu>Preprocessor>Element Type>Add/Edit/Delete命令，出现Element Types对话框，如图8.10所示，单击Add按钮，出现Library of Element Types对话框，在Library of Element Types列表框中选择Structural Solid中的Quad 8node 82单元，如图8.11所示，单击OK按钮关闭该对话框，单击Element Types对话框上的Options按钮，弹出PLANE 82 element type options对话框，在Element behavior K3 项下拉菜单中选择PLANE strs w/thk，如图8.12所示，单击OK按钮关闭该对话框，单击Close按钮关闭Element Types对话框。

ANSYS上机实验报告实验三：板壳的有限元分析班级：姓名：学号：一、实验题目图示正方形平板，承受垂直于板面的均布载荷作用P=20KN/m*m，板厚t=0.1m，平板外缘各边采用固定约束方式，材料选用低碳钢，弹性模量E=210GPa，u=0.33。

二、实验过程1、确定所采用的单位制：N，m，Pa。

2、问题类型：板壳问题。

3、利用ANSYS构造实体模型：1/4模型（正对称）和整体分析。

4、网格划分1）、定义材料属性：Main Menu: Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX: 210e9, PRXY: 0.33 →OK2）、定义单元类型：Main Menu: Preprocessor →Element Type→Add/Edit/Delete →Add →select Shell Elastic 8node 63 →OK (back to Element Types window)3）、定义实常数（厚度）：Main Menu: Preprocessor →Real Constants… →Add… →select Type 1→OK→input TK（I）: 0.1 ，TK（J）: 0.1 ，TK（K）: 0.1 ，TK（L）: 0.1 →OK→Close (the Real Constants Window)4）、划分网格：在size element edge length （单元边长值）处输入0.25、加载及求解。

加载（整体）过程：Main Menu: Solution →Define Loads →Apply →Structural →Displacement →On Nodes →拾取四个边线→OK →select Lab2:ALL OFF →OKMain Menu: Solution →Define Loads →Apply →Structural →pressure →On Areas →拾取面→OK →Value: 20000→OK加载（1/4）过程：Main Menu: Solution →Define Loads →Apply →Structural →Displacement →On Nodes →拾取右边线和上边线→OK →select Lab2:ALL OFF →Apply→拾取左边线→OK→select Lab2:UX→Apply→拾取下边线→OK→select Lab2:UY→OKMain Menu: Solution →Define Loads →Apply →Structural →pressure →On Areas →拾取面→OK →Value: 20000→OK求解：Main Menu: Solution →Solve →Current LS →OK(to close the solve Current Load Step window) →OK6、分析变形、位移和应力状况并抓图。

开孔薄板有限元分析报告一、有限元分析的目的通过对两种模型（一个上边开口的和另一个上下两边开口的模型）的静力分析，比较与其对应的理论解的不同，了解有限元仿真软件与理论计算存在的，进一步熟悉workbench求解有限元问题的一般步骤。

二、实体建模（两个模型）建立如下所示的模型，其中，边长300mm，宽80mm，厚5mm，边缘为半径是10mm的半孔。

上边开口的实体模型（模型A）上下两边开口的模型（模型B）模型A模型B模型采用的单元类型模型A：1 3862 LID186 (20 Node Quadratic Hexahedron)2 3858 SOLID186 (20 Node Quadratic Wedge)3 112 CONTA174 (Quadratic Quadrilateral Contact)4 112 TARGE170 (Quadratic Quadrilateral Target)5 64 SURF154 (3D Quadratic Quadrilateral)模型B：1 3856 SOLID186 (20 Node Quadratic Hexahedron)2 3866 SOLID186 (20 Node Quadratic Wedge)3 100 CONTA174 (Quadratic Quadrilateral Contact)4 100 TARGE170 (Quadratic Quadrilateral Target)5 64 SURF154 (3D Quadratic Quadrilateral)2.载荷与约束的施加方法（绘图表示并说明）；两模型施加的载荷与约束相同约束：单击static structural，选择长方体的左侧面，鼠标右键选择“insert>fixed support”载荷：选择长方体的左侧面，鼠标右键选择“insert>force”，大小为50N。

二、带孔平板的有限元分析1：问题描述图所示为一个有中心圆孔的薄板，薄板厚度t=0.01m，薄板弹性模量E=210000N/cm2,泊松比μ=0.3，p=100N/cm,ρ=2.7g/cm3此问题为平面应力问题，用有限元求解出带孔平板的应力集中问题，并与弹性力学的精确解进行比较。

2：求解步骤第一步：建立工作文件名和工作标题（1）选择Utility Menu—File—Change Jobname命令,出现Change Jobname对话框。

在Enter new jobname 输入栏中输入工作文件名plate，单击Ok按钮关闭该对话框。

（2）选择Utility Menu—File—Change Tile命令，出现Change Tile对话框，在输入栏中输入Stress analysis in a sheet，单击Ok按钮关闭该对话框。

第二步：设置计算类型选择Main Menu—Preference—Structural-Ok命令.第三步：选择单元类型选择Main Menu—Preprocessor—Element Type—Add/Edit/Delete命令,出现Element Type对话框，选择Solid-Quad 4node 42—Ok命令,再在Element Type对话框中选择Options—K3:Plane Strs w/thk/—Ok—Close命令.第四步：定义材料参数选择Main Menu—Preprocessor—Material Props—Material Models—双击Structural—双击Linear—双击Elastic—双击Isotropic命令，出现如下对话框填写Ex：2.1e5，PRXY：0.3；选择Ok命令。

第五步：定义实常数以确定平面问题的厚度选择Main Menu—Preprocessor—Real condtants—Add/Edit/Delete—Add—Type1—Ok命令，出现以下对话框，在Real condtant Set No中填写1，在THK中填写1，选择Ok—Close命令.第六步：创建几何模型1：生成平面方板选择Main Menu—Preprocessor—Modeling—Creating—Areas—Rectangle—By 2 Corners—Wp X:0, Wp Y:0,Width:100,Height:100—Ok1：生成圆孔平面选择Main Menu—Preprocessor—Modeling—Creating—Areas—Circle—Solid Circle—Wp X:50, Wp Y:50,Radis:5—Ok2：生成带孔方板选择Main Menu—Preprocessor—Modeling—Operate—Booleana—Subtract—Areas，鼠标点击方板1—Ok，在Multi-Entities窗口点击Ok，在Subtract Areas窗口点击Ok.. 鼠标点击圆孔2—Ok, 在Multi-Entities窗口点击Ok，在Subtract Areas窗口点击Ok.出现如下图1：第七步：网格划分选择Main Menu—Preprocessor—Meshing—MeshTool命令，在MeshTool窗口点击Size Controls下的Globle:Set—NDIV:29—Ok, 在MeshTool窗口点击Mesh—Pick all—Close命令。