|	Remember to set the directory 
*change_directory C:\Users\manto\documents\201819_esami\paom\20190523_paom_lab\upload
*save_as_model C:\Users\manto\documents\201819_esami\paom\20190523_paom_lab\upload\mesh_geom_mat_link.mud yes
|----------------------------------
| Evaluation of the modal responce of a coil with a hollow circular cross-section
| Its geometry is:
| Mean radius: Rm = 20 mm
| Total coils: n = 4,5
| Coil pitch: p = 15 mm
| Wire outer diameter: de = 12 mm 
| Wire inner diameter: di = 6 mm
|
| Material: Titanium: 
| E = 110 GPa, nu = 0.3,  
| rho = 4.7 e-9 tonn/mm^3
|----------------------------------
| 	MESH GENERATION START
|----------------------------------
| Define the points located at: 
| - the system origin O (0,0,0)
| - axial extention of the coil along the Z-axis that it is equal to L = p * n
|----------------------------------
*add_points
0
0
0
0
0
4.5*15
|----------------------------------
| Define the position of the coil mean radius by a node, ID 1.
|----------------------------------
*add_nodes
20
0
0
|----------------------------------
| The coil will be modelled by 1D elements,
| its axis is coherent with the Z-axis of the model global system, 
| therefore by imposing a roto-translation expantion with the center O
| the coil is modelled.
| 32 divisions are required in the circumferential direction (Z rotation)
| 15/32 number of circumferential divisions must be applied along the coil pitch (Z traslation)
| number of repetitions: total coil * number of circumferential direction
|----------------------------------
*expand_reset
*set_expand_point
0
0
0
*set_expand_rotation z 360/32
*set_expand_translation z 15/32
*set_expand_repetitions 32*4.5
*expand_remove
*expand_nodes
 1
# | End of List
|----------------------------------
| Remove the unused node by the SWEEP menu
|----------------------------------
*set_sweep_tolerance
0.0001
*sweep_nodes
all_existing
|----------------------------------
| The command prompt evidences this message: 
| Deleting 1 duplicate nodes!
| Deleting 0 collapsed elements!
|----------------------------------
|----------------------------------
| 	MESH GENERATION FINISH
|----------------------------------
|----------------------------------
| 	GEOMETRIC PROPERTIES START
|----------------------------------
*new_geometry 
*geometry_type mech_three_beam_gen
*geometry_name de12_di6
*geometry_option section:circular
*geometry_param radius
(12+6)/4
*geometry_param wall_thick
(12-6)/2
*geometry_param radius
*add_geometry_elements
all_existing
|----------------------------------
| 	3D plot of the geometric properties 
|   defined previously and applied to the line elements 
|----------------------------------
*identify_geometries 
*regen
*draw_beams_3d on
*draw_offsets on
*elements_solid
*regenerate
*draw_beams_3d off
*regenerate
*identify_none 
*regen
|----------------------------------
| 	GEOMETRIC PROPERTIES FINISH
|----------------------------------
|----------------------------------
| 	MATERIAL START
|----------------------------------
*new_mater standard 
*mater_option general:state:solid
*mater_name titanium
*mater_param general:mass_density 4.7e-9
*mater_param structural:youngs_modulus
110000
*mater_param structural:poissons_ratio
0.3
*add_mater_elements
all_existing
|----------------------------------
| 	MATERIAL FINISH
|----------------------------------
|----------------------------------
| 	LINK START
|----------------------------------
| I create two nodes at the coil axis at z=0
| and at the other free extremity z = 67.5
|----------------------------------
*add_nodes
0
0
0
0
0
4.5*15
|----------------------------------
| 	RBE2 at Z=0
|----------------------------------
*new_rbe2
*rbe2_name rbe2_z0
*rbe2_ret_node
 291
*add_rbe2_tied_nodes
 1
# | End of List
*rbe2_tied_dof 1
*rbe2_tied_dof 2
*rbe2_tied_dof 3
*rbe2_tied_dof 4
*rbe2_tied_dof 5
*rbe2_tied_dof 6
|----------------------------------
| 	RBE2 at Z=67.5
|----------------------------------
*new_rbe2
*rbe2_name rbe2_z67.5
*rbe2_ret_node
 292
*add_rbe2_tied_nodes
 290
# | End of List
*rbe2_tied_dof 1
*rbe2_tied_dof 2
*rbe2_tied_dof 3
*rbe2_tied_dof 4
*rbe2_tied_dof 5
*rbe2_tied_dof 6
|----------------------------------
| 	LINK FINISH
|---------------------------------
*save_model