Chawit Kaewnuratchadasorn, Jiaji Wang, Chul-Woo Kim, Xiaowei Deng
{"title":"Geometry physics neural operator solver for solid mechanics","authors":"Chawit Kaewnuratchadasorn, Jiaji Wang, Chul-Woo Kim, Xiaowei Deng","doi":"10.1111/mice.13405","DOIUrl":null,"url":null,"abstract":"This study developed Geometry Physics neural Operator (GPO), a novel solver framework to approximate the partial differential equation (PDE) solutions for solid mechanics problems with irregular geometry and achieved a significant speedup in simulation time compared to numerical solvers. GPO leverages a weak form of PDEs based on the principle of least work, incorporates geometry information, and imposes exact Dirichlet boundary conditions within the network architecture to attain accurate and efficient modeling. This study focuses on applying GPO to model the behaviors of complicated bodies without any guided solutions or labeled training data. GPO adopts a modified Fourier neural operator as the backbone to achieve significantly improved convergence speed and to learn the complicated solution field of solid mechanics problems. Numerical experiments involved a two-dimensional plane with a hole and a three-dimensional building structure with Dirichlet boundary constraints. The results indicate that the geometry layer and exact boundary constraints in GPO significantly contribute to the convergence accuracy and speed, outperforming the previous benchmark in simulations of irregular geometry. The comparison results also showed that GPO can converge to solution fields faster than a commercial numerical solver in the structural examples. Furthermore, GPO demonstrates stronger performance than the solvers when the mesh size is smaller, and it achieves over 3<span data-altimg=\"/cms/asset/df3c163a-e7c6-4e48-9813-d4de41ec9066/mice13405-math-0001.png\"></span><mjx-container ctxtmenu_counter=\"93\" ctxtmenu_oldtabindex=\"1\" jax=\"CHTML\" role=\"application\" sre-explorer- style=\"font-size: 103%; position: relative;\" tabindex=\"0\"><mjx-math aria-hidden=\"true\" location=\"graphic/mice13405-math-0001.png\"><mjx-semantics><mjx-mo data-semantic- data-semantic-role=\"unknown\" data-semantic-speech=\"times\" data-semantic-type=\"operator\"><mjx-c></mjx-c></mjx-mo></mjx-semantics></mjx-math><mjx-assistive-mml display=\"inline\" unselectable=\"on\"><math altimg=\"urn:x-wiley:10939687:media:mice13405:mice13405-math-0001\" display=\"inline\" location=\"graphic/mice13405-math-0001.png\" xmlns=\"http://www.w3.org/1998/Math/MathML\"><semantics><mo data-semantic-=\"\" data-semantic-role=\"unknown\" data-semantic-speech=\"times\" data-semantic-type=\"operator\">×</mo>$\\times$</annotation></semantics></math></mjx-assistive-mml></mjx-container> and 2<span data-altimg=\"/cms/asset/67f630b9-5c45-4ee8-a436-0fef5e94744b/mice13405-math-0002.png\"></span><mjx-container ctxtmenu_counter=\"94\" ctxtmenu_oldtabindex=\"1\" jax=\"CHTML\" role=\"application\" sre-explorer- style=\"font-size: 103%; position: relative;\" tabindex=\"0\"><mjx-math aria-hidden=\"true\" location=\"graphic/mice13405-math-0002.png\"><mjx-semantics><mjx-mo data-semantic- data-semantic-role=\"unknown\" data-semantic-speech=\"times\" data-semantic-type=\"operator\"><mjx-c></mjx-c></mjx-mo></mjx-semantics></mjx-math><mjx-assistive-mml display=\"inline\" unselectable=\"on\"><math altimg=\"urn:x-wiley:10939687:media:mice13405:mice13405-math-0002\" display=\"inline\" location=\"graphic/mice13405-math-0002.png\" xmlns=\"http://www.w3.org/1998/Math/MathML\"><semantics><mo data-semantic-=\"\" data-semantic-role=\"unknown\" data-semantic-speech=\"times\" data-semantic-type=\"operator\">×</mo>$\\times$</annotation></semantics></math></mjx-assistive-mml></mjx-container> speedup for a large degree of freedom in the two-dimensional and three-dimensional examples, respectively. The limitations of nonlinearity and complicated structures are further discussed for prospective developments. The remarkable results suggest the potential modeling applications of large-scale infrastructures.","PeriodicalId":156,"journal":{"name":"Computer-Aided Civil and Infrastructure Engineering","volume":"73 1","pages":""},"PeriodicalIF":8.5000,"publicationDate":"2025-01-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Computer-Aided Civil and Infrastructure Engineering","FirstCategoryId":"5","ListUrlMain":"https://doi.org/10.1111/mice.13405","RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}
引用次数: 0
Abstract
This study developed Geometry Physics neural Operator (GPO), a novel solver framework to approximate the partial differential equation (PDE) solutions for solid mechanics problems with irregular geometry and achieved a significant speedup in simulation time compared to numerical solvers. GPO leverages a weak form of PDEs based on the principle of least work, incorporates geometry information, and imposes exact Dirichlet boundary conditions within the network architecture to attain accurate and efficient modeling. This study focuses on applying GPO to model the behaviors of complicated bodies without any guided solutions or labeled training data. GPO adopts a modified Fourier neural operator as the backbone to achieve significantly improved convergence speed and to learn the complicated solution field of solid mechanics problems. Numerical experiments involved a two-dimensional plane with a hole and a three-dimensional building structure with Dirichlet boundary constraints. The results indicate that the geometry layer and exact boundary constraints in GPO significantly contribute to the convergence accuracy and speed, outperforming the previous benchmark in simulations of irregular geometry. The comparison results also showed that GPO can converge to solution fields faster than a commercial numerical solver in the structural examples. Furthermore, GPO demonstrates stronger performance than the solvers when the mesh size is smaller, and it achieves over 3 and 2 speedup for a large degree of freedom in the two-dimensional and three-dimensional examples, respectively. The limitations of nonlinearity and complicated structures are further discussed for prospective developments. The remarkable results suggest the potential modeling applications of large-scale infrastructures.
期刊介绍:
Computer-Aided Civil and Infrastructure Engineering stands as a scholarly, peer-reviewed archival journal, serving as a vital link between advancements in computer technology and civil and infrastructure engineering. The journal serves as a distinctive platform for the publication of original articles, spotlighting novel computational techniques and inventive applications of computers. Specifically, it concentrates on recent progress in computer and information technologies, fostering the development and application of emerging computing paradigms.
Encompassing a broad scope, the journal addresses bridge, construction, environmental, highway, geotechnical, structural, transportation, and water resources engineering. It extends its reach to the management of infrastructure systems, covering domains such as highways, bridges, pavements, airports, and utilities. The journal delves into areas like artificial intelligence, cognitive modeling, concurrent engineering, database management, distributed computing, evolutionary computing, fuzzy logic, genetic algorithms, geometric modeling, internet-based technologies, knowledge discovery and engineering, machine learning, mobile computing, multimedia technologies, networking, neural network computing, optimization and search, parallel processing, robotics, smart structures, software engineering, virtual reality, and visualization techniques.