实验一坝体的有限元建模及应力应变分析一、实验目的:1、掌握ANSYS软件基本的几何形体构造方法、网格划分方法、边界条件施加方法及各种载荷施加方法。2、熟悉有限元建模、求解及结果分析步骤和方法。3、能利用ANSYS软件对结构进行有限元分析。二、实验设备:微机,ANSYS软件。三、实验内容:计算分析模型如图所示,分析坝体的应力、应变。1m5m0.55m四、实验步骤:1进入ANSYS程序→ANSYS→changetheworkingdirectoryintoyours→inputInitialjobname:dam2设置计算类型ANSYSMainMenu:Preferences→selectStructural→OK3选择单元类型ANSYSMainMenu:Preprocessor→ElementType→Add/Edit/Delete→Add→selectSolidQuad4node42→OK(backtoElementTypeswindow)→Options…→selectK3:PlaneStrain→OK→Close(theElementTypewindow)4定义材料参数ANSYSMainMenu:Preprocessor→MaterialProps→MaterialModels→Structural→Linear→Elastic→Isotropic→inputEX:2.1e11,PRXY:0.3→OK5生成几何模型生成特征点ANSYSMainMenu:Preprocessor→Modeling→Create→Keypoints→InActiveCS→依次输入四个点的坐标:input:1(0,0),2(1,0),3(1,5),4(0.45,5)→OK生成坝体截面ANSYSMainMenu:Preprocessor→Modeling→Create→Areas→Arbitrary→ThroughKPS→依次连接四个特征点,1(0,0),2(1,0),3(1,5),4(0.45,5)→OK生成坝体截面如图一图一6网格划分ANSYSMainMenu:Preprocessor→Meshing→MeshTool→(SizeControls)lines:Set→依次拾取两条横边:OK→inputNDIV:15→Apply→依次拾取两条纵边:OK→inputNDIV:40→OK→(backtothemeshtoolwindow)Mesh:Areas,Shape:Quad,Mapped→Mesh→PickAll(inPickingMenu)→Close(theMeshToolwindow)图二7模型施加约束分别给下底边和竖直的纵边施加x和y方向的约束ANSYSMainMenu:Solution→DefineLoads→Apply→Structural→Displacement→Onlines→pickthelines→OK→selectLab2:UX,UY→OK给斜边施加x方向的分布载荷ANSYS命令菜单栏:Parameters→Functions→Define/Edit→1)在下方的下拉列表框内选择x,作为设置的变量;2)在Result窗口中出现{X},写入所施加的载荷函数:1000*{X};3)FileSave(文件扩展名:func)→返回:Parameters→Functions→Readfromfile:将需要的.func文件打开,任给一个参数名,它表示随之将施加的载荷→OK→ANSYSMainMenu:Solution→DefineLoads→Apply→Structural→Pressure→OnLines→拾取斜边;OK→在下拉列表框中,选择:Existingtable→OK→选择需要的载荷参数名→OK图三8分析计算ANSYSMainMenu:Solution→Solve→CurrentLS→OK(toclosethesolveCurrentLoadStepwindow)→OK9结果显示(1)变形情况ANSYSMainMenu:GeneralPostproc→PlotResults→DeformedShape…→selectDef+Undeformed→OK(backtoPlotResultswindow)如图四;图四(2)位移情况GeneralPostproc→PlotResults→ContourPlot→NodalSolu…→select:DOFsolution,UX,UY,Def+Undeformed,如图五、六、七所示图五X方向位移图图六Y方向位移图图七总的位移图(3)压力结果图GeneralPostproc→PlotResults→ContourPlot→NodalSolu→Stress,SX,SY,SZ,Def+Undeformed→OK如图八、九、十所示:图八X方向应力图图九Y方向应力图图十Z方向的应力图10结果分析由以上实验的变形和位移的云图可知,坝体最上边的变形量最明显,最大,主要是因为最上边所受载荷最大,故最大位移量发生在坝体的顶部。由应力云图可知,X轴最大应力发生在坝体最上边的最右边,而Y轴最大应力发生在坝底的两个拐角处:Z轴方向所受应力与X轴方向相同,故应加固坝底拐角和坝顶最外拐角的强度,从而能提高坝体的强度。实验二平板孔的有限元建模及应力应变分析如图所示,使用ANSYS分析平面带孔平板,分析在均布载荷作用下板内的应力分布。已知条件:F=20N/mm,L=200mm,b=100mm,圆孔半径r=20,圆心坐标为(100,50),E=200Gpa。板的左端固定。图1-2带孔平板模型1.问题描述实例类型:ANSYS结构分析。分析类型:线性静力分析。单元类型:PLANE82ANSYS功能示例:实体建模包括基本的建模操作,布尔运算和网格细化;施加均布载荷;显示变形后形状和应力等值线图、单元信息列表;基本的结果验证技巧。ANSYS帮助文件:在ANSYSStructuralAnalysisGuide了解StructuralStaticAnalysis分析知识,在ANSYSElementsReference部分了解Plane82单元的详细资料。2.建立有限元模型1).建立工作目录并添加标题以Interactive方式进入ANSYS,File菜单中设置工作文件名为Plane、标题为plane。2).定义材料属性其操作如下:GUI:PreProcessorMaterialPropsMaterialmodelsStructuralLinearElasticIsotropic在弹出对话框中键入EX=200000(单位Mpa),PRXY=0.3。3)定义单元实常数操作如下:选择单元对于任何分析,必须在单元类型库中选择一种或者多种合适的单元类型。单元类型决定了附加的自由度(位移,转角、温度等)。许多单元还需要设置一些单元选项,比如单元特性和假设。单元结果的打印输出选项等,对于本问题选择Plane82单元。选择单元得操作如下:GUI:PreProcessorMenuElementTypeAdd/Edit/Delete选择Plane82,单击OK按钮。GUI:PreProcessorMenuElementTypeAdd/Edit/DeleteOptions在K3项后面下拉菜单中选“Planestrsw/thk“。单元厚度为20mm,定义单元厚度操作如下:GUI:PreProcessorMenuRealConstantsAdd/Edit/DeleteAdd在弹出的对话框中THK后文本框中键入材料厚度值20。4).创建实体模型(1)创建关键点其操作如下:GUI:PreProcessorModelingCreatekeypiontinactiveCS1(0,0,0),2(0.1,0,0),3(0.2,0,0),4(0.2,0.05,0),5(0.2,0.1,0),6(0.1,0.1,0),7(0,0.1,0),8(0,0.05,0),9(0.08,0.05,0),10(0.1,0.03,0),11(0.12,0.05,0),12(0.1,0.07,0),13(0.1,0.05,0)(2)创建线其操作如下:GUI:PreProcessorModelingCreatelineslinesstraightline;(3)创建面其操作如下:GUI:PreProcessorModelingCreateAerasArbitarybylines(4)下面通过布尔“粘”操作生成圆孔,其操作如下:GUI:ProcessorModelingOperateBooleansglueAreas单击pickall,单击OK按钮,布尔操作完毕之后,实体模型为带孔平板。如图所示图1平板模型5).划分网格如图1可知将平板均分为四块,每一块分网相同,所以选择小块长边输入NDIV值40,短边为10,与圆弧相连长边NDIV值为8,短边为5,且圆弧上值设为5,同理,另外三块相同区域设置相同NDIV值,划分网格结果如图所示,其操作如下:GUI:PreProcessorMeshingmeshtoolmapmeshpickall如图2、图2平板划分网格图3施加载荷并求解1).定义约束由已知得,需要固定(Fix)板左边线,即需要约束线上节点所有自由度(AllDOFs),其操作如下:GUI:SolutionDefineLoadsApplyStructuralDisplacementOnLines弹出ApplyU,ROTonLines对话框。选择板左侧边线,在Lab2栏选AllDOF。单击Apply按钮。2).施加载荷在板右侧边施加均布载荷,载荷大小为20/20=1Mpa,施加载荷操作如下:GUI:SolutionDefineLoadsApplyStructuralPressureOnLines弹出ApplyPRESonLines对话框,键入载荷值-1(由于载荷方向离开板,为拉力所以为负值,反之为正),单击OK按钮。如图3图3平板施加载荷、约束3.求解GUI:SolutionSolveCurrentLS4.查看分析结果下面查看分析结果,对于静力分析主要是模型位移及等效应力等值线图或者节点结果数据列表。1).显示模型变形图其操作如下:GUI:GeneralPostprocPlotResultsDeformedShape如图4图4平板变形图2).显示位移等值线分布图其操作如下:GUI:GeneralPostprocPlotResultsNodalSolutionDoFSolutionDisplacementVectorsum如图5图5平板位移云图3).显示等效应力等值线图其操作如下:GUI:GeneralPostprocPlotResultsContourPlotNodalSoluStressVonMises如图6图6平板应力云图5结果分析板的最大变形量为0.124X10-5m,从ansys分析结果图中可以看出,位移最大的位置为平板圆孔右端,即发生最大变形量的位置为平板的最右段,即力作用的地方,左端为固定端,故最左端的位移为0。所受的最大应力为3.794pa,同时,由图发现最大应力出现在孔的上下两顶点,与实际情况相同,可知分析正确。实验总结:1、建模:熟悉基本的建模操作.2、施加载荷和求解:掌握定义对称约束和施加均布载荷的操作步骤,可以根据结构及其载荷的分布判断是否可以通过其对称性适当地简化模型。3、查看分析结果:掌握显示变形图形和应力等值线图的操作,掌握动画显示参数变化的操作步骤,能够显示变形动画;了解基本的验证技巧,特别是通过网格细化并进行结果对比来验证网格密度是否合理。附件I:实验一log文件/BATCH/COM,ANSYSRELEASE12.0.1UP2009041514:11:2606/22/2012/FILNAME,dam,1/TITL