变形体接触问题实例变形体接触问题实例| Version : Mentat 3.3*alias action*procedure_pause 0|*exec_p gen_image.proctop()*action | figure 7.0 - dummy*action | figure 7.1 - dummy*action | figure 7.2 - dummy| step 1|push(meshgen)push(coord_system)*grid_u_domain0 8*grid_u_spacing1*grid_v_domain0 5*grid_v_spacing0.5*set_grid on*fill_viewpop(1)push(plot_screen)*set_point_labels onpop(1)| create a surface corresponding to the tube| length=8.0, inner diam=4.0, outer diam=4.3|*add_surfacespoint( 0.000000000000e+00, 4.000000000000e+00, 0.000000000000e+00)point( 8.000000000000e+00, 4.000000000000e+00, 0.000000000000e+00)point( 8.000000000000e+00, 4.500000000000e+00, 0.000000000000e+00)point( 0.000000000000e+00, 4.500000000000e+00, 0.000000000000e+00)push(move)*set_move_translations0 -0.2 0*move_points3 4# | End of List*action | figure 7.3pop(1)push(coord_system)*system_cylindricalpop(1)*add_points5 15 016 15 0*set_grid off*fill_viewpush(coord_system)*system_rectangularpop(1)push(duplicate)*set_duplicate_translations0 .6 0*duplicate_points56# | End of Listpop(1)*add_surfaces5687*action | figure 7.4push(move)*set_move_translations0 1.25 0*move_surfaces2# | End of List*action | figure 7.5pop(1)| mesh the tube 8x3 elements| mesh the tool 14x6 elements|push(convert)*set_convert_divisions8 3*convert_surfaces1# | End of List*set_convert_divisions14 6*convert_surfaces2# | End of Listpush(plot_screen)*set_points off*set_point_labels off*set_surfaces off*regenerate*action | figure 7.6| select the marc element type|pop(1)pop(1)pop(1)push(jobs)push(element_types)popup(mech_eltype_axisym_solid_pm)*action | figure 7.7*element_type 10popdown(mech_eltype_axisym_solid_pm)all_existingpop(1)pop(1)| fix the far let side of the tube in axial direction|push(apply)push(apply_mechanical)*apply_type fixed_displacementset($usersubfv,FORCDT)popup(fixed_disp_pm)*apply_dof xpopdown(fixed_disp_pm)*add_apply_nodes1 10 19 28# | End of List*action | figure 7.8pop(1)pop(1)| create an artificial connection between the tube| and the tool by means of a spring|push(links)set($class,spring)set($notclass1,tie)set($notclass2,servo)push(springs_dashpots)*zoom_box*zoom_box(1,0.439909,0.375318,0.540816,0.496183)*link_class spring *link_value stiffness10e3*link_class spring *spring_node 0131*link_class spring *spring_dof 01*link_class spring *spring_node 19*link_class spring *spring_dof 11*action | figure 7.9pop(1)pop(1)| step 2| | define and assign material properties to both tube and tool|*fill_viewpush(mat_props)*material_type mechanical:isotropicpopup(mat_isotropic)*material_value isotropic:youngs_modulus30e6*material_value isotropic:poissons_ratio.3*material_type plasticitypopup(iso_plastic_popmenu)*material_value plasticity:yield_stress3.6e4popdown(iso_plastic_popmenu)popdown(mat_isotropic)*add_material_elements1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24# | End of List*identify_materials *regen*action | figure 7.10*new_material*material_type mechanical:isotropicpopup(mat_isotropic)*material_value isotropic:youngs_modulus40e6*material_value isotropic:poissons_ratio.3*material_type plasticitypopup(iso_plastic_popmenu)*material_value plasticity:yield_stress6.0e4popdown(iso_plastic_popmenu)popdown(mat_isotropic)*add_material_elements25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 5051 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 7677 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101102 103 104 105 106 107 108# | End of List*action | figure 7.11pop(1)| step 3| assign the elements of tube and tool to contact bodies|push(contact)push(contact_bodies)*contact_body_nametube*contact_deformablepopup(contact_deform_pm)popdown(contact_deform_pm)*add_contact_body_elements1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24# | End of List*new_contact_body*contact_body_nametool*contact_deformablepopup(contact_deform_pm)popdown(contact_deform_pm)*add_contact_body_elements25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 5051 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 7677 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101102 103 104 105 106 107 108# | End of List*identify_contact *regen*action | figure 7.12*identify_none *regenpop(1)pop(1)| step 4 | define the time dependent loading of the flaring tool| and apply this to the element edges at the far right| of the flaring tool|push(apply)push(apply_mechanical)push(table)*set_table_typepopup(table_type_pm)list(table_type_list,0)timepopdown(table_type_pm)*set_table_xmax87*set_table_xstep87*set_table_ymax2400*table_add0 09 90039 240087 0*table_filled*action | figure 7.13*show_modelpop(1)*new_apply*apply_type edge_loadset($usersubfv,FORCEM)popup(apply_edge_load_popmenu)*apply_value p*apply_table p0set($dof,p0)popup(apply_select_table_popmenu)list(select_table_list,0)table1po