DESIGN METHOD BASED ON THE CONCEPT OF EMERGENCE AND ITS APPLICATION

Similar documents
GENERATION OF DIVERSE DYNAMIC BEHAVIORS USING THE EMERGENT DESIGN SYSTEM

Optimal Shape and Topology of Structure Searched by Ants Foraging Behavior

ME 176 Final Exam, Fall 1995

Probabilistic Engineering Design

ME 176 Final Exam, Fall 1997

Stress and Strain State in the Segmental Linings during the Excavation by Tunnel Boring Machines

DEVIL PHYSICS THE BADDEST CLASS ON CAMPUS AP PHYSICS

Section 6: PRISMATIC BEAMS. Beam Theory

An Elastic Contact Problem with Normal Compliance and Memory Term

VISUAL PHYSICS ONLINE RECTLINEAR MOTION: UNIFORM ACCELERATION

Module 1. Energy Methods in Structural Analysis

Solution to 1-D consolidation of non-homogeneous soft clay *

6.4 A cylindrical specimen of a titanium alloy having an elastic modulus of 107 GPa ( psi) and

Development and Validation of a Two-Phase Model for Reinforced Soil by Considering Nonlinear Behavior of Matrix

Chemical softening and hardening of geomaterials in situ

ENGINEERING COUNCIL CERTIFICATE LEVEL ENGINEERING SCIENCE C103 TUTORIAL 4 - STRESS AND STRAIN

Biomaterial Scaffolds

Evaluation of pressure and bulk modulus for alkali halides under high pressure and temperature using different EOS

MATHEMATICAL MODELLING OF HUMAN LIMB IN MATLAB

OPTIMIZATION OF FLOWS AND ANALYSIS OF EQUILIBRIA IN TELECOMMUNICATION NETWORKS

CONSTRUCTION PROCESS NUMERICAL SIMULATION AND SEISMIC ASSESSMENT OF MALLORCA CATHEDRAL

Gravitational Poynting theorem: interaction of gravitation and electromagnetism

ULTIMATE BEARING CAPACITY OF FOOTING ON SANDY SOIL AGAINST COMBINED LOAD OF VERTICAL, HORIZONTAL AND MOMENT LOADS

Boundary Element Method Calculation of Moment Transfer Parameters in Slab-Column Connections

MODELLING BAMBOO AS A FUNCTIONALLY GRADED MATERIAL

Stresses in Sedimentary Strata, including Coals, and the Effects of Fluid Withdrawal on Effective Stress and Permeability

Optimal Mass Distribution Prediction for Human Proximal Femur with Bi-modulus Property

Illustrative exercises for the lectures

Automatic Behavior Generation of Artificial Fish Using Elasticity Momentum

FUZZY FINITE ELEMENT METHOD AND ITS APPLICATION

IN several centuries, human has tried to understand all

4. A Physical Model for an Electron with Angular Momentum. An Electron in a Bohr Orbit. The Quantum Magnet Resulting from Orbital Motion.

Semi-implicit Treatment of the Hall Effect in NIMROD Simulations

3D topology optimization using cellular automata

Web Resource: Ideal Gas Simulation. Kinetic Theory of Gases. Ideal Gas. Ideal Gas Assumptions

Nonexistence of Plane Symmetric Micro Model in Presence of Cosmological Constant in Barber s Modified Theory of General Relativity

Fig. 1. Circular fiber and interphase between the fiber and the matrix.

Parameters Identification of Equivalent Circuit Diagrams for Li-Ion Batteries

MODELLING FLUID FLOW IN STRESS SENSITIVE PETROLEUM RESERVOIR CONSIDERING FAULT REACTIVATION PROBLEM

Biomechanical model of the shear stress distribution in the femoral neck

Topology Optimization of Compliant Mechanism with Geometrical Advantage

Microscopic Momentum Balance Equation (Navier-Stokes)

Progress In Electromagnetics Research Symposium 2005, Hangzhou, China, August

arxiv: v1 [physics.comp-ph] 17 Jan 2014

Mechanics of Biomaterials

Classical Mechanics NEWTONIAN SYSTEM OF PARTICLES MISN NEWTONIAN SYSTEM OF PARTICLES by C. P. Frahm

Determination of Mechanic Resistance of Osseous Element through Finite Element Modeling

Towards Green Distributed Storage Systems

A. Idesman. Keywords: time integration, spurious oscillations, numerical dispersion

TOPOLOGY STRUCTURAL OPTIMIZATION USING A HYBRID OF GA AND ESO METHODS

Investigation on Ring Valve Motion and Impact Stress in Reciprocating Compressors

A matrix Method for Interval Hermite Curve Segmentation O. Ismail, Senior Member, IEEE

FREEWAY WEAVING. Highway Capacity Manual 2000 CHAPTER 24 CONTENTS EXHIBITS

Archetype-Blending Multiscale Continuum Method

Lecture 5 Charge Density & Differential Charge. Sections: 2.3, 2.4, 2.5 Homework: See homework file

On the Linear Threshold Model for Diffusion of Innovations in Multiplex Social Networks

BAYESIAN PREMIUM AND ASYMPTOTIC APPROXIMATION FOR A CLASS OF LOSS FUNCTIONS

A014 Uncertainty Analysis of Velocity to Resistivity Transforms for Near Field Exploration

Topology Optimization Mathematics for Design

A Guidance Law for a Mobile Robot for Coverage Applications: A Limited Information Approach

On the Nature of Coherent Turbulent Structures in Channel Bends: Burst-Sweep Orientations in Three-Dimensional Flow Fields

ENGINEERING COUNCIL DYNAMICS OF MECHANICAL SYSTEMS D225 TUTORIAL 3 CENTRIPETAL FORCE

Effect of Fractal Dimension of Fine Aggregates on the Concrete Chloride Resistance

Transmission lines using a distributed equivalent circuit

A study of jacking force for a curved pipejacking

Chapter 14 Thermal Physics: A Microscopic View

Sample Average Approximation for Stochastic Empty Container Repositioning

Strain accumulation in sand due to cyclic loading: drained cyclic tests with triaxial extension

FOCUS ON CONCEPTS Section 7.1 The Impulse Momentum Theorem

Hasna BenSaid Faculté des sciences de Gafsa Campus Universitaire Sidi Ahmed Zarroug 2112 Gafsa, Tunisie

SUPPLEMENTARY MATERIAL. Authors: Alan A. Stocker (1) and Eero P. Simoncelli (2)

Theory of Network Communication

Post Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method

E : Ground-penetrating radar (GPR)

DO PHYSICS ONLINE. WEB activity: Use the web to find out more about: Aristotle, Copernicus, Kepler, Galileo and Newton.

Scalar multiplication and algebraic direction of a vector

Finite Element Simulation of Bar-Plate Friction Welded Joints Steel Product Subjected to Impact Loading

Noise constrained least mean absolute third algorithm

THE ANALYSIS OF THE CONVECTIVE-CONDUCTIVE HEAT TRANSFER IN THE BUILDING CONSTRUCTIONS. Zbynek Svoboda

Access from the University of Nottingham repository:

On Calculation of Turbulent Free-surface Flows with Particle Method

A possible mechanism to explain wave-particle duality L D HOWE No current affiliation PACS Numbers: r, w, k

The number of marks is given in brackets [ ] at the end of each question or part question. The total number of marks for this paper is 72.

1. A pure shear deformation is shown. The volume is unchanged. What is the strain tensor.

STANDARD SAMPLE. Reduced section " Diameter. Diameter. 2" Gauge length. Radius

v v Downloaded 01/11/16 to Redistribution subject to SEG license or copyright; see Terms of Use at

Peak Strain and Displacement Sensors for Structural Health Monitoring

Non-linear consolidation of soil with vertical and horizontal drainage under time-dependent loading

An Optimal Split-Plot Design for Performing a Mixture-Process Experiment

A spectral Turán theorem

Solid Mechanics Homework Answers

On resilience of distributed routing in networks under cascade dynamics

MATLAB SYMBOLIC COMPUTATION FOR THE STEADY STATE MODELING OF SYMMETRICALLY LOADED SELF EXCITED INDUCTION GENERATOR. Gurung K., Freere P.

Numerical Modelling of Dynamic Earth Force Transmission to Underground Structures

Discrete Element Modelling of a Reinforced Concrete Structure

Real Gas Thermodynamics. and the isentropic behavior of substances. P. Nederstigt

MATERIAL MECHANICS, SE2126 COMPUTER LAB 4 MICRO MECHANICS. E E v E E E E E v E E + + = m f f. f f

The International Association for the Properties of Water and Steam

Blast design parameters and their impact on rock fragmentation

Kinetic plasma description

Transcription:

DESIGN METHOD BASED ON THE CONCEPT OF EMERGENCE AND ITS APPLICATION Koichiro Sato¹ and Yoshiyuki Matsuoka 2 ¹ Graduate School of Science and Technology, Keio Uniersity, Yokohama, Japan Koichiro_Sato@a7.keio.jp, matsuoka@mech.keio.ac.jp 2 Faculty of Science and Technology, Keio Uniersity, Yokohama, Japan matsuoka@mech.keio.ac.jp ABSTRACT: In the early process of design, dierse ideas of design must be obtained from global solution search under the constraint of unclear design conditions. It is difficult to derie an adequate design solution to a design object with complicated conception because the current method of deriing dierse solutions depends solely on designer's experience and knowledge. Therefore, it is thought that, with a method of deriing dierse solutions globally, it leads to the further improement of design enironment. In this research, we propose the system based on emergent design that deries dierse solutions that can be applied to the early process of design, by using the concept of emergence that explains the origin and eolution of the organism. Moreoer, we aim to indicate the possibility of deriing dierse solutions using the emergent design system, by applying this system to artifact design and comparing this system with the conentional optimization method. Keywords: Emergence, dierse design solutions, early process of design 1

1. INTRODUCTION The design process can be diided roughly into two processes: the early processes of design corresponding to conceptual design and basic design, and the late processes of design corresponding to detailed design. In the early process of design, dierse ideas of design must be obtained from global solution search under the constraint of unclear design conditions. On the other hand, in the late process of design, a unique solution is obtained from local solution search under the constraint of clear design conditions. In conentional engineering design, a unique solution is obtained from local solution search because it uses the optimization methods. This method is not applicable to the early process of design obtaining dierse solutions while it can apply to the late process of design. In the early process of design, it is difficult to derie an adequate design solution to a design object with complicated conception because the current method of deriing current dierse solutions depends solely on designer's experience and knowledge. Therefore, it is thought that, with a method of deriing dierse solutions globally, it leads to the further improement of design enironment. In this research, we propose the emergent design system that deries dierse solutions that can be applied to the early process of design, by using the concept of emergence that explains the origin and eolution of the organism. Moreoer, we aim to indicate the possibility of deriing dierse solutions using the emergent design system, by applying this system to artifact design and comparing this system with the conentional optimization method. 2. EMERGENCE AND EMERGENT DESIGN METHOD 2. 1. ANALOGY BETWEEN EMERGENCE PROCESS AND DESIGN PROCESS In the natural world, arious organism species exist in the same enironment. In the fields of biology and ecology, scientists began to regard that arious organism species were created in the process of emergence. In the field of system engineering, emergence is defined as obtaining new functions, structures, or actions by-two-way dynamic processes (Kitamura S). The first process (bottom-up process) inoles the generation of global order using local interactions among elements that behae independently. The second process (top-down process) inoles the binding behaiors of elements using the global order (Fig.1). 2

Global Order New function, structure, and action Bottom-up Top-down Elements behaing independently Local interaction Figure 1: Concept of emergence. The process of design can be expressed as follows. First, in the early process of design, an ealuation is gien to the design proposal, which is obtained from designers intuition. The design proposal is dismissed when the ealuation does not meet the standard, and the design proposal is adopted when the ealuation fulfils the standard. Then the design proposal obtained in the early process is optimized in the late process of design, and the design solution is deried. In comparison with this design process and the emergent process, the following similarities can be found. First, the process of generating the design proposal with the ealuation meeting a certain standard is similar to the bottom-up process that generating the entire feature by autonomous components interaction in the emergence. Second, the process of optimizing detailed part of design proposal is similar to the top-down process that binds components by the entire feature in the emergence. These similarities indicate the possibility of applying the concept of emergence to the design. 2. 2. CONCEPT OF EMERGENT DESIGN METHOD The design method based on emergence is called the emergent design method. This method has two processes; the bottom-up process and the top-down process of emergence. In the bottom-up process, a design proposal is obtained when dierse ideas of design deried under unclear design conditions satisfy the standards set by the designer. In the top-down process, the design 3

ideas deried by the bottom-up process are optimized. As a result, it is possible to derie dierse design solutions maintaining the diersity of the design idea obtained in the bottom-up process. 3. EMERGENT DESIGN SYSTEM 3. 1. BASIC STRUCTURE OF EMERGENT DESIGN SYSTEM Figure 2 shows the emergent design system based on the emergent design method. The emergent design system consists of bottom-up process and top-down process. In the bottom-up process, dierse solution candidates that meet the low standard set by the designer are deried, and in the top-down process, to satisfy the constraint condition, those solution candidates are optimized. This system deries dierse solutions by going through these two processes. BOTTOM-UP PROCESS TOP-DOWN PROCESS Input condition Form generation by Cellular automata Ealuation Candidate for design solutions Homogenization method Calculate elastic constants Finite element analysis Update density Ealuation End Cellular automata in n p b i w i e i d d max e d in n ( 1 k) ( 0 1) = k + k = b i w n ie i i= 1 p = max ( d d ) e d :Input ector of cellular automata :Neighborhood information ector :Positional information ector k :Composite ratio :Indicates the existence or nonexistence of element(1 or 0) :Coefficient of ector direction :Unit ector of the direction from the neighboring element to the object element :Distance between the apex and the object element :Unit ector of the direction from 26 :Distance between the apex and the most distant cell from the apex p Design solutions Figure 2: Emergent design system. 4

3. 2. BOTTOM-UP PROCESS In bottom-up process, the forms are generated self-organizationally by using the cellular automata (CA). The diersity of an organism is noted, and rules referring to two properties for dierse organism morphogenesis, induction and apical dominance, are input ector to the CA(Fig.3). Firstly, we describe induction. An organism is formed by interaction between neighborhood cells. Neighborhood cells affect a cell and change it into a cell exhibiting different features. This property is called induction. The first input is defined as neighborhood information ector n, which is expressed by Eq. (1). n = 26 i= 1 biw i e i (1) Where i is the surrounding element number, b i indicates the existence or non-existence of element (1 or 0), w i is the coefficient of ector direction, and e i is the unit ector of the direction to the object element. Secondly, we describe apical dominance. In the deelopmental process, a certain tissue dominates; for example, the buds of plants and heads of animals. Such tissue is called an apex, and dominant action by the apex is called apical dominance. The second input is defined as positional information ector p, which is expressed by Eq. (2) = ( d max d) p e d (2) Induction Neighborhood Information Vector Apical Dominance Apex Dominant Tissue Positional Information Vector w 6 e i 1 2 3 4 5 6 7 8 9 d max -d V p d d max Apex n w 1 e i Form - Generation Space Interaction between a cell and neighborhood cells Dominant action by the apex (a) Induction (b) Apical dominance Figure 3: Modelling of input ector. 5

Where d max is the distance between the apex and the most distant cell from the apex, d is the distance between the apex and the object element, and e d is the unit ector of the direction to the object element. Moreoer, the composite ratio k is set, and input ector in is defined as expressed by Eq. (3). = k + ( 1 k) in n p (3) Dierse forms can be generated by changing the composite ratio k. Thus, dierse design ideas are generated self-organizationally in the bottom-up process. 3. 3. TOP-DOWN PROCESS In the top-down process, dierse solutions are deried by optimizing dierse ideas of design generated with the bottom-up process respectiely. In this research, we execute topology optimization based on the homogenization method. The homogenization method is explained as follows. The design domain is assumed to be composed of microstructures with periodicity, and the homogenized material constant is calculated from this relation between microstructures and macrostructures. In the topology optimization based on the homogenization method, the entire design domain is assumed to be composed of porous microstructures with regular periodicity, and the macro elastic constant homogenized from porous microstructures is calculated by the homogenization method. In addition, the state equation of the domain is soled by the finite element method using the elasticity tensor of each homogenized part, and repeatedly calculated until satisfying the objectie function and the constraint condition (Kikuchi N). As a result of such topology optimization, the density of the design domain increases as the hole of porous microstructures narrows in the load part of the stress, and the density of the design domain decreases as the hole of porous microstructures extends in the no-load part of the stress. Hence, the optimal material distribution in the design domain can be deried. 6

4. APPLICATION OF EMERGENT DESIGN SYSTEM TO ARTIFACT DESIGN In this research, when applying the emergent design system to a concrete artifact design, we selected an artificial hip, an artifact that is difficult to conceie for designers. 4. 1. INPUT CONDITIONS IN BOTTOM-UP PROCESS The input conditions for the proposed system execution are: the form generation space that shows generation range of element, the initial seed element position that is position of element generation beginning, a size of element, the position of apex, the element-generation number that is the number of updating forms, and the composite ratio k. The form generation space and an initial seed element position were decided as shown in Figure 4. The form generation space was decided as follows. First, the femoral model like Figure 4(b) was made referring to an aerage size of thighbone of Japanese men shown in Figure 4(a) ( Baba H). The position of the ball when the artificial hip was inserted from the center of femoral head as shown from this model in Figure 4(c) was decided, and the form generation space like Figure 4(d) was decided by basing the lesser trochanter on the position and the femoral diameter, etc. So as not to decrease the flexibility of the generation form, the initial seed elements were set as contact side with the ball shown in Figure 4(d) as the elements that existed in all the design solutions. The element size 59 86 74 Greater trochanter Head of femur Center point Contact side with ball 32 31 128 44 38 46 24 20 30 26 412 Lesser trochanter 102 (a) Aerage dimension of Japanese man femur (Unit:mm) 20 (Unit:mm) (b) Femoral model (c) Form generation space decision (d) Form generation space Figure 4: Decision of form generation space Artificial Hip. 7

was set to the length of 2mm on a side in consideration of the number of elements, and the position of apex was set to the center point of the ball supporting the load. The elementgeneration number was set to 150, and the forms were generated respectiely as the composite ratio k was changed from 0.1 to 1.0 at 0.1 interals, and 50 solution candidates were generated. 4. 2. INPUT CONDITIONS IN TOP-DOWN PROCESS In top-down process, the femoral model without femoral head was set on the assumption that the stem is inserted inside of human body. The greater trochanter was kept in this femoral model, and the cancellous bone and the cortical bone were set as shown in Figure 5. The solution candidates obtained in a bottom-up process were inserted in this model, and optimized. Figure 5 shows the mechanics condition. In this research, the load to the femoral head and the abductor muscle were assumed at the one foot geostationary upright positioning, and the compressie load of 1500N was set from the perpendicular line to the inside of 13 at the stem head top part, the tension load of 1000N was set from the perpendicular line to the inside of 20 at the greater trochanter. Moreoer, the constraint condition is a complete restraint on the thighbone bottom side ( Sakamoto J. and Oda J). In addition, the material property alue used in the top down process is indicated in Table 1. In this research, the titanium alloy (Ti-6Al-4V) with a remarkable affinity for human body was used as the material of the stem ( Sakamaki T). The objectie function was Greater trochanter 1000N 13 1500N 20 Cancellous bone Table 1: Young's modulus and poisson's ratio. Cortical bone Cortical bone Young's modulus 15GPa Poisson's ratio 0.3 Cancellous bone 1GPa 0.3 Titanium alloy 126GPa 0.27 Figure 5: Loding and fixing conditions. 8

assumed to be aerage compliance (Eq. (4)), and the restriction condition was assumed to be a olume (Eq. (5)). minimize Φ = subject to 1 ε D εdω H (4) ρdω V (5) 2 T Ω Ω Here, Ω is the design domain, ε is the strain, DH is stress-strain matrix, ρ is dimensionless amount of the alue from 0 to 1 that shows whether the element exists, and V is olume. As for expression (4), the strain energy consered in the whole when the external force acts is shown, and forms with high stiffness are deried by minimizing the objectie function (Kikuchi N). 4. 3. RESULT AND CONSIDERATION Figure 6 and 7 shows the form generation process and examples of design candidates in bottomup process. It was confirmed that dierse forms, forms with different length of stem, forms with different topology, skeleton forms, to hae arious features were generated from these results. Figure 8 shows design candidates in the bottom-up process, the solution in the top-down process, and idea sketches. From Figure 8, it was confirmed that the structure to satisfy the mechanics condition to hae dierse topology had been generated. From these results, dierse forms were deried under the same initial condition in the bottom-up process. It was confirmed that dierse Figure 6: Form generation process. Figure 7: Examples of design candidates by bottom-up process. 9

Solution candidate by bottom-up process Solution by top-down process Idea sketch Figure 8: Examples of design solutions. solutions were deried by optimizing each form in the top-down process, which was not possible with the established optimization methods that only derie a unique solution depending on initial form. 5. EFFECTIVENESS ANALYSIS OF PROPOSED SYSTEM The form generation space filled with elements is set as initial form, and figure 9 shows a form to optimize the initial form by the topology optimization that uses the homogenization method. The solution shown in Figure 9 deried without going through the bottom-up process by the topology optimization that uses the homogenization method is defined as a conentional solution. Moreoer, the solution deried through the bottom-up process and the top down process was Optimization Figure 9: Solution by conentional optimization method. 10

700 600 Single load transmission Solutions by this system Conentional solution Maximum equialent stress [MPa] 500 400 300 200 (b) Structure different from conentional solution 100 (a) Frame structure Forked load transmission 0 0.05 0.10 0.15 0.20 0.25 Mass [kg] Figure 10: Comparison between mass and maximum equialent stress. defined as the solution by this system. Figure 10 shows the result of comparing 20 solutions deried by this system with the conentional solution about the mass and maximum equialent stress. The stem of Figure 10(a) is lightened to about 25% compared with the conentional solution because the form is like a skeleton that has a number of holes. Howeer, there is a part where high stress is shown in the structure of skeleton form. Moreoer, the form in the upper part of the stem of Figure 10(b) transmits the load from the ball by a different route from the conentional solution, and the form under the stem has improed the stiffness by uniting forked form at the bottom. From these results, it was confirmed that dierse solutions were deried; for example, a solution that is lightened maintaining the strength equal to the conentional solution and a solution that is lightened more than half weight compared to the conentional solution though the maximum equialent stress indicates high alue. Thus, the solutions, different from the conentional solution, were deried with reducing the weight and maintaining the strength, and the possibility of deriing dierse solutions using this system was indicated. 11

6. CONCLUSION In this research, we proposed the emergent design system that deried dierse solutions and was applicable to the early process of design by using the concept of emergence, and applied to the artificial hip design. The possibility of deriing arious solutions from this system was indicated by comparing the solutions deried from this system with the solution deried from conentional optimization method. In addition, this work is supported in part by Grant in Aid for the 21 st Century Center of Excellence for System Design: Paradigm Shift from Intelligence to Life from Ministry of Education, Culture, Sport, and Technology in Japan. REFERENCES: Kitamura S.,Emergent Systems-Toward a New Paradigm for Artificial Systems. Toward System Theory of Function Emergence. Society of Instrument and Control Engineers, 1996, 35(7), 492-495. Kikuchi N. Optimum design theory by homogenization method. Bulletin of the Japan Society for Industrial and Appleied Mathematics, 1993,3(1), 2-26 Baba H. Anthropology 13 ols. And 3 additional ols.: Osteometry, 1991 (Yuzankaku Inc, Tokyo). Sakamoto J. and Oda J. A Study on Deelopment of the Total Hip Prosthesis Design Fitted for Japanese Patients (2nd Report, Shape Accuracy Estimation of the Femur CAD Model and Improement in the Canal Fill Ratio of Prosthesis). Japan Society Mechanical Engineering, 2004, 70(697), 1186-1192. Sakamaki T. The Recent Adance of the Artificial Hip Joint. Journal of the Keio Medical Society, 1997, 74(5), 293-302. 12