evolutionary structural optimisation
DESCRIPTION
Evolutionary Structural Optimisation. KKT Conditions for Topology Optimisation. KKT Conditions (cont’d). KKT Conditions (cont’d). Strain energy density should be constant throughout the design domain This condition is true if strain energy density is evenly distributed in a design. - PowerPoint PPT PresentationTRANSCRIPT
![Page 1: Evolutionary Structural Optimisation](https://reader035.vdocument.in/reader035/viewer/2022062811/568160e2550346895dd0132c/html5/thumbnails/1.jpg)
1
Evolutionary Structural Optimisation
![Page 2: Evolutionary Structural Optimisation](https://reader035.vdocument.in/reader035/viewer/2022062811/568160e2550346895dd0132c/html5/thumbnails/2.jpg)
2
kkuuvukukuL
ukukvukuukukuuL
ukukku
ukuukukuuC
kuVCLfkuVVCL
fkuVV
ukuC
T
ee
ee
T
e
eee
e
T
e
T
e
T
e
eee
e
T
e
T
e
T
e
eeee
eee
te
212
21
21
201
0
2
where,
0
function, Lagrangian00toSubject
Minimise
KuuT
KKT Conditions for Topology Optimisation
![Page 3: Evolutionary Structural Optimisation](https://reader035.vdocument.in/reader035/viewer/2022062811/568160e2550346895dd0132c/html5/thumbnails/3.jpg)
3
KKT Conditions (cont’d)
€
∂L∂ρ e
= uT ∂k∂ρ e
u+ λ 2∂k∂ρ e
u+ λ1ve + ∂u∂ρ e
2uT k + λ 2k( )
Since λ 2 is arbitary, we select λ 2 to eliminate the ∂u∂ρ e
term, i.e. λ 2 = −2uT
∴ ∂L∂ρ e
= uT ∂k∂ρ e
u− 2uT ∂k∂ρ e
u+ λ1ve = −uT ∂k∂ρ e
u+ λ1ve = 0
Let the strain energy of a solid element (ρ e =1) be se ,se= uek0ueand,
uT ∂k∂ρ e
u = ∂ρ∂ρ e
uTk0u =se
Therefore,seλ1ve
=1
![Page 4: Evolutionary Structural Optimisation](https://reader035.vdocument.in/reader035/viewer/2022062811/568160e2550346895dd0132c/html5/thumbnails/4.jpg)
4
KKT Conditions (cont’d)
Strain energy density should be constant throughout the design domain
This condition is true if strain energy density is evenly distributed in a design.
Similar to fully-stressed design.
Need to compute strain energy density Finite Element Analysis
€
seλ1ve
=1
![Page 5: Evolutionary Structural Optimisation](https://reader035.vdocument.in/reader035/viewer/2022062811/568160e2550346895dd0132c/html5/thumbnails/5.jpg)
5
Evolutionary Structural Optimisation (ESO)
Fully-stressed design – von Mises stress as design sensitivity.
Total strain energy = hydrostatic + deviatoric (deviatoric component usually dominant in most continuum)
Von Mises stress represents the deviatoric component of strain energy.
Removes low stress material and adds material around high stress regions descent method
Design variables: finite elements (binary discrete) High computational cost. Other design requirements can been incorporated by
replacing von Mises stress with other design sensitivities – 0th order method.
![Page 6: Evolutionary Structural Optimisation](https://reader035.vdocument.in/reader035/viewer/2022062811/568160e2550346895dd0132c/html5/thumbnails/6.jpg)
6
ESO Algorithm
1. Define the maximum design domain, loads and boundary conditions.
2. Define evolutionary rate, ER, e.g. ER = 0.01.3. Discretise the design domain by generating finite element
mesh. 4. Finite element analysis.5. Remove low stress elements,
6. Continue removing material until a fully stressed design is achieved
7. Examine the evolutionary history and select an optimum topology that satisfy all the design criteria.€
σ e ≤ ER × 1+σ min( )
![Page 7: Evolutionary Structural Optimisation](https://reader035.vdocument.in/reader035/viewer/2022062811/568160e2550346895dd0132c/html5/thumbnails/7.jpg)
7
Cherry
Initial design domain
Fixed
Gravitational Load
ESO solution
![Page 8: Evolutionary Structural Optimisation](https://reader035.vdocument.in/reader035/viewer/2022062811/568160e2550346895dd0132c/html5/thumbnails/8.jpg)
8
Michell Structure Solution
![Page 9: Evolutionary Structural Optimisation](https://reader035.vdocument.in/reader035/viewer/2022062811/568160e2550346895dd0132c/html5/thumbnails/9.jpg)
9
ESO: Michell Structure
![Page 10: Evolutionary Structural Optimisation](https://reader035.vdocument.in/reader035/viewer/2022062811/568160e2550346895dd0132c/html5/thumbnails/10.jpg)
10
ESO: Long Cantilevered Beam
![Page 11: Evolutionary Structural Optimisation](https://reader035.vdocument.in/reader035/viewer/2022062811/568160e2550346895dd0132c/html5/thumbnails/11.jpg)
2.5D Optimisation Reducing thickness relative to sensitivity values rather than
removing/adding the whole thickness
11
Mesh Size 14436, less the 268 elements removed from mouth of spanner
Load Case 1: 2N/mm
Load Case 2: 2N/mm
Roller support
Non-Design Domain
![Page 12: Evolutionary Structural Optimisation](https://reader035.vdocument.in/reader035/viewer/2022062811/568160e2550346895dd0132c/html5/thumbnails/12.jpg)
Spanner
12
![Page 13: Evolutionary Structural Optimisation](https://reader035.vdocument.in/reader035/viewer/2022062811/568160e2550346895dd0132c/html5/thumbnails/13.jpg)
Thermoelastic problems Both temperature and mechanical loadings FE Heat Analysis to determine the temperature distribution Thermoelastic FEA to determine stress distribution due to
temperature Then ESO using these stress values
477
720
24
Design Domain
P UniformTemperature
![Page 14: Evolutionary Structural Optimisation](https://reader035.vdocument.in/reader035/viewer/2022062811/568160e2550346895dd0132c/html5/thumbnails/14.jpg)
Plate with clamped sides and central load
14
T = 7CT = 5C
T = 0C T = 3C
![Page 15: Evolutionary Structural Optimisation](https://reader035.vdocument.in/reader035/viewer/2022062811/568160e2550346895dd0132c/html5/thumbnails/15.jpg)
Group ESO Group a set of finite elements Modification is applied to the entire set Applicable to configuration optimisation
15
![Page 16: Evolutionary Structural Optimisation](https://reader035.vdocument.in/reader035/viewer/2022062811/568160e2550346895dd0132c/html5/thumbnails/16.jpg)
Example: Aircraft Spoiler
16
![Page 17: Evolutionary Structural Optimisation](https://reader035.vdocument.in/reader035/viewer/2022062811/568160e2550346895dd0132c/html5/thumbnails/17.jpg)
Example: Optimum Spoiler Configuration
17
![Page 18: Evolutionary Structural Optimisation](https://reader035.vdocument.in/reader035/viewer/2022062811/568160e2550346895dd0132c/html5/thumbnails/18.jpg)
Multiple Criteria Using weighted average of sensitivities as removal/addition
criteria
18
AR 1.5Mesh 45 x 30
P
P
Maximise first mode frequency & Minimise mean compliance
![Page 19: Evolutionary Structural Optimisation](https://reader035.vdocument.in/reader035/viewer/2022062811/568160e2550346895dd0132c/html5/thumbnails/19.jpg)
Optimum Solutions (70% volume)
19
wstiff:wfreq = 1.0:0.0
wstiff:wfreq = 0.7:0.3 wstiff:wfreq = 0.5:0.5
wstiff:wfreq = 0.0:1.0wstiff:wfreq = 0.3:0.7
![Page 20: Evolutionary Structural Optimisation](https://reader035.vdocument.in/reader035/viewer/2022062811/568160e2550346895dd0132c/html5/thumbnails/20.jpg)
20
Chequerboard Formation
Numerical instability due to discretisation. Closely linked to mesh dependency.
Piecewise linear displacement field vs. piecewise constant design update
![Page 21: Evolutionary Structural Optimisation](https://reader035.vdocument.in/reader035/viewer/2022062811/568160e2550346895dd0132c/html5/thumbnails/21.jpg)
Smooth boundary: Level-set function Topology optimisation based on moving smooth boundary Smooth boundary is represented by level-set function Level-set function is good at merging boundaries and
guarantees realistic structures Artificially high sensitivities at nodes are reduced, and
piecewise linear update numerically more stable Manipulate implicitly through
21
0,, yxyx
yxyxyxyxyxyx
,0,,0,,0,
http://en.wikipedia.org/wiki/File:Level_set_method.jpg
![Page 22: Evolutionary Structural Optimisation](https://reader035.vdocument.in/reader035/viewer/2022062811/568160e2550346895dd0132c/html5/thumbnails/22.jpg)
Topology Optimisation using Level-Set Function
Design update is achieved by moving the boundary points based on their sensitivities
Normal velocity of the boundary points are proportional to the sensitivities (ESO concept)• Move inwards to remove material if sensitivities are low• Move outwards to add material if sensitivities are high
Move limit is usually imposed (within an element size) to ensure stability of algorithm
Holes are usually inserted where sensitivities are low (often by using topological derivatives, proportional to strain energy)
Iteration continued until near constant strain energy/stress is reached.
22
![Page 23: Evolutionary Structural Optimisation](https://reader035.vdocument.in/reader035/viewer/2022062811/568160e2550346895dd0132c/html5/thumbnails/23.jpg)
Numerical Examples
23