# This is an example of an indata file for the Impact Finite Element Program
#
# It is a simple mix of rod and solid elements which are connected to each other
# The model is fixed in the upper left corner and is pulled by a constant
# force in the lower right corner.
#
# Sometimes this combination can be tricky and need a firm hand on the timestep
# control. In this case, it has been set to a fixed value. Autostep does not
# suffice.
#
# Please note that you can view the results with GID while the program is
# running so there is no need to wait for the completion of this simulation.
nodes
1 x = 0.0 y = 0.0 z = 0.0
2 x = 5.0 y = 0.0 z = 0.0
3 x = 10.0 y = 0.0 z = 0.0
4 x = 15.0 y = 0.0 z = 0.0
5 x = 20.0 y = 0.0 z = 0.0
6 x = 25.0 y = 0.0 z = 0.0
7 x = 40.0 y = 0.0 z = 0.0
8 x = 45.0 y = 0.0 z = 0.0
9 x = 50.0 y = 0.0 z = 0.0
10 x = 55.0 y = 0.0 z = 0.0
11 x = 70.0 y = 0.0 z = 0.0
12 x = 75.0 y = 0.0 z = 0.0
13 x = 80.0 y = 0.0 z = 0.0
14 x = 85.0 y = 0.0 z = 0.0
15 x = 90.0 y = 0.0 z = 0.0
16 x = 95.0 y = 0.0 z = 0.0
17 x = 100.0 y = 0.0 z = 0.0
18 x = 105.0 y = 0.0 z = 0.0
19 x = 115.0 y = 0.0 z = 0.0
20 x = 120.0 y = 0.0 z = 0.0 load = end_load
21 x = 0.0 y = 5.0 z = 0.0
22 x = 5.0 y = 5.0 z = 0.0
23 x = 10.0 y = 5.0 z = 0.0
24 x = 15.0 y = 5.0 z = 0.0
25 x = 20.0 y = 5.0 z = 0.0
26 x = 25.0 y = 5.0 z = 0.0
27 x = 40.0 y = 5.0 z = 0.0
28 x = 45.0 y = 5.0 z = 0.0
29 x = 50.0 y = 5.0 z = 0.0
30 x = 55.0 y = 5.0 z = 0.0
31 x = 70.0 y = 5.0 z = 0.0
32 x = 75.0 y = 5.0 z = 0.0
33 x = 80.0 y = 5.0 z = 0.0
34 x = 85.0 y = 5.0 z = 0.0
35 x = 90.0 y = 5.0 z = 0.0
36 x = 95.0 y = 5.0 z = 0.0
37 x = 100.0 y = 5.0 z = 0.0
38 x = 105.0 y = 5.0 z = 0.0
39 x = 115.0 y = 5.0 z = 0.0
40 x = 120.0 y = 5.0 z = 0.0
41 x = 5.0 y = 10.0 z = 0.0
42 x = 10.0 y = 10.0 z = 0.0
43 x = 20.0 y = 10.0 z = 0.0
44 x = 25.0 y = 10.0 z = 0.0
45 x = 40.0 y = 10.0 z = 0.0
46 x = 45.0 y = 10.0 z = 0.0
47 x = 50.0 y = 10.0 z = 0.0
48 x = 55.0 y = 10.0 z = 0.0
49 x = 60.0 y = 10.0 z = 0.0
50 x = 65.0 y = 10.0 z = 0.0
51 x = 70.0 y = 10.0 z = 0.0
52 x = 75.0 y = 10.0 z = 0.0
53 x = 80.0 y = 10.0 z = 0.0
54 x = 85.0 y = 10.0 z = 0.0
55 x = 90.0 y = 10.0 z = 0.0
56 x = 95.0 y = 10.0 z = 0.0
57 x = 115.0 y = 10.0 z = 0.0
58 x = 120.0 y = 10.0 z = 0.0
59 x = 5.0 y = 15.0 z = 0.0
60 x = 10.0 y = 15.0 z = 0.0
61 x = 20.0 y = 15.0 z = 0.0
62 x = 25.0 y = 15.0 z = 0.0
63 x = 30.0 y = 15.0 z = 0.0
64 x = 35.0 y = 15.0 z = 0.0
65 x = 40.0 y = 15.0 z = 0.0
66 x = 45.0 y = 15.0 z = 0.0
67 x = 50.0 y = 15.0 z = 0.0
68 x = 55.0 y = 15.0 z = 0.0
69 x = 60.0 y = 15.0 z = 0.0
70 x = 65.0 y = 15.0 z = 0.0
71 x = 70.0 y = 15.0 z = 0.0
72 x = 75.0 y = 15.0 z = 0.0
73 x = 80.0 y = 15.0 z = 0.0
74 x = 85.0 y = 15.0 z = 0.0
75 x = 90.0 y = 15.0 z = 0.0
76 x = 95.0 y = 15.0 z = 0.0
77 x = 115.0 y = 15.0 z = 0.0
78 x = 120.0 y = 15.0 z = 0.0
79 x = 5.0 y = 20.0 z = 0.0
80 x = 10.0 y = 20.0 z = 0.0
81 x = 20.0 y = 20.0 z = 0.0
82 x = 25.0 y = 20.0 z = 0.0
83 x = 30.0 y = 20.0 z = 0.0
84 x = 35.0 y = 20.0 z = 0.0
85 x = 40.0 y = 20.0 z = 0.0
86 x = 45.0 y = 20.0 z = 0.0
87 x = 50.0 y = 20.0 z = 0.0
88 x = 55.0 y = 20.0 z = 0.0
89 x = 60.0 y = 20.0 z = 0.0
90 x = 65.0 y = 20.0 z = 0.0
91 x = 70.0 y = 20.0 z = 0.0
92 x = 75.0 y = 20.0 z = 0.0
93 x = 80.0 y = 20.0 z = 0.0
94 x = 85.0 y = 20.0 z = 0.0
95 x = 90.0 y = 20.0 z = 0.0
96 x = 95.0 y = 20.0 z = 0.0
97 x = 100.0 y = 20.0 z = 0.0
98 x = 105.0 y = 20.0 z = 0.0
99 x = 110.0 y = 20.0 z = 0.0
100 x = 115.0 y = 20.0 z = 0.0
101 x = 120.0 y = 20.0 z = 0.0
102 x = 125.0 y = 20.0 z = 0.0
103 x = 5.0 y = 25.0 z = 0.0
104 x = 10.0 y = 25.0 z = 0.0
105 x = 20.0 y = 25.0 z = 0.0
106 x = 25.0 y = 25.0 z = 0.0
107 x = 30.0 y = 25.0 z = 0.0
108 x = 35.0 y = 25.0 z = 0.0
109 x = 40.0 y = 25.0 z = 0.0
110 x = 45.0 y = 25.0 z = 0.0
111 x = 50.0 y = 25.0 z = 0.0
112 x = 55.0 y = 25.0 z = 0.0
113 x = 60.0 y = 25.0 z = 0.0
114 x = 65.0 y = 25.0 z = 0.0
115 x = 70.0 y = 25.0 z = 0.0
116 x = 75.0 y = 25.0 z = 0.0
117 x = 80.0 y = 25.0 z = 0.0
118 x = 85.0 y = 25.0 z = 0.0
119 x = 90.0 y = 25.0 z = 0.0
120 x = 95.0 y = 25.0 z = 0.0
121 x = 100.0 y = 25.0 z = 0.0
122 x = 105.0 y = 25.0 z = 0.0
123 x = 110.0 y = 25.0 z = 0.0
124 x = 115.0 y = 25.0 z = 0.0
125 x = 120.0 y = 25.0 z = 0.0
126 x = 125.0 y = 25.0 z = 0.0
127 x = 0.0 y = 30.0 z = 0.0
128 x = 5.0 y = 30.0 z = 0.0
129 x = 10.0 y = 30.0 z = 0.0
130 x = 15.0 y = 30.0 z = 0.0
131 x = 0.0 y = 35.0 z = 0.0
132 x = 5.0 y = 35.0 z = 0.0
133 x = 10.0 y = 35.0 z = 0.0
134 x = 15.0 y = 35.0 z = 0.0
1001 x = 0.0 y = 0.0 z = -5.0
1002 x = 5.0 y = 0.0 z = -5.0
1003 x = 10.0 y = 0.0 z = -5.0
1004 x = 15.0 y = 0.0 z = -5.0
1005 x = 20.0 y = 0.0 z = -5.0
1006 x = 25.0 y = 0.0 z = -5.0
1007 x = 40.0 y = 0.0 z = -5.0
1008 x = 45.0 y = 0.0 z = -5.0
1009 x = 50.0 y = 0.0 z = -5.0
1010 x = 55.0 y = 0.0 z = -5.0
1011 x = 70.0 y = 0.0 z = -5.0
1012 x = 75.0 y = 0.0 z = -5.0
1013 x = 80.0 y = 0.0 z = -5.0
1014 x = 85.0 y = 0.0 z = -5.0
1015 x = 90.0 y = 0.0 z = -5.0
1016 x = 95.0 y = 0.0 z = -5.0
1017 x = 100.0 y = 0.0 z = -5.0
1018 x = 105.0 y = 0.0 z = -5.0
1019 x = 115.0 y = 0.0 z = -5.0
1020 x = 120.0 y = 0.0 z = -5.0
1021 x = 0.0 y = 5.0 z = -5.0
1022 x = 5.0 y = 5.0 z = -5.0
1023 x = 10.0 y = 5.0 z = -5.0
1024 x = 15.0 y = 5.0 z = -5.0
1025 x = 20.0 y = 5.0 z = -5.0
1026 x = 25.0 y = 5.0 z = -5.0
1027 x = 40.0 y = 5.0 z = -5.0
1028 x = 45.0 y = 5.0 z = -5.0
1029 x = 50.0 y = 5.0 z = -5.0
1030 x = 55.0 y = 5.0 z = -5.0
1031 x = 70.0 y = 5.0 z = -5.0
1032 x = 75.0 y = 5.0 z = -5.0
1033 x = 80.0 y = 5.0 z = -5.0
1034 x = 85.0 y = 5.0 z = -5.0
1035 x = 90.0 y = 5.0 z = -5.0
1036 x = 95.0 y = 5.0 z = -5.0
1037 x = 100.0 y = 5.0 z = -5.0
1038 x = 105.0 y = 5.0 z = -5.0
1039 x = 115.0 y = 5.0 z = -5.0
1040 x = 120.0 y = 5.0 z = -5.0
1041 x = 5.0 y = 10.0 z = -5.0
1042 x = 10.0 y = 10.0 z = -5.0
1043 x = 20.0 y = 10.0 z = -5.0
1044 x = 25.0 y = 10.0 z = -5.0
1045 x = 40.0 y = 10.0 z = -5.0
1046 x = 45.0 y = 10.0 z = -5.0
1047 x = 50.0 y = 10.0 z = -5.0
1048 x = 55.0 y = 10.0 z = -5.0
1049 x = 60.0 y = 10.0 z = -5.0
1050 x = 65.0 y = 10.0 z = -5.0
1051 x = 70.0 y = 10.0 z = -5.0
1052 x = 75.0 y = 10.0 z = -5.0
1053 x = 80.0 y = 10.0 z = -5.0
1054 x = 85.0 y = 10.0 z = -5.0
1055 x = 90.0 y = 10.0 z = -5.0
1056 x = 95.0 y = 10.0 z = -5.0
1057 x = 115.0 y = 10.0 z = -5.0
1058 x = 120.0 y = 10.0 z = -5.0
1059 x = 5.0 y = 15.0 z = -5.0
1060 x = 10.0 y = 15.0 z = -5.0
1061 x = 20.0 y = 15.0 z = -5.0
1062 x = 25.0 y = 15.0 z = -5.0
1063 x = 30.0 y = 15.0 z = -5.0
1064 x = 35.0 y = 15.0 z = -5.0
1065 x = 40.0 y = 15.0 z = -5.0
1066 x = 45.0 y = 15.0 z = -5.0
1067 x = 50.0 y = 15.0 z = -5.0
1068 x = 55.0 y = 15.0 z = -5.0
1069 x = 60.0 y = 15.0 z = -5.0
1070 x = 65.0 y = 15.0 z = -5.0
1071 x = 70.0 y = 15.0 z = -5.0
1072 x = 75.0 y = 15.0 z = -5.0
1073 x = 80.0 y = 15.0 z = -5.0
1074 x = 85.0 y = 15.0 z = -5.0
1075 x = 90.0 y = 15.0 z = -5.0
1076 x = 95.0 y = 15.0 z = -5.0
1077 x = 115.0 y = 15.0 z = -5.0
1078 x = 120.0 y = 15.0 z = -5.0
1079 x = 5.0 y = 20.0 z = -5.0
1080 x = 10.0 y = 20.0 z = -5.0
1081 x = 20.0 y = 20.0 z = -5.0
1082 x = 25.0 y = 20.0 z = -5.0
1083 x = 30.0 y = 20.0 z = -5.0
1084 x = 35.0 y = 20.0 z = -5.0
1085 x = 40.0 y = 20.0 z = -5.0
1086 x = 45.0 y = 20.0 z = -5.0
1087 x = 50.0 y = 20.0 z = -5.0
1088 x = 55.0 y = 20.0 z = -5.0
1089 x = 60.0 y = 20.0 z = -5.0
1090 x = 65.0 y = 20.0 z = -5.0
1091 x = 70.0 y = 20.0 z = -5.0
1092 x = 75.0 y = 20.0 z = -5.0
1093 x = 80.0 y = 20.0 z = -5.0
1094 x = 85.0 y = 20.0 z = -5.0
1095 x = 90.0 y = 20.0 z = -5.0
1096 x = 95.0 y = 20.0 z = -5.0
1097 x = 100.0 y = 20.0 z = -5.0
1098 x = 105.0 y = 20.0 z = -5.0
1099 x = 110.0 y = 20.0 z = -5.0
1100 x = 115.0 y = 20.0 z = -5.0
1101 x = 120.0 y = 20.0 z = -5.0
1102 x = 125.0 y = 20.0 z = -5.0
1103 x = 5.0 y = 25.0 z = -5.0
1104 x = 10.0 y = 25.0 z = -5.0
1105 x = 20.0 y = 25.0 z = -5.0
1106 x = 25.0 y = 25.0 z = -5.0
1107 x = 30.0 y = 25.0 z = -5.0
1108 x = 35.0 y = 25.0 z = -5.0
1109 x = 40.0 y = 25.0 z = -5.0
1110 x = 45.0 y = 25.0 z = -5.0
1111 x = 50.0 y = 25.0 z = -5.0
1112 x = 55.0 y = 25.0 z = -5.0
1113 x = 60.0 y = 25.0 z = -5.0
1114 x = 65.0 y = 25.0 z = -5.0
1115 x = 70.0 y = 25.0 z = -5.0
1116 x = 75.0 y = 25.0 z = -5.0
1117 x = 80.0 y = 25.0 z = -5.0
1118 x = 85.0 y = 25.0 z = -5.0
1119 x = 90.0 y = 25.0 z = -5.0
1120 x = 95.0 y = 25.0 z = -5.0
1121 x = 100.0 y = 25.0 z = -5.0
1122 x = 105.0 y = 25.0 z = -5.0
1123 x = 110.0 y = 25.0 z = -5.0
1124 x = 115.0 y = 25.0 z = -5.0
1125 x = 120.0 y = 25.0 z = -5.0
1126 x = 125.0 y = 25.0 z = -5.0
1127 x = 0.0 y = 30.0 z = -5.0
1128 x = 5.0 y = 30.0 z = -5.0
1129 x = 10.0 y = 30.0 z = -5.0
1130 x = 15.0 y = 30.0 z = -5.0
1131 x = 0.0 y = 35.0 z = -5.0 constraint = fixed
1132 x = 5.0 y = 35.0 z = -5.0
1133 x = 10.0 y = 35.0 z = -5.0
1134 x = 15.0 y = 35.0 z = -5.0
# nip stands for number of intergrationpoints and can also be 1
# but then you will experience numerical instability (hourglassing)
elements of type solid_iso_6
1 nodes = [1,2,1002,1001,21,22,1022,1021] nip = 8 material = e_steel
2 nodes = [2,3,1003,1002,22,23,1023,1022] nip = 8 material = e_steel
3 nodes = [3,4,1004,1003,23,24,1024,1023] nip = 8 material = e_steel
4 nodes = [22,23,1023,1022,41,42,1042,1041] nip = 8 material = e_steel
5 nodes = [41,42,1042,1041,59,60,1060,1059] nip = 8 material = e_steel
6 nodes = [59,60,1060,1059,79,80,1080,1079] nip = 8 material = e_steel
7 nodes = [79,80,1080,1079,103,104,1104,1103] nip = 8 material = e_steel
8 nodes = [103,104,1104,1103,128,129,1129,1128] nip = 8 material = e_steel
9 nodes = [127,128,1128,1127,131,132,1132,1131] nip = 8 material = e_steel
10 nodes = [128,129,1129,1128,132,133,1133,1132] nip = 8 material = e_steel
11 nodes = [129,130,1130,1129,133,134,1134,1133] nip = 8 material = e_steel
12 nodes = [5,6,1006,1005,25,26,1026,1025] nip = 8 material = e_steel
13 nodes = [25,26,1026,1025,43,44,1044,1043] nip = 8 material = e_steel
14 nodes = [43,44,1044,1043,61,62,1062,1061] nip = 8 material = e_steel
15 nodes = [61,62,1062,1061,81,82,1082,1081] nip = 8 material = e_steel
16 nodes = [81,82,1082,1081,105,106,1106,1105] nip = 8 material = e_steel
17 nodes = [82,83,1083,1082,106,107,1107,1106] nip = 8 material = e_steel
18 nodes = [83,84,1084,1083,107,108,1108,1107] nip = 8 material = e_steel
19 nodes = [63,64,1064,1063,83,84,1084,1083] nip = 8 material = e_steel
20 nodes = [84,85,1085,1084,108,109,1109,1108] nip = 8 material = e_steel
21 nodes = [85,86,1086,1085,109,110,1110,1109] nip = 8 material = e_steel
22 nodes = [65,66,1066,1065,85,86,1086,1085] nip = 8 material = e_steel
23 nodes = [45,46,1046,1045,65,66,1066,1065] nip = 8 material = e_steel
24 nodes = [27,28,1028,1027,45,46,1046,1045] nip = 8 material = e_steel
25 nodes = [7,8,1008,1007,27,28,1028,1027] nip = 8 material = e_steel
26 nodes = [9,10,1010,1009,29,30,1030,1029] nip = 8 material = e_steel
27 nodes = [29,30,1030,1029,47,48,1048,1047] nip = 8 material = e_steel
28 nodes = [47,48,1048,1047,67,68,1068,1067] nip = 8 material = e_steel
29 nodes = [67,68,1068,1067,87,88,1088,1087] nip = 8 material = e_steel
30 nodes = [87,88,1088,1087,111,112,1112,1111] nip = 8 material = e_steel
31 nodes = [88,89,1089,1088,112,113,1113,1112] nip = 8 material = e_steel
32 nodes = [89,90,1090,1089,113,114,1114,1113] nip = 8 material = e_steel
33 nodes = [69,70,1070,1069,89,90,1090,1089] nip = 8 material = e_steel
34 nodes = [49,50,1050,1049,69,70,1070,1069] nip = 8 material = e_steel
35 nodes = [48,49,1049,1048,68,69,1069,1068] nip = 8 material = e_steel
36 nodes = [11,12,1012,1011,31,32,1032,1031] nip = 8 material = e_steel
37 nodes = [31,32,1032,1031,51,52,1052,1051] nip = 8 material = e_steel
38 nodes = [51,52,1052,1051,71,72,1072,1071] nip = 8 material = e_steel
39 nodes = [71,72,1072,1071,91,92,1092,1091] nip = 8 material = e_steel
40 nodes = [91,92,1092,1091,115,116,1116,1115] nip = 8 material = e_steel
41 nodes = [92,93,1093,1092,116,117,1117,1116] nip = 8 material = e_steel
42 nodes = [93,94,1094,1093,117,118,1118,1117] nip = 8 material = e_steel
43 nodes = [73,74,1074,1073,93,94,1094,1093] nip = 8 material = e_steel
44 nodes = [53,54,1054,1053,73,74,1074,1073] nip = 8 material = e_steel
45 nodes = [52,53,1053,1052,72,73,1073,1072] nip = 8 material = e_steel
46 nodes = [33,34,1034,1033,53,54,1054,1053] nip = 8 material = e_steel
47 nodes = [13,14,1014,1013,33,34,1034,1033] nip = 8 material = e_steel
48 nodes = [15,16,1016,1015,35,36,1036,1035] nip = 8 material = e_steel
49 nodes = [16,17,1017,1016,36,37,1037,1036] nip = 8 material = e_steel
50 nodes = [17,18,1018,1017,37,38,1038,1037] nip = 8 material = e_steel
51 nodes = [35,36,1036,1035,55,56,1056,1055] nip = 8 material = e_steel
52 nodes = [55,56,1056,1055,75,76,1076,1075] nip = 8 material = e_steel
53 nodes = [75,76,1076,1075,95,96,1096,1095] nip = 8 material = e_steel
54 nodes = [95,96,1096,1095,119,120,1120,1119] nip = 8 material = e_steel
55 nodes = [96,97,1097,1096,120,121,1121,1120] nip = 8 material = e_steel
56 nodes = [97,98,1098,1097,121,122,1122,1121] nip = 8 material = e_steel
57 nodes = [19,20,1020,1019,39,40,1040,1039] nip = 8 material = e_steel
58 nodes = [39,40,1040,1039,57,58,1058,1057] nip = 8 material = e_steel
59 nodes = [57,58,1058,1057,77,78,1078,1077] nip = 8 material = e_steel
60 nodes = [77,78,1078,1077,100,101,1101,1100] nip = 8 material = e_steel
61 nodes = [99,100,1100,1099,123,124,1124,1123] nip = 8 material = e_steel
62 nodes = [100,101,1101,1100,124,125,1125,1124] nip = 8 material = e_steel
63 nodes = [101,102,1102,1101,125,126,1126,1125] nip = 8 material = e_steel
# D is the diameter of the rod. In this case, it is in mm.
elements of type rod_2
64 nodes = [4,5] D = 5.0 material = e_steel
65 nodes = [8,9] D = 5.0 material = e_steel
66 nodes = [114,115] D = 5.0 material = e_steel
67 nodes = [14,15] D = 5.0 material = e_steel
68 nodes = [122,123] D = 5.0 material = e_steel
constraints of type boundary_condition
fixed Vx = 0 Vy = 0 Vz = 0 Ax = 0 Ay = 0 Az = 0
materials of type elastoplastic
e_steel rho = 0.0000078 E = 210 nu = 0.3 yield_stress = 0.180 EP = 0.1
# There is also an elastic material available as follows:
# material of type elastic
# e_steel rho = 0.0000078 E = 210 nu = 0.3
loads
end_load Fx = 0.1 Fy = 0 Fz = 0.1 Mx = 0 My = 0 Mz = 0
trackers of type rodforce
1 element = [64] filename = rod_force_64.trk
2 element = [65] filename = rod_force_65.trk
3 element = [66] filename = rod_force_66.trk
4 element = [67] filename = rod_force_67.trk
5 element = [68] filename = rod_force_68.trk
controls
run from 0 to 20 step 0.0001
print every 0.1 step
print tracker every 0.01 step