TOCHNOG User's manual - a free explicit/implicit FE program

Dennis Roddeman

• Table of Contents
• Basic information
  • Introduction
  • Legal Issues
  • Distribution files and installation
  • How to perform a calculation and how to get started
  • Plotting programs
  • Space discretization, time discretization
  • Program capabilities
  • Files used by Tochnog
• Equations
  • Convection and diffusion of heat
  • Material deformation and flow
  • Contact analysis
  • Ground water flow
  • Wave equation
• Examples
  • Example 1: Backward facing step (fluid flow)
  • Example 2: Confined compression of fluid filled porous material
  • Example 3: Plasticity in plate with circular hole
  • Example 4: Two-dimensional convection and diffusion
  • Example 5: Inverse modeling of indenter test
  • Example 6: Impact of metal plate, local mesh refinement
  • Example 7: Global mesh refinement over polynomial domain
  • Example 8: Residual remeshing near heat spot
  • Example 9: 2D Impact of metal plate, rebuilding of mesh after large deformations.
  • Example 10: Shear of low-tension material (cracking)
  • Example 11: Propagation of a disturbance in the wave equation.
  • Example 12: Contact frictional heat generation.
  • Example 13: Automatically embedded tendon with local mesh refinement
  • Example 14: Continuous metal forming.
  • Example 15: Generation of two holes in 2D mesh.
  • Example 16: Bearing capacity of foundation.
  • Example 17: Hertz contact problem.
  • Example 18: Thermally induced stresses in plate.
  • Example 19: Nonlocal plasticity in softening bar.
  • Example 20: Three dimensional local mesh refinement near heat source.
  • Example 22: Constant volume cylindrical compression with Modified Cam Clay.
  • Example 23: Automatic local mesh refinement near shear band.
  • Example 24: 3D Impact of metal plate, rebuilding of mesh after large deformations.
• Input file, general remarks
• Input file, initialization part
  • echo switch (first record of initialization part)
  • number_of_space_dimensions number_of_space_dimensions
    (second record of initialization part)
  • derivatives (third record of initialization part, if specified)
  • condif_temperature
  • groundflow_velocity
  • groundflow_pressure
  • materi_damage
  • materi_density
  • materi_diffusion
  • materi_displacement
  • materi_history_variables number_of_variables
  • materi_maxwell_stress number_of_chains
  • materi_plasti_f
  • materi_plasti_f_nonlocal
  • materi_plasti_kappa
  • materi_plasti_rho
  • materi_strainenergy
  • materi_strain_elasti
  • materi_strain_intergranular
  • materi_strain_plasti
  • materi_strain_total
  • materi_stress
  • materi_velocity
  • materi_velocity_integrated
  • materi_void_fraction
  • materi_work
  • residue
  • wave_scalar
  • wave_fscalar
  • end_initia (last record of initialization part)
• Input file, data part
  • area_element_group index geometry_entity_item geometry_entity_index element_group
  • area_element_group_sequence_element index name
  • area_element_group_sequence_elementgroup index group_0 group_1 ...
  • area_element_group_sequence_geometry index geometry_entity_item geometry_entity_index
  • area_element_group_sequence_time index time_0 time_1 ...
  • area_node_dataitem index geometry_entity_item geometry_entity_index data_item_name
  • area_node_dataitem_double index value_0 value_1 ...
  • area_node_dataitem_integer index value_0 value_1 ...
  • bounda_force index node_range unknown_0 unknown_1 ...
  • bounda_sine index start_time freq_0 amp_0 freq_1 amp_1 ...
  • bounda_time index time load time load ...
  • bounda_time_file index switch
  • bounda_time_user index switch
  • bounda_unknown index node_range unknown_0 unknown_1 ...
  • change_dataitem index data_item_name data_item_index data_item_number operat
  • change_dataitem_time index time value ...
  • change_dataitem_discrete index switch
  • change_dataitem_time_user index switch
  • check_combination index switch
  • condif_convection index a_c T_r
  • condif_convection_geometry index node_0 node_1 node_2 ...
  • condif_radiation index a_r T_r
  • condif_radiation_geometry index node_0 node_1 node_2 ...
  • contact_friction friction
  • contact_geometry index geometry_entity_item geometry_entity_index
  • contact_geometry_switch index switch
  • contact_heatgeneration factor
  • contact_penalty_pressure pressure_penalty
  • contact_penalty_temperature temperature_penalty
  • contact_penalty_velocity velocity_penalty
  • control_data_delete index data_item_name index_range
  • control_data_put index data_item_name index_range
  • control_data_put_double index ...
  • control_data_put_double_node index ...
  • control_data_put_integer index ...
  • control_distribute index distribution_type data_item_name data_item_number ...
  • control_distribute_values index distribution_value distribution_value
  • control_eigen index type_of_eigen_analysis number_of_modes
  • control_eigen_scale index largest_component
  • control_eigen_values index eigen_value_0 eigen_value_1 ...
  • control_mesh_adjust_geometry index geometry_entity_item_0
    geometry_entity_index_0 geometry_entity_item_1
    geometry_entity_index_1
  • control_mesh_delete_geometry index geometry_entity_item geometry_entity_index
  • control_mesh_delete_geometry_element index element_name_0 element_name_0 ...
  • control_mesh_delete_geometry_elementgroup index element_group_0 element_group_1 ...
  • control_mesh_delete_geometry_factor index factor_0 factor_1 ...
  • control_mesh_delete_geometry_movenodes index switch
  • control_mesh_delete_small index eps
  • control_mesh_generate_beam index element_group geometry_entity_item geometry_entity_index
  • control_mesh_generate_contactspring index element_group geometry_entity_item geometry_entity_index
  • control_mesh_generate_contactspring_element index element_0 element_1
  • control_mesh_generate_spring1 index element_group geometry_entity_item geometry_entity_index
  • control_mesh_generate_spring2 index element_group geometry_entity_item geometry_entity_index
  • control_mesh_generate_truss index element_group geometry_entity_item geometry_entity_index
  • control_mesh_generate_trussbeam index element_group geometry_entity_item geometry_entity_index
  • control_mesh_generate_truss_beam_loose index switch
  • control_mesh_generate_truss_beam_macro index macro_0 macro_1 ...
  • control_mesh_macro index macro_item element_group
    n ...
  • control_mesh_macro_element index element_type
  • control_mesh_macro_parameters index x y ...
  • control_mesh_macro_set_node_boundary index switch
  • control_mesh_merge index switch
  • control_mesh_merge_epscoord index epsilon
  • control_mesh_merge_macro index macro_0 macro_1 ...
  • control_mesh_merge_not index geometry_entity_item geometry_entity_index
  • control_mesh_new_mesh index delta
  • control_mesh_new_mesh_element index element_type
  • control_mesh_refine_globally index refinement_type
  • control_mesh_refine_globally_geometry index geometry_entity_item
    geometry_entity_index
  • control_mesh_refine_locally index percentage
  • control_mesh_refine_locally_geometry index geometry_entity_item
    geometry_entity_index
  • control_mesh_refine_locally_not index geometry_entity_0 geometry_entity_index_0
  • control_mesh_refine_locally_only index geometry_entity_0 geometry_entity_index_0
  • control_mesh_refine_locally_unknown index unknown
  • control_mesh_remesh index switch
  • control_mesh_remesh_factor index geometry_factor residue_factor
  • control_mesh_renumber index lowest_element lowest_node
  • control_mesh_split index switch
  • control_mesh_split_only index geometry_entity geometry_entity_index
  • control_options_convection index switch
  • control_options_inertia index switch
  • control_options_relaxation relax_0 relax_1 ...
  • control_options_skip_plasticity index switch
  • control_options_solver index solver_type
  • control_options_solver_petsc_ksptype index ksptype
  • control_options_solver_petsc_pctype index pctype
  • control_print index data_item_name_0 data_item_name_1 ...
  • control_print_database index switch
  • control_print_data_versus_data index data_item_name_0 index_0 number_0
    data_item_name_1 index_1 number_1 ...
  • control_print_dx index switch
  • control_print_element index data_item_name
  • control_print_filter index print_filter_index_0 print_filter_index_1 ...
  • control_print_gid index switch
  • control_print_gmv index switch
  • control_print_history index data_item_name_0 data_item_index_0 number_0 ...
  • control_print_matlab index switch
  • control_print_plotmtv index switch
  • control_print_tecplot index switch
  • control_print_unknowns index switch
  • control_print_unknownsrhside index switch
  • control_print_vtk index switch
  • control_relaxation_condif_temperature index group_condif_temperature_relaxation
  • control_relaxation_groundflow_pressure index group_groundflow_pressure_relaxation
  • control_relaxation_materi_velocity index group_materi_velocity_relaxation
  • control_repeat index number_of_repeats control_index
  • control_repeat_until_item index control_index repeat_criterium data_item_name
    data_item_index data_item_number
  • control_repeat_until_tolerance index tolerance
  • control_restart index switch
  • control_timestep index step_size time_increment step_size time_increment ...
  • control_timestep_iterations index number_of_iterations
  • control_timestep_iterations_automatic index ratio_criterium maximum_timestep
  • control_timestep_iterations_automatic_stop index switch
  • control_timestep_multiplier index multiplier
  • control_unknown_freeze index unknown_0 unknown_1 ...
  • control_unknown_reset_geometry index geometry_item_name geometry_item_index
  • control_unknown_reset_unknown index unknown_0 unknown_1 ...
  • control_unknown_reset_value index value
  • dependency_diagram index unknown_value_0 ...data_item_value_0 ...
  • dependency_item index data_item element_group unknown n
  • dof_label unknown_0 unknown_1 ...
  • element index element_name node_0 node_1 node_2 ...
  • element_beam_direction index dir_x dir_y dir_z
  • element_beam_moment index force_x_first_node force_y_first_node moment_z_first_node force_x_second_node force_y_second_node moment_z_second_node ...
  • element_group index element_group
  • element_contactspring_direction index dirN_x dirN_y dirN_z dirT1_x dirT1_y dirT1_z dirT2_x dirT2_y dirT2_z
  • element_contactspring_force index force_N force_T1 force_T1
  • element_spring_force index force
  • element_tendon_direction index dir_x dir_y dir_z dir_x dir_y dir_z ...
  • element_tendon_number index tendon_0 tendon_1 ...
  • element_tendon_strain index eps_0 eps_1 ...
  • element_tendon_stress index sig_0 sig_1 ...
  • element_tendon_volume index volume_0 volume_1 ...
  • element_truss_direction index dir_x dir_y dir_z
  • element_truss_force index force
  • force_element_edge index force_0 force_1 ...
  • force_element_edge_factor index a₀ a₁ ...a_n
  • force_element_edge_geometry index node_0 node_1 node_2 ...
  • force_element_edge_sine index start_time freq_0 amp_0 freq_1 amp_1 ...
  • force_element_edge_time index time load time load ...
  • force_element_edge_normal index force
  • force_element_edge_normal_factor index a₀ a₁ ...a_n
  • force_element_edge_normal_geometry index geometry_item_name geometry_item_index
  • force_element_edge_normal_sine index start_time freq_0 amp_0 freq_1 amp_1 ...
  • force_element_edge_normal_time index time load time load ...
  • force_element_edge_water index density_water g dir_x dir_y dir_z
  • force_element_edge_water_geometry index geometry_item_name geometry_item_index
  • force_element_edge_water_time index time load time load ...
  • force_element_volume index force_0 force_1 ...
  • force_element_volume_factor index a₀ a₁ ...a_n
  • force_element_volume_geometry index element_0 element_1 element_2 ...
  • force_element_volume_sine index start_time freq_0 amp_0 freq_1 amp_1 ...
  • force_element_volume_time index time load time load ...
  • force_gravity g_x g_y g_z
  • force_gravity_time time load time load ...
  • geometry_bounda_factor index factor_0 ...
  • geometry_bounda_sine_x index a b
  • geometry_bounda_sine_y index a b
  • geometry_bounda_sine_z index a b
  • geometry_brick index x_c y_c z_c l_x l_y l_z tolerance
  • geometry_circle index x_c y_c radius tolerance
  • geometry_circle_segment index x_c y_c radius side_x side_y tolerance
  • geometry_cylinder index x_0 y_0 z_0 x_1 y_1 z_1 radius tolerance
  • geometry_cylinder_segment index x_0 y_0 z_0 x_1 y_1 z_1 radius side_x side_y side_z tolerance
  • geometry_ellipse index x_c y_c a b tolerance
  • geometry_line index x_0 y_0 z_0 x_1 y_1 z_1 radius
  • geometry_point index x y z radius
  • geometry_polynomial index a₀ a₁ ...a_n x_0 x_1 y_0 y_1 tolerance
  • geometry_quadrilateral index x_0 y_0 z_0 x_1 y_1 z_1 x_2 y_2 z_2
    x_3 y_3 z_3 tolerance
  • geometry_set index geometry_entity_0 geometry_entity_index_0 geometry_entity_1
    geometry_entity_index_1 ...
  • geometry_sphere index x_c y_c z_c radius tolerance
  • geometry_sphere_segment index x_c y_c z_c radius side_x side_y side_z tolerance
  • geometry_triangle index x_0 y_0 z_0 x_1 y_1 z_1 x_2 y_2 z_2 tolerance
  • groundflow_density r
  • groundflow_phreaticlevel ...
  • groundflow_phreaticlevel_n nx ny
  • groundflow_phreaticlevel_bounda method
  • group_axisymmetric index switch
  • group_beam_inertia index I
  • group_beam_memory index memory_type
  • group_beam_plane index first_direction second_direction
  • group_beam_young index E
  • group_condif_absorption index a
  • group_condif_capacity index C
  • group_condif_conductivity index k
  • group_condif_density index density
  • group_condif_flow index beta₁   beta₂   beta₃
  • group_contactspring_cohesion index c
  • group_contactspring_direction index dir_x   dir_y   dir_z
  • group_contactspring_memory index memory_type
  • group_contactspring_stiffness index k_N k_T
  • group_contactspring_friction index f
  • group_contactspring_friction_automatic index switch
  • group_groundflow_capacity index C
  • group_groundflow_materidivergence index switch
  • group_groundflow_permeability index pe_x pe_y pe_z
  • group_integration_method index method
  • group_integration_points index type
  • group_materi_damage_mazars index e₀ a_t   b_t   a_c   b_c   b
  • group_materi_damping index d
  • group_materi_density index density
  • group_materi_density_groundflow index density_wet density_dry
  • group_materi_elasti_camclay_g index G
  • group_materi_elasti_camclay_poisson index n
  • group_materi_elasti_compressibility index co
  • group_materi_elasti_lade index B R l
  • group_materi_elasti_poisson index poisson
  • group_materi_elasti_transverse_isotropy index dir_0 dir_1 dir_2 C_aaaa C_bbbb
    C_aabb C_abab C_bcbc
  • group_materi_elasti_volumetric_poisson index n
  • group_materi_elasti_volumetric_young_order index n
  • group_materi_elasti_volumetric_young_values index e_0   s_0   e_1   s₁   ...
  • group_materi_elasti_young index E
  • group_materi_elasti_young_polynomial index E₀   E₁ ...
  • group_materi_elasti_young_power index p₀ E₀ a
  • group_materi_expansion_linear index a
  • group_materi_expansion_volume index b
  • group_materi_failure_cruching index threshold delete_time
  • group_materi_failure_damage index threshold delete_time
  • group_materi_failure_plasti_kappa index threshold delete_time
  • group_materi_failure_rupture index threshold delete_time
  • group_materi_failure_voidfraction index threshold delete_time
  • group_materi_hyper_besseling index K₁ K₂ a
  • group_materi_hyper_blatz_ko index G b
  • group_materi_hyper_mooney_rivlin index K₁ K₂
  • group_materi_hyper_neohookean index K₁
  • group_materi_hyper_reducedpolynomial index K₁ K₂ ...
  • group_materi_hyper_volumetric_linear index K
  • group_materi_hyper_volumetric_murnaghan index K b
  • group_materi_hyper_volumetric_ogden index K b
  • group_materi_hyper_volumetric_polynomial index K_0   K_1   ...
  • group_materi_hyper_volumetric_simotaylor index K
  • group_materi_maxwell_chain index E_0 t_0 ...E_n-1 t_n-1
  • group_materi_membrane index switch
  • group_materi_memory index memory_type
  • group_materi_plasti_boundary index group_0 group_1
  • group_materi_plasti_boundary_factor index factor
  • group_materi_plasti_camclay index M k l
  • group_materi_plasti_cap index c f a R e_v^p p_b ...
  • group_materi_plasti_compression index sigy
  • group_materi_plasti_diprisco index g b_f b_p c_p t_p q_c q_e x_c x_e b_f⁰ C
  • group_materi_plasti_druckprag index phi c phiflow
  • group_materi_plasti_druckprag_tensioncutoff index switch
  • group_materi_plasti_gurson index sigy q1 q2 q3
  • group_materi_plasti_hlc index sigy m21 m22 m23 m1
  • group_materi_plasti_heatgeneration factor
  • group_materi_plasti_hypo_cohesion index c
  • group_materi_plasti_hypo_intergranularstrain index R m_R m_T b_r c
  • group_materi_plasti_hypo_pressuredependentvoidratio index switch
  • group_materi_plasti_hypo_wolffersdorff index j h_s n e_c0 e_d0 e_i0 alpha beta
  • group_materi_plasti_kinematic_hardening index a
  • group_materi_plasti_matsuokanakai index phi c phiflow
  • group_materi_plasti_mohrcoul index phi c phiflow
  • group_materi_plasti_mohrcoul_softening index phi_0 c_0 phiflow_0 phi_1 c_1 phiflow_1 kappa_crit
  • group_materi_plasti_mohrcoul_tensioncutoff index switch
  • group_materi_plasti_tension index sigy
  • group_materi_plasti_user index switch
  • group_materi_plasti_visco_always index switch
  • group_materi_plasti_visco_exponential index g a
  • group_materi_plasti_visco_power index h p f^ref
  • group_materi_plasti_vonmises index s_y0
  • group_materi_plasti_vonmises_nadai index C k₀ n
  • group_materi_stokes index switch
  • group_materi_viscosity index n
  • group_materi_viscosity_user index switch
  • group_materi_viscosity_bingham index s_y g m
  • group_materi_viscosity_exponential index n₀ m
  • group_materi_viscosity_heatgeneration switch
  • group_spring_direction index dir_x   dir_y   dir_z
  • group_spring_plasti index F_y
  • group_spring_stiffness index k
  • group_time index birth death
  • group_truss_area index A
  • group_truss_density index r
  • group_truss_memory index memory_type
  • group_truss_rope index switch
  • group_truss_plasti index s_y
  • group_truss_young index E
  • group_type index type_name_0 type_name_1 ...
  • group_user_umat index switch
  • group_user_data index ...
  • group_volume_factor index factor
  • group_wave_speed_of_sound index c
  • inverse_iterations number_of_iterations
  • inverse_iteration_number index number
  • inverse_parameter index data_item_name data_item_index data_item_number
  • inverse_parameter_limits index lower_value higher_value
  • inverse_parameter_step index size
  • inverse_parameter_variation index variation
  • inverse_target index data_item_name data_item_index data_item_number
  • inverse_target_data index target_value penalty
  • inverse_target_timestep index timestep_index
  • materi_diffusion_minimum c_minimum
  • node index coord_0 coord_1 coord_2
  • node_boundary index switch
  • node_damping index damping_x damping_y damping_z
  • node_dof index dof_0 dof_1 ...
  • node_dof_calcul index ...
  • node_dof_start_refined index dof_0 dof_1 ...
  • node_eigen index unknown_0 unknown_1 ...unknown_0 unknown_1 ...
  • node_mass index mass_x mass_y mass_z
  • node_rhside index rhside_0 rhside_1 ...
  • node_start_refined index coord_0 coord_1 coord_2
  • node_stiffness index stiffness_x stiffness_y stiffness_z
  • options_convection switch
  • options_elementloop switch
  • options_inertia switch
  • options_matrix_group switch
  • options_matrix_length matrix_length
  • options_mesh specifier_x specifier_y specifier_z
  • options_nonlocal nonlocal_radius
  • options_processors index nproc
  • options_relaxation relax_0 relax_1 ...
  • options_residuefactor factor_0 factor_1 ...
  • options_skip_plasticity switch
  • options_solver solver
  • options_stabilization switch
  • post_calcul unknown operat ...
  • post_error_item index data_item_name data_item_index number
  • post_error_result index percentage
  • post_integrate index data_item_name data_item_index data_item_number ...
  • post_global index name_0 name_1 ...
  • post_integrate_result index result
  • post_line index x_0 y_0 z_0 x_1 y_1 z_1
  • post_line_operat index operat
  • post_line_dof index dof_0 dof_1 ...
  • post_line_dof_calcul ...
  • post_line_moment index switch
  • post_line_n index n
  • post_node index data_item operat geometry_entity geometry_entity_index
  • post_node_result index result_0 result_1 ...
  • post_node_rhside_fixed value_0 value_1 ...
  • post_node_rhside_free value_0 value_1 ...
  • post_node_rhside_ratio ratio
  • post_node_rhside_ratio_unknowntypes unknowntype_0 ...
  • post_point index x y z
  • post_point_dof index dof_0 dof_1 ...
  • post_point_dof_calcul ...
  • post_quadrilateral index x_0 y_0 z_0 x_1 y_1 z_1 x_2 y_2 z_2 x_3 y_3 z_3
  • post_quadrilateral_dof index dof_0 dof_1 ...
  • post_quadrilateral_dof_calcul ...
  • post_quadrilateral_n index n
  • print_arithmetic switch
  • print_control switch
  • print_define switch
  • print_failure switch
  • print_filter index data_item_name data_item_index number_0 number_1 ...
  • print_lastdatabase switch
  • print_where switch
  • slide_geometry index geometry_entity geometry_entity_index
  • slide_friction index friction_coefficient
  • slide_penalty index penalty_coefficient
  • target_item index data_item_name data_item_index number
  • target_value index value tolerance
  • tendon index x_0 y_0 z_0 x_1 y_1 z_1 tendon_area
  • tendon_elasti index stiffness
  • tendon_expansion index a
  • tendon_plasti index yield_stress
  • tendon_stress index initial_stress
  • tendon_stress_time index time_0 factor_0 time_1 factor_1 ...
  • time_calculation cpu_time_in_seconds
  • time_current current_time
  • volume_factor a₀ a₁ ...a_n
  • end_data (last record of data part)
• Runtime file
• Interaction analyzes and advanced analyzes
  • Fluid-structure interaction
  • Consolidation analysis: ground water flow in deforming solid
  • Heat transport in ground water flow
  • Heat transport in materials
• Practical tips
  • Continuing an analysis
  • Saving CPU time.
  • Saving computer memory.
  • To model geometrically linear material
  • When you find yourself in times of troubles
• User supplied subroutines
• References
• GNU GENERAL PUBLIC LICENSE

This document was translated from L^AT_EX by H^EV^EA and H^AC^HA.