Computing within Materials: Self-Adaptive Materials and Self-organizing Agents

From Macroscale to Microscale Computing

Stefan Bosse
University of Bremen, Dept. Mathematics & Computer Science
University of Koblenz, Fac. Computer Science
11.4.2018

Overview

Material Informatics: Material-integrated Computing


Smart Materials

Smart Materials: Fusion of Sensorial- and Adaptive Materials with Information Processing

Sensorial Material

A material or structure with integrated sensing and data processing (ICT)

Adaptive Material

A material or structure with integrated sensing, data processing and actuation that can control and change material or structure properties

figadapStructMcEvoyMcEvoy, 2015 [1]

Smart Materials

In our understanding, a smart material provides the following major features:

  1. Perception using various kinds of sensors, e.g., measuring of strain, displacement, temperature, pressure, forces;

  2. Changing of local material and structure properties by actuators, e.g., stiffness or damping variation;

  3. Integrated Information and Communication Technologies (ICT);

  4. Distributed Approach: Local sensor processing and actuator control - Global cooperation and coordination.

  5. Robustness and Self-Organization

Smart Materials

figsmartMat


Fig. 1. (Left) General architecture of a Sensorial- and Adaptive Material = Smart Material (Right) Functional Decomposition: Sensing, Acting, Processing ↔ Data + Instruction Streams

Material Informatics: Computing within Materials

  • Traditionally computation is separated from sensing and control

  • Smart Materials poses the tight coupling of computation, communication, sensing, and control with loosely coupled nano computers

  • Algorithmic scaling and distribution are required

    figscaling

Computing Power and Efficiency

  • A normalized computing efficiency of a computer (considering only the data processing unit) can be defined by
(1)
\[\varepsilon  = \frac{C}{{AP}}
\]
A
Chip Area in square millimetres
C
Computing Power in (Integer) Mega Instructions Per Second (MIPS)
P
Electrical Power in Watt
  • The computing efficiency can be used to compare different computers and devices, i.e., giving a scaling factor:
    s=εx/εy

Technologies

Existing “Nano” Computers

  • Smart Dust Millions of loosely coupled nano computers, e.g.,
    • embedded in materials
    • scattered on surfaces
    • dispersed in liquids, foils, ..
    • about 10mm3 volume

Micro Mote M3

figm3

ELM System

figelm

Comparison of Computers

Micro Mote (M3 ) ELM System Atmel Tiny 20 Freescale KL03 ARM Cortex Smart Phone
Processor Arm Cortex M0 C8051F990 (SL) AVR Arm Cotex M0+ Arm Cortex A9
Clock 740kHz max. 32kHz - 48MHz 1GHz
CPU Chip Area 0.1mm2 9mm2 1mm2 4mm2 7mm2/ROM
Sensors Temperature - - - Temp, Light, Sound, Accel., Press., Magn.
Communication 900MHz radio, optical optical electrical - 3G/4G, WLAN, USB, Bluetooth, NFC
Harvester, Battery Solar cell, Thin film Solar cell, Coin - - -
Power Consumption 70mW / CPU 160mW / CPU 20mW 3mW @ 48MHz 100mW avg.,
Manufacturing 180nm CMOS - - - 40nm CMOS
Package Wire bonded Silicon Stack PCB Single Chip Single Chip
Computing Eff. ε 150 0.02 0.6 4.0 0.53

Optimization and Adaptation


Physical Model

  • A component composed of a Smart Material consists of:

    • Mass (body) nodes or material regions
    • Interconnection elements with parameterizable properties
Multi-body Physics Model

Solid body mass nodes with springs connecting nodes:

  • A spring is an actuator with two optimization target variables:
    Stiffness s, Damping d

  • Each spring is a strain sensor delivering the sensor value σ for the computation of the observation variable

figmbp

Optimization Goals

Reduction of global and local stress, strain, or forces of arbitrariely shaped components under varying load situations



Control and Optimization Cycle

  1. Perception using sensors

  2. Comparison of local and global observation variables

  3. Modification of material parameters by actuators


Algorithms

Global

  • A global observation variable x is used to compute a correction of the target variable s

  • The correction function kx uses the ratio of the local and the global observable

do with x  {ε,σ,U}
 X:=0; n do X := X + xn
 X := X/||
 n do
  rn := kx(xn/X, sn) 
  sn := sn * rn
until |Err| < Err0

Segment

  • Network is partitioned in segments
  • Observation variable is computed for each segment
  • Target variable is computed for each segment
do with x  {ε,σ,U}
 Si𝕊 do
  Xs,i:=0;
  nSi do Xs,i := Xs,i + xn
  Xs,i := Xs,i / |Si|
  nSi do
    rn := kx(xn/Xs,i, sn)
    sn := sn * rn
until |Err| < Err0

Neighbour

  • Neighbour node negotiation
  • Observation variable is limited to node boundary
  • Swapping of increments of target variable
do with x  {ε,σ,U}
 {ni,nj  N | ij 
  |pos(ni)-pos(nj)|=1} do
  if xi/xj < 1 
     si-Δs > s0 
     sj+Δs < s1 then
    si := si - Δs
    sj := sj + Δs
  end if
always

Multi-Agent Systems


Sensor Data Distribution

  • Each node agent sends its sensor values to neighbour nodes within a range of radius 1

  • This completes the set of sensor each node requires

  • Remote signals are used for sensor distribution using &Delta-distance routing

figROIaction1


Fig. 2. Sensor Distribution by neighbour nodes using using remote signals (or tuples)

Distributed Observable Computation

  • Global observable computation in a large-scale distributed network is expensive and difficult (failures)

    • Random walk and directed diffusion can be used to approximate a global observable
  • Segmented approach requires network segmentation

    • Without central instance difficult;

    • Instead a floating segment window is placed around each node

    • Each node has observable from all direct neighbours (North, South, West, East, Up, Down)

    • A chained distribution of data is used in each segment (N nodes N segments!)

    • Observable values with distance r=1 and r=2 are collected by each node to compute region observable

Distributed Observable Computation

Distance r=1

Observable from direct neighbours

Distance r=2

Direct neighbours deliver also values form their neighbours (opposite to request direction)

figROIaction2


Fig. 3. Distribution of neighbour observable values for region accumulation using remote signals (or tuples)

Distributed Adaptation

  • The global and the segment algorithm need no further distribute coordination
  • After the region observable (global or segment) is computed the actuators (springs) of each node can be modified basing on ratio of the local and the region observable value
  • Neighbour negotiation approach do not require a region observable

figROIAction3


Fig. 4. Negotiation is used between two neighbour nodes to achieve a stiffness reconfiguration

Simulation

Multi-domain Simulation Framework

figmultisimFlow


Fig. 5. Hybrid simulation environment with Abaqus, Matlab, and SEJAM2

Simulation Example

  • Device under Test: Plate (8x5x3 nodes), large hole, external load

MAS World

  • Event-based agent behaviour activates sensor processing, distribution, and adaptation only if something changes

Physical World

Simulation Results

  • Global and segment optimization achieves 40% decrease of total strain and maximum strain energy of mass elements using a linear correction function

  • Neighbour negotiation is simple but not as efficient as segment approach

figresultsAdapt

Conclusions and Outlook

  • Smart Materials poses the tight coupling of computation, communication, sensing, and control with loosely coupled nano computers

  • Algorithmic scaling and distribution are required

  • Distributed information processing paradigma: Multi-agent Systems

  • Multi-domain simulation enables the development and evaluation of different optimization strategies for smart adaptive materials and structures

  • The SEJAM simulator enables the simulation and analysis of coupled physical and computational systems

  • Global and segment optimization achieves 40% decrease of total strain and maximum strain energy of mass elements using a linear correction function

References

[1] M. A. McEvoy and N. Correll, “Thermoplastic variable stiffness composites with embedded, networked sensing, actuation, and control,” Journal of Composite Materials, vol. 49, no. 15, 2015.