案例：道路工程问题的计算机仿真.docVIP

题号 回 答 内 容 得分 第一题：平面断裂分析（二维断裂问题应力强度因子计算） 问题描述 构件长200mm，高50mm，厚度5mm，距下端中心线左20mm处有一长15mm竖向垂直裂缝，弹性模量2000MP,泊松比0.3,上端中心受竖向向下集中荷载大小为F=1000N,简化为平面应力/应变问题，计算应力强度因子。 ANSYS主要计算过程 定义单元类型 本题简化为平面应力问题，非对称，选用六节点三角形单元PLANE183模拟加载过程，选择Plane strs w/thk，在real constant set number1 ,for PLANE183，THK取5。 图1 Preferences for CUI Filtering 图2 Element types 定义材料性能和参数：main menu|preprocessor|material props|material models-define material behavior Ex=2000,prxy=0.3 图3 定义材料参数 创立几何模型、划分网格 创建关键节点KEYPOINT1（100,0），2（100,50），3（0,50），4（-50,50），5（-100,0），6（-20,0），7（-20,15），8（-20,0）。关键节点6,8重合，但是必须定义为两个点。 通过关键节点创造面：面1由关键点34567组成，面2由关键点12378组成。 图4 关键节点 图5创建面 图6 网格划分 对7号节点选择为尖端，设置如下： 图7 KSCON命令 图8 划分网格 通过PlotCtrls--HardCopy--TO FILE 将网格划分后的结果存为JPG格式。 加载求解 本题加载考虑比较简单，简支梁，关键节点5处xy约束，关键节点1y方向约束，节点3处受竖直向下集中力大小为1000N。 图9 荷载和约束 利用KCALC计算应力强度因子时x轴必须平行于断纹面，y轴垂直于断纹面，坐标系必须为激活的模型坐标系和结果坐标系。 求解的时候自定义局部坐标系11：local ,11,0,0,0,0,90（将整体坐标系绕z轴旋转90°）得到新的坐标系11，将工作坐标系改到11，csys，11 rsys,11 ,模型坐标系和节点坐标系均为11。 图10 局部坐标系下的模型 选择Main Menu|Solution|Solve Current LS,单击OK，ANSYS开始求解计算。 查看求解结果 等效应力等值线图 Main Menu|General Plot Results|Contour Plot|Nodal Solu,在Contour Nodal Solution Data对话框 Item to be Contoured 列表框中选择Nodal Solution|stress|von Mises stress,单击OK。 图11 等效应力等值线图 应力强度因子 图12、13 节点编号 图14 PATH 图15KCSYM Model type 选择全裂纹 定义裂纹面路径：Main Menu|General Postproc|Path Operations|Define Path|by Nodes,拾取路径上的点，其中裂纹尖端为第一点，整个裂纹面四个附加点，两个点沿一个裂纹面，另两个点沿另一个裂纹面。所选点x坐标应相同，即x轴和裂纹面平行。 最终计算结果： *** CALCULATE MIXED-MODE STRESS INTENSITY FACTORS **** ASSUME PLANE STRESS CONDITIONS ASSUME A FULL-CRACK MODEL (USE 5 NODES) EXTRAPOLATION PATH IS DEFINED BY NODES: 158 165 162 164 161 WITH NODE 158 AS THE CRACK-TIP NODE USE MATERIAL PROPERTIES FOR MATERIAL NUMBER 1 EX = 2000.0 NUXY = 0.30000 AT TEMP = 0.0000 **** KI = 0.62191 , KII = 0.28328 , KIII = 0.0000 **** 三、个人做题小结 通过这个例子我对ANSYS软件使用有了一定的认识，从最开始的畏惧，一点一点摸索，终于形成这份成果，虽然模型很粗糙，设定也简化了不少，将梁简化为类似平板，但对我而言，