有限元上机实验报告
ABAQUS有限元上机报告 南理工
有限元上机实验报告[——以Abaqus软件进行的有限元分析]汪健强1008320139实验1——平面问题应力集中分析目的要求：掌握平面问题的有限元分析方法和对称性问题建模的方法。
通过简单力学分析，可以知道本实验问题属于平面应力问题，基于结构和载荷的对称性，可以只取模型的1/4进行分析。
用8节点四边形单元分析X=0截面σx的分布规律和最大值，计算圆孔边的应力集中系数，并与理论解对比。
一、实验过程概述：1、启动ABAQUS/CAE2、创建部件3、创建材料和截面属性4、定义装配件5、设置分析步6、定义边界条件和载荷7、划分网格8、提交分析作业9、后处理10、退出ABAQUS/CAE二、实验结果：（1）边界受力图（1）X方向应力分量σx应力云图：（2）左边界直线与圆弧边交点的σx值为： 2.96714 MPa；（2）左右对称面上的σx曲线：三、实验内容分析：a)模型全局σx应力分布：σx应力集中分布于中心圆孔与x、y轴相交的地方，且与x轴相交处应力为负，与y轴相交处应力为正；沿圆周向周围，σx迅速减小；沿y 方向的σx应力大于沿x方向的σx应力。
b)应力集中系数为 2.92975，小于理论值3.0。
误差来源：有限元分析方法是将结构离散化，网格划分得越稀疏，计算出的结果就越偏离理论值。
分的越密集，结果越接近与理论值。
四、实验小结与体会：通过本次实验，对理论课所学有限元基本方法有了一个更加直观、深入的理解。
通过对Abaqus软件三个步骤：前处理、分析计算、后处理的操作，了解了这款软件的基本应用和它对有限元的一些很好的应用。
试验中，遇到诸多问题，仔细思考，加之请教老师，逐一解决，确实很有收获。
更增加了对有限元的认识，和对其功能之强大有了更深的理解。
实验二平面问题有限元解的收敛性一、实验目的和要求：（1）在ABAQUS软件中用有限元法探索整个梁上σx和σy的分布规律。
（2）计算梁底边中点正应力σx的最大值；对单元网格逐步加密，把σx的计算值与理论解对比，考察有限元解的收敛性。
有限元法基础及应用实验上机报告实例
一，实验描述：1、本作业属于哪类问题在本作业中，根据板的结构特点和受力情况，确定该问题属于平面对称应力问题，定义分析类型为静力学分析。
2、本文采用如何的单位制本题中，长度单位为m，故为方便起见，采用力的单位为N，压强的单位为pa，时间的单位为s，质量的单位为kg。
3、单元类型：对单元描述；材料；实常数单元类型为选取shell Elastic 4node63单元材料：弹性模量为2.1e11pa，泊松比为0.33。
由题意确定时常数，即厚度为0.1m4、划分网格。
网格划分设置。
单元数，节点数。
网格化分设置：设置单元边长值为0.1m，指定单元形状为Areas。
5、加载描述（1）对整体模型，首先对四周固定端加载位移约束，设定其位移值为0；然后，对面施加面载荷，设置面载荷为20000N；接着将施加在实物体上的载荷转换到有限元模型上，并显示施加在有限元模型上的载荷。
加载完后，对该有限元模型进行求解。
（2）对四分之一模型，首先对两边固定端加载位移约束，设定其位移值为0；然后，对两坐标轴所在的边施加对称载荷，最后对面施加面载荷，设置面载荷为20000N；接着将施加在实物体上的载荷转换到有限元模型上，并显示施加在有限元模型上的载荷。
加载完后，对该有限元模型进行求解。
6、后处理：最大MISIS应力和最大位移的位置和大小。
绘制结构的应力和变形图。
二，实验步骤（一）绘制整体实体模型1，在ANSYS中构造实体模型，如下图所示2，根据结构特点及所受载荷地情况，选取shell Elastic 4node63单元，设置材料常数：弹性模量E=2.1e11，u=0.33，单元边值为0.1m对其进行网格剖分，网格划分图如下：3，正确施加载荷和边界条件，结果如下：边界条件施加载荷20000N/m求解以后4，绘制平板的应力和变形图，并给出最大应力和变形的位置及大小：应力图应变图从图中可以看出其中实体边界中点位置的应力最大为0.583e+07pa，最大变形在中间圆弧的位置，0.144e-3m（二）绘制四分之一的实体模型1，在ANSYS中构造实体模型，如下图所示2,根据结构特点及所受载荷地情况，选取shell Elastic 4node63单元，设置材料常数：弹性模量E=2.1e11，u=0.33，单元边值为0.1m对其进行网格剖分，网格划分图如下：3，正确施加载荷和边界条件，注意此处有一个对称载荷的加载，结果如下：4，绘制平板的应力和变形图，并给出最大应力和变形的位置及大小：应力图应变图三，实验小结这次ANSYS上机实验课是使用shell中的Elastic 4node63单元，这使我对ANSYS软件中的单元有了更深的认识，同时对平面问题的静力分析的基本思路和操作步骤更加熟悉。
有限元上机实验报告
有限元ANSYS实验报告
学校：华北水利水电大学
学院：机械学院
专业：机械设计制造及其自动化
姓名：
学号：2010
指导老师：纪占玲
（一）带孔板壳模型静力分析一、新建文件
二、预处理，选择材料模型类型等。
三、建模
四、划分单元
五、施加约束、载荷
六、求解
七、查看结果
（二）内六角扳手静力分析
问题：
一个截面宽度为10mm的内六角扳手，在手柄的顶部施加扭矩为100N，然后在相同的部位施加垂直向下的力20N，分析在两种荷载作用下扳手的应力分布。
尺寸如下：截面宽度10mm、形状为正六边形、手柄长20cm、杆长7.5cm，倒角半径1cm、弹性模量2.1×10¹¹Pa，泊松比0.3 。
一、新建文建，预处理和上面一样，把不同的模型类型选择如下：
二、建模
三、划分单元网格，并生成实体模型
四、施加约束、载荷
五、查看结果
（三）其它练习实例。
有限元上机报告
有限元上机报告模板一、实验目的题目：图示折板上端固定，右侧受力F=1000N，该力均匀分布在边缘各节点上；板厚t=2mm，材料选用低碳钢，弹性模量E=210GPa，u=0.33。
此题属于平面应力问题，采用的单元类型为：Solid Quad-8node单位制：力（N）、弹性模量（MPa）、长度（mm）。
二、实体建模的方法1.定义矩形：Main Menu>Preprocessor>Modeling>Create>Areas>Rectangle>By Dimensions。
输入：X1=0，X2=30，Y1=0，Y2=60。
点击对话框上的Apply按钮创建第一个矩形。
输入X1=30，X2=60，Y1=0，Y2=30。
点击对话框上的OK按钮创建第二个矩形。
2. 面图元的加运算：Main Menu>Preprocessor>Modeling>Operate>Booleans>Add>Areas。
点击对话框上的Pick All按钮，所有的面被加在一起。
点击ＳＡＶＥ＿ＤＢ保存数据。
三、单元类型1. 此题属于平面应力问题，采用Solid单元的Plane82。
2. 材料属性：弹性模量为21e4MPa，泊松比为0.33。
3. 实常数：板厚2mm。
四、划分网格以边长为2mm的正方形划分单元，总共有675个单元，736个节点。
具体操作如下：1. Main Menu >Preprocessor>Meshing>Mesh Tool.2. 在Mesh Tool面板上点击Size Controls模块内的Global Set 按钮，在弹出的对话框中输入Element edge length为2，按ok完成。
3. 在Mesh下拉框中选Areas。
4. 点击Mesh按钮网格划分，点击Pick All按钮，点击Close关闭Mesh Tool 面板。
南京理工大学2016年有限元上机实验报告(ABAQUS)
点线性等参元（完全积分 Quad，Linear；减缩积分 Quad，linear，Reduced integration；非协调模式 Quad，Linear，Incompatible modes）和 8 节点二次等参 元（Quad，Quadratic） 。
7 创建并提交分析。 ○ 8 查看结果并分析。 ○
4 计算结果分析讨论与结论
4.1 粗网格下梁中部应力分量和上下边法向应力对比
1 理论解： ○
X 方向正应力由下式计算：
已知 q=1N/mm2 ，h=160mm，L=1000mm， ymax
h 代入上式得： 2
3
x max
6 106 1 1 1 3 0 0.08 106 4 29.497MPa 3 0.16 4 2 4 5
分别应用 3 节点三角形单元、4 节点线性等参元（完全积分、减缩积分、 非协调模式） 、8 节点二次等参元完全积分进行下列各项数值实验：1）用粗网 格求解梁中部应力分量 x 最大值和上下边法向应力分量，并通过精确解对采用 不同单元的 x 计算精度进行对比分析；2）对粗网格下梁中部铅直（y 向）位移 进行对比分析；3）通过多次网格加密，对比试验 3 节点三角形单元和 8 节点二 次等参元的收敛速度。总结出研究结论，撰写实验报告。
2 3 节点三角形单元计算结果： ○
x 的应力云图
梁中部应力分量 x 变化曲线
上边法向应力分量
4
下边法向应力分量
梁中部应力分量 x 最大值为 17.03Mpa。 梁上边法向应力分量最大值为-1.3428Mpa 梁下边法向应力分量最大值为 0.3428Mpa
3 4 节点线性等参单元完全积分： ○
-0.130665
有限元上机实验报告
有限元上机实验报告结构数值分析与程序设计上机实验院系 :土木工程与力学学院专业:土木工程班级:姓名:学号 :指导教师:1、调试教材 P26-30 程序 FEM1。
1.1 、输入数据文件为 :6,4,12,6,1.0E0,0.0,1.0,0.0,13,1,25,2,43,2,56,3,50.0,2.00.0,1.01.0,1.00.0,0.01.0,0.02.0,0.01,3,7,8,10,121.2 、输出数据文件为 :NN NE ND NFIX E ANU T GM NTYPE6 4 12 60.1000E+01 0.000 1.0000.0000E+00 1 NODE X-LOAD Y-LOAD1 0.000000E+00 -0.100000E+012 0.000000E+00 0.000000E+003 0.000000E+00 0.000000E+004 0.000000E+00 0.000000E+005 0.000000E+00 0.000000E+006 0.000000E+00 0.000000E+00NODE X-DISP Y-DISP1 -0.879121E-15 -0.325275E+012 0.879121E-16 -0.125275E+013 -0.879121E-01 -0.373626E+004 0.117216E-15 -0.835165E-155 0.175824E+00 -0.293040E-156 0.175824E+00 0.263736E-15ELEMENT X-STR Y-STR XY-STR1 -0.879121E-01 -0.200000E+01 0.439560E+002 0.175824E+00 -0.125275E+01 0.256410E-153 -0.879121E-01 -0.373626E+00 0.307692E+004 0.000000E+00 -0.373626E+00 -0.131868E+00 2、修改 FEM1,计算 P31例 2-2 。
上机报告--有限元分析
有限元分析上机报告上海电机学院机械学院目录Project1 桁架结构静力有限元分析 (1)Project2 梁结构静力有限元分析 (3)Project3 平面结构静力有限元分析 (5)Project4 实体结构静力有限元分析 (6)有限元分析上机报告Project1 桁架结构静力有限元分析上机地点上机日期指导教师班级学号姓名一、上机目的二、主要仪器设备三、上机结果1、根据任务要求提交图片或者数据；2、提交.log文件；3、验证有限元分析结果。
四、体会与建议有限元分析上机报告Project2 梁结构静力有限元分析上机地点上机日期指导教师班级学号姓名一、上机目的二、主要仪器设备三、上机结果1、根据任务要求提交图片或者数据；2、提交.log文件；3、验证有限元分析结果。
四、体会与建议有限元分析上机报告Project3 平面结构静力有限元分析上机地点上机日期指导教师班级学号姓名一、上机目的二、主要仪器设备三、上机结果1、根据任务要求提交图片或者数据；2、提交.log文件。
四、体会与建议有限元分析上机报告Project4 实体结构静力有限元分析上机地点上机日期指导教师班级学号姓名一、上机目的二、主要仪器设备三、上机结果1、根据任务要求提交图片或者数据；2、提交.log文件。
四、体会与建议出师表两汉：诸葛亮先帝创业未半而中道崩殂，今天下三分，益州疲弊，此诚危急存亡之秋也。
然侍卫之臣不懈于内，忠志之士忘身于外者，盖追先帝之殊遇，欲报之于陛下也。
诚宜开张圣听，以光先帝遗德，恢弘志士之气，不宜妄自菲薄，引喻失义，以塞忠谏之路也。
宫中府中，俱为一体；陟罚臧否，不宜异同。
若有作奸犯科及为忠善者，宜付有司论其刑赏，以昭陛下平明之理；不宜偏私，使内外异法也。
侍中、侍郎郭攸之、费祎、董允等，此皆良实，志虑忠纯，是以先帝简拔以遗陛下：愚以为宫中之事，事无大小，悉以咨之，然后施行，必能裨补阙漏，有所广益。
将军向宠，性行淑均，晓畅军事，试用于昔日，先帝称之曰“能”，是以众议举宠为督：愚以为营中之事，悉以咨之，必能使行阵和睦，优劣得所。
有限元上机实验报告(董妍)
有限元实验报告
董妍-学号 5070309005
1000000，同时时间步长的控制采取根据温度自适应的方式控制步长。 将 Max # increments 设为 1000000，Initial Time Step 设为 0.01，Finish when exceed 设为 900， Max Temperature Change Allowed 设为 10。 前两个参数是在定义为自适应步长控制 时需要给出的希望完成给定时间长度内分析所需的最大时间增量步数和建议初始时间 步长。第三个是指定当所有节点的温度都高于这一指定的基准温度时，程序结束运行。 而最后一个参数是指最大允许的温度改变值——改变这个值可以控制计算结果的精度。
2 / 13
Hale Waihona Puke 有限元实验报告董妍-学号 5070309005
elements 为 0 则说明操作成功。 如果为 200， 则点击 Flip Elements 下的 All Selected 翻转 所有单元后再次检查。 在 SWEEP 命令下点击 ALL，Unused PTS 和 NODES，扫除一些重叠的损坏的点。 再次进入 CHECK 下点击 Cross Elements 察看，如果重叠元素为 0 则操作成功。 进入 RENUMBER 命令下查看目前的节点数和单元数，然后点击 Renumber All 为所有 的节点和单元重新编号以便后续操作。这样，一个有限元模型就建立好了（如图） 。
2 2 3
为了比较不同参数下对内部温度场分布的影响情况， 我计划在实验过程中将内部换热系数、材料导热系数、材 料比热以及模具尺寸作一些变化，考察这几个参数的变化 对温度场分布的影响。参数的选择会在“参数设计”部份 详细说明。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
有限元法基础及应用上机报告南京理工大学2015年12月上机实验一1 实验题目设计一个采用减缩积分线性四边形等参元的有限元模型，通过数值试验来研究网格密度、位移约束条件与总刚度矩阵奇异性、沙漏扩展、求解精度的关系，并验证采用减缩积分时保证总刚度矩阵非奇异的必要条件。
总结出你的研究结论，撰写实验报告。
2 实验目的通过实验来研究减缩积分方案中网格密度和位移约束条件对总体刚度矩阵奇异性和求解精度的影响，以此加深对有限元减缩积分的理解，和对减缩积分中保证总体刚度矩阵非奇异性的认识。
3建模概述先保持位移约束条件不变，研究网格密度对总体刚度矩阵奇异性和求解精度的影响，并验证采用减缩积分时保证总刚度矩阵非奇异的必要条件。
如下图1所示，建立一个简支和链杆的约束条件，然后不断增加网格密度，通过ABAQUS 来计算位移和应力的变化规律。
简支（两个约束）链杆（一个约束）积分点（3个独立关系式）节点（两个自由度）4 计算结果分析讨论与结论 1）1*1单元四边形减缩积分实验载荷 布种/单元应力云图2）2*1单元四边形减缩积分实验载荷 单元应力云图3）4*4单元四边形减缩积分实验载荷布种单元应力云图结果分析5 实验体会与小结单元刚度矩阵的特征：(1)对称性(2)奇异性(3)主元恒正K相同(4)平面图形相似、弹性矩阵D、厚度t相同的单元，eK的分块子矩阵按结点号排列，每一子矩阵代表一个结点，占两行两(5)e列，其位置与结点位置对应。
整体刚度矩阵的特征：(1)对称性(2)奇异性(3)主元恒正(4)稀疏性(5)非零元素呈带状分布。
[K]的物理意义是任意给定结构的结点位移所得到的结构结点力总体上满足力和力矩的平衡。
为消除[K]的奇异性，需要引入边界条件，至少需给出能限制刚体位移的约束条件。
对于一个给定形式的单元，如果采用精确积分，则插值函数中所有项次在|J|=常数的条件下能被精确积分，并能保证刚度矩阵的非奇异性。
如果采用减缩积分，因为插值函数中只有完全多项式的项次能被精确积分，因此需要进行刚度矩阵非奇异必要条件的检查。
上机实验二1 实验题目图示一个简支梁平面应力模型。
梁截面为矩形，高度h=160mm ，长度L=1000mm ，厚度t=10mm 。
上边承受均布压力q =1N/mm 2，材料E=206GPa ，μ=0.29。
X 方向正应力弹性力学理论解为：)534()4(622223-+-=h y h y q y x L h q x σ分别应用3节点三角形单元、4节点线性等参元（完全积分、减缩积分、非协调模式）、8节点二次等参元完全积分进行下列数值实验：1）用较粗单元网格求解梁中部应力分量σx 的最大值和上下边法向应力分量，并对各单元计算精度进行比较分析；2）对粗网格下梁中部最大位移进行对比和分析；3）通过网格加密对比试验3节点三角形单元和8节点二次等参元的收敛速度。
总结出研究结论，撰写实验报告。
2 实验目的（1）通过ABAQUS 软件，用有限元法分析整个梁上的σ_x 的分布规律，讨论σ_y 的有限元解与材料力学解的区别；（2）用有限元法求梁顶边和底边中点正应力σx 的最大值；（3）比较3节点三角形单元、4节点线性等参元（完全积分、减缩积分、非协调模式）、8节点二次等参元完全积分求得的参数与理论解的接近程度； （4）逐步加密单元网格，把有限元法求得的σx 值与理论值进行对比，考察解的收敛性；（5）针对以上力学模型，对比分析3节点三角形平面单元和8节点四边形平面单元的求解精度和收敛性；（6）绘制σx 的误差——计算次数曲线，并进行分析说明。
3 建模概述（1）启用ABAQUS/CAE 程序；（2）创建部件（Module ：Part ），选择2D Planar ，Approximate size=2000，绘制长度为1000，高度为200的矩形；（3）创建材料和截面属性（Module：Property），弹性模量为E=210000 MPa，泊松比为0.3；（4）定义装配件（Module：Assembly），选择Independent；（5）设置分析步（Module：Step）；（6）定义边界条件和载荷（Module：Load），将左右两边界在中点分成两段，在梁的上侧面施加1Mpa的载荷，选取左边界中点设置U1=0和U2=0，选取右边界中点设置U2=0；（7）划分网格（Module：Mesh），全局尺寸设为50，采用3节点线性平面应力三角形单元CPS3划分网格；（8）提交分析作业（Module：Job）；（9）后处理（Module：Visualization）：①显示的应力云图；②查询底边中点的最大值，与理论最大应力值比较考察其精度；③绘制底边上各点的应力的曲线；（10）细化网格验证收敛性，重设单元尺寸为20和10，给出和的应力云图，考察收敛于理论解的程度；（11）高阶单元分析，将单元类型改为CPS8，单元尺寸分别取为100、50、20：①绘制和的云图；②查询底边中点的最大值，对比（10）和（11）中的结果，体会高精度有限元网格的精度和收敛性；③查询模型中从上至下的值。
4 计算结果分析讨论与结论（1）尺寸大小为50的3节点三角形单元计算云图如下：1）x方向应力云图梁中部x方向应力分量最大值为：19.57512）y 方向应力云图上边界xσ方向应力分布曲线：上边界法向最大应力：-1.54372下边界xσ方向应力分布曲线：下边界法向最大应力：0.5437243）y方向位移云图梁中部y方向最大位移：-0.157743（2）尺寸大小为50的四节点完全积分计算云图如下：1）x方向应力云图梁中部x方向应力分量最大值为：28.6935 2）y方向应力云图上边界x σ方向应力分布曲线：上边界法向最大应力：-3.1434下边界xσ方向应力分布曲线：下边界法向最大应力：2.14343）位移云图梁中部y 方向最大位移：-0.188218（3）尺寸大小为50的四节点减缩积分1）x 方向应力云图梁中部x 方向应力分量最大值为：23.33512）y 方向应力云图上边界x方向应力分布曲线：上边界法向最大应力：-1.04348下边界x 方向应力分布曲线：下边界法向最大应力：0.04348213）位移云图梁中部y 方向最大位移：-0.240197（4）尺寸大小为50的四节点非协调模式积分1）x 方向应力云图梁中部x 方向应力分量最大值为：29.27662）y 方向应力云图上边界xσ方向应力分布曲线：上边界法向最大应力：-0.900612下边界xσ方向应力分布曲线：下边界法向最大应力：-0.09938853）位移云图梁中部y方向最大位移：-0.196154（5）尺寸大小为50的八节点二次等参元完全积分1）x方向应力云图梁中部x方向应力分量最大值为：29.5572 2）y方向应力云图上边界x σ方向应力分布曲线：上边界法向最大应力：-1.04859下边界xσ方向应力分布曲线：下边界法向最大应力：0.04859323）位移云图梁中部y 方向最大位移：-0.198043（6）尺寸大小为20的3节点三角形单元1）x 方向应力云图梁中部x 方向应力分量最大值为：25.68172）y 方向应力云图上边界x方向应力分布曲线：上边界法向最大应力：-1.35765下边界x 方向应力分布曲线：下边界法向最大应力：0.3576523）位移云图梁中部y 方向最大位移：-0.187793（7）尺寸大小为10的3节点三角形单元1）x 方向应力云图梁中部x 方向应力分量最大值为：27.90662）y 方向应力云图上边界xσ方向应力分布曲线：上边界法向最大应力：-1.19271下边界xσ方向应力分布曲线：下边界法向最大应力：0.1927053）位移云图梁中部y方向最大位移：-0.195956（8）尺寸大小为20的8节点二次等参元1）x方向应力云图梁中部x方向应力分量最大值为：29.5085 2）y方向应力云图上边界x σ方向应力分布曲线：上边界法向最大应力：-1.01358下边界xσ方向应力分布曲线：下边界法向最大应力：0.01358293）位移云图梁中部y 方向最大位移：-0.199315（9）尺寸大小为10的8节点四边形等参单元1）x 方向应力云图梁中部x 方向应力分量最大值为：27.59812）y 方向应力云图上边界x方向应力分布曲线：上边界法向最大应力：-1.00365下边界x 方向应力分布曲线：下边界法向最大应力：0.003645013）位移云图梁中部y 方向最大位移：-0.200394单元尺寸为50的不同类型单元有限元分析结果：相同类型单元不同网格密度下的有限元分析对比：在同等网格密度下计算水平方向应力分量时，二次单元的计算精度要优于线性单元。
在单元阶数相同的情况下，三角形单元计算精度次于四节点单元。
此外，减缩积分的线性四节点单元因为存在沙漏模式，其应力的计算精度要比同为线性四节点单元的完全积分单元和非协调模式的单元低。
在同等网格密度下计算法向应力分量时，减缩积分的线性四节点单元的计算精度最高，其次依次是八节点二次单元、非协调模式的线性四节点单元、三角形单元、完全积分的线性四节点单元。
用八节点二次单元和线性四节点非协调单元计算应力是比较合适的方案。
在同等网格密度下计算位移，八节点二次单元的精度最高，非协调模式的线性四节点单元的精度次之，这两种单元的误差都在5%以内，而其他单元的计算精度都偏大。
5 实验体会与小结通过本次实验，我基本掌握了使用ABAQUS软件建模求解平面应力问题，学会了查询应力以及位移分布情况的方法。
了解了单元尺寸变化和单元类型不同对结果精度的影响。
单元尺寸越小，精度越高，收敛性越好。
高精度单元收敛速度快，运算精度高。
但要根据实际情况选择合理的单元，提高计算效率。
上机实验三1 实验题目一个矩形平板，长1000mm,宽100mm ，厚度5mm 。
材料的E=200GPa ，0.3=μ， 3-6/107.82mm Kg ⨯=ρ。
板的一对短边简支。
进行下列计算分析并撰写实验报告。
1）在相同的粗网格（厚度方向1层单元）下，分别用8节点线性六面体全积分等参元、8节点六面体非协调元、20节点二次六面体等参元计算其前六阶自由振动频率（Hz ），对计算结果用表格作汇总对比和分析（不要附图）。
2）采用适当网格密度的8节点六面体非协调元，用隐式直接积分法对该平板进行瞬态响应研究。
板的上表面受对称三角形脉冲均布压力，最大值0.05MPa 。
要求：（1）通过计算获得平板下表面中心位置沿长边方向的正应力分量响应最大值随载荷脉宽（0.2毫秒、1.0毫秒、6毫秒、100毫秒、1秒、静载）变化的曲线；（2）板的厚度改50mm ，研究上述规律的变化，并进行归纳和对比分析讨论。