Mesh Orientation Research Articles

An Unstructured Mesh Coordinate Transformation‐Based FDTD Method

ABSTRACTWe propose a novel unstructured mesh finite‐difference time‐domain (FDTD) method for solving electromagnetics problems with complicated geometries. The method, which solves the TE‐mode reduced form of Maxwell's equations, can handle both material interfaces and anisotropic material. Using the transformation optics principle, which describes how fields and material tensors change under coordinate transformations, we locally transform each cell in the mesh to a reference unit‐square computational domain where the usual FDTD update is performed. This comes at a cost: employing unstructured grids and coordinate transformations requires more complicated data structures, a mesh orientation process, and potentially introduces an anisotropic material tensor at every mesh cell. Nonetheless, we find that the method maintains the same desirable properties of the classic FDTD method (explicit, divergence‐free B‐field, nondissapative, and second‐order accuracy) while also gaining conforming material interfaces and boundaries in complicated geometry. Even further, we prove that the method is stable under a Courant condition and a fairly nonrestrictive mesh condition, hence defeating the late‐time stability issue plaguing prior nonorthogonal FDTD methods. To verify the method, we conduct convergence studies on three electromagnetic cavity problems with known exact solutions. For these numerical studies, we find that the method maintains second‐order convergence and stability.

The Influence of Thickness and Mesh Orientation on The Impact Resistance Performance of Geogrid Reinforced Oil Palm Shell Concrete

The Influence of Thickness and Mesh Orientation on The Impact Resistance Performance of Geogrid Reinforced Oil Palm Shell Concrete

Effectiveness of Different Configurations of Ferrocement Retrofitting for Seismic Protection of Confined Masonry: A Numerical Study

A ferrocement layer, which consists of a wire mesh and cement mortar, is a popular retrofitting method for existing structural elements, particularly wall or slab panels. This paper presents a study on the effectiveness of different configurations of ferrocement for seismic retrofitting of confined masonry through finite element analysis. The masonry panel was modeled using expanded brick-unit elements, where the element was expanded in size by as much as half of the mortar thickness, and an interacting zero-thickness interface was applied to mimic the elastic-plastic and damage behavior during tension, shear, and compression. The concrete damage plasticity (CDP) model was used to model the confining reinforced concrete frame and overlay mortar in the ferrocement layer, and the reinforcing bars and wire mesh were modeled using elastic-plastic behavior. In the present numerical study, nine models were subjected to cyclic and pushover shear test simulations, considering the effects of the number of ferrocement layers and the wire mesh orientation. The volumetric ratio of the wire mesh to the masonry (ρwm) ranged from 0.48% to 1.92%, whereas the ratio of the mortar overlay to the masonry (ρmo) varies from 10.42% to 41.66%. Based on the increase in the lateral strength, the model with the largest volume of the ferrocement layer exhibited the largest increase in strength. However, the most cost-effective retrofitting configuration was presented by model DS-1-45, in which a single layer of ferrocement was applied on both sides of the wall using 45° of wire mesh orientation. The DS-1-45 model provided a lateral strength increase of more than 6 times compared to the original unreinforced model. Doi: 10.28991/CEJ-2024-010-09-02 Full Text: PDF

Rapid particle generation from an STL file and related issues in the application of material point methods to complex objects

In this paper, we focus on three issues related to applications of material point methods (MPMs) to objects with complex geometries. They are material point generation, compatibility of material points with a mesh, and sensitivity to mesh orientation. An efficient method of generating material points from a stereolithography (STL) file is introduced. This material point generation method is independent of the mesh used in MPM calculations. The compatibility between the material points and the mesh is then studied. We also show that the original MPM and the dual domain material point (DDMP) method are sensitive to mesh orientation. These issues are related to the calculation of the internal force and are concerns of the MPMs. They become more prominent when MPMs are applied to complex geometries. Our numerical results show that the recently developed local stress difference (LSD) algorithm (Perez et al. in J Comp Phys 498:112681, 2024) can be used to effectively address them.

Oscillation-free implicit pressure explicit concentration discontinuous Galerkin methods for compressible miscible displacements with applications in viscous fingering

The system of compressible miscible displacements is widely adopted to model surfactant flooding in enhanced oil recovery techniques, where a low-viscosity fluid is injected underground to replace the high-viscosity oil. When the mobility ratio of the injected fluid to oil is high, the waterflood front tends to be unstable and exhibits a finger-like growth pattern, known as viscous fingering. Due to its unstable nature, the viscous fingering phenomenon is sensitive to mesh orientation and numerical discretization. Therefore, high-order numerical methods are preferable to reduce numerical artifacts and mesh dependence. In this paper, we propose a high-order discontinuous Galerkin method for the coupled nonlinear system of compressible miscible displacements to simulate the viscous fingering fluid instability in porous media. We adopt the IMplicit Pressure Explicit Concentration time marching approach based on implicit-explicit Runge-Kutta methods to achieve high-order temporal accuracy. Additionally, we introduce an oscillation-free damping term to control the spurious oscillations encountered in the waterflood front due to the large gradient of saturation. We have conducted ample numerical tests in two space dimensions to demonstrate the effectiveness and robustness of the proposed schemes in recovering viscous fingering.

A fibre direction informed continuum damage model with enriched kinematics for modelling matrix cracking in composites

Continuum Damage Mechanics models are widely used for simulation of ply-level damage such as matrix cracking in composites. Although their ease of implementation within the finite element framework is appealing, several limitations exist such as a smeared representation of sharp cracks, mesh orientation bias of crack growth, inaccurate energy dissipation due to approximate definition of characteristic length etc. In this work, a new ‘semi-discrete’ continuum damage model is proposed for matrix cracking, that includes sharp crack kinematics. The mesh orientation bias is eliminated by a fibre-directed growth driving criterion. The crack characteristic length is exactly computed due to an accurate geometric representation of the crack surface. The model also addresses a well-known problem of spurious stress transfer across crack faces when deformation becomes large. Through several challenging benchmark problems, the model is shown to perform nearly as good as discrete crack models while being simpler and more robust.

Experimental investigation on the effects of stainless-steel mesh reinforcing layers on low-velocity impact response of hybrid thermoplastic glass fiber composites

This study aims to assess the hybridization effect on the perforation threshold of Low-Velocity Impact (LVI) in thermoplastic glass composite laminates, incorporating layers of resin-impregnated stainless-steel mesh. Reinforcing methodologies such as hybridization are recently being adopted as a practical approach to increasing the energy-absorbing capacity of polymer composites. In the current paper, a multi-step hot press lamination method has been employed to fabricate the hybrid composite laminates strengthened with stainless-steel mesh layers. Several stacking sequences, metal mesh wire sizes, orientation and position relative to the impactor have been examined under various LVI energies. It was revealed that the LVI penetration energy was increased for the thermoplastic-based composite laminates reinforced with stainless-steel mesh layers. Furthermore, the LVI penetration energy threshold was significantly influenced by the metal mesh wire size, orientation and stacking sequence. Finally, the backlight method capability was assessed to detect the after-impact interlaminar damages.

Structural behavior of RC one-way slabs strengthened with ferrocement and FRP composites

Existing research lacks focus on the structural behavior of reinforced concrete (RC) one-way slabs reinforced with ferrocement, and there's a notable absence of comparative studies between fiber-reinforced polymers (FRPs) and ferrocement strengthening methods for these slabs. Given the cost-effectiveness and widespread use of ferrocement in structural reinforcement, this study aimed to address these gaps through an experimental program. Three sizes of wire mesh, categorized as Type-I (small), Type-II (medium), and Type-III (large), were employed in this study. Chemical or mechanical anchors were used to attach ferrocement jackets. Moreover, 6, 12, or 18 anchors were used to assess the effect of anchor configuration. The goal was to enhance the structural performance of slabs and compare them with slabs reinforced using FRP jackets. The study focused on preventing debonding of the strengthening layers, employing either mechanical or chemical anchors. All slabs exhibited ductile failure with flexural cracks. The peak load and ultimate deflection were enhanced by up to 49.00% and 109.07%, respectively, by the application of ferrocement jackets, whereas the dissipated energy was increased by up to 174.00 %. Notably, the use of chemical anchors demonstrated a superior ability to delay debonding and enhance ductility compared to mechanical anchors. Slabs reinforced with glass FRP (GFRP) showed delayed debonding relative to carbon FRP (CFRP) reinforced slabs, indicating the superior performance of chemical anchors with GFRP layers. Moreover, the type and size of wire mesh significantly influenced performance, with small and medium-sized mesh configurations enhancing ductility, while large-sized mesh exhibited relatively earlier debonding. The orientation of the wire mesh also played a crucial role, with parallel orientation to the longitudinal axis of slabs yielding better performance.

Mechanical properties of laminate materials based on polylactic acid and polyvinyl chloride meshes as reinforcement

The 3D printing parameters are known to have a significant impact on manufactured parts, and the layered morphology of these parts makes mechanical design analysis for engineering applications difficult. In this work, the tensile strengths and microhardness of 3D printed polylactic acid (PLA) specimens with different orientations and numbers of individual layers of mesh material (polyvinyl chloride – PVC) were investigated as a laminate composite. Composite specimens were obtained using 3D printing via fused deposition modelling (FDM). Moreover, the influence of printing parameters (i.e. infill density and layer height) and the number and orientation of reinforced meshes on the mechanical response was investigated. Fracture strength of PLA/PVC laminate composites ranges from 31.30 MPa (3 PVC mesh layers; mesh height position: 25 %│50 %│75 %; infill density: 60 %; PVC mesh orientation: 90°│45°│90°; layer height: 0.2 mm) to 18.62 MPa (without PVC mesh; infill density: 30 %; layer height: 0.1 mm) demonstrating a significant impact of the number of the PVC mesh layers, infill density of PLA and layer height on the final mechanical parameters of printing PLA/PVC elements. The surface hardness at the micro load level showed that the number of reinforcement layers affects the microhardness value, as well as material filling and mesh orientation. The specimen with the following parameters gave the best results: layer height: 0.2 mm; 3 PVC mesh layers; infill density: 60 %; PVC mesh orientation: 90°│45°│90°. The average hardness values for one layer and three layers of mesh were in accordance with tensile test results.

Efficient modelling of progressive damage due to quasi-static indentation on multidirectional laminates by a mesh orientation independent kinematically enriched continuum damage model

This work proposes a computationally efficient high-fidelity modelling framework for analysis of damage in composites subjected to quasi-static indentation. A ply-by-ply modelling approach is adopted, and a new semi-discrete continuum damage model is used for intra-ply cracking that allows for cracks to grow independent of the ply mesh pattern. This feature greatly simplifies the meshing effort since the requirement of a ply-oriented mesh and imposition of tie constraints at mesh-mismatched ply interfaces is eliminated. The model is further enhanced to realise kinematic interactions between the cracked and uncracked material domains at the constitutive level. Inter-ply delamination is simulated using cohesive elements. Through a challenging problem of static indentation on a quasi-isotropic laminate, it is shown that the model can capture solution-dependent multiple discrete ply cracks and detailed crack-delamination interaction, but with reduced computational time and complexity, compared to a reference model with ply-oriented mesh and pre-seeded cracks at known locations.

Effect of Textile Reinforced Mortar components on the mechanical behavior of historical masonry structures under cyclic loading

Masonry structures without reinforcement are generally considered unsuitable in seismically active areas due to their poor performance as a bracing system.Historical masonry is particularly affected by this issue. The need to find solutions to this problem has attracted the interest of the scientific community to preserve this cultural heritage.Textile reinforced mortar (TRM) has proven his efficiency to reinforce masonry structures and various researches have investigated this concept of reinforcement in the case of diagonal compression tests. Despite this, several parameters must be studied in order to optimize this reinforcement technique in boundary conditions close to seismic situation.The study presented here aimed at assessing the effect of varying the mesh size of the textile and the orientation of its threads in strengthening historical masonry constructions with testing conditions similar to those of an earthquake (cyclic tests). It also aimed to assess the influence of varying lime mortar on the reinforcement performance by using two different limes in the TRM reinforcement.The principal additional knowledge of the present work is to point out the effect of mesh orientation and mesh opening in terms of dissipated energy and in terms of delaying the debonding time of TRM (a longer time for the multidirectional EQ grid) which is the main issue in the case of TRM structural reinforcement.An increase in peak load (+110 % for RW-EQ-L2) and in the dissipated energy (+200 % for RW-EQ-L2) were hence observed during the tests. These results provide valuable information about masonry walls strengthening and repair for the scientific community.

Smooth Lagrangian Crack Band Model Based on Spress-Sprain Relation and Lagrange Multiplier Constraint of Displacement Gradient

Abstract A preceding 2023 study argued that the resistance of a heterogeneous material to the curvature of the displacement field is the most physically realistic localization limiter for softening damage. The curvature was characterized by the second gradient of the displacement vector field, which includes the material rotation gradient, and was named the “sprain” tensor, while the term “spress” is here proposed as the force variable work-conjugate to “sprain.” The partial derivatives of the associated sprain energy density yielded in the preceeding study, sets of curvature resisting self-equilibrated nodal sprain forces. However, the fact that the sprain forces had to be applied on the adjacent nodes of a finite element greatly complicated the programming and extended the simulation time in a commercial code such as abaqus by almost two orders of magnitude. In the present model, Smooth Lagrangian Crack Band Model (slCBM), these computational obstacles are here overcome by using finite elements with linear shape functions for both the displacement vector and for an approximate displacement gradient tensor. A crucial feature is that the nodal values of the approximate gradient tensor are shared by adjacent finite elements. The actual displacement gradient tensor calculated from the nodal displacement vectors is constrained to the approximate displacement gradient tensor by means of a Lagrange multiplier tensor, either one for each element or one for each node. The gradient tensor of the approximate gradient tensor then represents the approximate third-order displacement curvature tensor, or Hessian of the displacement field. Importantly, the Lagrange multiplier behaves as an externally applied generalized moment density that, similar to gravity, does not affect the total strain-plus-sprain energy density of material. The Helmholtz free energy of the finite element and its associated stiffness matrix are formulated and implemented in a user’s element of abaqus. The conditions of stationary values of the total free energy of the structure with respect to the nodal degrees-of-freedom yield the set of equilibrium equations of the structure for each loading step. One- and two-dimensional examples of crack growth in fracture specimens are given. It is demonstrated that the simulation results of the three-point bend test are independent of the orientation of a regular square mesh, capture the width variation of the crack band, the damage strain profile across the band, and converge as the finite element mesh is refined.

Isogeometric smooth crack-band model (isCBM) using spress–sprain relations adapted to microplane theory

Recent optimum fitting of the data from eleven different types of distinctive fracture tests of quasibrittle materials, including the gap tests and the size effect tests, revealed that the crack band model (CBM) with the M7 damage constitutive law outperforms all the other existing computational fracture models for concrete. This is attributed to physically correct boundary conditions, especially those at the crack faces,to the crack front of finite width which allows using a triaxial tensorial damage constitutive law, and to choosing for this law a realistic material behavior description, the M7. Nevertheless the CBM has limitations—the width of the crack band cannot be varied, the statistically smooth damage distribution across the band cannot be resolved, and the propagation direction is biased by the mesh orientation when a regular mesh is used. These limitations have been overcome by the smooth crack band model (sCBM), which uses unconventional homogenization of damage mechanics to yield an additional Helmholtz free energy density, called the “sprain energy”. This energy is a function of the curvature tensor (or a hessian) of the vectorial displacement field, called the “sprain curvature tensor” (or sprain tensor for brevity), which is a sum of the strain gradient tensor and, importantly, the rotation gradient tensor. But, in a previous study, difficulties arose in computational implementation—the use of constant-strain elements required the curvature to be resisted by sprain forces that represented the derivatives of sprain energy with respect to nodal displacements and had to be applied on the nodes of adjacent elements. The need for adjacent nodes is here avoided by a variational derivation of the Isogeometric Smooth Crack-Band Model (isCBM). The microplane damage constitutive model M7 for concrete is adopted, which is made possible by establishing the relations of the sprain stress (or spress, for brevity) and sprain on microplanes of arbitrary orientations. Using Isogeometric Analysis (IGA) based on non-uniform rational B-Splines (NURBS), the excessive displacement curvature and its resistance (spress) can be described by the shape functions with C1-continuity. Computational results show that the orientation bias of regular meshes is eliminated while keeping the constitutive law intact, which is important for compatibility with the existing material laws in general. In addition, comparisons with experimental data reveal that isCBM with M7 can reproduce the correspondence between the size of the fracture process zone and the fracture energy, as indicated by gap tests. Finally, the isCBM does not compromise the strain-path-dependence of material response.

Comparative verification of hydro-mechanical fracture behavior: Task G of international research project DECOVALEX–2023

Numerical simulations become a necessity when experimental approaches cannot cover the required physical and time scale of interest. One of such area is a simulation of long-term host rock behaviors for nuclear waste disposal and simulation tools involved in the assessment must go through rigorous validation tests. The DECOVALEX project (Development of COupled models and their VAlidation against EXperiments) is dedicated to this purpose by international participants.aawww.decovalex.org. This work is part of the ongoing phase DECOVALEX–2023 (D–2023, Task G) particularly aiming to simulate fracture behaviors under various conditions. Here, we cross-verified a variety of numerical methods including continuous and discontinuous approaches against four benchmark exercises with emphasis on numerical accuracy and parameterization of the various numerical approaches. The systematic inter-comparisons of test cases highlight advantages and disadvantages of the different numerical models. Numerical details on discretization effects (e.g. mesh density and orientation) and domain size were investigated in detail for practical applications. It became evident that meticulous attention to mesh resolution and domain size is imperative for achieving accurate numerical simulations, even for static cracks. Moreover, when comparing numerical methods to closed-form solutions for static cracks, all models successfully reproduced the maximum crack opening but encountered challenges near the crack tips. Finally, the paper discusses how to convert between and therefore compare parameters of various numerical approaches. Our benchmark studies reveal that each model necessitates a distinct number of parameters, even in simple scenarios like static crack aperture benchmarks. It is generally more practical to employ fewer parameters to mitigate model over-parameterization and enhance experimental feasibility.

Comparison of local and nonlocal regularization approaches in Eulerian-based finite element analyses

Finite element (FE) results might suffer from significant mesh dependency, especially when modeling shear bands in strain-softening soils. On the other hand, the shear bands in the field are generally thin. The computational costs might dramatically increase when modeling such thin shear bands, especially for large-scale problems, such as progressive landslides. To overcome these issues, FE analyses are generally performed with larger shear band thicknesses by adopting ‘softening scaling’ rules based on local and nonlocal (averaged) strains at the shear bands. The present study compares the performances of local and nonlocal approaches, with a specific focus on softening scaling and large deformation. Analyses are performed using a Eulerian-based FE program, which can model large strains in the shear band without numerical issues related to mesh distortion. Implementing strain-softening behavior in local and nonlocal methods, two idealized cases are simulated: (i) biaxial compression test, (ii) slope failure due to upslope surface loading. FE simulations show that the macroscopic response (i.e., load–displacement behavior) can be modeled using both local and nonlocal regularization techniques. Shear band thickness increases with the progress of shearing. In local analysis, mesh orientation has a considerable effect on shear band thickness. The existence of neighboring shear bands could affect nonlocal strain calculation and the simulation results.

Fabrication and characterization of laminated natural fibers and SS303 wire mesh reinforced epoxy-based hybrid composite

In today’s world, researchers are actively seeking the most suitable materials to meet the growing industrial demands. These materials need to possess specific characteristics such as lightweight nature, high strength, excellent mechanical properties, and environmental friendliness. Due to these requirements, relying solely on mono fiber composite materials is insufficient. Therefore, it has become necessary to incorporate different reinforcements into polymer matrices to achieve the desired mechanical properties. This study focuses on the fabrication of a novel hybrid composite using jute fibers, luffa fibers, palmyra leaf fibers, and AISI 303 (0.5 mm) wire mesh combined with epoxy resin. The fabrication process involved utilizing vacuum bagging techniques. To assess the mechanical properties of the hybrid composites, two wire mesh orientations (45° and 90°) and a total of 12 stacking sequences of jute/wire mesh/luffa/palmyra (J/W/L/P) were selected. The mechanical characterization of the hybrid composites included tensile, impact, flexural, and interlaminar strength tests. Furthermore, a dynamic mechanical analyzer (DMA-8000) was employed to investigate the elastic behavior, including storage modulus (E′), loss modulus ( E″), and damping factor (tan δ). Additionally, the fracture surface's microstructure was examined using a scanning electron microscope (SEM). Among the various stacking sequences and fiber orientations, the 45° fiber orientation and the stacking sequences JWLP, JLWP, JWPL, JPWL, PWJL, and PJWL exhibited superior ultimate tensile strength. Notably, the PWJL stacking sequence displayed the highest average tensile strength (74.83 MPa), flexural strength (131.60 MPa), percentage elongation (1.35%), and interlaminar strength (1.47 kN). Moreover, dynamic mechanical analysis revealed that the hybrid composite sample with 90° wire mesh orientation exhibited a peak energy absorption of 0.820 at the transition region, specifically at a frequency of 0.5 Hz.

Delving Deep into Pixel Alignment Feature for Accurate Multi-View Human Mesh Recovery

Regression-based methods have shown high efficiency and effectiveness for multi-view human mesh recovery. The key components of a typical regressor lie in the feature extraction of input views and the fusion of multi-view features. In this paper, we present Pixel-aligned Feedback Fusion (PaFF) for accurate yet efficient human mesh recovery from multi-view images. PaFF is an iterative regression framework that performs feature extraction and fusion alternately. At each iteration, PaFF extracts pixel-aligned feedback features from each input view according to the reprojection of the current estimation and fuses them together with respect to each vertex of the downsampled mesh. In this way, our regressor can not only perceive the misalignment status of each view from the feedback features but also correct the mesh parameters more effectively based on the feature fusion on mesh vertices. Additionally, our regressor disentangles the global orientation and translation of the body mesh from the estimation of mesh parameters such that the camera parameters of input views can be better utilized in the regression process. The efficacy of our method is validated in the Human3.6M dataset via comprehensive ablation experiments, where PaFF achieves 33.02 MPJPE and brings significant improvements over the previous best solutions by more than 29%. The project page with code and video results can be found at https://kairobo.github.io/PaFF/.

Large deformation and brittle failure calculated using the dual-domain material point method

The dual domain material point (DDMP) method is explored as a candidate to be implemented in a general purpose code to perform simulations of materials with complex geometry that undergo large history-dependent deformation and failure. To test its candidacy, we study its mesh convergence, its sensitivity to mesh orientation, and its ability to handle softening and failure of a material. Simulations of large deformation and simulations of mechanical failure are performed using both DDMP and the material point method (MPM). When cell-crossing of material points is not an issue and when there are a sufficient number of material points in each computation cell, the numerical error decreases with the square of the cell size as expected for both MPM and DDMP. DDMP has reduced error compared with MPM when there are many instances of material points crossing cell boundaries due to the continuous nature of the modified gradient of the shape functions. Simulations of a specimen under tension are also performed where the background mesh is aligned and misaligned with the tension direction. MPM displays a significant mesh-dependent stress field, DDMP shows negligible mesh dependency. Despite a mesh orientation-dependent stress field from MPM, the critical tension and failure mode from both MPM and DDMP calculations have negligible mesh dependency when using a non-local failure model. If only the failure mode is important (i.e., local stresses are unimportant), MPM with a non-local failure model is a suitable method for modeling failure with small deformations. However, if local stresses are also important or if there are large deformations with many cell-crossings before failure, DDMP should be the method that is used. A needed improvement for DDMP is identified from our numerical simulations.

Sliding or Stumbling on the Staircase: Numerics of Ocean Circulation Along Piecewise‐Constant Coastlines

AbstractCoastlines in most ocean general circulation models are piecewise constant. Accurate representation of boundary currents along staircase‐like coastlines is a long‐standing issue in ocean modeling. Pioneering work by Adcroft and Marshall (1998, https://doi.org/10.3402/tellusa.v50i1.14514) revealed that artificial indentation of model coastlines, obtained by rotating the numerical mesh within an idealized square basin, generates a spurious form drag that slows down the circulation. Here, we revisit this problem and show how this spurious drag may be eliminated. First, we find that physical convergence to spatial resolution (i.e., the main characteristics of the flow are insensitive to the increase of the mesh resolution) allows simulations to become independent of the mesh orientation. An advection scheme with a wider stencil also reduces sensitivity to mesh orientation from coarser resolution. Second, we show that indented coastlines behave as straight and slippery shores when a true mirror boundary condition on the flow is imposed. This finding applies to both symmetric and rotational‐divergence formulations of the stress tensor, and to both flux and vector‐invariant forms of the equations. Finally, we demonstrate that the detachment of a vortex flowing past an outgoing corner of the coastline is missed with a free‐slip (zero vorticity) condition at the corner. These results provide guidance for a better numerical treatment of coastlines (and isobaths) in ocean general circulation models.

Automated and data-driven plate computation for presurgical cleft lip and palate treatment

PurposePresurgical orthopedic plates are widely used for the treatment of cleft lip and palate, which is the most common craniofacial birth defect. For the traditional plate fabrication, an impression is taken under airway-endangering conditions, which recent digital alternatives overcome via intraoral scanners. However, these alternatives demand proficiency in 3D modeling software in addition to the generally required clinical knowledge of plate design.MethodsWe address these limitations with a data-driven and fully automated digital pipeline, endowed with a graphical user interface. The pipeline adopts a deep learning model to landmark raw intraoral scans of arbitrary mesh topology and orientation, which guides the nonrigid surface registration subsequently employed to segment the scans. The plates that are individually fit to these segmented scans are 3D-printable and offer optional customization.ResultsWith the distance to the alveolar ridges closely centered around the targeted 0.1 mm, our pipeline computes tightly fitting plates in less than 3 min. The plates were approved in 12 out of 12 cases by two cleft care professionals in a printed-model-based evaluation. Moreover, since the pipeline was implemented in clinical routine in two hospitals, 19 patients have been undergoing treatment utilizing our automated designs.ConclusionThe results demonstrate that our automated pipeline meets the high precision requirements of the medical setting employed in cleft lip and palate care while substantially reducing the design time and required clinical expertise, which could facilitate access to this presurgical treatment, especially in low-income countries.