| Open the reference model | |Keep in mind to modify the current directory |*open_model c:\users\manto\documents\201819_esami\paom\20190509_paom_lab\piede_cdm2_geom2d3d_set.mfd |*reset_view |*fill_view |----------------------------------- | PLOT curves |----------------------------------- *set_curve_tolerance 0.01 |----------------------------------- |----------------------------------- | START MESH CON-ROD |----------------------------------- |----------------------------------- | AUTOMESH: Curve divisions |----------------------------------- *set_curve_div_type_fix *set_curve_div_type_fix_avgl *set_curve_div_avgl 1 *set_curve_div_rest_evn *set_curve_div_applyrest_cvs *apply_curve_divisions conrod all_selected |----------------------------------- | AUTOMESH: Curve divisions | REFINEMENT |----------------------------------- |----------------------------------- | Outer circumferential edge at small end |----------------------------------- *select_clear_curves *set_curve_div_type_fix_ndiv *set_curve_div_num 43 *set_curve_div_rest_off *apply_curve_divisions 21854 # | End of List |----------------------------------- | Middle circumferential edges at small end |----------------------------------- *select_clear_curves *set_curve_div_type_fix_ndiv *set_curve_div_num 32 *set_curve_div_rest_off *apply_curve_divisions 21862 21818 # | End of List |----------------------------------- | Vertical edge at small end |----------------------------------- *select_clear_curves *set_curve_div_type_fix_ndiv *set_curve_div_num 3 *set_curve_div_rest_off *apply_curve_divisions 21866 # | End of List *select_clear_curves |----------------------------------- | AUTOMESH: Curve divisions | planar mesh Quadrilateral 4-nodes elements | conrod shank |----------------------------------- *af_planar_quadmesh 21854 12618 21849 21871 21837 12712 21863 21869 21818 21862 21866 # | End of List *select_clear_curves |----------------------------------- | STORE and PLOT the elements in a collector | named shank_elems |----------------------------------- *select_clear_elements *select_elements all_existing *store_elements shank_elems all_selected *identify_sets *regen *identify_none *regen *select_clear_elements *invisible_set shank_elems |----------------------------------- | AUTOMESH: Curve divisions | planar mesh Quadrilateral 4-nodes elements | conrod pocket |----------------------------------- *af_planar_quadmesh 12712 21834 21837 21863 21870 # | End of List *select_clear_curves |----------------------------------- | STORE and PLOT the elements in a collector | named pocket_elems |----------------------------------- *select_clear_elements *select_elements all_visible *store_elements pocket_elems all_selected *select_clear_elements *identify_sets *regen *visible_all_sets *identify_none *regen |----------------------------------- | FROM QUAD4 TO QUAD8 | CHANGE ELEMENTS CLASS |----------------------------------- *change_elements_quadratic *change_elements_class *set_change_class quad8 *change_elements_class all_existing |----------------------------------- | CREATE NODES: | - AT THE SMALL END INNER RADIUS(ID 2459) | - AT MIDDLE SMALL END (ID 2460) | - AT OUTER SMALL END (ID 48) |----------------------------------- *add_nodes 0 23/2 0 *add_nodes 0 25.003/2 0 |----------------------------------- | CREATE NODE AT THE SMALL END INNER RADIUS | create a 1D quadratic element (line3) | 1 node: inner node (ID 2459) | 2 node: outer node (ID 48) | 3 node: middle node (ID 2460) |----------------------------------- *set_element_class line3 *add_elements 2459 48 2460 |----------------------------------- | select the element to be subdivided by class (line3) |----------------------------------- *select_clear_elements *select_elements *select_elements_class line3 |----------------------------------- | MESH GENERATION: subdivide | The line3 element will be subdivided in three. |----------------------------------- *subdivide_reset *sub_divisions 3 1 1 *subdivide_elements all_selected |----------------------------------- | select the elements previously subdivided by class (line3) |----------------------------------- *select_clear_elements *select_elements *select_elements_class line3 |----------------------------------- | MESH GENERATION: expand | From 1D quadratic elements (line3) | to 2D-planar quadratic elements (quad8). |----------------------------------- *expand_reset *set_expand_repetitions 64 *set_expand_rotation z -180/64 *expand_elements all_selected |----------------------------------- | STORE and PLOT the elements in a collector | named small_end_elems |----------------------------------- *visible_all_sets *identify_sets *regen *invisible_set shank_elems *invisible_set pocket_elems *select_clear_elements *select_elements all_visible *store_elements small_end_elems all_selected *select_clear_elements *visible_all_sets |----------------------------------- | MESH GENERATION: sweep nodes |----------------------------------- *set_sweep_tolerance 0.05 *sweep_nodes all_existing |----------------------------------- | FINISH MESH CON-ROD |----------------------------------- |----------------------------------- |----------------------------------- | START MESH BUSH |----------------------------------- | outer node ID: 3251 | inner node ID: 3252 | middle node ID: 3253 |----------------------------------- *add_nodes -50 23/2 0 -50 21/2 0 -50 22/2 0 |----------------------------------- | create a 1D quadratic element (line3) | 1 node: outer node (ID 3251) | 2 node: inner node (ID 3252) | 3 node: middle node (ID 3253) |----------------------------------- *set_element_class line3 *add_elements 3251 3252 3253 |----------------------------------- | select the element to be subdivided by class (line3) |----------------------------------- *select_clear_elements *select_elements *select_elements_class line3 |----------------------------------- | MESH GENERATION: subdivide | The line3 element will be subdivided in two. |----------------------------------- *subdivide_reset *sub_divisions 2 1 1 *subdivide_elements all_selected |----------------------------------- | select the elements previously subdivided by class (line3) |----------------------------------- *select_clear_elements *select_elements *select_elements_class line3 |----------------------------------- | MESH GENERATION: expand | From 1D quadratic elements (line3) | to 2D-planar quadratic elements (quad8). |----------------------------------- *expand_reset *set_expand_point -50 0 0 *set_expand_repetitions 64 *set_expand_rotation z -180/64 *expand_elements all_selected |----------------------------------- | STORE and PLOT the elements in a collector | named bush_elems |----------------------------------- *visible_all_sets *identify_sets *regen *invisible_set shank_elems *invisible_set pocket_elems *invisible_set small_end_elems *select_clear_elements *select_elements all_visible *store_elements bush_elems all_selected *select_clear_elements *visible_all_sets |----------------------------------- | FINISH MESH BUSH |----------------------------------- |----------------------------------- | START MESH GUDGEON PIN |----------------------------------- |----------------------------------- | Outer gudgeon pin curve |----------------------------------- *select_clear_curves *set_curve_div_type_fix_ndiv *set_curve_div_num 72 *set_curve_div_rest_off *apply_curve_divisions 21879 # | End of List *select_clear_curves |----------------------------------- | Inner and vertical gudgeon pin curves |----------------------------------- *set_curve_div_type_fix *set_curve_div_type_fix_avgl *set_curve_div_avgl 0.5 *set_curve_div_rest_evn *set_curve_div_applyrest_cvs *apply_curve_divisions 21876 21877 21878 all_selected |----------------------------------- | AUTOMESH: Curve divisions | planar mesh Quadrilateral 4-nodes elements | gudgeon pin |----------------------------------- *af_planar_quadmesh gudgeon_pin *select_clear_curves |----------------------------------- | STORE and PLOT the elements in a collector | named gudgeon_pin_elems |----------------------------------- *invisible_set bush_elems *invisible_set pocket_elems *invisible_set shank_elems *invisible_set small_end_elems *select_elements all_visible *store_elements gudgeon_pin_elems all_selected *select_clear_elements *visible_all_sets |----------------------------------- | FROM QUAD4 TO QUAD8 | CHANGE ELEMENTS CLASS |----------------------------------- *change_elements_class *change_elements_quadratic *set_change_class quad8 *change_elements_class all_visible |----------------------------------- | FINISH MESH GUDGEON PIN |----------------------------------- |----------------------------------- | FINISH MESH |----------------------------------- |----------------------------------- | START GEOMETRIC PROPERTIES |----------------------------------- *identify_geometries *regen |----------------------------------- | SHANK ELEMS GEOMETRIC PROPERTIES |----------------------------------- *new_geometry *geometry_type mech_planar_pstress *geometry_name *geometry_name shank_elems_18mm *geometry_param norm_to_plane_thick 18 *add_geometry_elements shank_elems |----------------------------------- | POCKET ELEMS GEOMETRIC PROPERTIES |----------------------------------- *new_geometry *geometry_type mech_planar_pstress *geometry_name pocket_elems_4mm *geometry_param norm_to_plane_thick 4 *add_geometry_elements pocket_elems |----------------------------------- | SMALL END ELEMS GEOMETRIC PROPERTIES |----------------------------------- *new_geometry *geometry_type mech_planar_pstress *geometry_name small_end_elems_20mm *geometry_param norm_to_plane_thick 20 *add_geometry_elements small_end_elems |----------------------------------- | BUSH ELEMS GEOMETRIC PROPERTIES |----------------------------------- *new_geometry *geometry_type mech_planar_pstress *geometry_name bush *geometry_name bush_elems_20mm *geometry_param norm_to_plane_thick 20 *add_geometry_elements bush_elems |----------------------------------- | GUDGEON PIN ELEMS GEOMETRIC PROPERTIES |----------------------------------- *new_geometry *geometry_type mech_planar_pstrain *geometry_name gudgeon_pin_elems_60mm *geometry_param norm_to_plane_thick 60 *add_geometry_elements gudgeon_pin_elems *identify_none *regen |----------------------------------- | FINISH GEOMETRIC PROPERTIES |----------------------------------- |----------------------------------- | START MATERIALS |----------------------------------- *new_mater standard *mater_option general:state:solid *mater_name steel *mater_option structural:type:elast_plast_iso *mater_param structural:youngs_modulus 210000 *mater_param structural:poissons_ratio 0.3 *add_mater_elements all_existing *identify_materials *regen *identify_none *regen |----------------------------------- | FINISH MATERIALS |----------------------------------- |----------------------------------- | COMPONENTS POSITIONING START |----------------------------------- | PLOT THE BUSH AND GUDGEON PIN | ELEMENTS AND CURVES |----------------------------------- *invisible_all_sets *visible_set bush *visible_set bush_elems *visible_set gudgeon_pin *visible_set gudgeon_pin_elems |----------------------------------- | MESH GENERATION: MOVE (STEP 1) | BOTH BUSH AND GUDGEON PIN |----------------------------------- *move_reset *set_move_translation x 50 *set_move_translation z 50 *move_combined all_visible *invisible_set bush *invisible_set bush_elems |----------------------------------- | MESH GENERATION: MOVE (STEP 2) | GUDGEON PIN ONLY |----------------------------------- *move_combined all_visible *visible_all_sets |----------------------------------- | MESH GENERATION: SWEEP |----------------------------------- *set_sweep_tolerance 0.0001 *sweep_all |----------------------------------- | Command prompt | Deleting 8 duplicate nodes! | Deleting 0 collapsed elements! | Deleting 0 duplicate elements! | Deleting 228 duplicate points! |----------------------------------- |----------------------------------- | COMPONENTS POSITIONING FINISH |----------------------------------- |----------------------------------- | CONTACT START |----------------------------------- |----------------------------------- | CONTACT BODY | BODY 1--> Small_end_cb : Deformable Discrete |----------------------------------- *new_cbody mesh *contact_option state:solid *contact_option skip_structural:off *contact_body_name small_end_cb *contact_option defo_desc:discrete *add_contact_body_elements small_end_elems |----------------------------------- | CONTACT BODY | BODY 2--> Bush_cb : Deformable Analytical |----------------------------------- *new_cbody mesh *contact_option state:solid *contact_option skip_structural:off *contact_body_name bush_cb *contact_option defo_desc:analytical *add_contact_body_elements bush_elems |----------------------------------- | CONTACT BODY | BODY 3--> Gudgeon_pin_cb: Deformable Discrete |----------------------------------- *new_cbody mesh *contact_option state:solid *contact_option skip_structural:off *contact_body_name gudgeon_pin_cb *contact_option defo_desc:discrete *add_contact_body_elements gudgeon_pin_elems |----------------------------------- *identify_contact *regen *identify_none *regen |----------------------------------- | EXCLUDE DISCONTINUITY NODES FROM THE | ANALYTICAL CONTACT BODY e.g. bush |----------------------------------- *edit_contact_body bush_cb *contact_option discontinuity:manual *add_cbody_dc_nodes 3251 3252 3903 3387 # | End of List |----------------------------------- | CONTACT INTERACTION | MESHED (DEFORMABLE) vs MESHED (DEFORMABLE) |----------------------------------- *new_interact mesh:mesh *interact_option state_1:solid *interact_option state_2:solid *interact_name unilateral_frictionless *interact_option dist_tol:redefined *interact_param dist_tol 0.001 |----------------------------------- | CONTACT TABLE |----------------------------------- *new_contact_table *contact_table_name ctable_interference_fit |----------------------------------- | contact between the bush and the conrod |----------------------------------- *ctable_entry small_end_cb bush_cb *contact_table_option small_end_cb bush_cb contact:on *prog_string ctable:old_interact unilateral_frictionless *ctable_entry_interact small_end_cb bush_cb *contact_table_option small_end_cb bush_cb detection:first_second *contact_table_option small_end_cb bush_cb intfer_fit:on *contact_table_property small_end_cb bush_cb closure 0.03 |----------------------------------- | Table to modulate the interference_fit |----------------------------------- *new_md_table 1 1 *table_name table1_interference_fit *set_md_table_type 1 time *set_md_table_max_v 1 2 *set_md_table_method_data_points *table_add 0 0 1 1 2 1 *ctable_entry small_end_cb bush_cb *contact_table_property_table small_end_cb bush_cb closure table1_interference_fit |----------------------------------- | contact between the bush and the gudgeon_pin |----------------------------------- *ctable_entry gudgeon_pin_cb bush_cb *contact_table_option bush_cb gudgeon_pin_cb contact:on *edit_interact unilateral_frictionless *prog_string ctable:old_interact unilateral_frictionless *ctable_entry_interact bush_cb gudgeon_pin_cb *contact_table_option bush_cb gudgeon_pin_cb detection:second_first *ctable_entry gudgeon_pin_cb bush_cb *contact_table_option bush_cb gudgeon_pin_cb contact:off |----------------------------------- | LINK: SPRING | Node ID 3252: bush | Node ID 3988: gudgeon_pin | Used to delete the rigid motion of this component | that is not involved at the interference fit stage |----------------------------------- *new_link *link_class spring *spring_option spring_type:fixed_dof *link_class spring *spring_dof 0 2 *link_class spring *spring_dof 1 2 *link_class spring *spring_param stiffness 1 *link_class spring *spring_node 0 3252 *link_class spring *spring_node 1 3988 |----------------------------------- | BCs START |----------------------------------- |----------------------------------- | Select the nodes at the symmetry plane YZ | at X=0 |----------------------------------- *select_clear *select_filter_none *select_method_box *select_mode_and *select_nodes -1000 0.01 -1000 1000 -1000 1000 |----------------------------------- | BCs symm_yz_nx |----------------------------------- *new_apply *apply_type fixed_displacement *apply_name symmyz_nx *apply_dof x *apply_dof_value x *add_apply_nodes all_selected *select_clear |----------------------------------- | SHANK FIXED AT THE FREE EXTREMITY |----------------------------------- |----------------------------------- | Select the nodes at the shank free edge | at Y=0 |----------------------------------- *select_clear *select_filter_none *select_method_box *select_mode_and *select_nodes -1000 1000 -1000 -51.49 -1000 1000 |----------------------------------- | BCs farfield |----------------------------------- *new_apply *apply_type fixed_displacement *apply_name farfield *apply_dof y *apply_dof_value y *add_apply_nodes all_selected *select_clear |----------------------------------- | START LOADCASE |----------------------------------- *new_loadcase *loadcase_type struc:static *loadcase_name lcase1_interference *loadcase_ctable ctable_interference_fit *loadcase_option converge:displacements *loadcase_value displacement 0.01 *loadcase_value nsteps 4 |----------------------------------- | FINISH LOADCASE |----------------------------------- |----------------------------------- | START JOBS |----------------------------------- *new_job *job_class structural *job_name job1_interference *add_job_loadcases lcase1_interference *job_contact_table ctable_interference_fit *add_post_tensor stress *add_post_var von_mises *job_option nod_quantities:manual *add_post_nodal_quantity Displacement *add_post_nodal_quantity Ext_Force *add_post_nodal_quantity Rea_Force *add_post_nodal_quantity Cont_Nor_Strss *add_post_nodal_quantity Cont_Nor_Force *add_post_nodal_quantity Cont_Status *add_post_nodal_quantity Near_Cont_Dist |----------------------------------- | FINISH JOBS |----------------------------------- |----------------------------------- | CHECK AND SAVE MODEL |----------------------------------- *save_model *check_job *update_job |----------------------------------- | Found 1 errors and 1 warnings! |----------------------------------- |----------------------------------- | MESH GENERATION: CHECK AND FLIP ELEMENTS |----------------------------------- *check_upside_down *flip_elements all_selected *check_upside_down |----------------------------------- | CHECK AND SAVE MODEL | Found 0 errors and 0 warnings! |----------------------------------- *save_model *check_job *update_job *submit_job 1 *monitor_job |----------------------------------- | OPEN THE RESULTS |-----------------------------------