Ansys 分析在弹性静力学问题中的应用.docVIP

Ansys 分析在弹性静力学问题中的应用 yihe2212@ 问题：矩形含孔薄板，边界受拉伸，求孔边最大应力。 尺寸：板长100mm ，宽50mm ，孔半径：10mm 边界受均匀拉伸：1Mpa 材料属性：弹性模量200Gpa ，泊松比 0.3 1．选定问题所属领域 命令 Preferencestructural 2．选择4节点的四边形单元 命令 Preference PreprocessorElement TypeAdd/Edit/DeleteAdd 3．选择模型材料特性（弹性模量、泊松比） 命令 Preference Preprocessor Material Props Material Models Structural Linear Elastic Isotropic 设置参数EX （弹性模量）为200GPa ，设置参数PRXY （泊松比）为0.3 4．几何建模 ? 关键点 PreferencePreprocessorModelingCreateKeypointsIn Active CS 或者使用命令K,1,0,0,0 依次做出以下5个Keypoints 1(0,0,0, 2(0.1,0,0, 3(0.1,0.05,0, 4(0,0.05,0, 5(0.05,0.025,0 ? 线 PreferencePreprocessorModelingCreateLinesLinesStraight Line 或者使用命令K,1,1,2（绘制关键点1到关键点2的线段），并将该线段编号为 1 ? 面 Main MenuPreprocessorModelingCreateAreasArbitraryBy Lines 或者使用命令AL,L1,L2,L3,L4,L5,L6,L7,L8 得到的几何图形如下所示： 5．划分网格 命令 Main Menu Preprocessor Meshing Size Cntrls ManualSize Lines Picked Lines L2,L4划分为50个单元，L1,L3划分为100个单元，L5,L6,L7,L8划分为50个单元，网格划分结果如下图所示： 6．施加约束 在矩形平板左侧边界施加约束 命令 Main Menu Preprocessor Loads Define Loads Apply Structural Displacement On Lines 7．施加载荷 矩形平板右侧边界受均匀拉伸，1Mpa 命令 Main Menu Preprocessor Loads Define Loads Apply Structural Pressure On Lines 8．求解 命令 Main MenuSolutionSolveCurrent LS 9．检查结果 命令 Main MenuGeneral PostprocPlot ResultsContour PlotNodal Solu 得到矩形平板的应变分布图，如下： 从图中可看出，应变分布从左往右依次增大。 再查看结果的应力分布图，结果如下： 从图中可以看出，危险截面位于圆孔的最顶端和最底端处。 yihe2212@ 圆孔孔边应力集中的经典解 设该矩形薄板只在两边受有均布拉力 q, 如图所示，则可得出该问题的基尔斯的 解答： q r2 q r2 r2 σ ρ = (1 ? 2 + cos 2? (1 ? 2 (1 ? 3 2 , 2 2 ρ ρ ρ q r2 q r4 σ ? = (1 + 2 ? cos 2? (1 + 3 4 , 2 2 ρ ρ τ ρ? = τ ?ρ 沿着孔边， ρ q r2 r2 = ? sin 2? (1 ? 2 (1 + 3 2 . 2 ρ ρ = r ，环向正应力是 σ ? = q(1 ? 2 cos 2? 它的几个重要数值如下表所示： ? σ? 0o 30 o 0 45 o q 60 o 2q 90 o 3q ?q 画出孔边的应力分布如下图所示： -6- yihe2212@ Ansys 分析结果与理论解结果对比分析 通过 List---Results---Nodal Solution 得到各个节点的压力分布表， 取其中几个列表 如下： NODE 301 302 352 402 SX -19784 3.77E+06 3.76E+06 -19516 SY -1.45E+06 55505 55505 -1.43E+06 SZ 0 0 0 0 SXY 173.42 408.76 966.25 172.36 SYZ 0 0 0 0 SXZ 0 0 0 0 其中，各节点对应如下图： 302 402 352 经典解所对应的结果： NODE 301 302 352 402 301 σ -1E+06