Supplementary Materials for

Similar documents
Strain Rate and Temperature Effects on the Nonlinear Behavior of Woven Composites

Supplementary Information for Design of Bending Multi-Layer Electroactive Polymer Actuators

Modeling and Analysis of Thermal Bimorph Using COMSOL

STUDY OF THERMAL DIFFUSIVITY IN HEAT-INSULATING MATERIALS

Long-term Life Prediction of CFRP Structures Based on MMF/ATM Method

THERMAL ENDURANCE OF UNREINFORCED UNSATURATED POLYESTERS AND VINYL ESTER RESINS = (1) ln = COMPOSITES & POLYCON 2009

Measuring Temperature with a Silicon Diode

DESIGN OF THE DIE PROFILE FOR THE INCREMENTAL RADIAL FORGING PROCESS *

EVALUATION OF THERMAL CONDUCTIVITY IN PITCH- BASED CARBON FIBER REINFORCED PLASTICS

Uniaxial compressive stress strain model for clay brick masonry

Modelling of damage in composite materials using interface elements

Thermomechanics of Shape Memory Polymers and Composites

EFFECT OF SURFACE ASPERITY TRUNCATION ON THERMAL CONTACT CONDUCTANCE

Numerical Modeling of Self-Compacting Mortar Flow Using Discrete Element Method

LONGITUDINAL EFFECTS AND FOCUSING IN SPACE-CHARGE DOMINATED BEAMS. John Richardson Harris

Melting Temperature of Lithium, Sodium and Potassium at High Pressures

Easy Evaluation Method of Self-Compactability of Self-Compacting Concrete

Optimum Design of Assembled Cavity Dies for Precision Forging Process

2.141 Modeling and Simulation of Dynamic Systems Assignment #2

The Thermal Conductivity Theory of Non-uniform Granular Flow and the Mechanism Analysis

Description of the nonlinear viscoelastic and thixotropic properties of human blood by the modified model of Maxwell-Gurevich-Rabinovich

Incorporating strain gradient effects in a multi-scale constitutive framework for nickel-base superalloys

A DESIGN GUIDE OF DOUBLE-LAYER CELLULAR CLADDINGS FOR BLAST ALLEVIATION

Crystallization of Supercooled Liquid Elements Induced by Superclusters Containing Magic Atom Numbers Abstract: Keywords: 1.

DETERMINATION OF ELASTIC CONSTANTS CORNU S METHOD

POST-PERFORATION LENGTH AND VELOCITY OF KE PROJECTILES WITH SINGLE OBLIQUE TARGETS

ASSESSMENT OF THE VOLUME CHANGE BEHAVIOUR

Multiscale Entropy Analysis: A New Method to Detect Determinism in a Time. Series. A. Sarkar and P. Barat. Variable Energy Cyclotron Centre

Reading from Young & Freedman: For this topic, read the introduction to chapter 25 and sections 25.1 to 25.3 & 25.6.

Weak Form Equation-Based Finite Element Modeling of. Viscoelastic Asphalt Mixtures 1

Chapter 8 Deflection. Structural Mechanics 2 Dept of Architecture

Spine Fin Efficiency A Three Sided Pyramidal Fin of Equilateral Triangular Cross-Sectional Area

Buckling Behavior of 3D Randomly Oriented CNT Reinforced Nanocomposite Plate

SHAPE IDENTIFICATION USING DISTRIBUTED STRAIN DATA FROM EMBEDDED OPTICAL FIBER SENSORS

Chapter 8. Lecture Notes Dr. Rakhmad Arief Siregar Kolej Universiti Kejuruteraan Utara Malaysia

Random Vibration Fatigue Analysis with LS-DYNA

2. Electric Current. E.M.F. of a cell is defined as the maximum potential difference between the two electrodes of the

Thermal Cracking Simulations of Aged Asphalt Pavements using Viscoelastic Functionally Graded Finite Elements

REDUCTION OF FINITE ELEMENT MODELS BY PARAMETER IDENTIFICATION

Electrical conductivity enhanced dielectric and piezoelectric properties of ferroelectric 0-3 composites

Determination of fiber/matrix interface debond growth parameters from cyclic loading of single fiber composites

Department of Building Science, Tsinghua University, Beijing , China

Chapter 1 Introduction and Kinetics of Particles

Density and structure of undercooled liquid titanium

On alternative methods of determining Radius of Curvature using Newton s Rings set up

Journal of Engineering Science and Technology Review 11 (1) (2018) Research Article

A Pulley System Apparatus for a Laboratory Experience in Dynamics

Chemistry Department Al-kharj, October Prince Sattam Bin Abdulaziz University First semester (1437/1438)

MECHANICS OF MATERIALS

Experimental Design For Model Discrimination And Precise Parameter Estimation In WDS Analysis

Photoluminescence of Monolayer MoS 2 on LaAlO 3 and SrTiO 3 substrates

Chapter 4: Hypothesis of Diffusion-Limited Growth

Lecture #8-3 Oscillations, Simple Harmonic Motion

Extension of CSRSM for the Parametric Study of the Face Stability of Pressurized Tunnels

8.1 Force Laws Hooke s Law

ID 1163 CURE MONITORING OF THERMOSETTING RESIN COMPOSITES BY LACOMTECH DIELECTROMETRY

Feasibility Study of a New Model for the Thermal Boundary Resistance at an Interface of Solid Thin Films

Engineering Solid Mechanics

Physics 6A. Stress, Strain and Elastic Deformations. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

In the session you will be divided into groups and perform four separate experiments:

Comparison Studies on Dynamic Packaging Properties of Corrugated Paperboard Pads

The temperature dependences of electromechanical properties of PLZT ceramics

Current, Resistance Electric current and current density

DIRECT NUMERICAL SIMULATION OF DAMAGE PROGRESSION IN LAMINATED COMPOSITE PLATES USING MULTI-SCALE MODELLING

SIMULATION OF THE HEATING STEP WITHIN THE THERMOFORMING PROCESS USING THE FINITE DIFFERENCE METHOD

EVAPORATION EFFECT IN NONLINEAR PENETRATION OF HIGH ENERGY BEAM DRILLING

Model Fitting. CURM Background Material, Fall 2014 Dr. Doreen De Leon

Journal of Chemical and Pharmaceutical Research, 2014, 6(2): Research Article

Feature Extraction Techniques

Research on coupling theory and its experiments between torque turbulence and bearing load of multi-support-rotor systems

International Workshop SMART MATERIALS AND STRUCTURES

USEFUL HINTS FOR SOLVING PHYSICS OLYMPIAD PROBLEMS. By: Ian Blokland, Augustana Campus, University of Alberta

1 (40) Gravitational Systems Two heavy spherical (radius 0.05R) objects are located at fixed positions along

Simulation of Geomechanical Behavior during SAGD Process using COMSOL Multiphysics

EFFECT OF SHEAR STRESS ON FATIGUE LIFE OF COMPOSITE LAMINATES

External Transverse Direct Current Magnetic Field Effect On Optical Emission Of a Non-Thermal Atmospheric Pressure Argon Plasma Jet

MIT Amorphous Materials

The Solution of One-Phase Inverse Stefan Problem. by Homotopy Analysis Method

IDENTIFICATION OF STABILITY OF CONCRETE TUNNEL LINING USING COUPLED MODELING

Title. Author(s)Izumida, Yuki; Okuda, Koji. CitationPhysical review E, 80(2): Issue Date Doc URL. Rights. Type.

821. Study on analysis method for deepwater TTR coupled vibration of parameter vibration and vortex-induced vibration

ANALYTICAL INVESTIGATION AND PARAMETRIC STUDY OF LATERAL IMPACT BEHAVIOR OF PRESSURIZED PIPELINES AND INFLUENCE OF INTERNAL PRESSURE

Experimental Based Substructuring Using a Craig-Bampton Transmission Simulator Model

Modeling of Photovoltaic System in the Large Ocean-going Ship Based on Simulink

In this chapter we will study sound waves and concentrate on the following topics:

MECHANICS OF MATERIALS Design of a Transmission Shaft

Benchmarking of lamina failure tests from WWFE-I and WWFE-II with a three parameter micromechanics based matrix failure theory

SD Numerical Simulation Technique for Hydrodynamic Flow Gas- Solids Mixing Mantilla Núñez, Irla* 1 S. De Vicente C. 2

CHAPTER 5 THEORETICAL MODELING OF ELASTIC PROPERTIES

STUDY ON THE QUANTIFICATION AND CORRELATION OF GRADATIONAL SEGREGATION IN ASPHALT MIXTURES

3D SELF HEATING MODELING FOR ELECTRO- THERMAL CHARACTERISATION OF SiGe HBTs

Faraday's Law Warm Up

Influence of the Flux Pinning on Critical Current in HTc Superconductors with Nano-Sized Defects

Research in Area of Longevity of Sylphon Scraies

An analytical relation between relaxation time spectrum and molecular weight distribution

DETERMINATION OF ADSORTION LAYERS ON SILICON SORPTION ARTIFACTS USING MASS COMPARISON

All Excuses must be taken to 233 Loomis before 4:15, Monday, April 30.

2. A crack which is oblique (Swedish sned ) with respect to the xy coordinate system is to be analysed. TMHL

Data-Driven Imaging in Anisotropic Media

Nonlinear Evolutionary Mechanisms of Instability of Plane-Shear Slope: Catastrophe, Bifurcation, Chaos and Physical Prediction

Transcription:

advances.scienceag.org/cgi/content/full/3/4/e160890/dc1 Suppleentary Materials for Direct 4D printing via active coposite aterials Zhen Ding, Chao Yuan, Xirui Peng, Tiejun Wang, H. Jerry Qi, Martin L. Dunn This PDF file includes: Published 1 April 017, Sci. Adv. 3, e160890 (017) DOI: 10.116/sciadv.160890 section S1. Theroechanical properties of polyers section S. Additional results section S3. Geoetry of printed saples section S4. Modeling and siulation fig. S1. Theroechanical properties of elastoer (TangoBlack+) and glassy polyer SMP (VeroClear) as deterined by DMA. fig. S. Contribution in overall curvature fro CTE isatch strain and built-in copressive strain. fig. S3. A lattice structure (diensions in ) printed in an open configuration and then deployed into a copact configuration upon heating. fig. S4. Design of the printed saple in Fig. 3A. fig. S5. Design of the printed saple in Fig. 3B. fig. S6. Design of the printed saple in Fig. 3C. fig. S7. Design of the printed saple in Fig. 4A. fig. S8. Design of the printed saple in Fig. 4B. fig. S9. Design of the printed saple in Fig. 4C. fig. S10. Design of the printed saple in fig. S3. table S1. List of paraeters for constitutive odel.

section S1. Theroechanical properties of polyers We easured the storage odulus, tan delta, and theral strain versus teperature of both the SMP and the elastoer. These results were obtained with a layer printing tie of 68s, but differ insignificantly fro results with different printing ties. Upon heating the storage odulus of Veroclear decreases slightly until ~40 o C above which it reduces draatically. In contrast, the decrease of Tangoblack+ storage odulus has basically ceased after 10 o C (fig. S1A). The peak of tan delta in fig. S1B indicates the glass transition teperature (Tg) ~6 C and ~8 o C for Veroclear and Tangoblack+, respectively. Tangoblack+ has an alost linear theral strain ~1% fro roo teperature (5 o C) to 80 o C, while it is ~0.8% for Veroclear and nonlinear indicating a phase transition in that teperature range. The isatch theral strain between Veroclear and Tangoblack+ is less than 0.% in the easured teperature range. The residual theral strains of both Veroclear and Tangblack+ are alost zero after they were cooled back to roo teperature in a heating-cooling cycle. fig. S1. Theroechanical properties of elastoer (TangoBlack+) and glassy polyer SMP (VeroClear) as deterined by DMA. (A) Storage odulus versus teperature; (B) Tan delta versus teperature; (C) Theral strain in one heating (solid lines) and cooling (dash lines) cycle. section S. Additional results We quantify the ratio between the built-in copressive strain and CTE isatch strain towards bending curvature. The CTE isatch also contributes to the overall bending, but to a lesser extent, especially with the longer layer printing tie.

fig. S. Contribution in overall curvature fro CTE isatch strain and built-in copressive strain. Various shape actuations can be achieved via our direct 4d printing approach by caully configuring aterials in structure. Opposite to Fig. 4A, where the printed collapsed lattice structure deploys into an open configuration upon heating, a printed open lattice structure here draatically collapses into a copact one with a structural copression of ~6% and a lateral contraction of ~6% upon heating. Finite eleent siulations (the pic at the right botto of fig. S) show good agreeent with a predicted copression of 64% and lateral contraction of 7%. fig. S3. A lattice structure (diensions in ) printed in an open configuration and then deployed into a copact configuration upon heating. The structure has a size of 180 119 7 and a SMP/elastoer bilayer thicknesses of 0.4/0.4.

section S3. Geoetry of printed saples The designs of the printed saples are presented below. In the figures, all diensions are in. For Fig. 1C, two-layer flat strips of 80 5 0.6 with equal Tangoblack+ and Veroclear thickness layers fabricated with a layer printing tie of 68 s and a total printing tie of 66 in. For Fig. 1D, two-layer flat strips of 60 5 0.6 with equal Tangoblack+ and Veroclear thickness layers fabricated with a layer printing tie of 58 s and a total printing tie of 56 in. For Fig. A, pure Tangoblack+ or Veroclear flat strips of 60 8 1. fabricated with a series of layer printing ties (10, 13, 5, 37, 49, 58, 68 and 74 s) and total printing ties (11, 16, 31, 46, 60, 75, 90 and 104 in). For Fig. B, two-layer flat strips of 60 5 0.6 with equal Tangoblack+ and Veroclear thickness layers fabricated with a series of layer printing ties (10, 13, 5, 37, 49, 58, 68 and 74 s) and total printing ties (10, 13, 3, 34, 45, 56, 66, and 77 in). For Fig. C, two-layer flat strips of 60 5 t with equal Tangoblack+ and Veroclear thickness layers, where t (0.6, 1.,.4, 4.8 ) is the total thickness fabricated with a series of layer printing ties (10, 13, 5, 37, 49, 58, 68 and 74 s). The total printing ties for strip of 60 5 0.6 under layer printing ties of 10, 13, 5, 37, 49, 58, 68 and 74 s, for exaple, are 10, 13, 3, 34, 45, 56, 66, and 77 in, respectively. The total printing ties of other strips are basically proportional to their thicknesses. For Fig. D, two-layer flat strips of 60 5 1., where the Tangoblack+ volue (thickness) fraction varies with a series of values (.5, 0, 40, 50, 60, 80, 90, 95%) fabricated with a layer printing tie of 68 s and a total printing tie of 89 in. fig. S4. Design of the printed saple in Fig. 3A. Two-layer flat strips of 150 10 0.7, where Tangoblack+ (15 10 0.5 ) and Veroclear (15 10 0. ) segents alternate along the length fabricated with a layer printing tie of 68 s and a total printing tie of 70 in.

fig. S5. Design of the printed saple in Fig. 3B. Two-layer flat strips of 18 8 0.7 with Veroclear/Tangoblack+ thicknesses of 0./0.5 and an SMP separator is oriented at 45 o fabricated with a layer printing tie of 68 s and a total printing tie of 70 in. fig. S6. Design of the printed saple in Fig. 3C. Two-layer flat rings of 10/68.6- outer/inner diaeter and 1.- thickness, where Veroclear/Tangoblack+ has thicknesses of 0.4/0.8 and alternates along the circuference, fabricated with a layer printing tie of 68 s and a total printing tie of 89 in.

fig. S7. Design of the printed saple in Fig. 4A. Closed lattice with an overall size of 180 5. 7. Veroclear/Tangoblack+ (0.35/0.35 ) bilayers are printed perpendicular to the build tray and fabricated with a layer printing tie of 68 s and a total printing tie of 310 in.

fig. S8. Design of the printed saple in Fig. 4B. Closed lattice with an overall size of 19 6.3 6. Veroclear/Tangoblack bilayers are printed perpendicular to the build tray, have a series of thicknesses of 0.8/0.8, 0.6/0.6, 0.3/0.9, and 0./0.6, and are fabricated with a layer printing tie of 68 s and a total printing tie of 73 in. fig. S9. Design of the printed saple in Fig. 4C. Flat star-shaped bilayer structure consists of 6 identical equilateral triangles, each of which has outer/inner side length of 60/44.41. The

Veroclear/Tangoblack+ thicknesses are 0.88/1.77 and the structure is fabricated with a layer printing tie of 68 s and a total printing tie of 144 in. For Fig. 4D, printed flower consisting of five layers of petals with a gap of 0.5 between each. The petals in different layers are siilar in shape and size, except for the outer radius (4, 48, and 51 ). Each petal is a flat bilayer structure, and has a SMP/elastoer thicknesses of 0.3/0.3. The axiu petal width is less than 4.5. The elastoer is Tango+, while SMP is the digital aterials/coposites fro Pureverowhite, Veroyellow, Veroagenta, Veroblack and Verocyan. The printing tie for petals fro the botto layer to top layer is 5, 37, 49, 58 and 74 s, respectively, and the total printing tie for the flower is 35 in. Unlike the other saples, the flower is printed on the Stratasys J750 ultiaterial printer which supports the ultiple color aterials. In fig. S1, flat strips of ~0 7 1 (Veroclear and Tangoblack+) for storage odulus/glass transition easureents and ~5 8 0.6 for theral strain easureents. These saples are fabricated with a layer printing tie of 68 s and the total printing tie of 8 in for the forer and 66 in for the latter.

fig. S10. Design of the printed saple in fig. S3. Open lattice with an overall size of 180 119 7. Veroclear/Tangoblack bilayers are printed perpendicular to the build tray, have a thicknesses of 0.4/0.4, and are fabricated with a layer printing tie of 68 s and a total printing tie of 310 in. section S4. Modeling and siulation In order to describe the shape change behavior of the printed structures and reveal its echanis, we build a odel that accounts for the teperature-dependent theroechanical behavior of SMP aterial and the experientally-deterined built-in copressive strain in the elastoer. Multibranch odel is used to describe the theroechanical property of SMP aterial, in which one equilibriu branch and several theroviscoelastic nonequilibriu branches are arranged in parallel (4). Maxwell eleents are used in the nonequilibriu branches to represent the stress relaxation behavior of the aterial, and the total stress of the SMP aterial can be expressed as

Where n t t eq e S dt S ES es Enon exp ds 0 s 1 s T eq E S is the elastic odulus of the equilibriu branch, (S1) E non and are the elastic odulus and teperature dependent relaxation tie of the th nonequilibriu branch, es is the echanical strain calculated by subtracting the theral strain ST fro total strain. Based on the tie teperature superposition principle (TTSP), can be deterined by using the R relaxation tie at erence teperature shift where a T a shift T (S) R T is the teperature dependent shifting factor. The shifting factors can be calculated by cobining the Willias Landel Ferry (WLF) equation and Arrhenius-type equation (4). When the teperature is higher than the erence teperature T, the shifting factor can be expressed by using WLF equation C1 T T shift log a T, T T C T T (S3) where C 1, C and T are the aterial paraeters to be characterized. If the teperature is lower than the erence teperature T, the shifting factor can be expressed by Arrhenius-type equation: (S1) shift AF 1 1 c ln a T, T T k T T (S4) where A, Fc and k are the aterial constant, configurational energy and Boltzann s constant respectively. The DMA test results in the Suppleentary Inforation were used to identify the above paraeters eq non R including E S, E,, C 1, C and AFc/k. The storage odulus at 80 can be considered as the eq equilibriu odulus E S as the relaxation tie at this teperature in each nonequilibriu branch is inial. For the ulti-branch linear odel, the teperature dependent storage odulus E l E T and loss factor tan can be respectively represented as s E T E eq s T, loss odulus n non E T (S5A) 1 1 T

T T n non E El T 11 (S5B) El T tan E T (S5C) where is the test frequency. Eploying the nonlinear regression (NLREG) ethod (4), s non R E,, C 1, C and AFc/k can be deterined by fitting the storage odulus and tan curves shown in fig. S1A and B. The fitted paraeters are listed in table S1. For the elastoer, it is seen fro fig. S1A that the storage odulus is alost teperature-independent above roo teperature. So for siplicity, a linear elastic odel is used to describe the rubber aterial where E elas is the elastic odulus and int y E e total e theral e pr ing, (S6) elas elas elas elas elas total e elas, theral e elas, e are the total strain, theral strain ( elas T ) pr int ing elas and printing bulid-in strain of the elastoer, respectively.

table S1. List of paraeters for constitutive odel. Branch non R E (Pa) (s) Branch non R E (Pa) (s) 1 1.49E+08.00E-08 15 1.4E+08 97.9366 1.E+08 4.7E-07 16 1.1E+08 8849.19 3 1.3E+08 5.47E-06 17 1.41E+08 849.0 4 1.47E+08 5.89E-05 18 8189706 594.7 5 1.66E+08 0.000547 19 5681974 7.9E+04 6 1.89E+08 0.00454 0 1478539 65350.3 7.08E+08 0.03439 1 8031731.13E+05 8.51E+08 0. 171558 5.37E+06 9 1.78E+08 1 3 4830405 000000 10 1.44E+08 3.5059 4 1197657 8.54E+07 11 1.51E+08 9.451896 5 138314.00E+07 1 1.63E+08 30.3741 6 18.666 3.61E+08 13 1.6E+08 100 7 537188.00E+09 14 1.51E+08 315.367 eq E S 6E+06 AF (Pa) c/k -4000 C1 17.4 T ( ) C 66.35 S (/ ) 1.7E-4 E elas (Pa) 0.6E+06 elas (/ ).3E-4