Parametric deep energy approach for elasticity accounting for strain gradient effects

被引:142
作者
Vien Minh Nguyen-Thanh [1 ]
Anitescu, Cosmin [2 ]
Alajlan, Naif [3 ]
Rabczuk, Timon [3 ]
Zhuang, Xiaoying [1 ]
机构
[1] Leibniz Univ Hannover, Dept Math & Phys, Chair Computat Sci & Simulat Technol, Hannover, Germany
[2] Bauhaus Univ Weimar, Inst Struct Mech, D-99423 Weimar, Germany
[3] King Saud Univ, Coll Comp & Informat Sci, Dept Comp Engn, Riyadh, Saudi Arabia
关键词
Neural networks (NN); Physics-informed neural networks (PINNs); Partial differential equations (PDEs); Deep energy method; Elasticity; Strain gradient elasticity; INFORMED NEURAL-NETWORKS; LEARNING FRAMEWORK; SURFACE-ENERGY; ALGORITHM;
D O I
10.1016/j.cma.2021.114096
中图分类号
T [工业技术];
学科分类号
08 ;
摘要
In this work, we present a Parametric Deep Energy Method (P-DEM) for elasticity problems accounting for strain gradient effects. The approach is based on physics-informed neural networks (PINNs) for the solution of the underlying potential energy. Therefore, a cost function related to the potential energy is subsequently minimized P-DEM does not need any classical discretization and requires only a definition of the potential energy, which simplifies the implementation. Instead of training the model in the physical space, we define a parametric/reference space similar to isoparametric finite elements, which is in our example a unit square. The inputs are naturally normalized preventing the vanishing gradient problem and leading to much faster convergence compared to the original DEM. Forward-backward mapping is established by means of NURBS basis functions. Another advantage of this approach is that Gauss quadrature can be employed to approximate the total potential energy, which is the loss function calculated in the parametric domain. Backpropagation available in PyTorch with automatic differentiation is performed to calculate the gradients of the loss function with respect to the weights and biases. Once the network is trained, a numerical solution can be obtained in the reference domain and then is mapped back to the physical domain. The performance of the method is demonstrated through various numerical benchmark problems in elasticity and compared to analytical solutions. We also consider strain gradient elasticity, which poses challenges to conventional finite elements due to the requirement for C-1 continuity. (C) 2021 Elsevier B.V. All rights reserved.
引用
收藏
页数:27
相关论文
共 58 条
  • [1] Altan B., 1997, J. Mech. Behav. Mater., V8, P231, DOI DOI 10.1515/JMBM.1997.8.3.231
  • [2] An isogeometric analysis computational platform for material transport simulation in complex neurite networks
    Li A.
    Chai X.
    Yang G.
    Zhang Y.J.
    [J]. MCB Molecular and Cellular Biomechanics, 2019, 16 (02): : 123 - 140
  • [3] Artificial Neural Network Methods for the Solution of Second Order Boundary Value Problems
    Anitescu, Cosmin
    Atroshchenko, Elena
    Alajlan, Naif
    Rabczuk, Timon
    [J]. CMC-COMPUTERS MATERIALS & CONTINUA, 2019, 59 (01): : 345 - 359
  • [4] [Anonymous], 2009, Isogeometric Analysis: Towards Unification of Computer Aided Design and Finite Element Analysis
  • [5] [Anonymous], 1995, Bifurcation analysis in geomechanics
  • [6] [Anonymous], 1960, FINITE ELEMENT METHO
  • [7] Bathe K.J., 1996, Finite Element Procedures in Engineering Analysis
  • [8] ELEMENT-FREE GALERKIN METHODS
    BELYTSCHKO, T
    LU, YY
    GU, L
    [J]. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, 1994, 37 (02) : 229 - 256
  • [9] A unified deep artificial neural network approach to partial differential equations in complex geometries
    Berg, Jens
    Nystrom, Kaj
    [J]. NEUROCOMPUTING, 2018, 317 : 28 - 41
  • [10] Bergstra J, 2012, J MACH LEARN RES, V13, P281