材料成型CAE综合实验课程设计.doc

学 号： 课 程 设 计题 目材料成型cae综合实验学 院材料科学与工程学院专 业材料成型及控制工程班 级姓 名指导教师2011年12月16日实验课程名称：材料成型cae综合实验一 实验项目名称桁架的结构分析实验成绩实 验 者专业班级组 别同 组 者实验日期2011年12月13日一部分：实验预习报告（包括实验目的、意义，实验基本原理与方法，主要仪器设备及耗材，实验方案与技术路线等）一）实验目的掌握用ansys进行有限元结构分析的基本过程。二）基本原理和方法结构分析是有限元分析的基础内容，也是材料成型过程模拟的主要内容。有限元法是一种离散化的数值计算方法。离散后的单元和单元之间只通过节点相联系，所有场变量（位移、应力、温度等）都通过节点进行计算。对于每个单元，选取适当的插值函数，使得在子域内部、子域分界面上以及子域与外界分界面上都满足一定的条件。然后把所有单元的方程都组装起来，就得到整个结构的方程组。求解方程组，就可以得到方程的近似解。有限元刚度分析法的步骤：1）建立几何模型；2）对几何模型进行离散化处理；3）将单元节点位移作为基本未知量；4）选择位移模式；5）确定单元应变与位移、应力与应变的关系；6）根据虚功原理建立单元中节点力与节点位移的关系；7）根据作用力等效原则将每个单元所受的载荷移置到该单元的节点上；8）将各单元的刚度方程叠加，组装成整体刚度方程；9）根据边界条件修改刚度方程，消除刚体位移；10）求解整体刚度方程，得到节点位移；11）根据相应方程求解应力和应变；12）利用计算机图形方式，将计算结果以变形网格、等值线、彩色云图、动画等方式进行显示与分析。从应用角度看（如用ansys分析软件），整个过程可分为：l 前处理：步骤1）、2）l 加载求解：步骤7）、9）、10）、11）l 后处理：步骤12）其余为软件自动进行。三）实验内容桁架结构如图所示，各节点等距离分布，1、2节点间距1m，集中力p=10kn，确定桁架各杆件内力。第二部分：实验过程记录（可加页）（包括实验原始数据记录，实验现象记录，实验过程发现的问题等）上机实验步骤（gui法，参考）建立文件：utility menufilechange directoryinto yoursokutility menuj fileobnametrussok定义单元类型：main menupreferencestrucrturalokmain menupreprocessorelement typeadd/edit/deleteaddlink 2d spar 1ok closemain menupreprocessorreal constants add/edit/deleteaddlink 1okarea: 0.1ok close定义材料属性：main menupreprocessormaterial propsmaterial modelsstrucrurallinearelastic isotropicex: 30e6, prxy: 0.3ok生成节点：main menupreprocessormodelingcreatenodesin active csnode 1: 0,0,0apply node 2: 1,0,0apply 同理：node 3: 2,0,0node 4: 3,0,0node 5: 4,0,0node 6: 3,1,0node 7: 2,1,0node 8: 1,1,0okutility menuplotctrlsnumberingnode node numbers onok生成单元：main menupreprocessormodelingcreateelementsauto numberedthru nodespick 1,2applypick 2,3apply同理：3,44,55,64,63,63,72,81,87,86,7ok施加位移约束和载荷：main menusolutiondefined loadsapplystructuraldisplacementon nodespick 1,3lab 2: uyapplypick 5lab 2: ux, uyokmain menusolutiondefined loadsapplystructuralforce/momenton nodespick 2,4,7lab: fy, value: -10000ok求解：main menusolutionsolvecurrent lsokclose后处理：绘制变形图：main menugeneral postprocessorplot resultsdeformed shapedef shape onlyok（对显示结果进行存储：utility menuplot ctrlscapture imagefilesave as）显示节点位移：main menugeneral postprocessorlist resultsnodal solutiondof solution displacement vector sunok显示节点受力：main menugeneral postprocessorlist resultsnodal loadsall stru forc fok（或根据讨论的需要进行显示）教师签字_第三部分 结果与讨论（可加页）一、实验结果分析（包括数据处理、实验现象分析、影响因素讨论、综合分析和结论等）二、小结、建议及体会三、思考题实验结果可以显示出构件的变形情况，各节点的位移，可以看出最大位移为dmx=0.35576；小结与体会通过本实验，我认识到通过采用ansys软件对桁架的分析，可以简易直观在了解到结构的变形情况。只要输入结构的形式，约束条件和受力情况就可以得出结构的变形情况，很方便也很直观。实验课程名称：材料成型cae综合实验二 实验项目名称自主设计构件的结构分析实验成绩实 验 者专业班级组 别同 组 者实验日期2011年12月13日一部分：实验预习报告（包括实验目的、意义，实验基本原理与方法，主要仪器设备及耗材，实验方案与技术路线等）一）实验目的通过自主设计的构件，熟悉用ansys进行有限元结构分析的过程。二）基本原理和方法结构分析是有限元分析的基础内容，也是材料成型过程模拟的主要内容。有限元法是一种离散化的数值计算方法。离散后的单元和单元之间只通过节点相联系，所有场变量（位移、应力、温度等）都通过节点进行计算。对于每个单元，选取适当的插值函数，使得在子域内部、子域分界面上以及子域与外界分界面上都满足一定的条件。然后把所有单元的方程都组装起来，就得到整个结构的方程组。求解方程组，就可以得到方程的近似解。有限元刚度分析法的步骤：1）建立几何模型；2）对几何模型进行离散化处理；3）将单元节点位移作为基本未知量；4）选择位移模式；5）确定单元应变与位移、应力与应变的关系；6）根据虚功原理建立单元中节点力与节点位移的关系；7）根据作用力等效原则将每个单元所受的载荷移置到该单元的节点上；8）将各单元的刚度方程叠加，组装成整体刚度方程；9）根据边界条件修改刚度方程，消除刚体位移；10）求解整体刚度方程，得到节点位移；11）根据相应方程求解应力和应变；12）利用计算机图形方式，将计算结果以变形网格、等值线、彩色云图、动画等方式进行显示与分析。从应用角度看（如用ansys分析软件），整个过程可分为：l 前处理：步骤1）、2）l 加载求解：步骤7）、9）、10）、11）l 后处理：步骤12）其余为软件自动进行。三）实验内容1）参照实验一，自主设计一个构件的结构及其受力和约束状况。2）用ansys进行有限元结构分析。3）对分析结果进行讨论。桁架结构如图所示，各节点等距离分布，1、2节点间距1m，集中力p=10kn，确定桁架各杆件内力。第二部分：实验过程记录（可加页）（包括实验原始数据记录，实验现象记录，实验过程发现的问题等）上机实验步骤（gui法，参考）建立文件：utility menufilechange directoryinto yoursokutility menuj fileobnametrussok定义单元类型：main menupreferencestrucrturalokmain menupreprocessorelement typeadd/edit/deleteaddlink 2d spar 1ok closemain menupreprocessorreal constants add/edit/deleteaddlink 1okarea: 0.1ok close定义材料属性：main menupreprocessormaterial propsmaterial modelsstrucrurallinearelastic isotropicex: 30e6, prxy: 0.3ok生成节点：main menupreprocessormodelingcreatenodesin active csnode 1: 0,0,0apply node 2: 1,0,0apply 同理：node 3: 2,0,0node 4: 1,1,0okutility menuplotctrlsnumberingnode node numbers onok生成单元：main menupreprocessormodelingcreateelementsauto numberedthru nodespick 1,2applypick 2,3apply同理：1,42,43,4ok施加位移约束和载荷：main menusolutiondefined loadsapplystructuraldisplacementon nodespick 1lab 2: uyapplypick 3lab 2: ux, uyokmain menusolutiondefined loadsapplystructuralforce/momenton nodespick 4lab: fy, value: -10000ok求解：main menusolutionsolvecurrent lsokclose后处理：绘制变形图：main menugeneral postprocessorplot resultsdeformed shapedef shape onlyok（对显示结果进行存储：utility menuplot ctrlscapture imagefilesave as）显示节点位移：main menugeneral postprocessorlist resultsnodal solutiondof solution displacement vector sunok显示节点受力：main menugeneral postprocessorlist resultsnodal loadsall stru forc fok（或根据讨论的需要进行显示）教师签字_第三部分 结果与讨论（可加页）一、实验结果分析（包括数据处理、实验现象分析、影响因素讨论、综合分析和结论等）二、小结、建议及体会三、思考题在3点约束两个自由度，在1点处约束一个自由度，在4点加载向下的力，得到如上图，最大形变量dmx=.0.06595；小结与体会通过本实验，我认识了ansys在结构受力分析中的应用，可以大大的节省手工计算时间，还可以更加直观地看到结构的变形情况。实验课程名称：材料成型cae综合实验三 实验项目名称圆管的温度场和热应力分析实验成绩实 验 者专业班级组 别同 组 者实验日期2011年12月13日一部分：实验预习报告（包括实验目的、意义，实验基本原理与方法，主要仪器设备及耗材，实验方案与技术路线等）一）实验目的掌握用ansys进行温度场分析和热应力分析的基本过程。二）基本原理和方法 温度场分析和热应力分析是材料成型过程分析的重要内容，对材料成型cae起着至关重要的作用。有限元法是一种离散化的数值计算方法。离散后的单元和单元之间只通过节点相联系，所有场变量（位移、应力、温度等）都通过节点进行计算。对于每个单元，选取适当的插值函数，使得在子域内部、子域分界面上以及子域与外界分界面上都满足一定的条件。然后把所有单元的方程都组装起来，就得到整个结构的方程组。求解方程组，就可以得到方程的近似解。用ansys软件进行有限元分析，整个过程可分为：l 前处理：建立几何模型；对几何模型进行离散化处理等。l 加载求解：根据作用力等效原则将每个单元所受的载荷移置到该单元的节点上；根据边界条件修改刚度方程，消除刚体位移；求解整体刚度方程，得到节点位移；根据相应方程求解应力和应变等。l 后处理：利用计算机图形方式，将计算结果以变形网格、等值线、彩色云图、动画等方式进行显示与分析等。本实验主要进行热分析和热力耦合分析，将温度场作为载荷施加给所分析的构件。三）实验内容有一截面为圆环形的输暖管道，内外半径分别为200mm、800mm，管道内水的温度为80，管道外表层温度为10，求管道内的应力分布。（假设管道内充满水）该问题为轴对称问题，沿管道横截面取宽为50mm的矩形截面（如图）为计算模型。第二部分：实验过程记录（可加页）（包括实验原始数据记录，实验现象记录，实验过程发现的问题等）上机实验步骤（gui法，参考）建立文件：utility menu filechange directoryinto yoursokutility menu filejobnamethermalstressok建模：main menupreferencestrucrtural and thermalokmain menupreprocessorelement typeadd/edit/deleteaddthermal solidquad 4 node 55ok element typeoptionsk3: axisymmetricokmain menupreprocessormaterial propstemperature unitcelsuisokmain menupreprocessormaterial propsmaterial modelsthermalconductivity isotropickxx: 1.2okmain menupreprocessormodelingcreate areasrectangleby dimensionsx1, x2: 0.2, 0.8. y1, y2: 0, 0.05ok划分网格：main menupreprocessormeshingmeshingtoolsize controls: lines setspick 1,3okndiv: 10ok back to “meshtool windows” size controls: lines setspick 2,4okndiv: 2ok back to “meshtool windows”shape: quad freemeshpick allokclose “meshtool windows”热分析：main menusolutionanalysis typenew analysissteady stateokutility menuselectentitieslines, by num/pick, from fullapplypick 4okutility menuselectentitiesnodes, attached to, lines,allokmain menusolutiondefine loadsapplythermaltemperatureon nodespick alloklab 2: temp, value: 80okutility menuselectentitieslines, by num/pick, from fullapplypick 2okutility menuselectentitiesnodes, attached to, lines,allokmain menusolutiondefine loadsapplythermaltemperatureon nodespick alloklab 2: temp, value: 10okutility menuselecteverythingmain menusolutionload step optsoutput ctrlsdb/results filebasic quantities, last substep, value of n: 1okmain menusolutionsolvecurrent lsokclosemain menugeneral postprocplot resultscontour plotnodal solunodal solutiondof solutionnodal temperaturedeformed shape onlyok （或根据讨论的需要进行显示，对显示结果进行存储：utility menuplot ctrlscapture imagefilesave as）热应力耦合分析：main menupreprocessorelement typeswitch elem typethermal to strucokmain menupreprocessorelement typeadd/edit/deletadd optionsk3: axisymmetric okmain menupreprocessormaterial propsmaterial modelsstructurallinearelastic isotropicex: 1.2e11, prxy: 0.3 back to “define mareial model behavior”windows: structuralthermal expansionsecant coefficientisotropicalpx: 1.3e-6okutility menuselectentitieslines, by num/pick, from fullapplypick 3okutility menuselectentitiesnodes, attached to, lines,allokmain menupreprocessorcoupling/ceqncouple dofspick alloknset: 7, lab: uy okutility menuselecteverythingmain menupreprocessorcoupling/ceqncouple dofspick 1, 14, 24oknset: 8, lab: uxokmain menupreprocessorcoupling/ceqncouple dofspick 2, 12, 13oknset: 9, lab: uxokmain menusolutionanalysis typenew analysisstaticokmain menusolutiondefine loadsapplystructuraldisplacementon nodespick 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11lab2: uy, value: 0okmain menusolutiondefine loadsapplystructuraltemperaurefrom thermal analysisfname name of results file: thermalstress.rthokmain menusolutionsolvecurrent lsokclosemain menugeneral postprocplot resultscontour plotnodal solunodal solution stressx-component of stressdeformed shape onlyok （或根据讨论的需要进行显示，对显示结果进行存储：utility menuplot ctrlscapture imagefilesave as）教师签字_第三部分 结果与讨论（可加页）一、实验结果分析（包括数据处理、实验现象分析、影响因素讨论、综合分析和结论等）二、小结、建议及体会三、思考题第一幅图是温度场等值线图，现实的是内外径之间的温度场的分布，从10到80的渐变过程，温度逐渐升高；第二幅图是应力场图，其最大位移是dmx=.325e-4,smn=-.248e-5,smx=1105;随着温度的增加，管由内至外得热应力也逐渐变大；第三幅图是应力的分布。小结与体会通过本次实验，我了解到采用ansys对构件的温度与热应力分析，可以显示出的图片通过颜色的渐变来反映出温度与应力的变化过程，还可以得出应力壮态。实验课程名称：材料成型cae综合实验四 实验项目名称自主设计构件的温度场和热应力分析实验成绩实 验 者专业班级组 别同 组 者实验日期2011年12月13日一部分：实验预习报告（包括实验目的、意义，实验基本原理与方法，主要仪器设备及耗材，实验方案与技术路线等）一）实验目的通过自主设计的构件，熟悉用ansys进行有限元温度场和热应力分析的过程。二）基本原理和方法温度场分析和热应力分析是材料成型过程分析的重要内容，对材料成型cae起着至关重要的作用。有限元法是一种离散化的数值计算方法。离散后的单元和单元之间只通过节点相联系，所有场变量（位移、应力、温度等）都通过节点进行计算。对于每个单元，选取适当的插值函数，使得在子域内部、子域分界面上以及子域与外界分界面上都满足一定的条件。然后把所有单元的方程都组装起来，就得到整个结构的方程组。求解方程组，就可以得到方程的近似解。用ansys软件进行有限元分析，整个过程可分为：l 前处理：建立几何模型；对几何模型进行离散化处理等。l 加载求解：根据作用力等效原则将每个单元所受的载荷移置到该单元的节点上；根据边界条件修改刚度方程，消除刚体位移；求解整体刚度方程，得到节点位移；根据相应方程求解应力和应变等。l 后处理：利用计算机图形方式，将计算结果以变形网格、等值线、彩色云图、动画等方式进行显示与分析等。本实验主要进行热分析和热力耦合分析，将温度场作为载荷施加给所分析的构件。三）实验内容1）参照实验三，自主设计一个构件的结构及其受热和约束状况。2）用ansys进行有限元温度场和热应力分析。3）对分析结果进行讨论。有一截面为圆环形的输暖管道，内外半径分别为500mm、2000mm，管道内水的温度为90，管道外表层温度为20，求管道内的应力分布。（假设管道内充满水）该问题为轴对称问题，沿管道横截面取宽为100mm的矩形截面（如图）为计算模型。第二部分：实验过程记录（可加页）（包括实验原始数据记录，实验现象记录，实验过程发现的问题等）上机实验步骤（gui法，参考）建立文件：utility menu filechange directoryinto yoursokutility menu filejobnamethermalstressok建模：main menupreferencestrucrtural and thermalokmain menupreprocessorelement typeadd/edit/deleteaddthermal solidquad 4 node 55ok element typeoptionsk3: axisymmetricokmain menupreprocessormaterial propstemperature unitcelsuisokmain menupreprocessormaterial propsmaterial modelsthermalconductivity isotropickxx: 1.2okmain menupreprocessormodelingcreate areasrectangleby dimensionsx1, x2: 0.2, 0.8 y1, y2: 0, 0.03ok划分网格：main menupreprocessormeshingmeshingtoolsize controls: lines setspick 1,3okndiv: 10ok back to “meshtool windows” size controls: lines setspick 2,4okndiv: 2ok back to “meshtool windows”shape: quad freemeshpick allokclose “meshtool windows”热分析：main menusolutionanalysis typenew analysissteady stateokutility menuselectentitieslines, by num/pick, from fullapplypick 4okutility menuselectentitiesnodes, attached to, lines,allokmain menusolutiondefine loadsapplythermaltemperatureon nodespick alloklab 2: temp, value: 90okutility menuselectentitieslines, by num/pick, from fullapplypick 2okutility menuselectentitiesnodes, attached to, lines,allokmain menusolutiondefine loadsapplythermaltemperatureon nodespick alloklab 2: temp, value: 20okutility menuselecteverythingmain menusolutionload step optsoutput ctrlsdb/results filebasic quantities, last substep, value of n: 1okmain menusolutionsolvecurrent lsokclosemain menugeneral postprocplot resultscontour plotnodal solunodal solutiondof solutionnodal temperaturedeformed shape onlyok （或根据讨论的需要进行显示，对显示结果进行存储：utility menuplot ctrlscapture imagefilesave as）热应力耦合分析：main menupreprocessorelement typeswitch elem typethermal to strucokmain menupreprocessorelement typeadd/edit/deletadd optionsk3: axisymmetric okmain menupreprocessormaterial propsmaterial modelsstructurallinearelastic isotropicex: 1.2e11, prxy: 0.3 back to “define mareial model behavior”windows: structuralthermal expansionsecant coefficientisotropicalpx: 1.3e-6okutility menuselectentitieslines, by num/pick, from fullapplypick 3okutility menuselectentitiesnodes, attached to, lines,allokmain menupreprocessorcoupling/ceqncouple dofspick alloknset: 7, lab: uy okutility menuselecteverythingmain menupreprocessorcoupling/ceqncouple dofspick 1, 14, 24oknset: 8, lab: uxokmain menupreprocessorcoupling/ceqncouple dofspick 2, 12, 13oknset: 9, lab: uxokmain menusolutionanalysis typenew analysisstaticokmain menusolutiondefine loadsapplystructuraldisplacementon nodespick 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11lab2: uy, value: 0okmain menusolutiondefine loadsapplystructuraltemperaurefrom thermal analysisfname name of results file: thermalstress.rthokmain menusolutionsolvecurrent lsokclosemain menugeneral postprocplot resultscontour plotnodal solunodal solution stressx-component of stressdeformed shape onlyok （或根据讨论的需要进行显示，对显示结果进行存储：utility menuplot ctrlscapture imagefilesave as）教师签字_第三部分 结果与讨论（可加页）一、实验结果分析（包括数据处理、实验现象分析、影响因素讨论、综合分析和结论等）二、小结、建议及体会三、思考题实验四同三的原理相同，只改变了内外径,温度和截面宽度的情况下，得出如图所示的结果图显示了温度由20到80在构件中的分布情况，很直观。小结与体会通过这次实验，我认识ansys软件在热分析中的应用，而且可以直观地看出温度的分布。实验课程名称：材料成型cae综合实验五 实验项目名称环形焊接接头冷却温度场和应力场分析实验成绩实 验 者专业班级组 别同 组 者实验日期2010年12月23日一部分：实验预习报告（包括实验目的、意义，实验基本原理与方法，主要仪器设备及耗材，实验方案与技术路线等）一）实验目的通过对焊接接头温度场和应力场的有限元模拟，学习用ansys对实际工程问题进行数值分析的过程。二）基本原理和方法对焊接接头应力及温度场分析是材料成型cae中较为复杂的问题。它涉及到热与结构耦合等问题。有限元法是一种离散化的数值计算方法。离散后的单元和单元之间只通过节点相联系，所有场变量（位移、应力、温度等）都通过节点进行计算。对于每个单元，选取适当的插值函数，使得在子域内部、子域分界面上以及子域与外界分界面上都满足一定的条件。然后把所有单元的方程都组装起来，就得到整个结构的方程组。求解方程组，就可以得到方程的近似解。用ansys软件进行有限元分析，整个过程可分为：l 前处理：建立几何模型；对几何模型进行离散化处理等。l 加载求解：根据作用力等效原则将每个单元所受的载荷移置到该单元的节点上；根据边界条件修改刚度方程，消除刚体位移；求解整体刚度方程，得到节点位移；根据相应方程求解应力和应变等。l 后处理：利用计算机图形方式，将计算结果以变形网格、等值线、彩色云图、动画等方式进行显示与分析等。三）实验内容某一圆环由环形钢板和铁板焊接而成，焊接材料为铜，如图为其纵截面的1/2。圆盘初始温度为800，将圆环放置于空气中进行冷却，周围空气为30，对流系数为120w/(m2)。求5min后圆环内部的温度场和应力场分布。（材料参数见表）各种材料的属性见下表：材料温度弹性模量 gpa屈服强度gpa切变模量gpa导热系数w/(m)线胀系数1/比热容j/(kg)密度kg/m3泊松比钢302061.4020.649.81.1010-546578400.32001921.3319.844.84001751.1518.339.46001530.9215.634.08001250.6811.229.0铜301030.910.33991.6610-538689300.3200990.859.8389400900.758.9379600790.627.5366800580.455.2352铁301181.0411.880.51.1810-545578700.32001091.0110.263.5400930.918.650.3600750.766.939.4800520.565.129.6第二部分：实验过程记录（可加页）（包括实验原始数据记录，实验现象记录，实验过程发现的问题等）上机实验步骤（gui法，参考）建立文件：utility menu filechange directoryinto yoursokutility menu filejobnameweldok定义单元类型：main menupreprocessorelement typeadd/edit/deleteaddcouple fieldvector quad 13ok element typeoptionsk1: ux uy temp az, k3:axisymmetricokclose设置钢的弹性模量：main menupreprocessormaterial propsmaterial modelsstrucrurallinearelastic isotropic 单击4次add temperaturetemperature: 30,200,400,600,800, ex: 依次钢的弹性模量（单位pa）：2.06e11, 1.92e11, 1.75e11, 1.53e11, 1.25e11，prxy: 依次0.3ok设置钢的线膨胀系数：in ”define material model behavior” windowsstructuralthermal expansionsecand coefficientisotropic alpx: 1.10e-5ok设置钢的密度：in ”define material model behavior” windowsstructuraldensity den: 依次7840ok设置钢的双线性材料属性：in ”define material model behavior” windowsstructuralnonlinearinelasticrate independentisotropic hardening plasticitymises plasticitybilinear单击4次add temperaturetemperature: 30,200,400,600,800, yield stss: 依次钢的屈服强度（单位pa）, tang mod: 依次钢的切变模量（单位pa）ok设置钢的热传导率：in ”define material model behavior” windowsthermalconductivityisotropic单击4次add temperaturetemperature: 30,200,400,600,800, kxx: 依次钢的导热系数ok设置钢的比热容：in ”define material model behavior” windowsthermalspecific heat c: 465ok设置第2种和第3种材料的属性：in ”define material model behavior” windowsnew modeldefined material id: 2ok按设置第1种材料（钢）属性的方法定义第2种材料（铜）的属性。同理：设置第3种材料（铁）的属性。生成有限元模型：main menupreprocessormodelingcreate keypoints in active cs输入关键点：1（1, 0, 0）；2（1.8, 0, 0）；3（2.6, 0, 0）；4（1, 0.2, 0）；5（1.7, 0.2, 0）；6（1.9, 0.2, 0）；7（2.6, 0.2, 0），每次输入一组数据后必须单击apply按钮确认。main menupreprocessormodelingcreateareaarbitrarythrough kpspick 1, 2, 5, 4apply2, 3, 7, 6okmain menupreprocessormodelingcreatelinesarcby end kps & radpick 5, 6ok2okrad radius of the arc: 0.2okmain menupreprocessormodelingcreateareaarbitraryby lines2, 8, 9okutility menuplotarea （按save_db按钮保存或save as）main menupreprocessormeshingsizecntrlmanuasizelinespicked ines1, 3, 5, 7ok ndiv: 16okmain menupreprocessormeshingmesh attributesdefault attibsmat: 1okmain menupreprocessormeshingmeshareafree1okmain menupreprocessormeshingmesh attributesdefault attibsmat: 2okmain menupreprocesso