isoparametric elements and solution techniques. advanced design for mechanical system - lec...
Post on 15-Jan-2016
216 views
TRANSCRIPT
![Page 1: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/1.jpg)
Isoparametric elements and solution techniques
![Page 2: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/2.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 2
![Page 3: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/3.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 3
= ½ d1-2Tk1-2d1-2 +
+ ½ d2-4Tk2-4d2-4 +….=
= ½ DTKD
![Page 4: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/4.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 4
R=KD
• gauss elimination
• computation time:
(n order of K, b bandwith)
![Page 5: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/5.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 5
recall: gauss elimination
![Page 6: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/6.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 6
rotations
![Page 7: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/7.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 7
![Page 8: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/8.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 8
isoparametric elements
isoparametric: same shape functions for both displacements and coordinates
![Page 9: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/9.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 9
computation of B
x = du / dX
• but u=u(, ), v=v (, )
![Page 10: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/10.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 10
![Page 11: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/11.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 11
• J11* and J12* are coefficients of the first row of J-1
and
![Page 12: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/12.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 12
![Page 13: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/13.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 13
gauss quadrature
![Page 14: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/14.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 14
![Page 15: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/15.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 15
![Page 16: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/16.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 16
![Page 17: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/17.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 17
![Page 18: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/18.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 19
![Page 19: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/19.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 20
no strain at the Gauss points
so no associated strain energy
![Page 20: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/20.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 21
The FE would have no resistance to loads that would activate these modes
Global K singularUsually such modes superposed to ‘right’ modes
![Page 21: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/21.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 22
![Page 22: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/22.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 23
calculated stress =EBd are accurate at Gauss points
![Page 23: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/23.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 24
• the locations of greatest accuracy are the same Gauss points that were used for integration of the stiffness matrix
![Page 24: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/24.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 25
![Page 25: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/25.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 26
Rayleigh-Ritz method
Guess a displacement set that is compatible and satisfies the
boundary conditions
![Page 26: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/26.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 27
• define the strain energy as function of displacement set
• define the work done by external loads• write the total energy as function of the
displacement set• minimize the total energy as function of
the displacement and find• simulataneous equations that are solved
to find displacements
![Page 27: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/27.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 28
= (d)
d / d d1 = 0d / d d2 = 0d / d d3 = 0d / d d4 = 0……d / d dn = 0
![Page 28: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/28.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 29
![Page 29: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/29.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 30
![Page 30: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/30.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 31
patch tests
• only for those who develops FE
![Page 31: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/31.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 32
substructures
![Page 32: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/32.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 33
• divide the FEmodel in more parts
• create a FE model of each substructure
• Assemble the reduced equations KD=R
• Solve equations
![Page 33: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/33.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 34
Simmetry
![Page 34: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/34.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 35
![Page 35: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/35.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 36
![Page 36: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/36.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 37
Constraints
CD – Q =0
C is a mxn matrix
m is the number of constraints
n is the number of d.o.f.
How to impose constraints on KD=R
![Page 37: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/37.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 38
way 1 – Lagrange multipliers
=[1 2 …. m]T
T [CD-Q]=0
= 1/2DTKD – DTR + T [CD-Q]
![Page 38: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/38.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 39
• remember
dAD / dD = AT
dDTA/ dD = A
![Page 39: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/39.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 40
example
![Page 40: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/40.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 41
![Page 41: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/41.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 42
way 2- penalty method
![Page 42: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/42.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 43
½tT t = ½ [(CD-Q)T(CD-Q)]== ½ [(CD-Q)T(CD- Q)]== ½ [(CD-Q)TCD- (CD-Q)T Q)]== ½ [(DTCTCD-QTCD-DTCTQ+QTQ)]= ½[·];d(½[·])/dD==½[2(CTC)-(QTC)T- CTQ]==½[2(CTC)-(C)TQ- CTQ]==½[2(CTC)-CT Q- CTQ]=
![Page 43: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/43.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 44
=½[2(CTC)-CT Q- CTQ]== CTC-CT Q
(= T)
![Page 44: Isoparametric elements and solution techniques. Advanced Design for Mechanical System - Lec 2008/10/092](https://reader035.vdocument.in/reader035/viewer/2022062423/56649d625503460f94a4437f/html5/thumbnails/44.jpg)
Advanced Design for Mechanical System - Lec 2008/10/09 45