×

A comparison of different methods for calculating tangent-stiffness matrices in a massively parallel computational peridynamics code. (English) Zbl 1423.74953

Summary: In order to maintain the quadratic convergence properties of Newton’s method in quasi-static nonlinear analysis of solid structures it is crucial to obtain accurate, algorithmically consistent tangent-stiffness matrices. A goal of the study described in this paper was to establish the suitability of an underexplored method for numerical computation of tangent-stiffness operators, referred to as “complex-step”, and compare the method with other techniques for numerical derivative calculation: automatic differentiation, forward finite difference, and central finite difference. The complex-step method was implemented in a massively parallel computational peridynamics code for the purpose of this comparison. The methods were compared through profiling of the code for accuracy, speed, efficiency, and parallel scalability. This research provides data that can serve as practical guide for code developers and analysts faced with choosing which method best suits the needs of their application code. Additionally, motivated by the reproducible research movement, all of the code, examples, and workflow to regenerate the data and figures in this paper are provided as open source.

MSC:

74S30 Other numerical methods in solid mechanics (MSC2010)
65D25 Numerical differentiation
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Belytschko, T.; Moran, B.; Liu, W. K., Nonlinear finite element analysis for continua and structures, vol.1, (1999), Wiley
[2] T. Young, M. Mohlenkamp, Introduction to numerical methods and matlab programming for engineers, Ohio University Department of Mathematics, OH, 2009.
[3] Simo, J.; Hughes, T., Computational inelasticity, (1998), Springer New York · Zbl 0934.74003
[4] M. Parks, D. Littlewood, J. Mitchell, S. Silling, Peridigm users’ guide, SANDIA REPORT 2012-7800, Sandia National Laboratories (September 2012).
[5] Heroux, M.; Bartlett, R.; Howle, V.; Hoekstra, R.; Hu, J.; Kolda, T.; Lehoucq, R.; Long, K.; Pawlowski, R.; Phipps, E.; Salinger, A.; Thornquist, H.; Tuminaro, R.; Willenbring, J.; Williams, A.; Stanley, K., An overview of the trilinos project, ACM Trans. Math. Software, 31, 3, 397-423, (2005) · Zbl 1136.65354
[6] Chapra, S.; Canale, R., Numerical methods for engineers, (2010), McGraw-Hill, Inc NY
[7] Lyness, J.; Moler, C., Numerical differentiation of analytic functions, SIAM J. Numer. Anal., 4, 2, 202-210, (1967) · Zbl 0155.48003
[8] Lyness, J., Differentiation formulas for analytic functions, Math. Comp., 22, 102, 352-362, (1968) · Zbl 0159.45002
[9] Squire, W.; Trapp, G., Using complex variables to estimate derivatives of real functions, SIAM Rev., 110-112, (1998) · Zbl 0913.65014
[10] J. Newman, K. Anderson, D. Whitfield, 1998. Multi-disciplinary sensitivity derivatives using complex variables, computational fluid dynamics laboratory, Mississippi State University, ms, mssu-coe-erc-98-08.
[11] Anderson, W. K.; Newman, J. C.; Whitfield, D. L.; Nielsen, E. J., Sensitivity analysis for Navier-Stokes equations on unstructured meshes using complex variables, AIAA J., 39, 1, 56-63, (2001)
[12] Perez-Foguet, A.; Rodriguez-Ferran, A.; Huerta, A., Numerical differentiation for local and global tangent operators in computational plasticity, Comput. Methods Appl. Mech. Engrg., 189, 1, 277-296, (2000) · Zbl 0961.74078
[13] Perez-Foguet, A.; Rodriguez-Ferran, A.; Huerta, A., Numerical differentiation for non-trivial consistent tangent matrices: an application to the mrs-lade model, Internat. J. Numer. Methods Engrg., 48, 159-184, (2000) · Zbl 1005.74079
[14] Martins, J.; Kroo, I.; Alonso, J., An automated method for sensitivity analysis using complex variables, AIAA Paper, 689, 2000, (2000)
[15] Martins, J. R.; Sturdza, P.; Alonso, J. J., The complex-step derivative approximation, ACM Trans. Math. Software (TOMS), 29, 3, 245-262, (2003) · Zbl 1072.65027
[16] Lai, K.-L.; Crassidis, J., Extensions of the first and second complex-step derivative approximations, J. Comput. Appl. Math., 219, 1, 276-293, (2008) · Zbl 1145.65016
[17] Al-Mohy, A. H.; Higham, N. J., The complex step approximation to the Fréchet derivative of a matrix function, Numer. Algorithm, 53, 1, 133-148, (2010) · Zbl 1188.65054
[18] Jin, W.; Dennis, B. H.; Wang, B. P., Improved sensitivity analysis using a complex variable semi-analytical method, Struct. Multidiscip. Optim., 41, 3, 433-439, (2010)
[19] Voorhees, A.; Millwater, H.; Bagley, R., Complex variable methods for shape sensitivity of finite element models, Finite Elem. Anal. Des., 47, 10, 1146-1156, (2011)
[20] H.R. Millwater, C. Quintana, J. Garza, D. Wagner, Application of the complex variable finite element method, ZFEM, to structural mechanics sensitivity analysis.
[21] D. Systèmes, Abaqus v6. 12 Documentation-ABAQUS analysis user’s manual, ABAQUS Inc 6.
[22] Nishida, E.; Foster, J.; Briseno, P., Constant-strain-rate testing of a G10 laminate composite through optimized kolsky bar pulse shaping techniques, J. Compos. Mater., 47, 23, 2895-2903, (2013)
[23] Abreu, R.; Stich, D.; Morales, J., On the generalization of the complex step method, J. Comput. Appl. Math., 241, 84-102, (2013) · Zbl 1258.65027
[24] Phipps, E.; Gay, D.; Bartlett, R., Sacado: automatic differentiation tools for C++ codes, tech. rep. SAND2009-7540C, (2009), Sandia National Laboratories: Computer Science Research Institute, presented at Trilinos User’s Group Meeting 2009 in Albuquerue, New Mexico
[25] Heroux, M.; Willenbring, J., Trilinos users guide, tech. rep. SAND2003-2952, (2003), Sandia National Laboratories
[26] Gay, D., Semiautomatic differentiation for efficient gradient computations, tech. rep. SAND2004-4688P, (2004), Sandia National Laboratories: Optimization and Uncertainty Estimation Department
[27] Griewank, A., On automatic differentiation, tech. rep. CRPC-TR89003, (1989), Rice University: Center for Research on Parallel Computation
[28] Morton, K. W., Numerical solution of partial differential equations: an introduction, (2005), Cambridge University Press · Zbl 1126.65077
[29] Hughes, T. J.; Pister, K. S., Consistent linearization in mechanics of solids and structures, Comput. Struct., 8, 3, 391-397, (1978) · Zbl 0377.73046
[30] Hughes, T. J.; Taylor, R. L., Unconditionally stable algorithms for quasi-static elasto/visco-plastic finite element analysis, Comput. Struct., 8, 2, 169-173, (1978) · Zbl 0365.73029
[31] Team, S. S.M., Adagio 4.18 user’s guide, tech. rep. SAND2010-6313, (2010), Sandia National Laboratories
[32] M. Rezaiee Pajand, J. Alamatian, The dynamic relaxation method using new formulation for fictitious mass and damping, Struct. Eng. Mech., 34. · Zbl 1257.74173
[33] Silling, S., Reformulation of elasticity theory for discontinuities and long-range forces, J. Mech. Phys. Solids, 48, 1, 175-209, (2000) · Zbl 0970.74030
[34] Silling, S.; Epton, M.; Weckner, O.; Xu, J.; Askari, E., Peridynamic states and constitutive modeling, J. Elasticity, 88, 151-184, (2007) · Zbl 1120.74003
[35] Silling, S.; Lehoucq, R., Peridynamic theory of solid mechanics, Adv. Appl. Mech., 44, 73-168, (2010)
[36] Littlewood, D. J., Simulation of dynamic fracture using peridynamics, finite element modeling, and contact, ASME Conf. Proc., 44465, 209-217, (2010)
[37] Ha, Y.; Bobaru, F., Studies of dynamic crack propagation and crack branching with peridynamics, Int. J. Fracture, 162, 1-2, (2010), 229-224 · Zbl 1425.74416
[38] Bobaru, F.; Duangpanya, M., A peridynamic formulation for transient heat conduction in bodies with evolving discontinuities, J. Comput. Phys., 231, 7, 2764-2785, (2012) · Zbl 1253.80002
[39] Katiyar, A.; Foster, J.; Ouchi, H.; Sharma, M., A peridynamic formulation of pressure driven convective fluid transport in porous media, J. Comput. Phys., 261, 209-229, (2014) · Zbl 1349.76819
[40] Seleson, P.; Parks, M. L.; Gunzburger, M.; Lehoucq, R. B., Peridynamics as an upscaling of molecular dynamics, Multiscale Model. Simul., 8, 1, 204-227, (2009) · Zbl 1375.82073
[41] Seleson, P.; Parks, M. L.; Gunzburger, M., Peridynamic state-based models and the embedded-atom model, Commun. Comput. Phys., 15, 1, 179-205, (2014) · Zbl 1373.82029
[42] Rahman, R.; Foster, J.; Haque, A., A multiscale modeling scheme based on peridynamic theory, Int. J. Multiscale Comput. Eng., 12, 3, 223-248, (2014)
[43] Team, C. D., Cubit 14.0 user documentation, tech. rep., (2013), Sandia National Laboratories
[44] Cormen, T. H.; Leiserson, C. E.; Rivest, R. L.; Stein, C., Introduction to algorithms, (2009), The MIT Press · Zbl 1187.68679
[45] K. Davis, Effective use if the intel compilers offload features. (Sep. 2013). URL http://software.intel.com/en-us/articles/effective-use-of-the-intel-compilers-offload-features.
[46] C. Bischof, A. Carle, P. Khademi, A. Mauer, P. Hovland, Adifor 2.0 user’s guide, Mathematics and Computer Science Division Report ANL/MCS, Argonne National Laboratory, Argonne, Illinois.
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.