1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670 671 672 673 674 675 676 677 678 679 680 681 682 683 684 685 686 687 688 689 690 691 692 693 694 695 696 697 698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738 739 740 741 742 743 744 745 746 747 748 749 750 751 752 753 754 755 756 757 758 759 760 761 762 763 764 765 766 767 768 769 770 771 772 773 774 775 776 777 778 779 780 781 782 783 784 785 786 787 788 789 790 791 792 793 794 795 796 797 798 799 800 801 802 803 804 805 806 807 808 809 810 811 812 813 814 815 816 817 818 819 820 821 822 823 824 825 826 827 828 829 830 831 832 833 834 835 836 837 838 839 840 841 842 843 844 845 846 847 848 849 850 851 852 853 854 855 856 857 858 859 860 861 862 863 864 865 866 867 868 869 870 871 872 873 874 875 876 877 878 879 880 881 882 883 884 885 886 887 888 889 890 891 892 893 894 895 896 897 898 899 900 901 902 903 904 905 906 907 908 909 910 911 912 913 914 915 916 917 918 919 920 921 922 923 924 925 926 927 928 929 930 931 932 933 934 935 936 937 938 939 940 941 942 943 944 945 946 947 948 949 950 951 952 953 954 955 956 957 958 959 960 961 962 963 964 965 966 967 968 969 970 971 972 973 974 975 976 977 978 979 980 981 982 983 984 985 986 987 988 989 990 991 992 993 994 995 996 997 998 999 1000 1001 1002 1003 1004 1005 1006 1007 1008 1009 1010 1011 1012 1013 1014 1015 1016 1017 1018 1019 1020 1021 1022 1023 1024 1025 1026 1027 1028 1029 1030 1031 1032 1033 1034 1035 1036 1037 1038 1039 1040 1041 1042 1043 1044 1045 1046 1047 1048 1049 1050 1051 1052 1053 1054 1055 1056 1057 1058 1059 1060 1061 1062 1063 1064 1065 1066 1067 1068 1069 1070 1071 1072 1073 1074 1075 1076 1077 1078 1079 1080 1081 1082 1083 1084 1085 1086 1087 1088 1089 1090 1091 1092 1093 1094 1095 1096 1097 1098 1099 1100 1101 1102 1103 1104 1105 1106 1107 1108 1109 1110 1111 1112 1113 1114 1115 1116 1117 1118 1119 1120 1121 1122 1123 1124 1125 1126 1127 1128 1129 1130 1131 1132 1133 1134 1135 1136 1137 1138 1139 1140 1141 1142 1143 1144 1145 1146 1147 1148 1149 1150 1151 1152 1153 1154 1155 1156 1157 1158 1159 1160 1161 1162 1163 1164 1165 1166 1167 1168 1169 1170 1171 1172 1173 1174 1175 1176 1177 1178 1179 1180 1181 1182 1183 1184 1185 1186 1187 1188 1189 1190 1191 1192 1193 1194 1195 1196 1197 1198 1199 1200 1201 1202 1203 1204 1205 1206 1207 1208 1209 1210 1211 1212 1213 1214 1215 1216 1217 1218 1219 1220 1221 1222 1223 1224 1225 1226 1227 1228 1229 1230 1231 1232 1233 1234 1235 1236 1237 1238 1239 1240 1241 1242 1243 1244 1245 1246 1247 1248 1249 1250 1251 1252 1253 1254 1255 1256 1257 1258 1259 1260 1261 1262 1263 1264 1265 1266 1267 1268 1269 1270 1271 1272 1273 1274 1275 1276 1277 1278 1279 1280 1281 1282 1283 1284 1285 1286 1287 1288 1289 1290 1291 1292 1293 1294 1295 1296 1297 1298 1299 1300 1301 1302 1303 1304 1305 1306 1307 1308 1309 1310 1311 1312 1313 1314 1315 1316 1317 1318 1319 1320 1321 1322 1323 1324 1325 1326 1327 1328 1329 1330 1331 1332 1333 1334 1335 1336 1337 1338 1339 1340 1341 1342 1343 1344 1345 1346 1347 1348 1349 1350 1351 1352 1353 1354 1355 1356 1357 1358 1359 1360 1361 1362 1363 1364 1365 1366 1367 1368 1369 1370 1371 1372 1373 1374 1375 1376 1377 1378 1379 1380 1381 1382 1383 1384 1385 1386 1387 1388 1389 1390 1391 1392 1393 1394 1395 1396 1397 1398 1399 1400 1401 1402 1403 1404 1405 1406 1407 1408 1409 1410 1411 1412 1413 1414 1415 1416 1417 1418 1419 1420 1421 1422 1423 1424 1425 1426 1427 1428 1429 1430 1431 1432 1433 1434 1435 1436 1437 1438 1439 1440 1441 1442 1443 1444 1445 1446 1447 1448 1449 1450 1451 1452 1453 1454 1455 1456 1457 1458 1459 1460 1461 1462 1463 1464 1465 1466 1467 1468 1469 1470 1471 1472 1473 1474 1475 1476 1477 1478 1479 1480 1481 1482 1483 1484 1485 1486 1487 1488 1489 1490 1491 1492 1493 1494 1495 1496 1497 1498 1499 1500 1501 1502 1503 1504 1505 1506 1507 1508 1509 1510 1511 1512 1513 1514 1515 1516 1517 1518 1519 1520 1521 1522 1523 1524 1525 1526 1527 1528 1529 1530 1531 1532 1533 1534 1535 1536 1537 1538 1539 1540 1541 1542 1543 1544 1545 1546 1547 1548 1549 1550 1551 1552 1553 1554 1555 1556 1557 1558 1559 1560 1561 1562 1563 1564 1565 1566 1567 1568 1569 1570 1571 1572 1573 1574 1575 1576 1577 1578 1579 1580 1581 1582 1583 1584 1585 1586 1587 1588 1589 1590 1591 1592 1593 1594 1595 1596 1597 1598 1599 1600 1601 1602 1603 1604 1605 1606 1607 1608 1609 1610 1611 1612 1613 1614 1615 1616 1617 1618 1619 1620 1621 1622 1623 1624 1625 1626 1627 1628 1629 1630 1631 1632 1633 1634 1635 1636 1637 1638 1639 1640 1641 1642 1643 1644 1645 1646 1647 1648 1649 1650 1651 1652 1653 1654 1655 1656 1657 1658 1659 1660 1661 1662 1663 1664 1665 1666 1667 1668 1669 1670 1671 1672 1673 1674 1675 1676 1677 1678 1679 1680 1681 1682 1683 1684 1685 1686 1687 1688 1689 1690 1691 1692 1693 1694 1695 1696 1697 1698 1699 1700 1701 1702 1703 1704 1705 1706 1707 1708 1709 1710 1711 1712 1713 1714 1715 1716 1717 1718 1719 1720 1721 1722 1723 1724 1725 1726 1727 1728 1729 1730 1731 1732 1733 1734 1735 1736 1737 1738 1739 1740 1741 1742 1743 1744 1745 1746 1747 1748 1749 1750 1751 1752 1753 1754 1755 1756 1757 1758 1759 1760 1761 1762 1763 1764 1765 1766 1767 1768 1769 1770 1771 1772 1773 1774 1775 1776 1777 1778 1779 1780 1781 1782 1783 1784 1785 1786 1787 1788 1789 1790 1791 1792 1793 1794 1795 1796 1797 1798 1799 1800 1801 1802 1803 1804 1805 1806 1807 1808 1809 1810 1811 1812 1813 1814 1815 1816 1817 1818 1819 1820 1821 1822 1823 1824 1825 1826 1827 1828 1829 1830 1831 1832 1833 1834 1835 1836 1837 1838 1839 1840 1841 1842 1843 1844 1845 1846 1847 1848 1849 1850 1851 1852 1853 1854 1855 1856 1857 1858 1859 1860 1861 1862 1863 1864 1865 1866 1867 1868 1869 1870 1871 1872 1873 1874 1875 1876 1877 1878 1879 1880 1881 1882 1883 1884 1885 1886 1887 1888 1889 1890 1891 1892 1893 1894 1895 1896 1897 1898 1899 1900 1901 1902 1903 1904 1905 1906 1907 1908 1909 1910 1911 1912 1913 1914 1915 1916 1917 1918 1919 1920 1921 1922 1923 1924 1925 1926 1927 1928 1929 1930 1931 1932 1933 1934 1935 1936 1937 1938 1939 1940 1941 1942 1943 1944 1945 1946 1947 1948 1949 1950 1951 1952 1953 1954 1955 1956 1957 1958 1959 1960 1961 1962 1963 1964 1965 1966 1967 1968 1969 1970 1971 1972 1973 1974 1975 1976 1977 1978 1979 1980 1981 1982 1983 1984 1985 1986 1987 1988 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 2016 2017 2018 2019 2020 2021 2022 2023 2024 2025 2026 2027 2028 2029 2030 2031 2032 2033 2034 2035 2036 2037 2038 2039 2040 2041 2042 2043 2044 2045 2046 2047 2048 2049 2050 2051 2052 2053 2054 2055 2056 2057 2058 2059 2060 2061 2062 2063 2064 2065 2066 2067 2068 2069 2070 2071 2072 2073 2074 2075 2076 2077 2078 2079 2080 2081 2082 2083 2084 2085 2086 2087 2088 2089 2090 2091 2092 2093 2094 2095 2096 2097 2098 2099 2100 2101 2102 2103 2104 2105 2106 2107 2108 2109 2110 2111 2112 2113 2114 2115 2116 2117 2118 2119 2120 2121 2122 2123 2124 2125 2126 2127 2128 2129 2130 2131 2132 2133 2134 2135 2136 2137 2138 2139 2140 2141 2142 2143 2144 2145 2146 2147 2148 2149 2150 2151 2152 2153 2154 2155 2156 2157 2158 2159 2160 2161 2162 2163 2164 2165 2166 2167 2168 2169 2170 2171 2172 2173 2174 2175 2176 2177 2178 2179 2180 2181 2182 2183 2184 2185 2186 2187 2188 2189 2190 2191 2192 2193 2194 2195 2196 2197 2198 2199 2200 2201 2202 2203 2204 2205 2206 2207 2208 2209 2210 2211 2212 2213 2214 2215 2216 2217 2218 2219 2220 2221 2222 2223 2224 2225 2226 2227 2228 2229 2230 2231 2232 2233 2234 2235 2236 2237 2238 2239 2240 2241 2242 2243 2244 2245 2246 2247 2248 2249 2250 2251 2252 2253 2254 2255 2256 2257 2258 2259 2260 2261 2262 2263 2264 2265 2266 2267 2268 2269 2270 2271 2272 2273 2274 2275 2276 2277 2278 2279 2280 2281 2282 2283 2284 2285 2286 2287 2288 2289 2290 2291 2292 2293 2294 2295 2296 2297 2298 2299 2300 2301 2302 2303 2304 2305 2306 2307 2308 2309 2310 2311 2312 2313 2314 2315 2316 2317 2318 2319 2320 2321 2322 2323 2324 2325 2326 2327 2328 2329 2330 2331 2332 2333 2334 2335 2336 2337 2338 2339 2340 2341 2342 2343 2344 2345 2346 2347 2348 2349 2350 2351 2352 2353 2354 2355 2356 2357 2358 2359 2360 2361 2362 2363 2364 2365 2366 2367 2368 2369 2370 2371 2372 2373 2374 2375 2376 2377 2378 2379 2380 2381 2382 2383 2384 2385 2386 2387 2388 2389 2390 2391 2392 2393 2394 2395 2396 2397 2398 2399 2400 2401 2402 2403 2404 2405 2406 2407 2408 2409 2410 2411 2412 2413 2414 2415 2416 2417 2418 2419 2420 2421 2422 2423 2424 2425 2426 2427 2428 2429 2430 2431 2432 2433 2434 2435 2436 2437 2438 2439 2440 2441 2442 2443 2444 2445 2446 2447 2448 2449 2450 2451 2452 2453 2454 2455 2456 2457 2458 2459 2460 2461 2462 2463 2464 2465 2466 2467 2468 2469 2470 2471 2472 2473 2474 2475 2476 2477 2478 2479 2480 2481 2482 2483 2484 2485 2486 2487 2488 2489 2490 2491 2492 2493 2494 2495 2496 2497 2498 2499 2500 2501 2502 2503 2504 2505 2506 2507 2508 2509 2510 2511 2512 2513 2514 2515 2516 2517 2518 2519 2520 2521 2522 2523 2524 2525 2526 2527 2528 2529 2530 2531 2532 2533 2534 2535 2536 2537 2538 2539 2540 2541 2542 2543 2544 2545 2546 2547 2548 2549 2550 2551 2552 2553 2554 2555 2556 2557 2558 2559 2560 2561 2562 2563 2564 2565 2566 2567 2568 2569 2570 2571 2572 2573 2574 2575 2576 2577 2578 2579 2580 2581 2582 2583 2584 2585 2586 2587 2588 2589 2590 2591 2592 2593 2594 2595 2596 2597 2598 2599 2600 2601 2602 2603 2604 2605 2606 2607 2608 2609 2610 2611 2612 2613 2614 2615 2616 2617 2618 2619 2620 2621 2622 2623 2624 2625 2626 2627 2628 2629 2630 2631 2632 2633 2634 2635 2636 2637 2638 2639 2640 2641 2642 2643 2644 2645 2646 2647 2648 2649 2650 2651 2652 2653 2654 2655 2656 2657 2658 2659 2660 2661 2662 2663 2664 2665 2666 2667 2668 2669 2670 2671 2672 2673 2674 2675 2676 2677 2678 2679 2680 2681 2682 2683 2684 2685 2686 2687 2688 2689 2690 2691 2692 2693 2694 2695 2696 2697 2698 2699 2700 2701 2702 2703 2704 2705 2706 2707 2708 2709 2710 2711 2712 2713 2714 2715 2716 2717 2718 2719 2720 2721 2722 2723 2724 2725 2726 2727 2728 2729 2730 2731 2732 2733 2734 2735

========
LOGBOOK Version 2.11
========
09 Nov 1998  removed mistake from anisotropic.f
 faster code for orthotropic materials in stress.f and
elemmatr_pr.f (replacement of orthotropic.f); to be
implemented in elemmatr_la.f too.
10 Nov 1998  definition of sets within *NODE and *ELEMENT.
11 Nov 1998  several steps within one calculation.
 no allocation of mast1 and mast2 for static calculations.
16 Nov 1998  small correction to orientations.
 correction to material transformation in kyrinput.
17 Nov 1998  small corrections for multiple step calculations.
23 Nov 1998  correction for second order frequency calculations.
 resetting the loads for perturbative steps
24 Nov 1998  check for second order static calculations OK.
26 Nov 1998  frequency calculations with prestress works! At
buckling load the frequency reduces to zero.
10 Dec 1998  corrected an error in kyrinput (call of nquickx2).
Jan 1999  implemented LANSO for eigenvalue calculations:
did not work; abandoned.
 implemented ARPACK for eigenvalue calculations:
problems under Linux, which do not occur on the
IRIX. TO BE ANALYZED!
 implemented UMFPACK for static calculations: for
small problems half as fast as PROFILE; for big
problems too much memory grabbed; abandoned.
 implemented SPOOLES for static calculations: works
fine; for big problems and nested dissection ordering
about 25% faster than PROFILE (18,000 equations).
Probably even faster for really big problems and low
fillin (> 50,000 equations and < 5 % fillin).
3 Feb 1999  on a SGI, SPOOLES is about 7 times faster than the
profile method for about 110,000 degrees of freedom
(900 MB RAM needed).
7 Feb 1999  started reorganization: (*) instead of (1), implicit
none instead of implicit real*8(ah,oz)
Mar 1999  ARPACK works on Linux when compiled with FC=f77f2c; reason
unclear; however, it is not necessarily faster than
LEVALH. Maybe calculating the largest eigenvalues of the
inverse matrix is faster?
27 Mar 1999  MODAL DYNAMIC works. Testing of buckling with ARPACK started.
8 April 1999  calculating frequencies in shiftinvert mode with ARPACK/
SPOOLES works and is significantly faster than the direct
mode (factor 4 or 5).
 BUCKLING works!
18 April 1999  checked shiftinvert mode for big problems. Up to a factor
of 90 faster than the direct method!
 refined output selection: U, NT, S or E.
21 April 1999  united the routine composing the eigenvalue matrices
with the routine constructing the static stiffness matrix
(elemmatr_ei.f and elemmatr_st.f > c_3d20.f)
 automation of the test routine (compare)
20 May 1999  reorganisation of the element topology information
 inclusion of C3D8 elements.
7 juni 1999  in previous releases 2nd order static calculations used
the large displacment stiffness only, 2nd order modal
calculations used the stress stiffness only. Now both
stiffnesses are used in any of these cases. What is
left to do is the use of the differential large deformation
matrix for buckling calculations (independent of iperturb),
in addition to the stress stiffness differential used now.
6 juli 1999  correction of small errors
7 juli 1999  removed lumping in modal analysis. Results match ABAQUS
results very well (deviation < 0.5 %).
10 juli 1999  started removing lumping from linear dynamic calculations.
eigenmodes should be normalized with respect to the mass
matrix. To do!
11 juli 1999  normalization with respect to the mass matrix done. Still
differences with ABAQUS for beta (damping coef) != 0. Maybe
there is still an error for the supercritical regime. To
check
 removed levalh and subprograms: default and only solver for
modal analysis is ARPACK
31 juli 1999  took care of inhomogeneous equations: introduce a dummy
node with fixed displacements.
Sept. 1999  started to work on geometrically nonlinear calculations
26 Oct. 1999  geometrically nonlinear option implemented and checked.
 in second order frequency calculations the material
temperatures from the previous nonperturbative step are
kept.
2 Dec. 1999  changed the extrapolated stress and strain values in the
midpoints of the 20node elements to the ABAQUS convention,
i.e. extrapolation is made within each element after which
the mean is taken of values in neighboring elements.
Dec. 1999  started the implementation of hyperelastic solids,
hyperfoam materials and deformation plasticity.
25 Dec. 1999  removed an error in mafillsm and mafillpr: introduced
a new variable idof3
JanFeb 2000  check of the hyperelastic capabilities; merges disp and
stress into one subroutine results to cover nonlinear
calculations with prescribed displacements.
7 Feb 2000  corrected an error in materialdata for mixed isotropic/
orthotropic/anisotropic materials
20 Feb 2000  accelerated hyperelastic calculations by a factor of 15
5 Mar 2000  deformation plasticity works (mechanical forces)!
9 Mar 2000  deformation plasticity subject to thermal forces works!
30 April 2000  C3D20R seems to work very well for thick and thin shells
(just one layer of brick elements)
 upgraded to SPOOLES.2.2.
19 May 2000  started implementing implicit and explicit integration
dynamics
June 2000  implementation nonlinear dynamics finished;
needs checking
 material data are determined in each integration point
separately instead of per element
3 July 2000  corrected an error occurring in nonlinear calculations
with mpc's
July 2000  changed the stressstrain curve to be given for deformation
plasticity from 2nd PiolaKirchhoff stress vs. Lagrangian
strain to true stress vs. Eulerian strain.
18 July 2000  included extrapolation of one increment to the next:
accelerates convergence
5 Aug 2000  started implementing incremental plasticity with
isotropic and/or kinematic hardening
Sept 2000  started checking implicit and explicit dynamics
18 Sept 2000  corrected the allocation of jq: 3*nk+1 (needed in
mastruct)
 implicit and explicit dynamics seem to work!
However, the calculation of the initial acceleration due
to a sudden load change at the start of a step, requiring
the solution of a linear equation set with the full mass
matrix causes problems: the solution of the system is not
plausible (1.e24). The problem improves drastically by
adding some stiffness to the mass.
1 Oct 2000  the only environment variable left is the jobname
all other file names are derived from the jobname:
.inp,.dat,.frd,.eig and .sta
 changed CALCULIX into CalculiX
 stored increment information in the .sta file
15 Oct 2000  started to replace the FORTRAN input files by C input
files; purpose: automatic reallocation if the problem
data size exceed preset values
16 Oct 2000  postponed the replacement of FORTRAN by C to a later date;
6 Nov 2000  isotropic incremental plasticity works for beam
example!
16 Nov 2000  implemented interpolation between hardening curves at
different temperatures
9 Dec 2000  accelerated nonlinear calculations for given
displacements or temperatures: first iteration in
a new increment must be purely linear elastic
23 Dec 2000  put version 0.9 on our webpage
28 Dec 2000  started to implement creep; user routines option provided
for plastic hardening curves and the creep law.
17 Jan 2001  creep (viscoplasticity) works
25 Jan 2001  looked into iterative solvers (ITPACKnspcg and
pcgsolverTUMuenchen)
4 Feb 2001  abandoned iterative solvers: worked for small
examples (1,000 DOFs) but not for intermediate ones
(20,000 DOFs)
5 Feb 2001  started to implement cyclic symmetry
6 Feb 2001  fixed format is not longer supported. Use free format.
16 Feb 2001  invested some more time in pcgsolver. Works well
with Cholesky preconditioning. Convergence slow for
2D problems such as plates. For large, compact 3D problems
faster than SPOOLES. Pcgsolver can be selected for
linear problems (test phase).
27 Feb 2001  started to implement C3D8R, C3D20 and C3D10 elements.
10 Mar 2001  changed the Nested Dissection Ordering in SPOOLES to
MultiSection Ordering. Seems to be faster for medium
size problems (150,000 DOF).
 implemention of C3D8R, C3D20 and C3D10 elements
successfully finished.
26 Mar 2001  introduced preprocessor directives to cope with different
C  FORTRAN interface conventions on different machines
(underscore or not)
27 Mar 2001  changes input and output format from record length 80 to
record length 132
 started to code transformations (rectangular and
cylindrical) (*TRANSFORM).
21 Apr 2001  introduced an iterative procedure to improve the
results in a *BUCKLE procedure: the buckling factor
is linked to the value of sigma (ARPACK)
 transformations work
 started to implement the iterative procedure for
nonlinear calculations
22 Apr 2001  iterative procedure for nonlinear calculations works
30 Apr 2001  collapsed elements work now! (changed the add*f subroutines)
 corrected some storage inconsistensies in CalculiX.c; should
solve the problems encountered when several steps occur in
the same input deck.
2 May 2001  corrected an error leading to wrong results for
element types with weight different from 1 and anisotropic
material behavior (subroutines orthonl.f and anisonl.f:
forgot to multiply with weight)
14 May 2001  changed the Sloan renumbering routines to allow skylines
larger than 4byte integers (>2147483647).
17 May 2001  started to change mastruct from FORTRAN to C in order to
be able to reallocate the # of nonzero's in the matrix
21 May 2001  automatic reallocation of nonzero's works
26 May 2001  started automatic allocation of input data
(subroutine allocation.f)
5 June 2001  auomatic allocation works;
 started implementation of cyclic symmetry conditions
for frequency calculations
23 June 2001  corrected an error in the residual stress ( sign)
19 July 2001  started to work on tetrahedral meshing of point
clouds
9 Aug 2001  tetrahedral meshing of a box enclosing a point
cloud works
14 Sep 2001  corrected an error in file incplas.f:
"c9=c6*umb*3.d0" replaced by "c9=c6*3.d0"
18 Sep 2001  started change from updated Lagrangian to total
Lagrangian for incremental plasticity
7 Oct 2001  concluded change to total Lagrangian formulation
 started introduction of field xstiff and xstate
in preparation for umat routine
25 Oct 2001  wrote first user material subroutine
29 Oct 2001  umat (material) subroutine seems to work
24 Nov 2001  included stress at start of increment in umat
 new investigation of cyclic symmetry: frequencies
are correct, eigenmodes too, except at the
dependent boundary
25 Nov 2001  cyclic symmetry seems to work!
29 Nov 2001  correction for MPC's with only one term (= zero
SPC)
9 Dec 2001  mapping of the results for one sector to other sectors
in a cyclic symmetry calculation works.
7 Jan 2002  started coding rigid body motion
 worked on tensiononly umat routine
8 Jan 2002  corrected an error in the C translation in mastructcs.c
11 Jan 2002  corrected an integer internal overflow in arpackcs.c
30 Jan 2002  finished coding rigid body motion
 cyclic symmetry in conjunction with other MPC's works
7 Feb 2002  fixed some minor errors related to *FREQUENCY and
*BUCKLING
March 2002  started to work on the theory manual.
2 Apr 2002  corrected an error in frdcyc: the coordinates have
to be duplicated also for kode>1
10 Apr 2002  corrected an error in boundaries: deletion of SPCs
in local coordinate systems with OP=NEW did not work
properly
16 Apr 2002  continued to work on rigid body motion and nonlinear
mpc's in general
27 Apr 2002  *RIGID BODY works
 with RF only EXTERNAL forces are obtained. If a MPC
connects two nodes of the structure, the force
is internal and cannot be obtained using RF.
5 June 2002  allowed for nested *INCLUDE statements up to three
levels deep
12 June 2002  took out the normalization of the displacements in
arpackbu.c (buckling).
6 July 2002  changed getnewline: the input is freed from blanks
and changed to upper case. Thus, input is more
flexible.
17 July 2002  same name can be used for node sets and element sets
internally, a "N" or "E" is appended to distinguish
them
18 July 2002  sets can be defined using previously defined sets
 abbreviating u_calloc with NNEW
24 July 2002  a set can be defined using multiple *NSET or *ELSET
cards
 abbreviating realloc with RENEW
29 July 2002  corrected a bug in boundaries.f which occurred when
SPC's in transformed coordinates were removed with
OP=NEW
6 Sept 2002  corrected wrong file names in some of the *WARNING
and *ERROR messages
 initialized all allocation size variables in
CalculiX.c
16 Sept 2002  started to code umat_elastic_fiber.f to model fiber
reinforced hyperelastic materials
26 Sept 2002  umat_elastic_fiber.f seems to work
28 Sept 2002  changed the names of the stress tensors
1 Oct 2002  started to code the energy calculation
3 Oct 2002  the energy calculation seems to work
 implemented the local orientation option in the
umat routine
5 Oct 2002  got rid of the environment variable; ccx is started
with i flag for the input file (without .inp)
9 Oct 2002  started to change cascade: solving for the dependent
DOFs in the MPCs using SPOOLES
29 Oct 2002  finished Cversion of cascade. The decascading is
performed by calling SPOOLES to solve the
nonsymmetric system of equations.
3 Nov 2002  execution of renumber.f is removed for SPOOLES
4 Nov 2002  started to treat nonlinear MPC's more generally:
variable number of MPC terms between iterations
must be taken into account
8 Nov 2002  corrected an error in cycsymmods.f
16 Nov 2002  started to work on STRAIGHT, PLANE and MEAN ROTATION
nonlinear MPCs
21 Nov 2002  started to work on plane strain, plane stress,
axisymmetric, shell and beam elements
Dec 2002  started to work on single crystal plasticity
15 Jan 2003  STRAIGHT, PLANE and MEAN ROTATION MPCs seem to work
18 Jan 2003  *AMPLITUDE can also be used for linear static
calculations
 linear MPCs can have more than 9 terms
26 Jan 2003  updating the User's Manual and the test example set
 finishing the tests on plane strain, plane stress,
axisymmetric, shell and beam elements
 the use of SPOOLES in cascade leads to significant
longer run times. Original method is restored and
translated into C. SPOOLES maybe useful for MPC's
with a large radius of influence such as
incompressibility.
 fixed the connection between solid elements and 1D
or 2D elements.
8 Feb 2003  the internal state variables can be stored in the
.dat and .frd file
 improved the boundary conditions for plane strain,
plane stress and axisymmetric elements.
9 Feb 2003  single crystal umat seems to work
(caveat: for ithermal=0 the field eth(1..6) is not
defined)
13 Feb 2003  corrected some mistakes in cycsymmods.f
27 Feb 2003  changed the meaning of OP=NEW
 changed the effect of output options in multistep
analyses
01 Mar 2003  corrected an error in the energy calculation and
incremental plasticity: calculation of irreversible
quantities must start from the values at the start
of the increment and not rely on intermediate values
04 Mar 2003  solved the problem with beammr (definition of neq
in CalculiX.c for icascade!=0)
02 Apr 2003  started to work on wedges and 4node tets
 corrected the implementation of the alpha method
(starting acceleration zero); convergence is
accelerated
09 Apr 2003  changed integration scheme for the mass matrix in
dynamic calculations with discontinuous forces
10 Apr 2003  simplified shape functions for 20node hexa elements
13 Apr 2003  changed extrapolation coefficients for C3D20
 replaced alp=.2215 by alp=.2917 for explicit dynamic
calculations with 20node elements (e_c3d.f)
 important change in ikmpc and ikboun: 6 DOFs are
assigned to each node instead of 3 DOFs
27 Apr 2003  introduced for frequency calculations a stiffness
contribution due to centrifugal forces
01 May 2003  started working on stiffness contribution of
distributed surface loads
06 May 2003  stiffness contribution of distributed surface
loads seems to work
10 May 2003  introduced rotational DOFs for shells and beams
26 May 2003  added the parameter TIME=TOTAL TIME to the
*AMPLITUDE keyword card
04 June 2003  introduced bending moment and torque for 1d and 2d
structures
 allowed for seven DOFs in ikboun, ikforc and ikmpc:
DOF zero is reserved for the temperature
 tet4, wedge6 and wedge15 seem to work
08 June 2003  MPC's can contain rotational degrees of freedom now
11 June 2003  started to code thermal calculations
14 June 2003  changed the first dimension of v,vold,vini,veold,
accold,veini,accini,fn,nactdof,vt and fnt to four in
order to accommodate temperature and thermal flux
23 June 2003  introduced the logical parameters mass, stiffness,
buckling and rhs to decide which entities to build.
25 June 2003  changed xload(i) to xload(1..2,i), similarly for
nelemload and iamload to accommodate both convection/
radiation coefficients and environmental temperatures
for heat transfer calculations
17 Juli 2003  started forced convection and cavity radiation
31 Juli 2003  introduced the keyword card *CONTROLS to control
convergence and allow for linear calculations with
1d and 2d elements.
2 Aug 2003  allowed for linear calculations with 1D and 2D
elements (beams, shells..) by setting the
convergence criteria to 1.d+30.
3 Aug 2003  allowed for linearization of *MPC constraints and
*RIGID BODY constraints by adapting the convergence
criteria to 1.d+30.
31 Aug 2003  change integration point numbering for C3D15
(conform to ABAQUS)
5 Oct 2003  reintroduced damping in direct dynamic calculations;
improves performance without deteriorating the
quality of the results
7 Oct 2003  started to work on a restart file
29 Oct 2003  corrected a bug in the calculation of the
distributed load stiffness
 finished the restart capability
29 Nov 2003  reorganization of the output in the .dat file:
output is grouped per set (node set/element set)
 implementation of whole element output: ELSE and
EVOL, TOTALS=YES and TOTALS=ONLY
03 Dec 2003  read and write files for a RESTART are now
different: extension .rin for RESTART,READ and
.rout for RESTART,WRITE
 there is no default any more for *EL FILE and
*NODE FILE
09 Dec 2003  corrected a mistake in restartshort.f
23 Dec 2003  finished the theory manual!
26 Dec 2003  maxlenmpc is stored in the restart file
 removed a possible division through zero in
arpackcs.c (cyclic symmetry)
5 Jan 2004  introduced a new field typeboun to classify the type
of boundary conditions. Only the BC's defined by
*BOUNDARY should be deleted for OP=NEW.
types: R=rigidbody, P=planempc, S=straightmpc
M=midplane, U=usermpc and B=boundary.
6 Jan 2004  prescribing the displacements of all DOFs is now
possible (neq=0)
12 Jan 2004  improved the convergence of STRAIGHT and PLANE MPC's
14 Jan 2004  changed the syntax of
*INITIAL CONDITIONS,TYPE=STRESS: the residual stress
tensor must be given in each integration point,
not just one tensor per element
 corrected an error in allocation.f
19 Jan 2004  made a correction for dynamic calculations in
allocation.f
 made a correction for restart calculations:
additional definitions of amplitudes and sets is
allowed.
20 Jan 2004  started to check thermal analysis
21 Jan 2004  made a correction in incplas.f for thermal
viscoplastic calculations (J_mech)
24 Jan 2004  corrected an error in cycsymmods.f
26 Jan 2004  coded the output of heat flux and heat generation
2 Feb 2004  introduced nenerold to take energy requests into
account in frequency steps with preload.
22 Feb 2004  started the user subroutines for heat flux, the film
coefficient and the emissivity
24 Feb 2004  the elastically anisotropic material model with von
Mises viscoplasticity is automatically selected as
soon as a *PLASTIC, *CYCLIC HARDENING or *CREEP card
is combined with a *ELASTIC,TYPE=ORTHO card.
26 Feb 2004  finished the user subroutines dflux.f, film.f,
radiate.f and dload.f
28 Feb 2004  introduced uniform and nonuniform bodygenerated
heat flux; seems to work.
2 Mar 2004  started changing the way the input is read: the
importance of the order of the cards is minimized
8 Mar 2004  made the step time and total time available in the
umat routines
9 Mar 2004  started to make changes to distinguish between the
mechanical and thermal part of the equation system
(e.g. neq is replace by neq[0] and neq[1]). This
is needed for instationary thermal calculations
13 Mar 2004  blank lines in the input file are disregarded
 new materials can be defined after a restart
 primary creep was included in the Norton creep law
(power of total time)
17 Mar 2004  started to work on reducing the effect of the order
of the keywords in the input deck
22 Mar 2004  changed the calculation of the strain in
perturbative frequency and buckling steps: the
strain is calculated about the deformed
configuration
25 Mar 2004  reduction of the effect of the order in the input
deck seems to work
30 Mar 2004  the input data for *FREQUENCY were changed: now, you
can restrict the eigenfrequencies to an interval by
specifying its lower and upper value.
31 Mar 2004  transient heat transfer calculations seem to work
21 April 2004  forced convection heat transfer works
23 April 2004  first calculations with cavity radiation
22 May 2004  first acoustic frequency calculations
 static step following heat transfer step works
10 June 2004  changed sideload from character*5 to character*20
allows userdefined name for userdefined loading
12 June 2004  started changes to allow for modal dynamic calculations
of the standard wave equation (e.g. in acoustics,
shallow waves etc.)
30 June 2004  forces can be summed and the sum printed in the
.dat file (TOTALS=YES or TOTALS=ONLY)
 modal dynamics of phenomena governed by the
Helmholtz equations seems to work
12 July 2004  changed the * format in internal reads for integers
into '(i40)'
22 July 2004  included the changes by Manfred Spraul enabling
multithreading with SPOOLES
3 Aug 2004  made the SPOOLES call modular (in order to easily
include TAUCS and other solvers).
8 Aug 2004  changed the * format in internal reads for reals
into '(f40)'
 included options to call TAUCS and the SGIsolver
(at compile time with DTAUCS and DSGI)
11 Aug 2004  started to work on the storage of results in local
coordinates
13 Aug 2004  storing results in local coordinates works.
(GLOBAL=YES or GLOBAL=NO after the *EL PRINT, *EL FILE
*NODE PRINT and *NODE FILE keyword cards).
4 Sept 2004  the value of jout is kept across the increments
unless a new value is defined
5 Sept 2004  removed an error: xload was sorted, xloadold was not
now xloadold is sorted as well (routine isortiddc)
8 Sept 2004  coupled temperaturedisplacement calculations work
14 Sept 2004  started to code 6noded triangular 2d elements
(they are expanded to 15node wedges)
 started to work on 1D and 2D elements for
thermal and thermomechanical calculations
16 Sept 2004  splitting of gen3delem.f in smaller subroutines started
17 Sept 2004  splitting of gen3delem.f finished
6 Oct 2004  6nodes 2d elements work
 1d and 2d elements for thermal calculations work
 started to work on axisymmetric elements for
cavity radiation
7 Oct 2004  cavity radiation for axisymmetric elements works
10 Oct 2004  in case of divergence the actual solution fields
and residual forces are stored in the frd file
14 Oct 2004  started an interface between the ABAQUS umat user
routine and CalculiX umat user routine.
19 Oct 2004  introduced mechanical strain to calculate the
energy density and for the ABAQUS umat user routine
25 Oct 2004  changed the syntax of *NODAL DAMPING to provide
compatibility with ABAQUS
 introduced the heat transfer elements DC3D4,DC3D6,
DC3D10,DC3D15 and DC3D20 for compatility with ABAQUS.
Internally, they are identical to C3D4, C3D6 etc.
30 Oct 2004  CalculiX checks length of set names, amplitude names..
to verify whether they do not exceed 20 characters
31 Oct 2004  removed alph from linel: everything is based on
eth now (thermal strain)
7 Nov 2004  started to work on maximum distance MPC
16 Nov 2004  removed an error in mastruct.c (loop should start from
0 instead of 1)
18 Nov 2004  introduced a field fmpc for the MPC force
24 Nov 2004  corrected an error in mastructcs.c (cf. 16 Nov).
25 Nov 2004  introduced a variable idiscon to mark a discontinuity;
if a discontinuity occurs the displacements at the
start of the next increment are not extrapolated
28 Nov 2004  introduced a new field irowsgi in routine sgi.c
8 Dec 2004  included tieset and ntie in the restart files
11 Dec 2004  started to code the gap MPC.
21 Dec 2004  corrected some small errors in dyna and dynsolv
23 Dec 2004  introduced sorted search for amplitudes (identamta.f)
4 Jan 2005  finalized the DIST and GAP MPC; made the GAP MPC
accessible through a GAPUNI element and *GAP card.
 adjusted the year in the copyright statement
5 Jan 2005  worked on dealing with ABAQUS umat routines for nonlinear
materials (umat_abaqusnl)
22 Jan 2005  started to implement the possibility to define several
volumetric forces within one structure
2 Feb 2005  removed an error in e_c3d.f (for lumping)
 allowed for smaller increments in dynamic explicit
calculations
 checked the size of force and displacement residuals
if too big, the increment size is reduced
13 Feb 2005  several volumetric forces within one structure work;
 generalized gravity works
18 Feb 2005  applied constraints to nodes on a cyclic symmetry
axis in static calculations
26 Feb 2005  material, orientation, amplitude and set names
are scheduled to be 80 characters long, textpart
is scheduled to be 132 characters long.
2 Mar 2005  finished extending names and textpart
3 Mar 2005  create user routines utemp and cflux; added field vold
in user routines dflux, film and radiate.
12 Mar 2005  introduced the parameter TIME DELAY to shift the time
within an amplitude
29 Mar 2005  started the inclusion of nonzero SPCs in modal dynamic
calculations (similar to *BASE MOTION in ABAQUS)
2 Apr 2005  changed dynsolv from FORTRAN to C
17 Apr 2005  inclusion of nonzero SPCs in modal dynamic calculations
works
19 Apr 2005  started to work on steady state dynamics (harmonic loading)
3 May 2005  steady state dynamics works (harmonic loading)
10 May 2005  started steady state dynamics for nonharmonic periodic
loading
16 May 2005  started *SENSITIVITY to determine eigenvalues of
geometrically slightly perturbed structures
30 May 2005  steady state dynamics for nonharmonic periodic loading
works
22 Juni 2005  started a more efficient storage of the boundary
stiffness coefficients (those stiffness coefficients
which correspond to SPC's; important for modal
dynamic calculations)
29 Juni 2005  changed application of force to axisymmetric elements:
now, the force is the one applied to the sector the angle
of which is defined on the *SOLID SECTION card, and not
the total force over 2*pi
 new storage works (fields neq and nzs have now a
length 3)
5 July 2005  created user subroutine massflowrate.f
6 July 2005  removed an error in CalculiX.c: if nam>0, iamload
must be sorted too in routine isortiddc.
10 July 2005  corrected an error in radflowload.f
 allowed for description in .frd file
16 July 2005  corrected an error in gen3dconnect.f
23 July 2005  simplified some code in results.f
 worked on 1d/2d output of 1d/2d elements
26 July 2005  worked on section forces
30 July 2005  1d/2d output and section forces work
 for axisymmetric elements *CLOAD, *CFLUX and
*MASS FLOW RATE are to be defined for the complete
circumference
31 July 2005  for 1d/2d elements NLGEOM is selected automatically
11 Aug 2005  corrected an error in map3dto1d2d.f
5 Sept 2005  corrected an error in solidsections (axisymmetric
elements)
7 Sept 2005  changed nonlingeo, radflowload and results: now the
gas temperatures are taken into account in the
convergence check (cam)
9 Sept 2005  same issue is on Sept 7: if gas temperatures are
calculated the "displacement" convergence check is
mandatory, no matter the size of the residual forces
13 Sept 2005  removed sensitivity.c and sensitivities.f (did not
bring any gain).
14 Sept 2005  changed renumber.f such that gaps in the node numbering
do not increase the execution time
21 Sept 2005  inserted temporarily the ZienckiewiczZhu error
estimator
23 Sept 2005  finished the coding of the ZienckiewiczZhu error
estimator for C3D20R elements
28 Sept 2005  put part of nonlingeo.c into checkconvergence.c
29 Sept 2005  put part of nonlingeo.c into calcresidual.c
rename residual.f into storeresidual.f
1 Oct 2005  extrapolation of previous results as start values
for the next increment is not done at the start
of a new step, else it is always done (nonlingeo.c)
2 Oct 2005  introduced user routine sigini.f for the
specification of initial stress fields
10 Oct 2005  started a major revision of the fluid elements for
forced convection purposes: now
they consist of three nodes, in the end nodes the
temperature is unknown, in the middle node the mass
flow rate is to be given (DOF 1) with *BOUNDARY
 revision of the convection equations in radflowload.f
to reach agreement with QTRAN
11 Oct 2005  finished forced convection changes
24 Oct 2005  started the generation of cyclic symmetry conditions
for dissimilar meshes
6 Nov 2005  cyclic symmetry conditions for dissimilar meshes work
14 Nov 2005  corrected an error in mastructcs.c (element types
C3D4, C3D6 and C3D15 were not taken into account)
16 Nov 2005  got rid of the NORENUMBER option: not needed any more
1 Dec 2005  corrected an error in e_c3d.f and linel.f
(initial shear stresses)
9 Dec 2005  corrected an error in some umat routines: +beta
should be beta
12 Dec 2005  started to work on modal dynamics and steady state
dynamics for cyclic symmetric structures
22 Jan 2006  adapted dyna.c (modal dynamics) for cyclic symmetry
calculations
30 Jan 2006  performed some changes to prepare for gas dynamics
calculations: instead of *SOLID SECTION, *GAS SECTION
should be used.
9 Feb 2006  corrected an error related to nam and nam_ in CalculiX.c
and calinput.f
27 Feb 2006  corrected an error in dyna.c and steadystate.c: t0,
t1old, t1 and iamt1 must be reallocated for cyclic
symmetric structures
28 Feb 2006  started the implementation of a thermal user material
5 March 2006  refined the cyclic symmetry conditions for
dissimilar meshes.
6 March 2006  starting a harmonized treatment of linear and
nonlinear materials.
7 March 2006  the profile solver is inactivated
12 March 2006  the harmonization of linear and nonlinear materials
is finished.
17 March 2006  allowing for more than one gravity load in an element
26 March 2006  heat conduction through the edges of a shell element
and the face of a plane stress element is possible
2 April 2006  started the extension of the forced convection
formulation into an aerodynamic network:
implementation of the orifice element
14 April 2006  introduced a Laplacetype method to find initial
pressures in an aerodynamic network
20 April 2006  started coding liquid networks
23 April 2006  calculations with liquid networks work
9 May 2006  allowed for imaginary gravity and centrifugal
loading (for steady state dynamics)
12 May 2006  speeded up all *ident*f files
25 May 2006  completed the discharge coefficient files for the
orifice, bleed tapping and preswirl nozzle
30 May 2006  change in tempload: gas nodes do not move during
the calculation
3 June 2006  started to code contact conditions
21 June 2006  corrected an error in mafillsm.f and
materialdata_th.f
11 Juli 2006  CYCLIC SYMMETRY and TYPE=NODE are required parameters
for the *TIE card and the *SURFACE card, respectively
26 Juli 2006  small corrections in radmatrix
 nodes with prescribed boundary conditions but not
belonging to elements are also stored in the frd file
30 Juli 2006  implemented time points at which output can be
requested
7 Sept 2006  one .onf file instead of several
 changed the order of the Gauss points for 4node
integration of faces of hexahedral elements
11 Sept 2006  speeded up cyclic symmetric thermal calculations
(took the symmetry of the integration points into
account in subroutine e_c3d_th.f)
27 Sept 2006  changed spooles.c to accommodate for nonsymmetric
systems of equations
28 Sept 2006  switched to spooles for the solution of the
nonsymmetric fluid network equation sets.
7 Oct 2006  switched back to dgesv (spooles is slower for small
systems)
 accelerated e_c3d_th.f
9 Oct 2006  incorporated some slight accelerations in e_c3d_th.f
14 Oct 2006  introduced a shape function routine specifically
for axisymmetric elements
15 Oct 2006  corrected an error in rubber.f
16 Oct 2006  accelerated axisymmetric calculations in results.f
17 Oct 2006  accelerated gen3dfrom2d.f
 contact with one element works
2 Nov 2006  introduced the parameter CAVITY on the *RADIATE
card
6 Nov 2006  accelerated axisymmetric calculations in mafillgas.f
and resultgas.f
 got rid of a segmentation fault in nonlingeo.c
7 Nov 2006  corrected a typing mistake for axisymmetric elements
in results.f
14 Nov 2006  changed the angle for axisymmetric elements to 2
degrees
17 Nov 2006  improved the connection with the user routine radiate
in radmatrix.f
23 Nov 2006  changed the argument list of attach.f
25 Nov 2006  changed 4node contact elements into 4 to 9node ones
3 Dec 2006  change did not improve convergence: 23 Nov state
reinstalled
 included the gas pipe element and a new convergence
strategy for fluids
10 Dec 2006  introduced *SURFACE INTERACTION and *SURFACE BEHAVIOR
13 Dec 2006  changed contact force into contact pressure (alpha is
multiplied with the area of the triangle)
17 Dec 2006  allowed for userdefined viewfactor inputdeck name
 use of shape20h_pl for plane stress and plane strain
elements
 extension for viewfactor file is .vwf
18 Dec 2006  only nonzero viewfactors are stored
 error estimator (Sascha Merz) included
 sorting of nelemload includes face number
21 Dec 2006  corrected an error in radflowload.c and nonlingeo.c
22 Dec 2006  made some changes in gen3dnor.f
26 Dec 2006  accelerated plane stress/strain and axisymmetric
calculations in e_c3d.f
6 Jan 2007  changed contact stiffness matrix slightly
7 Jan 2007  changed 4node contact elements into 5 to 9node ones:
second try (consistent stiffness matrix).
9 Jan 2007  5 to 9node contact elements seem to work.
14 Jan 2007  adapted surfaces.f in order to cover beam, shell and
2D elements
15 Jan 2007  corrected an error in e_c3d.f
16 Jan 2007  corrected an error in gaspipe.f
20 Jan 2007  started some changes to take 2D contact into account
(file gen3dsurf.f)
23 Jan 2007  corrected a mistake in the restart file
 treated Norton in the same way as the creep routine
(incplas.f)
31 Jan 2007  frequency calculations for 1D and 2D elements seem to
work.
9 Feb 2007  made a small change in frdcyc.c
22 Feb 2007  corrected an error in mafillsm.f
25 Feb 2007  made a change in map3dto1d2d.f, arpack.c and arpackcs.c
(rotational speed needed for stiffness matrix in
frequency calculations)
4 Mar 2007  introducted restrictor elements in the gas network
6 Mar 2007  made a change in arpack.c to allow for frequency
calculations after plastic calculations
 changes in the restrictor element files
11 Mar 2007  removed error occurring for empty node or element sets
20 Mar 2007  added parameter ADD to *CFLUX card
22 Mar 2007  updated umat_aniso_creep.f
24 Mar 2007  corrected an error in linear 1d/2d calculations:
created file fillrigidmpc.f
26 Mar 2007  removed vjj from subroutine incplas; modified the
influence of the thermal strain
3 April 2007  modified umat_aniso_creep.f
16 April 2007  replaced the rigid body formulation by a knot
formulation for 1d/2d elements: allows expansion
19 April 2007  corrected an error in allocation.f
20 April 2007  changed the size of field inpc to (long long) in order
to be able to read larger input files
21 April 2007  allowed for initial plastic strains
24 April 2007  made a correction in readinput.c
30 April 2007  started to make corrections for consecutive thermal
and mechanical calculations (especially for 1d/2d
elements)
2 May 2007  inclusion of branch fluid elements (start)
4 May 2007  made a correction in gen3dboun.f
15 May 2007  reduced the storage needs for the input deck:
introduced field ipoinpc to take variable length
lines into account
20 May 2007  updated the documention for fluid elements
of type branch
21 May 2007  allowed for linear rigid body calculations
23 May 2007  modified ithermal into a field to treat the BC's
for plane stress, plane strain and axisymmetric
elements in mixed *STATIC and *HEAT TRANSFER
calculations.
28 May 2007  created fully liquidstructure coupling for a one
dimensional flow in a flexible tube
30 May 2007  simplified materialdata_me.f
3 June 2007  started to split overall check of mechanical
and thermal convergence: qa,cam.. vectors of
length 2 instead of scalars
5 June 2007  corrected an error in couptempdisps.f
7 June 2007  simplified initialaccel.c
9 June 2007  simplified mafillsm.f
16 June 2007  transformed the logical mass into a vector: index 0
for mechanical calculations, index 1 for thermal
calculations
20 June 2007  started to implement spring elements
27 June 2007  linear and nonlinear springs seem to work
7 July 2007  started to code dashpots for linear dynamic calculations
10 July 2007  corrected a couple of small errors in gen3dnor.f.
26 July 2007  tried different differential equation solvers for
damped linear dynamics (dderkf, ddeabm and ddebdf)
30 July 2007  dashpots in combination with steady state dynamics
works
6 Aug 2007  new fluid elements: free and forced vortex, absolute
torelative elements and viceversa, Moehring elements
10 Sep 2007  closed the .inp, .dat and .sta file properly at the
end of the calculation
11 Sep 2007  changed ithermal into ithermal(2) in file gen3dnor.f
16 Sep 2007  introduced an error message in gen3dboun and gen3dforc
to cover the situation in which rotational loadings
or constrains are applied to a node in any but the
first step
24 Sep 2007  made a change in results.f (treatment of MPC's)
25 Sep 2007  made a change to nonlingeo.c and prediction.c:
mechanical results were not correctly saved during
a heat transfer step.
30 Sep 2007  resolved the problem from 16 Sep 2007: definition of
rotational BC's and moments in any but the first step
is now OK.
2 Oct 2007  small corrections to creeps.f, plastics.f, gen3dboun.f,
gen3dforc.f and checkconvergence.c
6 Oct 2007  made a change in steadystate.c: for cyclic symmetry
only ngraph sectors are stored
 allowed for a userdefined amplitude definition
10 Oct 2007  continued to work on frequencydependent dashpot
constants
15 Oct 2007  correction in the expansion of MPC's in expand.c
16 Oct 2007  replaced umat_aniso_creep.f and expand.c
for an anisotropic creep subroutine the creep
strain goes in and the stress comes out
20 Oct 2007  started to implement an uncoupled
temperaturedisplacement procedure
25 Oct 2007  changed the way the imaginary part of cloads and
dloads is stored
7 Nov 2007  improvements in frd output in case of divergence
 linear application of user defined loads for
static calculations
14 Nov 2007  included participation factors for steady state
dynamics calculations in the .dat file
19 Nov 2007  modified dyna.c in order to take contact
into account
21 Nov 2007  added the variable OUTPUT to *VIEWFACTOR,WRITE.
 P*NP can be used for *STATIC calculations as well
01 Dec 2007  started to combine contact with modal dynamic
calculations
15 Dec 2007  removed an error in rhs.f and mafillsm.f: force on
dependent node in MPC is redistributed among the
independent nodes
18 Dec 2007  thickness of plane strain elements is reduced if
there are axisymmetric elements in the model
20 Dec 2007  corrected an error in tempload
9 Jan 2008  started the implementation of computational fluid
dynamics
10 Jan 2008  changed the format for elements in the .dat file to
allow for element numbers exceeding 99999
 implemented engineering constants for orthotropic
elastic materials
14 Jan 2008  made a change in gen3dnor.f
22 Jan 2008  worked on the implementation of isochoric elements
(C3D20RH; was started years ago)
23 Jan 2008  made a correction in nonlingeo.c, mafillgas.f and
resultgas.f
30 Jan 2008  corrected the node order of 3node elements in frd.f
 replaced f20.10 by f20.0 in buckles.f, frequencies.f
and heattransfers.f
02 Feb 2008  finished a preliminary version of cfd (coding only);
05 Feb 2008  continued work on rotational rigid body motion with
linear dynamics
 continued work on incompressible elements
09 Feb 2008  expanded v,vold,vini,vr,vi and vt to contain 5
entries for the variables (one extra space for
static pressure)
10 Feb 2008  checking expansion; nearly OK, only segststate
still deviates
11 Feb 2008  expansion to 5 entries works
12 Feb 2008  the calculation of idof is changed to allow for 8 dofs:
temperature, 3 displacements or velocities, static
pressure, 3 angles or rotation (instead of 7: static
pressure is new)
 changed the internal dofs for rotations to 5,6 and 7.
dof 4 is for static pressure. The external dofs for
rotations remain 4,5 and 6, for static pressure it is 8.
For temperature the internal dof is 0, the external
dof is 11.
16 Feb 2008  calculated a default position tolerance in case none is
given by the user.
17 Feb 2008  removed an error in calinput.f and near2d.f
 continued checking of precfd.f
21 Feb 2008  removed an error in arpackcs.c
24 Feb 2008  debugged initialcfd.f
1 Mar 2008  started adaptive increment size in modal dynamic
calculations.
4 Mar 2008  modal dynamics with adaptive increment size seems to
run
5 Mar 2008  corrected an error in dyna.c
11 Mar 2008  corrected an error in selcycsymmods.f and mpcrem.f
 in modal dynamics and steady state dynamics calulations
any previous loading is removed
12 Mar 2008  isochoric elements seem to work
 .sta file is written for modal dynamics calculations
too
21 Mar 2008  included code for multistage MPC's
 continued computational fluid dynamics
27 Mar 2008  made some changes to multistages.f
 modified the header of user subroutine dflux.f
31 Mar 2008  continued debugging of 3D fluid dynamics
3 April 2008  modifications to dyna.c
10 April 2008  started an update of the gas network elements
15 April 2008  continued update of gas network elements
19 April 2008  finished update of gas network elements
 extended output for cyclic frequency calculations
(PHS,MAXU,MAXS)
23 April 2008  made some quality improvements in incplas.f
28 April 2008  continued computational fluid dynamics
ordered nelemface
 corrected an error in mafilltrhs.f
29 April 2008  update of dyna.c
7 May 2008  introduced contact damping and friction for
modal dynamic calculations
15 May 2008  different meaning of MAXU and MAXS in arpackcs.c
 correction in steadystate.c
27 May 2008  changed the order of the integration points in
gauss3d6
31 May 2008  coded a linearized version of the mean rotation MPC
3 June 2008  implemented pretension section to simulate
pretension in bolts (3D only).
5 June 2008  connected the linear equations solvers to compfluid.c
for liquids (pressure equation)
 improved error message for zero columns in the
stiffness matrix
12 June 2008  started to work on the improvement of the viewfactor
calculations
18 June 2008  corrections in expand.c and pretensionsections.c
21 June 2008  corrected an error in expand.c
 introduced the FIXED parameter on the *BOUNDARY card
23 June 2008  changed 1.d10 to 1.d5 as smallest allowable
coefficient in a MPC
28 June 2008  implemented STEADY STATE calculations for the
linear dynamics procedure
1 July 2008  coded an update of CYCLIC MPC's for nonlinear
calculations (nonlinmpc)
3 July 2008  corrected an error in envtemp.f
6 July 2008  made some changes in compfluid.c: stable time increment
is calculated in each iteration anew
9 July 2008  continued fluid dynamics: change of boundary conditions
has to be limited to guarantee stability.
13 July 2008  improved the time increment management for
computational fluid dynamics
19 July 2008  set the mechanical strain to zero at the end of a
frequency step
 started changes in the contact formulation: 1. adjust
in the first iteration of the first increment of the
first step; 2. zero force for too large distances
3. contact area division number
30 July 2008  made a change to surfaces.f
8 August 2008  corrected an error in printoutint.f
4 Sept 2008  coded 2nd order derivatives in shape3tri,shape6tri,
shape4q and shape8q
9 Sept 2008  supplemented subroutine springstiff with the derivatives
of xi and eta; contact convergence should improve
 corrected errors in nonlinmpc (CYCLIC) and bounadd
14 Sept 2008  made corrections in radiates.f, incplas.f and
usermaterials.f
20 Sept 2008  introduced CE for the equivalent creep strain
for anisotropic materials
24 Sept 2008  changed to another viewfactor calculation
 minor changes in radmatrix, e_c3d_th and cascade
11 Oct 2008  change of extrapolation of the mass flow to the
end nodes (fluidextrapolate.f)
 correction in umat_abaqus.f and umat_abaqusnl.f
 correction for section forces in map3dto1d2d.f
13 Oct 2008  the FREQUENCY and TIME POINTS parameters are
mutually exclusive for data storage keyword cards.
28 Oct 2008  modified the calculation of section forces to
accommodate for large deformations
7 Nov 2008  modified the output format for SDV in the .dat
and .frd file
9 Nov 2008  no correction of dependent node for nonmatching
cyclic symmetry conditions
 cyclic conditions are considered to be linear
(no update due to large deformation)
 changed interface for umatht and dload (for
lubrication problems)
13 Nov 2008  wrote user routines for coupled mechanical
lubrication problems
15 Nov 2008  made a change in springstiff: consistent tangent
matrix was not correct
17 Nov 2008  started to change dyna to take contact into account
30 Nov 2008  driven lid cavity (incompressibile fluid dynamics)
works!
 started the implementation of additional MPC's for
middle nodes belonging to dependent contact
surfaces
7 Dec 2008  additional MPC's for 3D and 2D elements seem to work
10 Dec 2008  changed PE into PEEQ and CE into CEEQ (print requests
or frd requests for the equivalent plastic/creep
strain)
 created printing output (.dat file) for gas networks:
static pressure (PS), total pressure (PT) and
mass flow (MF)
11 Dec 2008  created printing output (.dat file) for 3D fluid flow:
velocities (V) and static pressure (PT)
15 Dec 2008  no motion of dependent variables in cyclic
symmetric calculations (generatecycmpcs.f)
16 Dec 2008  the triangulation file for cyclic symmetry conditions
with dissimilar meshes takes the job name and
ending .tri
18 Dec 2008  corrected an error in shape8q (calculation of xsi) and
all other 2D shape functions
24 Jan 2009  got to work airfoil with near zero viscosity and zero
wall conditions (compressible fluid); smoothing is
not exactly the same as in Nithiarasu's program;
25 Jan 2009  got to work airfoil with Euler conditions (MPC's for
fluids)
2 Feb 2009  corrected an error in gen3dnor.f
4 Feb 2009  ZienkiewiczZhu error estimator now also works for
frequency, buckling, modal dynamics and steady
state dynamics calculations (stx instead of sti in
the call of out.f)
7 Feb 2009  coded shock viscosity smoothing
 included PARDISO as solver
10 Feb 2009  corrected an error in remastruct: mass[0]=1
21 Feb 2009  wrote a routine to calculate the lift and drag force
5 Mar 2009  calculation of the velocity through derivation
in dyna and dynacont
 replaced logical by integer in fluid calculations
7 Mar 2009  worked on turbulence (SST model)
8 Mar 2009  started work on pretension for 2D elements
10 Mar 2009  finished pretension for 2D elements
19 Mar 2009  introduced the variable pnewdt in umat_user.f
21 Mar 2009  coded direct modal damping
22 Mar 2009  continued to work on 3Dfluids and cyclic MPC's
4 April 2009  provided for BX, BY and BZ as volumetric loading
including nonuniform loading (BXNU...)
 included the velocity in the subroutine parameters
of dload.f
10 April 2009  started acceleration of dyna
19 April 2009  added some parameters to user routine dload
21 April 2009  coded direct damping for steady state dynamics
22 April 2009  PU and PHS are calculated for all sectors
requested
26 April 2009  corrected an error in incplas.f (residual stress)
2 May 2009  started tied contact
10 May 2009  tied contact works
 corrected an error in e_c3d_trhs.f
19 May 2009  started Mortar contact
15 June 2009  allowed for local coordinate systems for
fluid output
20 June 2009  changed the calculation of rf for nodes belonging
to MPC's and 1d or 2d elements
 made an update for gas networks (not finished)
 started debugging turbulence
23 June 2009  made a correction in results.f and pretensionsections.f
26 June 2009  added the cylic symmetry axis to the displacement
headers in the frd file
27 June 2009  provided two entries in field iperturb in order to
distinguish between material and geometrical
nonlinearities; nlgeom not automatically active for
*DYNAMIC and *COMBINED TEMPERATUREDISPLACEMENT
procedures
28 June 2009  started the multithreading version of compfluid
8 July 2009  update of mortar contact
10 July 2009  removed an error in radmatrix.f
 updated modal dynamic and contact
23 July 2009  accelerated the execution of dynacont and dfdbj
25 July 2009  change in nonlingeo: ielas=1 is only active in the first
iteration of the first increment if ilin=1 AND if the
step is not dynamic
3 August 2009  some small changes to Mortar contact
5 August 2009  added istep and iinc to the ABAQUS umat files
 user material does not automatically imply geometric
nonlinearity
11 August 2009  checked for SPC's in contact middle slave nodes in
routine gencontmpc.c
20 Sept 2009  updated Mortar contact
21 Sept 2009  made a correction in frequencies.f, buckles.f,
modaldynamics.f and steadystatedynamics.f
3 Oct 2009  accounted for 'NLGEOM=NO' for ABAQUS compatibility
 made some corrections in dynacont.c
11 Oct 2009  allowed for a selection of nodes for modal dynamics
calculations in order to speed up the calculations
12 Oct 2009  started to change mint_ into mi(2) in order to store
the maximum number of dofs per node
13 Oct 2009  continued (not finished yet)
24 Oct 2009  changed nactdof(1:3,*) into nactdof(1:mi(2),*)
25 Oct 2009  introduced the parameter TIME RESET on the *STATIC and
*HEAT TRANSFER keyword card
26 Oct 2009  in nonlinear calculations the nodes are listed where the
residual forces and variable changes are maximum
28 Oct 2009  corrected an error in springstiff.f
 inverted the order of working through the amplitudes:
last comes first (a redefined amplitude overrules
the earlier definitions)
31 Oct 2009  ran valgrind and started removing complaints by
valgrind
1 Nov 2009  started to use mi(2) for defining the size of field
v, vold... > smaller fields for thermal and
mechanical calculations
4 Nov 2009  mi(2) works for values >= 3.
 started introducing sets to limit output to the
frd file
10 Nov 2009  output to the frd file can be limited to a node set
14 Nov 2009  made substantial changes in dyna.c and dynacont.c
to speed up modal dynamic contact calculations
15 Nov 2009  xl2 and similar fields start now from 1, not 0.
22 Nov 2009  made substantial changes in dyna.c and subprograms
to speed up execution (force CHANGES are calculated)
28 Nov 2009  major speed up changes for modal dynamics finished
12 Dec 2009  continued speed optimization of modal dynamics
calculations
15 Dec 2009  implemented local system output for tensors in
frdtensor.f
9 Jan 2010  started coding liquid channels as 1D networks
12 Jan 2010  von Karman vortex street about a cylinder works
17 Jan 2010  coded linear pressureoverclosure
23 Jan 2010  changed once again output for MAXS in arpackcs
(quousque tandem abutere......)
 worked on friction in modal dynamics + contact
 started coding a different shock capturing
procedure
25 Jan 2010  made some changes in e_c3d_th.f (no displacement
information available for purely thermal calculations)
31 Jan 2010  expanded restrictors and branches to liquid
applications (1D networks).
4 Feb 2010  changed the way in which, for a given slave node, the
corresponding master triangle is found for spring contact
purposes
14 Feb 2010  made a change in results.f
21 Feb 2010  for contact springs the pressure should be multiplied by
the slave area; this major change was started today.
25 Feb 2010  changed the proximity check in linkdissimilar.f
7 Mar 2010  made corrections for ZZS output and ngraph>1
 coded slave area calculations for 2D contact applications
13 Mar 2010  coded the parameter ADJUST on the *CONTACT PAIR card
 got rid of INFACE SLIDING
15 Mar 2010  introduced local time stepping for steady state CFD
 cleaned up e_c3d_[v1,p,v2,t]rhs.f
 previous stable cfd version: 13 Mar 2010
22 Mar 2010  nodes with largest displacements etc. in 2d structures
are the original node numbers now, not the expanded
numbers
26 Mar 2010  corrected bodyforce.f
 defined the size of xstiff in dyna.c and steadystate.c
as (long long)
 corrected an error in frdvector.f
28 Mar 2010  local time stepping in steady state CFD calculations seems
to work
30 Mar 2010  made some changes in preparation for the Bishop 4node
shell element
11 Apr 2010  changed the dimension of field IL and IU in all sort
routines from 21 to 31.
 started to code friction for static contact calculations
(penalty)
18 Apr 2010  changed the definition of c0 for linear contact springs in
gencontelem.f
 changed the search for corresponding master triangles
in gentiedmpc to kneigh=1 and neighbor search
19 Apr 2010  plane stress, plane strain and axisymmetric elements
should be defined in the z=0 plane. This is now checked
in gen3delem.f
20 Apr 2010  started to compare frd files for selected test examples
(in addition to dat files)
25 Apr 2010  changed fields for networks in order to allow for
geometric data as unknowns
 default for *EL FILE is GLOBAL=YES
1 May 2010  for a gate valve the extent of closure can be solved for:
example for geometric unknowns in networks
17 May 2010  worked on liquid channels
 corrected an error in allocation.f
20 May 2010  made a correct in dyna.c (in case a dload subroutine has
been defined)
27 May 2010  improved the linear contact spring relationship
 worked on liquid channels
9 June 2010  modified the linear contact spring relationship slightly
14 June 2010  started debugging backwater calculations
16 June 2010  change the convergence rules for changing numbers of
contact spring elements (LARGE SLIDING)
30 June 2010  simple backwater calculations work (sluice gate/reservoir)
12 July 2010  weir calculations work
5 Aug 2010  made some small corrections
5 Sep 2010  2d results are calculated in each iteration so that they
are available in the user subroutines
9 Sep 2010  added a user subroutine for concentrated loads
10 Sep 2010  worked on gap conductance
11 Sep 2010  corrected an error in frdcyc (output of CONT)
13 Sep 2010  included the C3D8I and C3D8R formulation coded by
otto bernhardi
16 Sep 2010  gap conductance works (heat transfer through mechanical
contact)
23 Sep 2010  coded an acoustic material (pressure only)
 coded Coulomb friction for penalty contact
6 Oct 2010  changed MPC's for nodal forces on shells in
gen3dforc.f
9 Oct 2010  worked on contraction, enlargement, drop and step
for liquid channels
12 Oct 2010  contraction, enlargement, drop and step for liquid
channels seem to work
20 Oct 2010  corrected an error in dyna.c and steadystate.c
23 Oct 2010  changed dyna for cyclic symmetry: only the part of the
eigenvectors which is needed is expanded.
31 Oct 2010  implemented a slightly different scheme for compressible
CFDcalculations
5 Nov 2010  original CFDscheme proved to be superior
7 Nov 2010  started to work on the acceleration of steadystate
12 Nov 2010  finished acceleration of dyna and steadystate
15 Nov 2010  accelerated expand.c
 extended the arguments of cload.f
 coded MAXE as output variable (for frequency
calculations with cyclic symmetry)
20 Nov 2010  started acceleration of CFDcode
21 Nov 2010  corrected small errors in dyna.c, dynacontc. and dfdbj.c
13 Dec 2010  tried another variant of shock smoothing (exact
simulation of the second order derivative )
1 Jan 2011  calculating presgradient after applying the pressure
boundary conditions
3 Jan 2011  calculating the pressure switch as a truely second
derivative (normed such that max in field is 1)
5 Jan 2011  removed shocksmoothing based on local von Mises stress
20 Jan 2011  introduced basic smoothing of 0.1 for Euler calculations
 included Couette and Poiseuille examples
6 Feb 2011  changed applybounp.f (MPC's)
9 Feb 2011  started to work on MPC's for networks
12 Feb 2011  MPC's for networks work
15 Feb 2011  updated fridaforc.f
 added nonhomogeneous terms to network MPC's
28 Feb 2011  modified results (network output)
2 Mar 2011  started changes for contact:
 proportional reduction of penetration for
overclosures
 change friction (*CHANGE FRICTION)
 model change for contact pairs (*MODEL CHANGE)
8 Mar 2011  improved *PRETENSIONSECTION
 *FRICTION CHANGE works
 *CHANGE MODEL works
 proportional reduction of penetration for
overclosures works
9 Mar 2011  corrected an error in pretensionsections.f and all
restart routines (restartread.f, restartwrite.f and
skip.f)
14 Mar 2011  ordered nelemload as soon as loading is being read
from file
15 Mar 2011  the radiation matrix is recalculated only if the
emissivity or the viewfactors changed
24 Mar 2011  corrected an error in radmatrix.f
27 Mar 2011  changed storage of radiation matrix from full
storage into sparse storage
5 April 2011  wrote an interface to solve asymmetric matrices
with PARDISO
7 April 2011  parallelized the calculation of the viewfactor
matrix for radiation (cave: logicals in FORTRAN do
not seem to be threadsafe)
9 April 2011  got rid of iptri for radiation and introduced a
new entry in kontri > accelerates execution
19 April 2011  introduced initialchannel.f
20 April 2011  changed the allocation of xstiff to (long long)
3 May 2011  started to parallelize results.f
14 May 2011  parallelization of results.f finished
 changed a parameter in gencontelem.f
19 May 2011  in user subroutine cload fn is provided. To this
end the elements are determined in dyna.c for
which the stresses have to be determined.
24 May 2011  made subroutine frditeration.f available to write
the displacements in each iteration to file (for
debugging purposes)
7 June 2011  documented beam element examples
 started to work on complex frequencies
9 June 2011  corrected an error in mafillsmcs.f
 changed mi(2)into mi(*) in the integer declarations
13 June 2011  changed ielorien(*), ielmat(*) and thicke(2,*) into
ielorien(mi(3),*), ielmat(mi(3),*) and thicke(mi(3),*)
for composite applications
4 July 2011  made a small change in spooles.c (number of cpus
for parallel radiation applications)
5 July 2011  inserted a factor of 2. for Coriolis in e_c3d.f and
e_corio.f
6 July 2011  made some changes for complex frequency calculations
17 July 2011  change in .eigfile: for calculations without cyclic
symmetry the square of the eigenvalue (square of the
eigenfrequency) is stored
 Coriolis works
20 July 2011  took care of displacement and stress output for
composites
25 July 2011  started to change SPC and MPCtreatment for fluids
30 July 2011  introduced label PRF for *NODE FILE: size and angle
of reaction forces for cyclic symmetry frequency
calculations
 the square of the eigenfrequencies is stored in the
.eig file (the square can be negative  buckling)
1 Aug 2011  for cyclic symmetry calculations the imaginary part
of the strains is stored as well (if E is requested)
7 Aug 2011  introduced the mechanical strain ME as output variable
 corrected some mistakes in umat_abaqusnl.f
9 Aug 2011  corrected a mistake in CalculiX.c/allocation.f
(initialization of xstate/definition of nstate)
7 Sept 2011  corrected errors in mastructcs.c and mafillsmcs.f
(springs in connection with cyclic symmetry)
8 Sept 2011  created an option to store the displacements of the
last iterations (for debugging purposes in case of
divergence)
13 Sept 2011  introduced a check whether in mechanical calculations
with temperature loading a temperature was assigned to
all nodes belonging to elements
2 Oct 2011  introduced fields t0g and t1g to take temperature gradients
across shells into account
 continued automatic remeshing for quadratic elements in
contact
8 Oct 2011  corrected the calculation of the external faces and solid
boundary nodes in precfd.f
9 Oct 2011  change treatment of slave contact nodes in 2d plane stress/
strain and axisymmetric elements in gen3dsurf.f
11 oct 2011  reversed change from 25 July 2011 and 8 Oct 2011. Reason:
wrong results of convergent JefferyHamel flow
12 Oct 2011  improved automatic remeshing for C3D10 with contact
16 Oct 2011  introduced an error estimator (ERR) based on the extrapolation
within element belonging to one and the same node
25 Oct 2011  OUTPUT=2D or OUTPUT=3D is required in the presence of
shells or beams (or SECTION FORCES)
27 Oct 2011  accelerated CFD (stored dpress in e_c3d_prhs.f)
30 Oct 2011  reversed change of 25.10
 introduced the parameter CONTACT ELEMENTS on the
*NODE FILE and *EL FILE cards
31 Oct 2011  worked on turbulence models for CFD.
7 Nov 2011  in CFD calculations the solid surface and freestream
surface have to be declared as nodal sets (only
used for turbulence)
9 Nov 2011  started updating gas networks
13 Nov 2011  adjusted the time increment for CFD w.r.t. the
turbulent viscosity
 default for beams and shells is OUTPUT=3D, for plane
stress, plane strain and axisymmetric elements it is
OUTPUT=2D
14 Nov 2011  coded *CHANGE MATERIAL and *CHANGE PLASTIC in order to
change plastic material data within a step
 no creep in a *STATIC step (change of 52 into 51)
16 Nov 2011  changed the error estimator > standard deviation instead
of maxmin
24 Nov 2011  made a small change to shape8hu.f
 coded the interpolation of initial fields in integration
points for quadratic elements remeshed due to contact
 incorporated some changes for gas networks
26 Nov 2011  no creep in a *STATIC step. Use *VISCO instead
1 Dec 2011  replaced CCX_NPROC by OMP_NUM_THREADS
7 Dec 2011  introduced field emeini for umat user routines
10 Dec 2011  corrected some errors in remeshcontact.c and
remeshcontactel.f
17 Dec 2011  made some corrections for the storage of contact variables
in the .dat file for modal dynamic calculations
 changed k and ipoint into long long in matrixstorage.c
22 Dec 2011  corrections to modal dynamics + contact friction
 correction to multithreading
2 Jan 2012  worked on tensiononly materials
12 Jan 2012  changed calculation of initial temperatures in initialnet.f
16 Jan 2012  changed e11.4 into e13.6 in all printout*.f files
22 Jan 2012  started change in contact: consistent derivation of
the continuous normal
28 Jan 2012  corrected an error in checkconvergence.c and springstiff.f
5 Feb 2012  contact: changed continuous normal back to mathematically
correct normal
8 Feb 2012  corrected an error in nonlingeo.c (emeini)
19 Feb 2012  updated Mortar contact (without friction)
25 Feb 2012  changed the format of the .frd file into a mixture of
ASCII and binary information in case *NODE OUTPUT or
*ELEMENT OUTPUT is used.
4 Mar 2012  changed format for node and element number in the .dat
file to i10 (from i6)
10 Mar 2012  made some changes in createmddof.f and cycsymmods.f
12 Mar 2012  corrections in steadystate.c and rectcylvi.f
13 Mar 2012  N on the *CYCLIC SYMMETRY MODEL card is compared with the
angle based on the geometry.
 extrapolation of the network mass flow to the end nodes
is done in radflowload > can be used as boundary condition
in a subsequent thermomechanical structural calculation
(determination of the heat transfer coefficients)
18 Mar 2012  corrected an error in dyna.c (several steps with
restricted output > initial conditions in subsequent
step for ALL nodes)
23 Mar 2012  introduced a random error of 0.1 % in gen3dfrom1d.f in
order to get rid of hourglassing in beam elements.
26 Mar 2012  removed change from 23 March.
 started to code the expansion of linear shells, beams,
axisymmetric, plane stress and plane strain elements
10 Apr 2012  removed an error in nonlingeo.c
 additional option for calculating viewfactors only
19 Apr 2012  change in triangulate: now several faces belonging to one
and the same element may belong to the master surface
of a cyclic symmetry tie
22 Apr 2012  corrected an error in compfluid.c: instationary
incompressible flow works (Couette)
29 Apr 2012  introduced larger time step for the stability terms
for incompressible fluid calculations
1 May 2012  removed call to applybounv in between mafillv1rhs and
mafillv2rhs
5 May 2012  started to change the application of boundary conditions
for incompressible fluids
10 May 2012  made a change in e_c3d.f for contact springs:
important in the first iteration for given
displacements
11 May 2012  change in springstiff and springforc: Coulomb friction
is only taken into account for dynamic calculations
(linear and nonlinear) and geometrically nonlinear
static calculations
14 May 2012  changed applyboun.f to include boundary conditions for
the turbulent variables
23 May 2012  corrected small errors (for thermal calculations)
24 May 2012  accelerated some CFDroutines
29 May 2012  made a change in gen3dboun.f and flowresult.f
30 May 2012  updated Mortar contact
3 June 2012  CFD: MPC's are updated between step 1 and step 2:
necessary for cyclic symmetry calculations
14 June 2012  started to work on submodels
27 June 2012  continued work on submodels
5 July 2012  started work on automatic cyclic symmetry
conditions for CFD
 submodels seem to work
18 July 2012  step time is taken into account for linear calculations
(amplitude definitions)
 corrected an error in thermal contact
23 July 2012  changed the order of the last line in field ifacew1 and
ifacew2
28 July 2012  made some corrections in frd.c
 include modifications by Stefanie Asenbauer
4 Aug 2012  modified plastic routines (*VISCO)
 adjustment in e_c3d.f for submodes with nlgeom
7 Aug 2012  made a correction in remeshcontactel.f for elements with
more than one face belonging to a slave or master
contact surface
6 Sept 2012  corrected an error in remeshsurf.f
13 Sept 2012  inserted cgx reading routines for binary frd files
(submodel applications)
 allowed for contact calculations with linear strains
(NLGEOM is not automatically activated for contact)
 allowed for storage of the triangulation of the master
surfaces for cyclic symmetry, contact and tied contact
29 Sept 2012  corrected errors pointed out by valgrind while running
the test examples
7 Oct 2012  worked on mafillsmas (solution of asymmetric systems
of equations)
20 Oct 2012  changed c0=1.d6*dsqrt(area) into c0=1.d3*dsqrt(area)
in gencontelem.f
21 Oct 2012  worked on CFD + cyclic symmetry
22 Oct 2012  contact penetration is reduced proportionally to the
step time only in the first step (change in gencontelem.f)
28 Oct 2012  changed springstiff: introduced vini
3 Nov 2012  removed change from 28 Oct.
5 Nov 2012  changed cylindrical into cartesian coordinates in
interpolcycsymcfd.c
 changed 1.d5 into 1.d30 in bounadd.f
8 Nov 2012  continued to work on CFD and cyclic symmetry
12 Nov 2012  worked on CFD and cyclic symmetry
14 Nov 2012  added a "C" in sixth position for linear elements after
remeshing of quadratic elements at contact surfaces
 removed remeshing of quadratic elements at contact
surfaces in frdoutput
17 Nov 2012  made some important changes to springstiff.f (consistent
stiffness matrix: inclusion of normal force contribution
in the tangential stiffness)
25 Nov 2012  made some changes for the automatic generation of a
integer*8 version
29 Nov 2012  cleaning up files: systematic application of
implicit none; replacing numbers in FORTRAN subroutine
calls by variables (needed for the integer*4 >
integer*8 transition).
5 Dec 2012  checking master and slave facial surfaces for identical
faces (in allocont.f).
9 Dec 2012  started combined solution procedure (instead of
alternating) of gas and solid temperatures for generic
network elements
8 Jan 2013  worked on dashpot forces in steady state dynamics
calculations
17 Jan 2013  changed the defaults for penalty contact
19 Jan 2013  modification for axisymmetric and plane strain elements
in e_c3d_th.f, advecforc.f and advecstiff.f
28 Jan 2013  corrected the ASCII output from CProgramme (included
(float))
31 Jan 2013  changed map3dto1d2d for a correct interpolation of the
forces
 created networkinum.f
26 Feb 2013  for cyclic symmetric modal dynamics or steady state dynamics
calculations point forces in other sectors can be defined
in the global carthesian system or local cylindrical
system.
 modification in linkdissimilar for nonmultistage
calculations
 compressible laminar cfdcalculations work for selected
examples (linear Couette, Carterproblem)
28 Feb 2013  introduced the options REMOVE and REMOVE ALL on the
*EQUATION card
6 March 2013  started change for inclusion of userdefined pressure
overclosure relationship for contact
 change in restart (storage of plicon, plkcon ...)
9 March 2013  splitted updatecfd.f into updatecon.f and con2phys.f
13 March 2013  changed 4 into 14 in frd.c for PHAXX...
 corrections for tabular data underneath *SURFACE
BEHAVIOR
17 March 2013  introduced linear regression for CFDcalculations with
nonzero shock coefficient
21 March 2013  changed CLOAD,DLOAD,CFLUX,DFLUX: if redefined within same
step for same dof: added instead of replaced
 add INLET and OUTLET (network elements)
 FLUX no also works for nonCFD applications
24 March 2013  started to linearize contact for subsequent frequency
calculations (arpack.c)
25 March 2013  made some important changes to complexfreq.c
2 April 2013  made some corrections in complexfreq.c
 reduced CFD time scales to dtimef and dt(*)
8 April 2013  worked on arpack.c and arpackcs.c (taking into account
contact element stiffness; same as on March 24).
 selecting CONTACT ELEMENTS on *NODE FILE visualizes
the remeshing of the quadratic elements adjacent to
a contact surface
13 April 2013  removed an error in knotmpc.f and fillknotmpc.f (wrong
sign of field e)
 started differentation in dimensionality of knots
28 April 2013  removed an error for S6elements from gen3dforc.f
30 April 2013  .eigfile modified to accept complex eigenvalues
 complexfreq.c modified in order to store complex
eigenvalues and eigenvectors
5 May 2013  started to modify steadystate.c in order to handle
nonorthogonal eigenvectors and complex eigenfrequencies
6 May 2013  continued to work on complexfreq.c and steadystate.c
 pretension beam elements are removed after generating
the pretension MPC's
11 May 2013  coded shape functions for a 21node hexahedral element
(reduced from a 27node Lagrange hexahedral element)
and a 9node quadrilateral element; possible
applications: penalty contact
12 May 2013  changed .frd output requests V, S and HFL for
CFDcalculations into VF, SF and HFLF.
 changed the fluid frdfile ending into
jobnamef.frd
16 May 2013  started to make changes for the remeshing of C3D20
elements into C3D21..C3D26 elements at contact
18 May 2013  removed first change from 12. May 2013.
19 May 2013  changed names of (contact)spring and dashpot elements:
the number in the name is the number of master nodes
(i.e. one less than before; done to accommodate
9node master faces)
26 May 2013  coded 11noded up to 14nodes quadratic tetrahedral
elements (for contact)
27 May 2013  introduced parameter pnewdt in incplas.f
28 May 2013  introduced parameter pnewdt in umat_aniso_plas.f and
umat_aniso_creep.f
31 May 2013  made a change in shellsections.f (allowed for composite
shells with just one layer)
 started implementing 2x5x5 integration scheme for beams
(B32R)
4 June 2013  changed the 2x5x5 integration scheme from a Gauss
Simpson scheme into a GaussKronrod scheme; the latter
scheme contains the 2x2x2 ordinary Gauss scheme as
subset (only for rectangular crosssection beams)
 implemented weights 1/100 and 24/100 for 8node faces
in contact instead of 1/12 and 1/3
6 June 2013  corrected an error for the output of fluid quantities
20 June 2013  increased the maximum number of points on a plastic stress
strain curve from 20 to 200.
25 June 2013  switched back to a 2x2x2 integration scheme for
beams with a circular cross section
2 July 2013  introduced VF, PSF, TSF, MACH, TTF, PTF, SF and HFLF for
output of fluid calculations
3 July 2013  corrected an error in pretensionsections.f
 got fluid>structure interactions to work
7 July 2013  made some changes to steps.f, frequencies.f, arpack.c,
arpackcs.c and resultsmech.f in order to take previous
plasticity into account in a frequency calculation
22 July 2013  started to work on heat induction (electromagnetics).
29 July 2013  continued electromagnetics.
3 August 2013  corrected some errors in connection with frequency
calculations in combination with contact
8 Oct 2013  corrected an error in dyna.c
15 Oct 2013  started debugging inductive heating
26 Oct 2013  started properly deallocating all fields in CalculiX.c
30 Oct 2013  actualized frdread.c
31 Oct 2013  removed all finite element fluid subroutines; started
coding finite volume routines
4 Nov 2013  continued to work on electromagnetics
27 Nov 2013  started coding facetoface penalty contact
4 Dec 2013  worked on electromagnetics and facetoface penalty
16 Jan 2014  finalizing facetoface penalty contact
21 Jan 2014  worked on finite volume CFD (the CBSmethod was removed
from the code in November/December 2013: didn't work well
enough). The finite volume code is based on the book by
Ferziger and Peric.
27 Jan 2014  changed the calculation of inum for networksl
28 Jan 2014  made a change to readinput.c
 made some small changes to isort*f
29 Jan 2014  corrected an error in springstiff_f2f.f
30 Jan 2014  corrected two errors in advecforc.f
 continued to work on CFD
1 Feb 2014  got rid of field nnn (no separate renumbering any more;
renumbering is assumed to be done by the solver, e.g.
SPOOLES and PARDISO)
 coded output of the real and imaginary part (or,
alternatively the magnitude and phase) of
contact values in cyclic symmetry frequency calculations
4 Feb 2014  introduced islavsurfold and pslavsurfold
5 Feb 2014  introduced the parameter NO CHANGE on the *VIEWFACTOR
card (avoids the repeated LU decomposition of the
radiation matrix)
8 Feb 2014  optimized code of shape8h.f
 optimized code of springstiff_f2f.f
 temporarily removed the *dtime=0. statement in
checkconvergence.c: to be checked!
12 Feb 2014  updated rewriting for properties and facetoface
penalty contact entities
13 Feb 2014  changed meaning of ttime: now it is the total time at the
start of the step
16 Feb 2014  transferred info from areaslav(3..6,*) to pmastsurf(3..6,*)
18 Feb 2014  got rid of xstiff in subroutines dyna.c and
steadystate.c
23 Feb 2014  change inputinfo.f, inputwarning.f and inputerror.f
keyword card is included in message
24 Feb 2014  created one version for integer*4 and integer*8
introduced macros ITG and ITGFORMAT
8 March 2014  added a mechanism to check for negative pressures in
facetoface contact
29 March 2014  changed the definition of NNEW
8 April 2014  worked on beams with a pipe section
14 May 2014  finite volume CFD seems to work for incompressible
laminar flow without thermal effects; the pressure
boundary condition is given at one element center,
the convection scheme is purely upwind.
 beams with pipe section work
 removed some errors in multiple step shell calculations
and shell calculations with restart
15 May 2014  tried upwind difference with central difference as
deferred correction; is better than pure upwind
difference.
21 May 2014  text underneath *HEADING is stored at the top of the
frdfile
 for axisymmetric structures with output=2d the
mass flow is the integral over 360 deg.
 gamma differencing is ok, but not as good as upwind with
deferred central differencing
22 May 2014  for axisymmetric structures with output=2d the
force is the integral over 360 deg.
 removed an error in mafillsmcs.f and mafillsmcsas.f
16 June 2014  corrected an error in mafillsm.f
25 June 2014  created incplas_lin.f for geometrically linear
calculations (not linked yet)
30 June 2014  corrected an important error in incplas.f, incplas_lin.f,
umat_aniso_plas.f and umat_aniso_creep.f
10 July 2014  created a .cvg file containing convergence information
12 July 2014  corrected the contact energy for all but the
exponential pressureoverclosure relationship
14 July 2014  coded the participation factors and the relative effective
modal mass in arpack.c and arpackcs.c
 corrected an error in files springstiff*f and
springforc*f
15 July 2014  coded Rayleigh damping for implicit dynamics
 coded forces as boundary condition for submodels
26 July 2014  buoyancy seems to work
30 July 2014  wrote incplas_lin.f for small deformation plasticity
(no nlgeom).
 changed the convergence defaults for *dynamic
(more cutbacks allowed, smallest relative increment
= 1.e10)
4 Aug 2014  changed the reading order: moved *SURFACE on
position 6 before *TRANSFORM (because of
TRANSFORMF)
7 Aug 2014  removed an error in remastruct.c: fini and fextini should
not be set to zero in each iteration.
2 Sept 2014  removed an error in printoutnode
 removed usercload in dyna.c and createmdelem.f: fn will
be removed from cload.f (too time consuming)
4 Sept 2014  activated dashpots for *DYNAMIC calculations
 nmethod,nactdof,bcont and fn were removed from
subroutine cload.f
14 Sept 2014  switched arrays and columns of field vel for CFD
preparation to use dgmres
15 Sept 2014  switched to an GMRES iterative solver for the velocity and
temperature CFD equations
21 Sept 2014  introduced iaxial in all subroutines of flux.f: for the
networks the mass flow variable is the partial mass
flow (2 %), for the element equations, however, the
total mass flow is to be taken
11 Oct 2014  introduced openmp in routines op*f
 made some changes to the inductive heating routines
14 Oct 2014  cylindrical incompressible Couette flow works
15 Oct 2014  made a correction involving the variable iaxial
16 Oct 2014  made a correction to forcadd.f
22 Oct 2014  added again a small amount of stiffness to the mass
matrix to determine the initial acceleration (this
section was removed during a couple of weeks, however,
the stiffness as proven to be necessary)
23 Oct 2014  implemented SF,SVF and HFLF for fluids (.datfile and
.frdfile)
1 Nov 2014  incompressible laminar CFD with wedge elements seems to
work
4 Nov 2014  the coordinates transferred to the user subroutines in
tempload.f are the ones in the material frame
5 Nov 2014  made a change in incplas.f (dxitril)
9 Nov 2014  increased the number of pressure iterations for CFD from 1
to 2 (necessary for deformed meshes)
19 Nov 2014  removed an error related to time points using total
time accross several steps
 replaced stop by "stop 201" > can be used to detect
if the program stopped in an irregular way
20 Nov 2014  replaced "stop 201" by "call exit(201)"
25 Nov 2014  enhanced and corrected the mean rotation formulation
 did an important change to gen3delem/gen3dboun: rotations
are not expanded any more using expandable rigid bodies,
instead mean rotation MPC's are used
28 Nov 2014  made a change to bounadd.f
9 Dec 2014  electromagnetism without heat induction might work
 eliminated c3d20ri (isochoric) elements
11 Dec 2014  worked on heat transfer within electromagnetism
21 Dec 2014  worked on the documentation of electromagnetism
23 Dec 2014  took the stick as internal variable in
springforc_f2f and springstiff_f2f
8 Jan 2015  allowed for facial surfaces in cyclic symmetry definitions
21 Jan 2015  started a major change in facetoface penalty contact
28 Jan 2015  replaced mean rotation MPC's by expandable rigid bodies
for nonlinear dynamics: knots can take on arbitrary
rotational values, mean rotation MPC's are limited to
90 degrees
31 Jan 2015  changed the sign of the damping coefficient for
nodal diameters > number of sectors/2 in dyna.c
 continued changes for facetoface contact
(semismooth Newton)
4 Feb 2015  removed the change in facetoface penalty contact.
 removed an error in incplas.f.
8 Feb 2015  adapted *MODEL CHANGE to (de)activate elements
 code *CHANGE SURFACE INTERACTION
 made xstate and xstateini available in gencontelem_f2f
 with "contact elements" all contact elements are stored
in file contactelement.inp. They can be loaded into cgx
and visualized according to step, increment, attempt
and iteration
11 Feb 2015  changed interpolatestate.f
 reversed the first change from 31 Jan 2015
24 Feb 2015  introduced the possibility to take contact into account
in a linear elastic calculation (application: fracture
mechanics)
28 Feb 2015  introduced a line search for facetoface penalty.
3 Mar 2015  checked for loss of contact at the start of a new
increment in gencontelem_f2f
12 Mar 2015  deactivated nasym for tied pressureoverclosure
22 Mar 2015  coded membranes and trusses (gen3dnor.f)
3 April 2015  modified attach.f (limitation of the solutions in
the subdomains to the borders of the subdomain)
15 April 2015  introduced an aleatoric removal of contact elements
in case of a local minimum of facetoface contact
22 April 2015  replaced attach.f: smaller change in increment size
between the iterations
23 April 2015  used qa[2] to force a cut back in nonlingeo if the
linear network does not converge (cf. checkconvnet.c)
26 April 2015  started compressible laminar CFD with the method
of openCFD
20 May 2015  derived the contact traction for facetoface
penalty contact from the stresses
 improved the use of equations in networks
9 June 2015  worked on zero mass flow conditions in CFDapplications
 implemented the calculations of energies in dynamic
applications
25 June 2015  continued compressible laminar flow
2 July 2015  corrected an error in jouleheating.f
14 July 2015  changed gen3dboun: applied BC to middle nodes as
well (additional MPC)
 corrected use of edge loading EDNORx for shells
26 July 2015  changed +bp into bp in correctvfa.f
1 Sept 2015  the energy is now calculated for all (nonlinear)
dynamic calculations
3 Sept 2015  made a change in gen3dboun: fixing u(n_2)=u(n) for
translational SPC's
7 Sept 2015  introduced a field "flux" in CFD
13 Sept 2015  continued to work on friction heating
21 Sept 2015  introducted variable iit in spring*_f2f.f
26 Sept 2015  made a correction for the temperatures in a
pretension mechanical calculation
29 Sept 2015  allowed for several contact pairs underneath one and the
same *CONTACT PAIR card
10 Oct 2015  introduced an extra loop in order to obtain accurate
results for transient cfdcalculations
18 Oct 2015  changed mafillpcomp (a term in the change of density
was dropped => improved stability)
24 Oct 2015  for cutbacks with surfacetosurface penalty the
stiffnesses are reduced by 100 and increased to the
original values at convergence; the increment size
is not reduced for the first cutback.
25 Oct 2015  changed error estimator ERR into largest difference of
worst principal stress in the elements adjacent to the
node
28 Oct 2015  stable compressible fluid version; (subsonic,
transonic and supersonic bumps work)
4 Nov 2015  made a change in mafillsmmain and mafillsm1
5 Nov 2015  accelerated the execution of mastruct.c, mastructcs.c,
mastructem.c, mastructf.c
14 Nov 2015  introduced several codeaccelerations in mastruct.c
(especially for contact)
19 Nov 2015  posted version 2.9
20 Nov 2015  changed insert.c and mastruct.c
26 Nov 2015  accelerated attach.f and distattach.f
29 Nov 2015  supersonic flow works
30 Nov 2015  supersonic flow works (II), slightly different version
1 Dec 2015  parallellized the asymmetric stiffness contributions
5 Dec 2015  continued to work on compressible laminar flow
major improvement of stability
7 Dec 2015  no line search for dynamic contact
10 Dec 2015  modified explicit dynamics
19 Dec 2015  viscous compressible stationary flow seems to work
(finite volume method)
7 Jan 2016  changed ifaceq(9.. into ifaceq(8.. and
ifacet(7.. into ifacet(6.. (i.e. no 26node
hex or 14node tet anymore).
9 Jan 2016  removed an error in beamextscheme.f
13 Jan 2016  creep and viscoplasticity in combination with an
anisotropic material is addressed by material user
subroutines only.
 removed an error in compfluid for compressible
materials: correctvfa should never be called twice
without calling mafillpcomp in between (recalculation
of ap(i))
21 Jan 2016  removed the requirement for 1.d6 relative time
increments for implicit dynamics in combination with
contact (checkconvergence.c)
24 Jan 2016  introduced the gammamethod in mafillvcomp with
Central Difference replaced by 0.4*Upwind Difference +
0.6*Central Difference
27 Jan 2016  reset compressible CFD to UD (Upwind Difference)
 started to parallellize nearly every CFDroutine (openMP)
30 Jan 2016  worked on parallellization of mafillpcomp.f
31 Jan 2016  finished parallellization of compressible cfd
1 Feb 2016  worked on parallellization of incompressible cfd
2 Feb 2016  unified the calculation of material properties in
compfluid.c
3 Feb 2016  made some corrections to frictionheating
 parallellized mafillp.f
6 Feb 2016  changed in mafilltcomp from CD to UD (forgotten on
the 27th of Jan) => bumpsuperfineshallow results
are smooth again now
13 Feb 2016  worked on explicit dynamics (Carlo Monjaraz)
3 March 2016  updated implicit dynamics (Matteo Pacher)
 switched everything to mafillsmmain
6 March 2016  put version 2.10 in the internet
 started to use the negative range of nactdof to
point to SPC's and MPC's > avoids the extensive
use of nident
7 March 2016  corrected an error in getglobalresults.c and in
storecontactprop.f
17 March 2016  updated sensitivity.c and subroutines
 improved mafillem.f, mafilldm.f and mafillcorio.f
28 March 2016  started looking into periodicity for CFD calculations
30 March 2016  made a correction in initialnet.f
 started to clean up fluidconstants.f
 made a correction in mafillpcomp.f
5 April 2016  made some changes to calculate sensitivities
w.r.t. the nodes in the mesh
9 April 2016  continued to work on sensitivities
 changed the dimension from 100 to 60 for the element
matrices and vectors except for electromagnetics
13 April 2016  translational periodic fluid conditions seem to work
24 April 2016  worked on rotational periodic fluid conditions
28 April 2016  continued to work on sensitivities
10 May 2016  changed the convergence criteria for gas networks
11 May 2016  continued to work on the convergence criteria for
gas networks
17 May 2016  removed the forced units of rpm from the rotational
speed entries in the vortex and Moehring fluid type
elements
22 May 2016  continued to work on cyclic symmetry conditions for
CFD
16 June 2016  updated the User's Manual for gas network elements (theory)
 created .net  files with detailed gas network
information
21 June 2016  introduced box section general beam elements
 continued to improve gas networks
23 June 2016  removed the addition of 0.5 to fluid section properties
and changed int into nint
 worked on improving the initial conditions for compressible
networks
7 July 2016  corrected an error regarding *node print and *el print:
commands stay active across step boundaries
 dsload can be used to apply a pressure to a surface
24 July 2016  cyclic symmetry and compressible flow works
26 July 2016  made a correction in map3dto1d2d_v.f
 cyclic symmetry and incompressible flow works
28 July 2016  made some changes in gen3dfrom2d.f
 changed the wedge integration for the mass in dynamic
calculations from 18 to 9 integration points (some of the
18 integration poins lie on the border, which causes
problems in case of an expansion of a triangular axi
element lying on the rotation axis)
30 July 2016  copied v into vold and stx into sti after the 5th call
of results in nonlingeo.c
