Axiomatic Spatial Registries, Discrete Voxel Lattice Kinematics, and Multi-Asset Interoperability in Urban Digital Twins

1. Socio-Spatial Friction & Abstract

The transition from continuous architectural manufacturing to autonomous, macro-scale discrete robotic assembly faces profound computational, mechanical, and spatial friction. Expanding localized unit voxel placements to urban-scale infrastructural domains inherently risks severe mathematical drift, manifesting as dynamic instability under nonlinear geometric deformations, matrix inversion singularities associated with recursive topographical transformations, and catastrophic collision vectors within unconstrained multi-asset robotic swarms operating in dynamically expanding, occluded topological spaces.

To resolve these systemic bottlenecks, the architectural logic imposes a mathematically uncompromising Cartesian foundation anchored to a rigid global depth datum, systematically guaranteeing positive elevation indexing to prevent computational collapse during massive global coordinate transformations. This highly formalized spatial registry couples with a bidirectional digital twin utilizing static condensation protocols to compress complex voxel sub-assemblies into mathematical super-elements, enabling real-time structural finite element analysis and phase-measuring profilometry. Bipedal inchworm robotics maneuver through this rigorously deterministic environment, executing hybrid discrete-continuous path planning and payload distribution optimizations pre-validated by time-delayed deep reinforcement learning networks prior to physical execution.

The eventual evolutionary convergence of this microscopic discrete kinematics with macro-scale Earth observation telemetry points toward a paradigm of autonomously metabolizing urban infrastructures. Continuous satellite monitoring of fluvial geomorphology and thermodynamic degradation will natively trigger generative algorithmic mesh-to-voxel pipelines, automatically deploying decentralized robotic swarms to iteratively assemble, reconfigure, or dismantle complex structural topographies in direct, zero-latency response to planetary environmental stimuli without centralized human orchestration.

2. Systemic Invariants & Tokens

ENTITY: Modular_Inchworm_Lattice_Assembler LAYER: Civil_Infrastructure LOGIC: Hybrid_Discrete_Continuous_Path_Planning RESOURCE: Polylactic_Acid_Voxel_Lattices

3. Parametric Operational Envelopes

4. Axiomatic Foundations

• Universal Cartesian Coordinate Anchor: All spatial, geometric, architectural, and topological variables are bound strictly to a global $Z$-axis foundation datum defined at exactly $Z = -500\text{ m}$. This absolute baseline serves as the mathematical origin for the digital twin spatial registry and coordinate transformation operations, ensuring universally positive structural elevation parameters to prevent negative array indexing, negative $Z$-coordinate matrix inversion singularities, and computational friction across all spatial derivations and path-planning vectors.
• Deterministic Discrete Voxel Latticing: The structural framework rejects continuous manufacturing paradigms in favor of volumetric, digital voxel units that function as mathematically deterministic digital materials. These lattices are formalized locally via 3D finite element analysis—employing Euler-Bernoulli beam theory or Timoshenko beam models based on the slenderness ratio—and aggregate into macro-geometries that exhibit stretch-dominated mechanical behavior governed by specific scaling relationships between relative stiffness ($E/E_ { s }$) and relative density ($\rho/\rho_ { s }$).
• Kinematic and Locomotion Operational Modes: The Modular Inchworm Lattice Assembler (MILAbot) traversal is strictly governed by the Multi-Body Dynamics for Inching-Locomotion Caterpillar Robots (MBD-ILAR) framework. Locomotion dynamics are constrained by the discrete topological nature of the lattice graph, requiring joint operations to alternate deterministically between Double Locked (DL-Mode), Single Unlocked (SU-Mode), and Locked Unlocked (LU-Mode) across a cyclic six-step peristaltic sequence mapped to the global datum.
• Zero-Tolerance Bidirectional State Synchronization: Real-time digital twin telemetry imposes a non-negotiable structural tolerance threshold where the vector disparity between the commanded state matrix and the actualized state matrix ($\Delta S = ||S_ { c m d } - S_ { a c t }||$) cannot exceed 0.5 of the lattice pitch. Any deviation beyond this invariant boundaries forces an autonomous interruption of the execution thread to trigger reactive realignment sequences and halt cascading physical assembly failures.

5. Field Equations & Analytical Calculus

The mathematical operationalization of the discrete robotic assembly system and its multi-scale digital twin synchronization framework requires a continuous mapping between local microscopic structural mechanics and global macroscopic geospatial reference frames. To decouple computational transformations from negative array indexing and prevent matrix inversion singularities during recursive spatial transformations over vast topographies, a universal Cartesian coordinate anchor is formalized at a deep subterranean vertical datum.

Let the global coordinate system be denoted as $\mathcal{G} = {X, Y, Z}$ and an arbitrary local coordinate frame embedded within a voxel element, robotic end-effector, or photogrammetry node be denoted as $\mathcal{L} = {x, y, z}$. The absolute spatial registry mapping a local point $p_ { \text{local }}$ to its global coordinate position $P_ { \text{global }}$ relative to the base datum is governed by an orthogonal rotation matrix $R$ and a translation vector $v_ { t }$, where $v_ { t }$ incorporates a static offset mapping the absolute physical origin to the vertical Z-axis foundation datum established at $Z = -500\text{ m}$.

$$ P_ { \text{global }} = R \cdot p_ { \text{local }} + v_ { t } $$

The translation vector $v_ { t }$ enforcing arithmetic consistency and strict positivity across the global digital twin spatial registry is defined by the column matrix expression:

$$ v_ { t } = \begin{bmatrix} x_ { 0 } \ y_ { 0 } \ z_ { 0 } - (-500) \end{bmatrix}^ { T } $$

At the structural micro-scale, individual volumetric struts within the discrete voxel lattices are modeled via 3D finite element analysis (FEA) based on their slenderness ratio using Euler-Bernoulli or Timoshenko beam theories. The local stiffness matrix $K_ { \text{loc }}$ for a single 3D voxel strut is a symmetric $12 \times 12$ matrix partitioned into $6 \times 6$ sub-matrices. The primary sub-matrix $K_ { 1 1 }$ captures the axial, torsional, and bending degrees of freedom at the primary node, derived directly via Hamilton's principle:

$$ K_ { 1 1 } = \begin{bmatrix} \frac{EA}{L} & 0 & 0 & 0 & 0 & 0 \ 0 & \frac{12EI_ { z }}{L^ { 3 }} & 0 & 0 & 0 & \frac{6EI_ { z }}{L^ { 2 }} \ 0 & 0 & \frac{12EI_ { y }}{L^ { 3 }} & 0 & -\frac{6EI_ { y }}{L^ { 2 }} & 0 \ 0 & 0 & 0 & \frac{GJ}{L} & 0 & 0 \ 0 & 0 & -\frac{6EI_ { y }}{L^ { 2 }} & 0 & \frac{4EI_ { y }}{L} & 0 \ 0 & \frac{6EI_ { z }}{L^ { 2 }} & 0 & 0 & 0 & \frac{4EI_ { z }}{L} \end{bmatrix} $$

where $E$ represents the material's Young's modulus, $G$ is the shear modulus, $A$ denotes the cross-sectional area, $L$ is the precise length of the discrete strut, $I_ { y }$ and $I_ { z }$ are the principal moments of inertia about the local $y$ and $z$ axes, and $J$ represents the polar moment of inertia. The off-diagonal block elements ($K_ { 1 2 }$, $K_ { 2 1 }$) and the secondary block ($K_ { 2 2 }$) enforce equilibrium across the discrete element interfaces.

When structural boundaries experience extreme displacement or severe in-plane loading, linear elastic approximations fail. The dynamic equations of motion must append geometric nonlinearity terms ($K_ { g }$, representing the work done by in-plane loads) and the fully nonlinear stiffness matrix ($K_ { \text{nl }}$):

$$ M\ddot{u} + (K + K_ { g } + K_ { \text{nl }})u = f_ { \text{ext }} $$

Within this formulations, $M$ is the symmetric mass matrix, $u$ is the vector of nodal displacements, $\ddot{u}$ defines the nodal acceleration vector, and $f_ { \text{ext }}$ represents the external force vectors. To facilitate real-time computational execution over millions of degrees of freedom, the system executes static condensation via Guyan reduction to compress complex voxel sub-assemblies into lower-dimensional super-elements. Partitioning the static equilibrium equation into active degrees of freedom $u_ { a }$ and condensed degrees of freedom $u_ { c }$ yields:

$$ \begin{bmatrix} K_ { a a } & K_ { a c } \ K_ { c a } & K_ { c c } \end{bmatrix} \begin{bmatrix} u_ { a } \ u_ { c } \end{bmatrix} = \begin{bmatrix} F_ { a } \ F_ { c } \end{bmatrix} $$

Assuming zero external forces acting directly on the condensed nodes ($F_ { c } = 0$), the dependent degrees of freedom are isolated explicitly as $u_ { c } = -K_ { c c }^ { - 1 }K_ { c a }u_ { a }$. Substituting this constraint back into the active partition yields the condensed matrix equation $K_ { \text{reduced }}u_ { a } = F_ { a }$, where the highly reduced stiffness matrix is isolated by:

$$ K_ { \text{reduced }} = K_ { a a } - K_ { a c }K_ { c c }^ { - 1 }K_ { c a } $$

To assemble the macro-structure global system matrix $K_ { \text{sys }}$, each elemental condensed matrix must undergo a congruent transformation using direction cosines. The local-to-global rotational transformation block $\lambda$ maps the localized axes to the global framework based on the intervening angles $\theta_ { i , J }$:

$$ \lambda = \begin{bmatrix} \cos\theta_ { x X } & \cos\theta_ { x Y } & \cos\theta_ { x Z } \ \cos\theta_ { y X } & \cos\theta_ { y Y } & \cos\theta_ { y Z } \ \cos\theta_ { z X } & \cos\theta_ { z Y } & \cos\theta_ { z Z } \end{bmatrix} $$

The full $12 \times 12$ orthogonal block diagonal transformation matrix $T$ is constructed directly from $\lambda$:

$$ T = \begin{bmatrix} \lambda & 0 & 0 & 0 \ 0 & \lambda & 0 & 0 \ 0 & 0 & \lambda & 0 \ 0 & 0 & 0 & \lambda \end{bmatrix} $$

The global stiffness matrix for an individual voxel element, $K_ { \text{glob }}^ { e }$, is rigorously established via the congruent transformation:

$$ K_ { \text{glob }}^ { e } = T^ { T } K_ { \text{reduced }}^ { e } T $$

The global structural topology is then solved using linear superposition across a sparse degree-of-freedom array to form the system-level equilibrium equation $K_ { \text{sys }}U_ { \text{sys }} = F_ { \text{sys }}$. Real-time simulation of structural displacement, yield thresholds, and localized friction points is resolved via inversion: $U_ { \text{sys }} = K_ { \text{sys }}^ { - 1 }F_ { \text{sys }}$.

For high-compliance soft-robotic locomotion within viscous operational fields, continuous continuum dynamics replace rigid-body linkage assumptions. Incorporating a constant volume hydrostatic skeleton constraint, the transverse displacement $w(x,t)$ of the soft continuum string is modeled via the spring-damper field equation:

$$ \rho\frac{\partial^ { 2 } w}{\partial t^ { 2 }} + c\frac{\partial w}{\partial t} - T\frac{\partial^ { 2 } w}{\partial x^ { 2 }} = 0 $$

where $\rho$ defines the linear mass density, $c$ is the viscous damping coefficient, and $T$ indicates the internal structural tension.

Path planning across the discrete topological lattice maps these continuum joint constraints to an optimized discrete graph search. The cost optimization function guiding the swarm navigation relies on an $A^*$ formulation optimized via Manhattan distance metrics to correctly assess orthogonal traversal overheads:

$$ f(n) = g(n) + h(n) $$

where $g(n)$ compiles the exact accumulated path cost from the origin to the current voxel node $n$, and $h(n)$ provides the estimated Manhattan distance from node $n$ to the target destination coordinates.

When simulating macro-scale urban interactions and infrastructure wind-loading profiles, the digital twin framework incorporates aerodynamic damping constraints. For an assembled voxel structure simplified to a representative mass-spring system with stiffness $K$, mass $M$, and cross-sectional area $S$ exposed to an oncoming fluid flow velocity $U$, the aerodynamic drag force $F_ { d }$ and the power input to the system oscillating with velocity amplitude $u^*$ are defined by:

$$ F_ { d } = C_ { d } S \frac{1}{2}\rho_ { a } U^ { 2 } $$

$$ P_ { \text{input }} = \frac{1}{2} \text{Re}(F_ { d } \cdot u^*) $$

Here, $C_ { d }$ is the aerodynamic drag coefficient and $\rho_ { a }$ is the density of the air. When the driving structural oscillation occurs below the fundamental resonance frequency ($\omega < \sqrt{K/M}$), the power input yields a negative real value, mathematically proving the manifestation of critical flow damping within the spatial infrastructure lattice.

6. Algorithmic State Imperatives

1. Initialize Global Spatial Registry Coordinate Datum: Establish a universal Cartesian coordinate anchor at an absolute mathematical baseline origin of $(0,0,0)$ where $Z = -500\text{ m}$ relative to the computational datum, mathematically forcing all subsequent spatial derivations, path-planning vectors, and topographical flow models to maintain a universally positive structural elevation parameter to eliminate negative Z-coordinate matrix inversion singularities during recursive spatial transformations.
2. Execute Local Voxel Stiffness Sub-Matrix Reduction: Compute the 3D finite element analysis (FEA) local stiffness matrix ($K_ { l o c }$) for an 8-node complex voxel with 48 degrees of freedom based on Hamilton's principle using Euler-Bernoulli or Timoshenko beam models, and apply static condensation via Guyan reduction ($K_ { r e d u c e d } = K_ { a a } - K_ { a c }K_ { c c }^ { - 1 }K_ { c a }$) assuming zero external forces on condensed nodes ($F_ { c } = 0$) to isolate active degrees of freedom ($u_ { a }$) and output an equivalent 2-node, 12-degree-of-freedom 3D beam element profile.
3. Assemble Sparse Global System Stiffness Matrix: Map the condensed local stiffness matrix ($K_ { r e d u c e d }^ { e }$) of each individual voxel element into the global coordinate system ($\mathcal{G}$) via the congruent transformation $K_ { g l o b }^ { e } = T^ { T } K_ { r e d u c e d }^ { e } T$ using an orthogonal block diagonal transformation matrix ($T$) structured from directional cosines ($\lambda$), then perform linear superposition of all elemental matrices into a sparse array to resolve the governing static equilibrium equation $U_ { s y s } = K_ { s y s }^ { - 1 }F_ { s y s }$.
4. Trigger Peristaltic Gait Cycle State Transitions: Coordinate the bipedal locomotion of the under-actuated relative robotic swarm platform (MILAbot) across the discrete topological lattice graph by cyclically alternating joint operational modes between Double Locked (DL-Mode), Single Unlocked (SU-Mode), and Locked Unlocked (LU-Mode) through a strict six-step sequence governed by the locomotion operator $T(c_ { i })$ while continuously validating that the gripping force matches the moment balance constraint $M_ { r e s i s t } \ge \sum(m_ { j } \cdot g \cdot \Delta x_ { j })$ against overturning payload masses.
5. Enforce Asynchronous WebSocket Telemetry Error-Correction: Establish a zero-latency bidirectional digital twin synchronization loop by offloading path-planning and inverse kinematics calculations from the physical SAMD51 microcontrollers to an external high-performance processing cluster via asynchronous WebSocket Python middleware, streaming actual joint states ($S_ { a c t }$), evaluating vector disparity against commanded states ($\Delta S = ||S_ { c m d } - S_ { a c t }||$), and autonomously terminating execution threads to trigger reactive realignment sequences if $\Delta S$ exceeds $0.5$ of the lattice pitch.

7. Kinematic Validation Protocols

• Absolute Spatial Registry and Base Datum Calibration: Real-time validation of local coordinate transformation vectors mapped from individual voxel elements, robotic end-effectors, photogrammetry nodes, and urban infrastructure meshes back to the global coordinate origin. Telemetry scripts must continuously verify that all path-planning vectors, structural stiffness matrices, and topographical flow models possess a universally positive structural elevation parameter based on the reference baseline to eliminate negative vertical axis coordinate values during recursive spatial transformations.
• Voxel Lattice Structural FEA and Local Stiffness Matrix Integrity: Stress testing and structural monitoring of discrete voxel lattices under static and dynamic loading. Individual struts must be validated using 3D finite element analysis to enforce mathematical boundary consistency across multi-scale nodes, dynamically switching between Euler-Bernoulli beam models for high slenderness ratios and Timoshenko beam models to capture shear effects in thick-framed elements.
• Nonlinear Dynamics and Guyan Reduction Verification: Computational tracking of geometric nonlinearity during severe in-plane loading or extreme structural displacement. Monitoring routines must track the execution of specialized eigenvalue solvers handling asymmetric matrices for modal analysis while verifying that real-time static condensation algorithms successfully compress complex voxel sub-assemblies from 48 degrees of freedom down to 12-degree-of-freedom equivalent 3D beam elements without destabilizing boundary stiffness profiles.
• Global System Stiffness Assembly and Static Equilibrium Tracking: Verification of the linear superposition of elemental global stiffness matrices mapped into a sparse degree-of-freedom array. The digital twin must continuously solve the global governing static equilibrium equation via matrix inversion to monitor localized structural friction points, structural deflection under gravitational load, and yield threshold proximity.
• Compounded Voxel Block Mechanical and Volumetric Throughput Validation: Lifecycle tracking of hierarchical discrete lattice assembly algorithms aggregating unit voxels into self-aligning compounded blocks. Protocols must verify stretch-dominated mechanical scaling behaviors, the structural integrity of reversible screw-released snap-fit connectors during vertical robotic installation, and real-time volumetric throughput minimums during multi-axis load deflection testing.
• MILAbot Actuator and Gear Reduction Torque Monitoring: High-resolution telemetry tracking of the primary motive forces generated by brushless DC motors. Systems must continuously parse simple field-oriented control logic and instantaneous angular feedback from high-resolution magnetic encoders through dedicated microcontrollers to verify that the multi-stage planetary and spur gear reduction assemblies meet structural torque demands when cantilevering loaded voxel blocks.
• Differentiated Gripper Load Handling and Placement Error Evaluation: Physical validation of structural anchor grippers and payload grippers under reactive operational moments. Telemetry must confirm that high-torque servo-actuated structure grippers rigidly lock into lattice voids to withstand overturning moments, while low-mass micro servo-actuated payload grippers utilize compliant kinematic linkages to maintain mechanical alignment within structural placement error tolerances.
• Peristaltic Gait Cycle and Locomotion Operator Synchronization: Sequential verification of the multi-body dynamics for inching-locomotion caterpillar robot framework. The system must monitor a strict physical six-step locomotion sequence—incorporating anchoring, rotational joint lifting, central prismatic joint linear translation, descent, structural load transfer, and negative contraction—while validating that the real-time homogenous transformation matrices maintain rigid-body constraints and moment balances.
• Bionic Slope-Climbing Suction Stability Assessment: Static numerical simulation and physical telemetry verification of electrically actuated suction cups integrated into the robotic chassis. Environmental testing must confirm uniform anchoring force distribution and structural slip prevention across complex angular slopes and variable lattice geometries.
• Continuous Soft Continuum String Damping Calibration: Real-time monitoring of soft-robotic continuum segments moving through viscous operational environments. Validation scripts must evaluate time-varying stiffness and damping coefficients derived from the principle of virtual power, ensuring compliance with constant volume hydrostatic skeleton constraints and verifying that asymmetric friction transitions driven by off-centered irradiation achieve directional motion hysteresis.
• Hybrid Discrete-Continuous Path Planning and Collision Avoidance Verification: Algorithmic validation of the multi-degree-of-freedom inverse kinematics engine synchronized with 3D voxel space A* search routines. The navigation system must compute localized joint trajectories using Manhattan distance heuristics to prevent physical linkages from colliding with rapidly expanding lattice structures, cross-referencing continuous volume intersections at every intermediate state transition.
• Bidirectional Digital Twin State Telemetry and Realignment Interruption: Deterministic state monitoring of physical hardware streamed via asynchronous WebSocket protocols to high-performance processing clusters. The digital twin must compute the exact vector disparity between commanded and actualized state matrices, executing an autonomous thread interruption and reactive realignment sequence if the error vector magnitude exceeds the structural tolerance threshold.
• DRL-Driven Dynamic Photogrammetry Position Optimization: Machine learning reward function validation for the V-STARS photogrammetry system. The deep reinforcement learning agent must actively optimize camera layout positions within the simulation loop to maximize field-of-view coverage for 3D point cloud reconstruction, mitigating visual line-of-sight occlusions caused by expanding voxel structures while enforcing simulated spatial registry collision parameters.
• Phase-Measuring Profilometry and Non-Linear Least Squares Calibration: Virtual-to-physical sensor validation utilizing Fringe Projection Profilometry loops. The system must execute non-linear least squares optimization to minimize the mean difference between virtual gamma images and physical sensor data, optimizing interdependent variables including phase-shift counts, baseline camera-to-projector spacing, and fringe densities to prevent configuration conflicts.
• Multi-Asset Interoperability and Semantic Knowledge Graph Evaluation: Inter-device communication tracking executed via the Asset Administration Shell protocol across heterogeneous hardware platforms. Automated control mechanisms must continuously parse semantic knowledge graphs representing manufacturing components to dynamically determine optimal factory configurations and real-time capacity utilization without rigid control architectures.
• Self-Learning Virtual Commissioning and Time-Delayed Safety Deployment: Simulation batch validation of Bayesian optimization and Genetic Algorithms for robotic trajectory energy minimization. Learned path-planning behaviors must be processed through a strict time-delayed deployment buffering system to guarantee that physical implementations adhere entirely to pre-validated kinematic safety envelopes.
• Generative AI Voxelization and Geometric Constraint Verification: End-to-end validation of the natural language to physical realization pipeline. Voxelization algorithms must intercept continuous generative meshes to filter geometry through component connectivity matrices, vertical stacking rules, load path evaluations, overhang stability limits, and robot reachability constraints before compilation into build order sequences.
• Urban-Scale Geospatial Intelligence and Subterranean Topology Mapping: Coordination testing of macro-scale city digital twins integrating 3D object models, laser scans, and subterranean sandstone cavern datasets. Spatial registry verification must ensure all subterranean environments are seamlessly continuous with surface-level civil infrastructure simulation matrices.
• Gravitational Socio-Spatial Flow and Infrastructure Robustness Analysis: Large-scale computational validation of transportation flow, employment shift, and infrastructure load matrices across Middle-layer Super Output Areas. Models must ingest live IoT, traffic camera, and flood sensor telemetry to identify cascading spatial failures and friction points in energy and transport networks.
• Aerodynamic Damping and Carbon Sequestration Fluid Dynamics Simulation: Finite element analysis of constructed voxel infrastructure subjected to severe wind loading. The system must compute aerodynamic drag forces and power input to mass-spring structural oscillations to verify critical flow damping below fundamental resonance frequencies, balancing structural metrics against urban carbon sequestration datasets.
• Earth Observation Fluidic and Thermodynamic Closed-Loop Validation: Multi-spectral automated water mask validation using satellite imagery constellations. The system must extract satellite-derived river surface temperatures and discharge rates to map fluvial geomorphology changes, autonomously triggering generative design pipelines, voxel discretization, and asset orchestration swarms to assemble protective infrastructure.

8. Geometric Voxel Assembly Matrix

9. Node Registry Payload JSON

{
  "node_id": "CIRG-FND-SEC-0001",
  "silo_id": "FND-SEC",
  "registry_metadata": {
    "title": "Axiomatic Spatial Registries, Discrete Voxel Lattice Kinematics, and Multi-Asset Interoperability in Urban Digital Twins",
    "date": "2026-05-26"
  },
  "systemic_tokens": {
    "entity": "Modular_Inchworm_Lattice_Assembler",
    "layer": "Civil_Infrastructure",
    "logic": "Hybrid_Discrete_Continuous_Path_Planning",
    "resource": "Polylactic_Acid_Voxel_Lattices"
  },
  "spatial_registry": [
    {
      "primitive_id": "PRIM-FND-SEC-0001",
      "geometry": "Absolute Origin Base Anchor",
      "vectors": {
        "x": 0.00,
        "y": 0.00,
        "z": 0.00
      },
      "scale": "Baseline Mathematical Datum Anchor",
      "hex_color": "#000000"
    },
    {
      "primitive_id": "PRIM-FND-SEC-0002",
      "geometry": "$1\\times1\\times 1$ Voxel Unit Lattice",
      "vectors": {
        "x": 0.00,
        "y": 0.00,
        "z": 500.00
      },
      "scale": "$65\\text{ mm}$ pitch, Stretch-dominated linear yield",
      "hex_color": "#3A86FF"
    },
    {
      "primitive_id": "PRIM-FND-SEC-0003",
      "geometry": "$2\\times2\\times2$ Compounded Voxel Block",
      "vectors": {
        "x": 0.00,
        "y": 0.00,
        "z": 500.00
      },
      "scale": "$130\\text{ mm} \\times 130\\text{ mm}$, Capacity: 3445 N",
      "hex_color": "#FF006E"
    },
    {
      "primitive_id": "PRIM-FND-SEC-0004",
      "geometry": "$3\\times3\\times3$ Compounded Voxel Block",
      "vectors": {
        "x": 0.00,
        "y": 0.00,
        "z": 500.00
      },
      "scale": "$195\\text{ mm} \\times 195\\text{ mm}$, Capacity: 8712 N",
      "hex_color": "#8338EC"
    },
    {
      "primitive_id": "PRIM-FND-SEC-0005",
      "geometry": "$4\\times2\\times2$ Compounded Assembly Block",
      "vectors": {
        "x": 0.00,
        "y": 0.00,
        "z": 500.00
      },
      "scale": "Volumetric Throughput: $4,394,000\\text{ mm}^ { 3 }/\\text{min}$",
      "hex_color": "#FFBE0B"
    },
    {
      "primitive_id": "PRIM-FND-SEC-0006",
      "geometry": "MILAbot Structure Gripper Node",
      "vectors": {
        "x": 0.00,
        "y": 0.00,
        "z": 500.00
      },
      "scale": "$2\\times2\\times1$ grid engagement, FeeTech FS117",
      "hex_color": "#FB5607"
    },
    {
      "primitive_id": "PRIM-FND-SEC-0007",
      "geometry": "MILAbot Voxel Payload Gripper Linkage",
      "vectors": {
        "x": 0.00,
        "y": 0.00,
        "z": 500.00
      },
      "scale": "Petal compliant linkage, $\\pm0.5$ pitch tolerance",
      "hex_color": "#00F5D4"
    },
    {
      "primitive_id": "PRIM-FND-SEC-0008",
      "geometry": "Bionic Suction Cup Traversal Anchor",
      "vectors": {
        "x": 0.00,
        "y": 0.00,
        "z": 500.00
      },
      "scale": "Electrically actuated slope-climbing interface",
      "hex_color": "#70E000"
    },
    {
      "primitive_id": "PRIM-FND-SEC-0009",
      "geometry": "Urban Digital Twin Macro Mesh Node",
      "vectors": {
        "x": 0.00,
        "y": 0.00,
        "z": 500.00
      },
      "scale": "Nottingham City Geospatial Registry (73 sq km)",
      "hex_color": "#FF5733"
    },
    {
      "primitive_id": "PRIM-FND-SEC-0010",
      "geometry": "Subterranean Sandstone Cavern Element",
      "vectors": {
        "x": 0.00,
        "y": 0.00,
        "z": 0.00
      },
      "scale": "850 network nodes mapped at surface level datum",
      "hex_color": "#8B4513"
    },
    {
      "primitive_id": "PRIM-FND-SEC-0011",
      "geometry": "MSOA Gravitational Flow Tracker",
      "vectors": {
        "x": 0.00,
        "y": 0.00,
        "z": 500.00
      },
      "scale": "8,436 areas, 71 million distinct trip vectors",
      "hex_color": "#2ECC71"
    },
    {
      "primitive_id": "PRIM-FND-SEC-0012",
      "geometry": "PlanetScope CubeSat Water Mask Voxel",
      "vectors": {
        "x": 0.00,
        "y": 0.00,
        "z": 500.00
      },
      "scale": "River Twin framework Basin flow gauge",
      "hex_color": "#1ABC9C"
    }
  ]
}