# 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