MARC实例.doc

下面将圆柱面的一些节点附着到二个曲面的相交线VIEW MSHOW VIEW3 M M显示视图VIEW3ATTACH M ATTACH-NODES-INTERSECT Mcylinder sphere 2 11 12 13 32 3334 G M下图中圆柱左端节点即为要附着到二个曲面相交线上的节点，左边的节点编号供参考。左图为ATTACH操作后的结果，有一部分单元变成畸形，需要修正。下面将利用STRETCH功能对一些单元进行调整。PLOT M MORE MINVISIBLE-SURFACES M ALL-EXIST MREGENERATE M M M使几何曲面从显示区消失。STRETCH M NODES M13 22 10 G12 27 9 G11 26 8 G 2 7 1 G32 48 30 G34 49 46 43 31 G进入STRETCH子菜单。输入一系列要修正的节点。第1步 对第1列节点进行拉伸第2步 对第2列节点进行拉伸第3步 对第3列节点进行拉伸第4步 对第4列节点进行拉伸 第5步 对第5列节点进行拉伸第6步 对第6列节点进行拉伸经过一连串的拉伸调整，修正后的圆柱面网格如左图所示。下面将把半球面与圆柱面附着一起。MOVE M RESET M TRANSLATIONS M0 3 0 ELEMENTS M ALL-UNSELECTED M M先将圆柱面单元沿Y轴移动3。PLOT M MORE M VISIBLE-ELEMENTS MALL-SELECTED M M REGENERATE M RETURN M RETURN M将目前INVISIBLE的半球面单元变为VISIBLE。NODES-EDIT M84 0 71.0469 75.1 83 14.7362 72.875 70.1951 82 22.8594 78.3594 60.1217 79 24.9911 84.25 50.6664 74 23.3348 90.1406 41.0287 75 16.5643 95.8281 31.2592 76 -.890113 99.8906 25.0159 修改节点位置将84点移至13点处，使二点重合。一般用鼠标器操作，修改完成后，对话区将显示新位置的坐标。其它节点的修改与84点类似。SWEEP M ALL M MMOVE M TRANSLATIONS M0 3 0 ELEMENTS M ALL-UNSELECTED M M用SWEEP命令消除由EDIT引起的节点重复编号，将圆柱面移回原来位置。RENUMBER M ALL M M用RENUMBER对节点号单元号等重排，得到分析用的网格。左图为生成的网格图边界条件的施加网格生成后，要给模型必要的边界条件才能进行分析，本题主要是施加对称边界条件及内压。VIEW M SHOW ALL VIEWS M M SELECT MCLEAR SELECT M M显示所有的视图，消除对半球面单元的选择。（集名仍留在数据库中）PLOT M ELEMENTS-SOLID M ELEMENTS SIZE (%)90 M REGENERATE M FILL M ELEMENT SIZE (%)100 M ELEMENTS-WIREFRAME MDRAW-FACES M 以实体形式显示单元，并将单元收缩为原来的90，检查网格是否有错，确认无错后，恢复单元大小，重新用线框形式显示单元。 MORE MEDGE-OUTLINE M REGENERATE M以单元的边来显示模型的轮廓，以便后面所述的边界条件的施加。先施加左右对称引起的位移边界条件为便于对称位移边界条件的施加，对一部分圆弧上的节点使用局部坐标系。BOUNDARY CONDITIONS MMECHANICAL M进入边界条件子菜单，边界条件类型为MECHANICAL。 NAME M yz1 给边界条件取名。 FIXED DISPLACEMENT M边界条件为指定位移。 X-DISPLACEMENT M OK M约束X方向位移。 NODES-ADD M 10 13 16 19 22 25 28 31 34 37 40 43 46 49 61 66 71 82 87 92 97 102 107 112 117 G指定施加的节点编号，最后用G确认。 NEW M新的边界条件定义指示。 NAME M yz2 新的边界条件名。 TRANSFORMS M CYLINDRICAL M 0 0 0 1 0 0对于圆弧上的节点采用圆柱坐标系。由点（1，0，0）与（1，0，0）连线为圆柱的轴。 13 34 66 71 82 92 97 102 107 G FIXED DISPLACEMENT M采用此局部坐标系的节点指定。 Y-ROTATION M OK M绕局部坐标系Y轴的转动自由度被约束。 NODES-ADD M 13 34 66 71 82 92 97 102 107 G将约束条件施加于采用局部坐标系的节 点上。 NEW M NAME M yz3 新的边界条件。FIXED DISPLACEMENT M Z-ROTATION M OK M约束在整体坐标系绕Z轴转动自由度。 NODES-ADD M 61 87 G施加本约束的节点为圆弧二个端点。以下为上面三个边界条件施加的结果。接着将施加半球面下侧节点（xy平面上的节点）位移约束条件。NEW MNAME M xy1将定义新的边界条件，名为xy1。TRANSFORMS M CYLINDRICAL M 0 0 0 1 0 0 采用圆柱坐标系为局部坐标系。50 57 58 59 60 8384 85 86 G MFIXED DISPLACEMENT M指定除二个端点以外的在xy平面上的节点采用局部坐标系。 X-DISPLACEMENT M Y-DISPLACEMENT M Z-DISPLACEMENT M Y-ROTATION M OK M约束径向。周向、高度方向的位移，以及沿周向的转动。NODES-ADD M50 57 58 59 60 83 84 85 86 G指定施加上述约束的节点，这些节点也就是使用局部坐标系的节点。NEW MNAME M xy2将定义新的边界条件，名为xy2。FIXED DISPLACEMENT M Y-DISPLACEMENT M Z-DISPLACEMENT M X-ROTATION M OK M约束在整体坐标系下y、z方向的位移以及绕x轴的转动。NODES-ADD M61 87 G M将约束条件施加到二个端点上。ID BOUNDARY CONDS M将所有边界条件显示出来，确认。ID BOUNDARY CONDS M最后施加内压确认结束。MECHANICAL M将定义新的边界条件。 NEW M FACELOAD M PRESSURE M1 OK M FACES-ADD M ALL-EXIST M在所有的单元面上施加值为1的分布 载荷。PLOT M MORE M EDGES-FULL M显示所有单元的边。REGENERATE M M M显示结果见前页图所示。至此全部边界条件的定义、施加即告结束。材料特性的定义由于是简单的线弹性问题，只要定义杨氏模量及波松比即可。MAIN M MATERIAL PROPERTIES M ISOTROPIC M YOUNGS MODULUS M 21000 YOISSONS RATIO M 0.3 OK M ELEMENTS-ADD M ALL-EXIST M在MAIN菜单里，检取MATERIAL PROPERTIES即进入了材料特性定义子菜单，再选取ISOTROPIC，输入杨氏模量、波松比，然后指定到所有单元上。几何特性的定义对于本题为定义单元厚度MAIN MGEOMETRIC PROPERTIES M 3-D MSHELL MTHICKNESS M1.0 OK MELEMENTS-ADD M ALL-EXIST M进入GEOMETRY PROPERTIES子菜单对三维单元中的壳单元定义厚度，并指定到所有单元上。选取作业参数并提交MARC运算分析在本步中将完成作业运行所需的全部数据输入，并在MENTAT上提交MARC程序运行分析，由于是线弹性问题，因而载荷工况可以不定义。MAIN MJOBS M MECHANICAL MJOB RESULTS MSELECT TENSOR M STRESS MSELECT VARIABLES M VON-MISES MOK M进行JOBS子菜单，选择应力分析为本作业的分析类型，在JOB RESULTS子菜单中选择后处理要显示的内容，此题选择应力张量和Von Mises等效应力。JOB PARAMETERS M#SHELL/BEAMLAYERS M 3 OK MRETURN M此JOB PARAMETER子菜单中，选择SHELL/BEAMLAYER，用键盘输入沿壳单元厚度方向的积分点数，对于线弹性问题用3层已足够准确。 ELEMENT TYPES M3-D SHELLS M 75 OK MALL-EXIST M MRETURN M最后的单元类型选择。在各单元种类名中选取3D SHELL后图形区会弹出一个菜单，在菜单中选取75（厚壳单元）并指定到所有单元。 CHECK M最后检查模型。FILES M SAVE AS Mdome_tube yes M RUN MSUBMIT 1 MMONITOR M在提交MARC运行作业以前，别忘了存贮数据库，以免运行MARC程序后进行的后处理操作将数据库冲掉。最后提交作业给MARC，并对运行过程进行监视。后处理分析结束后，对分析结果进行后处理OK M MAIN MRESULTS M OPEN DEFAULT MDEFORMED SHAPE-DEF ONLY MPLOT STYLE-CONTOUR BANDS MPLOT QUANTITY MEQUIVALENT VON MISES STRESS MNEXT INC M在RUN菜单中的EXIT NUMBER右边出现3004表示分析正常结束，检取OK光钮。在MAIN菜单中检取RESULTS取入后处理子菜单。OPEN DEFAULT将刚在MENTAT上提交完成作业的后处理文件打开，DEFONLY表示仅显示变形图。显示方式为云图。显示内容为Von Mises等效应力。NEXT INC使文件进到增量步0，模型的变形图及Von Mises等效应力云图同时显示出来。CLOSE M显示结束，关闭后处理文件。左图为变形及Von Mises等效应力云图整个分析结束，本例题可供用户进行线弹性分析参考。在MENTAT中建模、递交分析和结果后处理的菜单命令记录在c文件中，内容如下：| Version : Mentat 3.2*new_modelyes*reset*system_spherical*add_nodes100 0 0100 0 75*fill_view*show_view 4*fill_view*set_element_class line2*add_elements 1 2*sub_divisions6 1 1*subdivide_elementsall_existing*sweep_all*set_expand_rotations0 0 18*set_expand_repetitions5*expand_elementsall_existing*fill_view*set_symmetry_normal0 -1 0*symmetry_elementsall_existing*sweep_nodes*sweep_all*renumber_all*remove_elements 19 20 15 14# | End of List*remove_nodes 31 32# | End of List*sweep_all*renumber_all*fill_view*store_elementselem_domeall_existing*save_as_model dome.mud yes*new_modelyes*system_cylindrical*add_nodes25 0 025 0 50*add_elements 1 2*sub_divisions6 1 1 *subdivide_elementsall_existing*sweep_all*set_expand_rotations0 0 30*set_expand_repetitions3*expand_elementsall_existing*set_symmetry_normal0 -1 0*symmetry_elementsall_existing*sweep_all*renumber_all*store_elementselem_tubeall_existing*fill_view*save_as_model tube.mud yes*system_rectangular*set_move_rotations90 0 0*set_move_translations0 130 50*move_elementsall_existing*fill_view*merge_models dome.mud*fill_view*set_surface_type sphere*add_surfaces0 0 0100*set_surface_type cylinder*add_surfaces0 0 500 130 502525*set_plot_divisions_med*regenerate*select_setselem_dome*invisible_elementsall_unselected*regenerate*attach_nodes_intersect 2 1 84 83 82 79 74 75 76# | End of List*visible_elements*visible_elementsall_existing*invisible_elementsall_selected*regenerate*show_view 3*fill_view*attach_nodes_intersect 2 1 2 11 12 13 32 33 34# | End of List*stretch_nodes132210# | End of List12279# | End of List11268# | End of List271# | End of List*stretch_nodes 32 47 41 29# | End of List 33 48 30# | End of List3449464331# | End of List*move_reset*set_move_translations0 3 0*move_elementsall_unselected*visible_elementsall_selected*regenerate*edit_nodes 71*edit_nodes 76 -8.901210221221e-001 9.989062500000e+001 2.501585133398e+001 75 1.654629050347e+001 9.582812500000e+001 3.125912834005e+001 74 2.333484221523e+001 9.014062500000e+001 4.102864899859e+001 79 2.499111518087e+001 8.425000000000e+001 5.066645481197e+001 82 2.285936809288e+001 7.835937500000e+001 6.012172368690e+001 83 1.473627923807e+001 7.287500000000e+001 7.019510025273e+001 84 0.000000000000e+000 7.104687500000e+001 7.500000000000e+001*sweep_all*renumber_all*fill_view*set_move_translations0 -3 0*move_elementsall_unselected*renumber_all*select_clear*show_view 4*reset_view*fill_view*elements_solid*regenerate*set_element_size 80*regenerate*set_element_size 100*regenerate*elements_wireframe*set_faces off*edges_outline*regenerate*set_surfaces off*regenerate*set_points off*regenerate*show_view 1*fill_view*apply_type fixed_displacement*apply_dof x*add_apply_nodes 10 13 16 19 22 25 28 31 34 37 40 43 46 49 61 66 71 82 87 92 97 102 107 112 117# | End of List*apply_nameyz1*new_apply*apply_nameyz2*apply_type fixed_displacement*transform_cylindrical0 0 01 0 0*show_view 4*transform_cylindrical0 0 01 0 0 97 102 107 112 117 92 66 82 13 34 71# | End of List*fill_view*apply_type fixed_displacement*apply_dof ry*add_apply_nodes 97 102 107 112 117 92 66 82 13 34 71# | End of List*fill_view*new_apply*apply_nameyz3*apply_type fixed_displacement*apply_value rz*apply_dof rz*add_apply_nodes 61 87# | End of List*new_apply*apply_namexy1*transform_cylindrical0 0 01 0 0 60 59 58 57 50 83 84 85 86# | End of List*apply_type fixed_displacement*apply_value x*apply_dof x*apply_dof y*apply_dof z*apply_dof ry*add_apply_nodes 60 59 58 57 50 83 85 86# | End of List*new_apply*apply_namexy2*apply_type fixed_dis