研究生有限元作業.docVIP

研究生作业 问题1 以悬臂梁为例，通过有限元计算，比较三角元，矩形元和8节点等参元的计算精度（分别用位移合应力作为标准，平面应力问题）。取悬臂梁长50mm，截面宽1mm，高2mm,E=210GPa,P=1N,泊松比0.3. 分析过程： 一：建立模型 Preferences中选择Structural ; PreprocessorElementTypeAdd/Edit/Delete中选择Solid—Quad 4node 42Apply, Solid—Quad 8node 183Apply, Solid—Triangle 6node 2OK默认分别为1，2，3种单元; 分别在OptionElement behavior 中选择Plane strs w/thk ; 在Real ConstantsAdd/Edit/Delete设定厚度为1mm 对应3种实常数; Material PropsMaterial ModelsStructuralLinearElasticIsotropic中EX=2.1e5 Mp, PRXY=0.3；ModelingCreateAreasRectangleBy Dimensions 输入X1=0,X2=50,Y1=0,Y2=2 SAVE　DB（之后可从该步直接调用） MsehingMeshToolElement Attributes中选择Global—Set出现对话框，选择单元类型及对应的常数，先选择第1种。（这里我们之前设置了3种，都可选择，注意对应关系） Size ControlsLinesSet选择矩形两条长边，在Element edge length 中设置0.5；重复同样操作，将宽作同样设置 MeshAreas;Shape--Quad,Free；点击Mesh,完成网格划分。 这样划分的节点数为505，单元数400. 二：加载 LoadsDefine LoadsApplyStructuralDisplacementOn Lines选择左边矩形宽，约束所有自由度，位移值为0。LoadsDefine LoadsApplyStructuralForce/MomentOn Keypoints这里PlotKeypoints,显示矩形的四个关键点，然后选中右上角的关键点，出现对话框。 选择方向FY，value为-1. SAVE　DB 三：理论计算及有限元结果 SolutionSolveCurrent LS General PostprocPlot ResultsDeformed Shape选择Def + undeformed 理论计算： 由材料力学悬臂梁挠度 该平面应力梁最左端表层处拉伸应力为： General PostprocPlot ResultsContour PlotNodal SolutionStress---von Mises stressOK 最大值处Y向位移相对误差 合应力最大值相对误差 以下为矩形8节点单元183所得Y向位移及合应力结果： 最大值处Y向位移相对误差 合应力最大值相对误差 图 1-5 8节点矩形单元Y向节点位移 图 1-6 8节点矩形单元节点合应力 采用6节点三角形单元，为了保证节点数仍为505，网格划分是单元变长为1mm,得到节点数505，单元数200的三角形网格。 图 1- 7 6节点三角形网格 最大值处Y向位移相对误差 合应力最大值相对误差 图 1-8 6节点三角形单元Y向节点位移 图 1-9 6节点三角形单元节点合应力 四：结果比较及分析 节点数相同的情况下结果对比 单元形状 单元数 节点数 最大Y向位移误差 最大合应力误差 6节点三角形 200 505 0.0071% 0.0507% 4节点矩形 400 505 0.0682% 0.5413% 8节点矩形 126 509 0.0235% 0.9800% 结论： 在节点数相同的情况下，6节点三角形单元的模拟相对误差最小，最精确。 问题2 一端固定的欧拉压杆，长度为l=25.4mm，截面宽度b=1mm和厚度h=0.51mm，弹性模量E=210GPa,端部挠度w与载荷P以及端部挠度w与轴向位移u的关系为 其中，考虑轴向线性变形，u的表达式修正为，用有限元计算杆的后屈曲路径，并与理论解进行比较，确定误差。 分析过程： 一 建立有限元模型 1.选择单元 Beam 2D plastic23,设置单元实常数mm3 截面积 mm2，梁高； 2．选择材料属性为结构线弹性各项同性材料，杨氏模量E=2.1e5,泊松比； 3.先创建起始1及末尾21两个节点，坐标分别为（0，0，0），（25.4，0，0），然后在之间填充19个节点；将1和2点创建成一个单元，然后COPY生成20个单