Topological Semantics for Lumped Parameter Systems Modeling
Randi Wang, Vadim Shapiro
Technical Report, Spatial Automation Laboratory, University of WisconsinMadison 

Multiscale shape–material modeling by composition
We propose a formal framework for modeling multiscale material structures by recursive composition of twoscale material structures. The framework comprises three components: (1) single scale shapematerial models, supported by single scale queries, to represent the geometry and spatial distribution of material property on each coarse and fine scales, (2) mechanisms to link the scales by establishing an explicit relationship between shapematerial properties at fine scale and material properties at the coarse scale, and (3) multiscale queries abstracting fundamental multiscale operations by recursive composition. While the first component is consistent with classical solid heterogeneous material modeling, the second component manifests itself as a pair of conceptually new upscaling and downscaling functions. We show that classical solid modeling queries, exemplified by point membership testing, distance computation, and material evaluation, generalize to the corresponding multiscale queries that support implicit representations of multiscale structures as a composition of distinct single scale solid material models. The concept of neighborhood is indispensable in all three components. The framework provides a formal and consistent extension of solid modeling framework that underlies most commercial systems in use today, encompasses the variety of different approaches to multiscale modeling, identifies open issues and research problems with existing twoscale modeling methods, and provides foundations for nextgeneration systems by identifying key objects, classes, representation schemes, and API queries. Xingchen Liu, Vadim Shapiro.
ComputerAided Design, Volume 102, September 2018, Pages 194203 

Shape Aware Quadratures
Vaidyanathan Thiagarajan, Vadim Shapiro
Preprint, December 12, 2017 

Error Analysis of Adaptively Weighted Numerical Integration
Vaidyanathan Thiagarajan, Vadim Shapiro
Preprint, August 28, 2017 

LinearTime Thermal Simulation of AsManufactured FDM Components
Yaqi Zhang, Vadim Shapiro
Journal of Manufacturing Science and Engineering, Vol. 140, July 2018 

Structural Analysis of Laminated Composites using FunctionBehaviorStructure Model and Virtual Material Method
Goldy Kumar
Dissertation, May 2017, University of Wisconsin–Madison 

SAMPLE – BASED MATERIAL STRUCTURE MODELING
Xingchen Liu
Dissertation, May 2017, University of Wisconsin–Madison 

Samplebased synthesis of twoscale structures with anisotropy
Additive manufacturing transforms material into threedimensional parts incrementally, layer by layer or path by path. Subject to the build direction and machine resolution, an additively manufactured part deviates from its design model in terms of both geometry and mechanical performance. In particular, the material inside the fabricated part often exhibits spatially varying material distribution (heterogeneity) and direction dependent behavior (anisotropy), indicating that the design model is no longer a suitable surrogate to consistently estimate the mechanical performance of the printed component. Xingchen Liu, Vadim Shapiro.
ComputerAided Design, Volume 90, May 2017, Pages 199209 

SampleBased Synthesis of Functionally Graded Material Structures
Spatial variation of material structures is a principal mechanism for creating and controlling spatially varying material properties in nature and engineering. While the spatially varying homogenized properties can be represented by scalar and vector fields on the macroscopic scale, explicit microscopic structures of constituent phases are required to facilitate the visualization, analysis and manufacturing of functionally graded material (FGM). The challenge of FGM structure modeling lies in the integration of these two scales. We propose to represent and control material properties of FGM at macroscale using the notion of material descriptors which include common geometric, statistical, and topological measures, such as volume fraction, correlation functions and Minkowski functionals. At microscale, the material structures are modeled as Markov random fields: we formulate the problem of design and (re)construction of FGM structure as a process of selecting neighborhoods from a reference FGM, based on target material descriptors fields. The effectiveness of the proposed method in generating a spatially varying structure of FGM with target properties is demonstrated by two examples: design of a graded bone structure and generating functionally graded lattice structures with target volume fraction fields. Xingchen Liu, Vadim Shapiro.
Journal of Computing and Information Science in Engineering, Volume 17(3), September 2017
A preliminary version appeared as technical paper in 2016 ASME International Design and Engineering Technical Conference, IDETC/CIE 2016, Charlotte, North Carolina


Homogenization of material properties in additively manufactured structures
Additive manufacturing transforms material into threedimensional parts incrementally, layer by layer or path by path. Subject to the build direction and machine resolution, an additively manufactured part deviates from its design model in terms of both geometry and mechanical performance. In particular, the material inside the fabricated part often exhibits spatially varying material distribution (heterogeneity) and direction dependent behavior (anisotropy), indicating that the design model is no longer a suitable surrogate to consistently estimate the mechanical performance of the printed component. Xingchen Liu, Vadim Shapiro.
ComputerAided Design, Volume 78, September 2016, Pages 7182 

Adaptively Weighted Numerical Integration in the Finite Cell Method
Vaidyanathan Thiagarajan, Vadim Shapiro
Comput. Methods Appl. Mech. Engrg. 311 (2016): 250279 

Effective Contact Measures
Mikola Lysenko, Vadim Shapiro
ComputerAided Design, 70 (2016) 134143 

The New Frontiers in Computational Modeling of Material Structures
William Regli, Jarek Rossignac, Vadim Shapiro, Vijay Srinivasan
ComputerAided Design 77 (2016): 012 

Efficient 3D analysis of laminate structures using ABDequivalent material models Laminate composites are widely used in automotive, aerospace, and increasingly in consumer industries, due to their reduced weight and superior structural properties. However, structural analysis of complex laminate structures remains challenging. 2D finite element methods based on plate/shell theories may be accurate and efficient, but they generally do not apply to the whole structure and require identification and preprocessing of the regions where the underlying assumptions hold. Fully automated structural analysis using solid 3D elements with sufficiently high order basis functions is possible in principle, but is rarely practiced due to the significant increase in the cost of computational integration over a large number of laminate plies. We propose a procedure to replace the original laminate by much simpler new virtual material models. These virtual material models, under the usual assumptions made in lamination theory, have the same constitutive relationship as the corresponding 2D plate model of the original laminate, but require only a small fraction of computational integration costs in 3D FEA. We describe implementation of 3D FEA using these material models in a meshfree system using second order Bspline basis functions. Finally, we demonstrate their validity by showing agreement between computed and known results for standard problems. Goldy Kumar, Vadim Shapiro.
Finite Elements in Analysis and Design, Vol. 106, 2015, pages 4155 

Random Heterogeneous Materials via Texture Synthesis
Computer models of random heterogeneous materials are becoming increasingly important in order to support the latest advances in material science, biomedical applications and manufacturing. Such models usually take the form of a microstructure whose geometry is reconstructed from a small material sample, an exemplar. A widely used traditional approach to material reconstruction relies on stochastic optimization to approximate the material descriptors of the exemplar, such as volume fraction and twopoint correlation functions. This approach is computationally intensive and is limited to certain types of isotropic materials. The current invention show that formulating material reconstruction as a Markov Random Field (MRF) texture synthesis leads to a number of advantages over the traditional optimization based approaches. These include improved computational efficiency, preservation of many material descriptors, including correlation functions and Minkowski functionals, ability to reconstruct anisotropic materials, and direct use of the grayscale material images and two dimensional crosssections. Quantifying the quality of reconstruction in terms of correlation functions as material descriptors suggests a systematic procedure for selecting a size of neighborhood, a key parameter in the texture synthesis procedure. We support our observations by experiments using implementation with Gaussian pyramid and periodic boundary conditions. Xingchen Liu, Vadim Shapiro.
Computational Material Science, Vol. 99, March 2015, pages 177189 

Geometric Interoperability Via Queries
The problem of geometric (model and system) interoperability is conceptualized as a nontrivial generalization of the problem of part interchangeability in mechanical assemblies. Interoperability subsumes the problems of geometric model quality, exchange, and interchangeability, as well as system integration. Until now, most of the interoperability proposals have been datacentric. Instead, we advocate a querycentric approach that can deliver interoperable solutions to many common geometric tasks in computer aided design and manufacturing, including model acquisition and exchange, metrology, and computer aided design/analysis integration. Christoph Hoffmann, Vadim Shapiro, Vijay Srinivasan.
ComputerAided Design, Volume 46, January 2014, Pages 148–159 

Linear Algebraic Representation For Topological Structures
With increased complexity of geometric data, topological models play an increasingly important role beyond boundary representations, assemblies, finite elements, image processing, and other traditional modeling applications. While many graph and indexbased data structures have been proposed, no standard representation has emerged as of now. Furthermore, such representations typically do not deal with representations of mappings and functions and do not scale to support parallel processing, open source, and clientbased architectures. We advocate that a proper mathematical model for all topological structures is a (co)chain complex: a sequence of (co)chain spaces and (co)boundary mappings. This in turn implies all topological structures may be represented by a collection of sparse matrices. We propose a Linear Algebraic Representation (LAR) scheme for mod 2 (co)chain complexes using CSR matrices and show that it supports a variety of topological computations using standard matrix algebra, without any overhead in space or running time. A full open source implementation of LAR is available and is being used for a variety of applications. Antonio DiCarlo, Alberto Paoluzzi, Vadim Shapiro.
ComputerAided Design, Volume 46, January 2014, Pages 269–274 

Reduced Material Model of Composite Laminates For 3D Finite Element Analysis
Laminate composites are widely used in automotive, aerospace, medical, and increasingly in consumer industries, due to their reduced weight, superior structural properties and costeffectiveness. However, structural analysis of complex laminate structures remains challenging. 2D finite element methods based on plate and shell theories may be accurate and efficient, but they generally do not apply to the whole structure, and require identification and preprocessing (dimensional reduction) of the regions where the underlying assumptions hold. Differences in and limitations of theories for thin/thick plates and shells further complicate modeling and simulation of composites. Fully automated structural analysis using 3D elements with sufficiently high order basis functions is possible in principle, but is rarely practiced due to significant increase in computational integration cost in the presence of large number of laminate plies. G. Kumar, V. Shapiro.
2012 ASME International Design and Engineering Technical Conference, IDETC/CIE 2014, Buffalo, New York 

Adaptively Weighted Numerical Integration over Arbitrary Domains
In adaptively weighted numerical integration, for a given set of quadrature nodes, order and domain of integration, the quadrature weights are obtained by solving a system of suitable moment fitting equations in least square sense. The moments in the moment equations are approximated over a simplified domain that is homeomorphic to the original domain, and then are corrected for the deviation from the original domain using shape sensitivity analysis. Using divergence theorem, the moments reduce to integrals over the boundary of the simplified domain. The proposed approach supports accurate and efficient computation of quadrature weights for integration of a priori unknown functions over arbitrary 2D and 3D solid domains. Experimental results (2D) indicate that adaptively weighted integration compares favorably with more traditional approaches. Because the adaptively weighted integration avoids excessive domain subdivision, it is useful in many applications and meshfree analysis in particular. V. Thiagarajan, V. Shapiro.
Computers & Mathematics with Applications, Volume 67, 2014. 

Solving Inverse Configuration Space Problems By Adaptive Sampling
Given two shapes in relative motion, an important class of inverse configuration problems are solved by determining relative configurations that maintain setinclusion relationships (noninterference, containment, or contact) between the shapes. This class of inverse problems includes the wellknown problem of constructing a configuration space obstacle, as well as many other problems in computational design such as sweep decomposition, accessibility analysis, and dynamic packaging.We show that solutions to such problems may be efficiently approximated directly in the 6D configuration space SE(3) of relative motions by adaptive sampling. The proposed method relies on a wellknown fact that the manifold of the group SE(3) is a Cartesian product of two manifold subgroups: the group of rotations SO(3) and the group of translations R3. This property allows generating desired configurations by combining samples that are generated in these subgroups independently and adaptively. We demonstrate the effectiveness of the proposed approach on several inverse problems including the problem of sweep decomposition that arises in reverse engineering applications. S. Nelaturi, V. Shapiro
Proceedings of 2012 ACM Symposium on Solid and Physical Modeling (SPM2012), Dijon, France, October 2012. 

Analysis of multimaterial bonded assemblies on a nonconforming mesh
Bonded multimaterial assemblies arise frequently in design, manufacturing, architecture, and materials design. It is a common wisdom that finite element analysis of such assemblies usually requires all components to be represented by compatible finite element meshes; application of meshfree methods in such situations is often considered problematic due to the need to impose additional interface conditions. Neither approach scales to deal with realistically complex models arising in many applications. We propose a simple extension of meshfree analysis on a nonconforming mesh for linear structural analysis of such multimaterial assemblies. The method is simple, can be implemented within most FEA packages and does not require either compatible meshing or complex interface boundary conditions. Our numerical experiments demonstrate that computed results are in good agreement with known analytical and computational results for well studied multimaterial bonded assemblies (lap and butt joints). We also demonstrate application of the proposed method to realistically complex assembly of a mounted sculpture that cannot be easily analysed by other methods. G. Kumar, V. Shapiro
2012 ASME International Design and Engineering Technical Conference, IDETC/CIE 2012, Chicago, IL 

Configuration Workspaces of SeriesParallel Mechanisms
The workspace of a mechanism is the set of positions and orientations that is reachable by its end effector. Workspaces have numerous applications, including motion planning, mechanism design, and manufacturing process planning, but their representation and computation is challenging due to high dimensionality and geometric/topological complexity. We propose a new formulation of the workspace computation problem for a large class of mechanisms represented by seriesparallel constraint graphs. A wide variety of allowable constraints include all lower pair, some higher pair, and noncollision constraints. We show that the workspace of such mechanisms may be computed by a constraint propagation algorithm. After the space of all rigid body motions is discretized, these operations can be efficiently implemented using the Fast Fourier Transform and a depth first search. In contrast to algebraic formulations, the proposed method assures that all configurations in the computed workspace not only satisfy pairwise constraints but can be reached without breaking and reassembling the mechanism. B. McCarthy, M. Lysenko, V. Shapiro.
2012 ASME International Design and Engineering Technical Conference, IDETC/CIE 2012, Chicago, IL 

Geometric Interoperability for Resilient Manufacturing
Resilient manufacturing systems require adaptable, trustable, and affordable solutions in modelbased engineering (MBE) and platformbased engineering (PBE). Limited or poor geometric interoperability of the software supporting manufacturing and other engineering activities within the product life cycle is becoming a barrier not only for MBE but also curtails the potential benefits of PBE. Using a querybased approachthat is informed by the highly successful systems strategy of serviceoriented architecture (SOA), these barriers can be overcome, suggesting the concepts of model interchangeability, interoperability, and integration.We illustrate the proposed methodology on several geometric interoperability tasks commonly arising in manufacturing. C. Hoffmann, V. Shapiro, V. Srinivasan
Technical Report, Spatial Automation Laboratory, University of WisconsinMadison. Also available as a Purdue University technical report CSD TR #11015. 

Fourier Collision Detection
In this paper, we investigate a new approach to narrowphase collision detection for rigid objects based on the Fourier transform. This new collision test scales with respect to accuracy, and we are able to rigorously establish an upper bound on the error of our test relating the Hausdorfi metric to the number of Fourier coeffients used. Because our new form of the collision test is also a smooth inequality, it can be used as a holonomic unilateral constraint in many applications, such as path planning, rigid body dynamics, nesting or tool placement, replacing the need for more adhoc normal/contact based constraint solvers. Moreover, we also show how this constraint can be directly differentiated via Fourier multipliers with only a constant factor overhead, which leads to a simple method for constructing a Jacobian for both normal forces and rotational torques. M. Lysenko
Technical Report, Spatial Automation Lab, University of WisconsinMadison 

Geometric Issues in Computer Aided Design/Computer Aided Engineering Integration
The longstanding goal of computer aided design (CAD)/computer aided engineering (CAE) integration demands seamless interfaces between geometric design and engineering analysis/simulation tasks. The key challenge to this integration stems from the distinct and often incompatible roles geometric representations play, respectively, in design and analysis. This paper critically examines and compares known meshbased and meshfree approaches to CAD/CAE integration, focusing on the basic tasks and components required for building fully integrated engineering applications. For each task, we identify the fundamental requirements and challenges and discuss how they may be met by known techniques and proposed solutions. V. Shapiro, I. Tsukanov, A. Grishin
Journal of Computing and Information Science in Engineering, Volume 11, June 2011 

Rapid Mapping and Exploration of Configuration Space
We describe a GPUbased computational platform for sixdimensional configuration mapping, which is the description of the configuration space of rigid motions in terms of collision and contact constraints. The platform supports a wide range of computations in design and manufacturing, including three and six dimensional configuration space obstacle computations, Minkowski sums and differences, packaging problems, and sweep computations. We demonstrate dramatic performance improvements in the special case of configuration space operations that determine interferencefree or containmentpreserving configurations between moving solids. Our approach treats such operations as convolutions in the six dimensional configuration space that are efficiently computed using the Fast Fourier Transform (FFT). The inherent parallelism of FFT algorithms facilitates a straightforward implementation of convolution on GPUs with existing and freely available libraries, making all such six dimensional configuration space computations practical, and often interactive. S. Nelaturi, M. Lysenko, V. Shapiro
2011 ASME International Design and Engineering Technical Conference, IDETC/CIE 2011, Washington DC 