fem01.fm

Similar documents
242 외국어로서의독일어제 41 집.,.,. III.3 Y,. (2016). (2016),,,.,..., II.,,,

통편집.hwp

¿Ü±¹¹ýÁ¦³»Áö09054)

<C8AFB0E6B9FDBFACB1B85F3236B1C75F33C8A32E687770>

<C7D5C0C7BEC8C0C75FC1A6BDC3BFCD5FC3A4B9ABC0C75FBDC2C0CE5FB9E9B0E6C0CF5FC3D6C1BEBABB2E687770>

730 Basel - Waldsut - Scaausen - Singen Hocreinban 730 Von Basel Bad B bis Scwörstadt Verbundtari Regio Verkersverbund Lörrac (RVL) Von Wer-Brennet bi

고등독일어브로셔내지_양도원

FI_2017_007


보험판매와 고객보호의 원칙

Mappe_holzimgarten.pdf

Kompetent.Kernig.Kreativ.tiv.

Interstage5 SOAP서비스 설정 가이드

산선생의 집입니다. 환영해요

SW¹é¼Ł-³¯°³Æ÷ÇÔÇ¥Áö2013

UML

untitled

23

Staatsministeriums für Arbeit und Soziales, Familie und Integration, 2015). 이번 글은 우리 주변에서 흔히 접할 수 있는 직업 중의 하나인 소매업체 판매원의 직업훈련 프 로그램을 소개한다. Ⅱ. 소매업체 판매원

example code are examined in this stage The low pressure pressurizer reactor trip module of the Plant Protection System was programmed as subject for

SLTION UI SP RILL INSRTS SIZ ITM MOL SRIPTION P MIN MX SRIS ~8 SP RILL INSRTS - M4 INW OHRINSTZ - M4 Ø7.86 (#) Ø4.3 (#8) 5 SRIS Y,Z,,~4 SP RILL INSRTS

ETL_project_best_practice1.ppt

Orcad Capture 9.x

415.2 Duisburg - Essen - Dortmund - Hamm (West) Regionalverkehr Duisburg - Kamen Verbundtari VRR Kamen - Hamm Verkehrsgemeinschat Ruhr-Lippe (VR

°ø¾÷-01V36pš

455 Köln - Wuppertal - Hagen Hamm Dortmund - - Münster 455 Köln - Solingen Hb Verbundtari VRS Solingen Hb - Unna Verbundtari VRR Unna - Hamm/Münster V


<31B1E8C1A4B7E6B9DAC1F6BCF6B1E8B9CCBCF72E687770>

4 CD Construct Special Model VI 2 nd Order Model VI 2 Note: Hands-on 1, 2 RC 1 RLC mass-spring-damper 2 2 ζ ω n (rad/sec) 2 ( ζ < 1), 1 (ζ = 1), ( ) 1

목 차 1. 국외출장목적 1 2. 국외출장개요 1 3. 예산집행계획 2 4. 출장국가 ᆞ 방문기관세부내용 3 5. 출장결과 6

untitled

untitled

436 8., {(x, y) R 2 : y = x, < x 1} (, 1] φ(t) = (t, t), (, 2] ψ(t) = (t/2, t/2), [1, ) σ(t) = (1/t, 1/t).. ψ φ, σ φ. (φ, I) φ(i) φ : I φ(i). 8.2 I =

2 지방계약연구제 2 호 Ⅱ. 기본원칙규정 : 행정절차법제 54 조 제 62 조 2) 제 54 조공법상계약의허용 상이한법규가존재하지않는한, 공법적영역에서의법률관계는계약을통해체결, 변 경, 해지될수있다 ( 공법상의계약 ). 특히, 관청은행정행위를발하는대신, 행정행위의 2

<C7D1B1B9B0E6C1A6BFACB1B8C7D0C8B828C0CCC1BEBFF85FC0CCBBF3B5B75FBDC5B1E2B9E9292E687770>

+변협사보 4월호

Interview einer Erzieherin / eines Erziehers Das Interview wird mit Personen aus verschiedenen Ländern geführt und jeweils eines aus Deutschland und e

No Slide Title

공학박사학위 논문 운영 중 터널확대 굴착시 지반거동 특성분석 및 프로텍터 설계 Ground Behavior Analysis and Protector Design during the Enlargement of a Tunnel in Operation 2011년 2월 인하대

가장죽 - 독일어 문법 훑기

2011³»ÁöÆíÁýÃÖÁ¾

Alle Menschen werden Brüder, ( 합창 알토, 테너, 바리톤 ) 시류가가차없이갈라놓은것을네마법이다시묶는구나. 네온화한날개가머무는곳에서모든사람이형제가되리라. (A, T, B Solo) Wem der große Wurf gelungen, Eines F

<C0CCBCBCBFB52DC1A4B4EBBFF82DBCAEBBE7B3EDB9AE2D D382E687770>

¿Ü±¹¹ýÁ¦³»Áö0904È£4)

11234(3)

목차 1. 제품 소개 특징 개요 Function table 기능 소개 Copy Compare Copy & Compare Erase

歯J PDF

도식악보활용\(민경훈\)

Oracle Apps Day_SEM

Kursbuc der Deutscen Ban Reine - Münster - Hamm 410 Reine - Hamm Verkersgemeinscat Münsterland (VGM) Verbundarausweise gel

s SINUMERIK 840C Service and User Manual DATA SAVING & LOADING & & /

슬라이드 1

<31325FB1E8B0E6BCBA2E687770>

LXR 설치 및 사용법.doc

untitled

초보자를 위한 ADO 21일 완성

PowerPoint 프레젠테이션

08-이용갑

LEICA C-LUX

Your title goes here

e hwp

13 Who am I? R&D, Product Development Manager / Smart Worker Visualization SW SW KAIST Software Engineering Computer Engineering 3

PWR PWR HDD HDD USB USB Quick Network Setup Guide xdsl/cable Modem PC DVR 1~3 1.. DVR DVR IP xdsl Cable xdsl Cable PC PC DDNS (

<3033C6AFC1FD5FC0CCB3B2C0CE2E687770>

OR MS와 응용-03장

歯Enet_목차_.PDF

Boström, Familienerbrecht und Testierfreiheit in Schweden und anderen skandinavischen Ländern, in: Familienerbrecht und Testierfreiheit im europäische

09´ëÇк°¼Ò°³-¼�¿ïÄ·2

정보기술응용학회 발표

강의10

PowerPoint 프레젠테이션

독일기악미학\(이경희\)

Microsoft PowerPoint - AC3.pptx

VOL /2 Technical SmartPlant Materials - Document Management SmartPlant Materials에서 기본적인 Document를 관리하고자 할 때 필요한 세팅, 파일 업로드 방법 그리고 Path Type인 Ph

歯엑셀모델링

MPLAB C18 C

< C7D0B3E2B5B520B4EBBCF6B4C920C7D8BCB3C1F620C1A632BFDCB1B9BEEE26C7D1B9AEBFB5BFAA5FB5B6C0CFBEEE20492E687770>

M&S. EDISON. e-science EDISON_CFD[2] EDISON. (geometry kernel software),,, EDISON_CFD.,.. 4. e-science EDISON M&S.., EDISON , b f. Find to min

歯 PDF


USB USB DV25 DV25 REC SRN-475S REC SRN-475S LAN POWER LAN POWER Quick Network Setup Guide xdsl/cable Modem PC DVR 1~3 1.. DVR DVR IP xdsl Cable xdsl C

untitled

C# Programming Guide - Types

10김묘선

<C0FCC8C4BCD2BCB3B0FAC7E3B9ABC1D6C0C7C0FBB9CCC0C7BDC42DBABBB9AE2E687770>

제17권 1호(통권 제38호).hwp

DBPIA-NURIMEDIA

75)..... ( ).,,,,, < 국문초록 >.....,, * ( ), ( ), ( )

solution map_....

14.091~100(328-하천방재).fm

VZ94-한글매뉴얼

(2) : :, α. α (3)., (3). α α (4) (4). (3). (1) (2) Antoine. (5) (6) 80, α =181.08kPa, =47.38kPa.. Figure 1.

27(5A)-15(5868).fm

12(4) 10.fm

제목을 입력하세요.

Microsoft PowerPoint - User Manual pptx

PowerPoint 프레젠테이션

피고인의 불리한 증인에 대한 반대신문이 허용되지 않은 것은 소송절차 전 체과정을 통하여 적정하지 않았다고 판단될 때에만 협약을 위반한 것이다. 법정에 출석하지 않은 증인의 수사단계 진술의 증거가치는 반대신문이 행해 지지 않은 것에 대하여 사법기관에 얼마나 책임이 있는지

2005 7

ODS-FM1

Transcription:

Überblick Methode der finiten Elemente Beispiele Software im Internet Schritte für die Programmierung

Methode der finiten Elemente Bisher Finite Differenzen: Diskretisierung des Ortes und der Zeit Iterationsverfahren Untersuchung der Stabilität und der Konvergenzgeschwindigkeit FD nicht direkt anwendbar unregelmässig geformtes Grundgebiet Variationsaufgaben Methode der finiten Elemente allgemeiner analytischer Ansatz für die Näherungslösung der partiellen Differentialgleichung Methode der gewichteten Residuen: Differentialgleichungen Rayleigh-Ritz Verfahren: Variationsprobleme

Methode der gewichteten Residuen Ziel: Finde zu einer partielle Differentialgleichung eine System von gewöhnlichen Differentialgleichungen, oder ein algebraisches Gleichungssystem (zeitunabhängige Probleme) u( x, t) zb: Temperatur, soll vom Ort x und der Zeit t abhängen Feldvariable u( x, t) Differentialoperator L und dem Quellterm f(x,t) ist Lösung einer Differentialgleichung mit dem Lux [ (, t) ] = f( x, t) Beispiel: Poissongleichung u = ρ

Residuum Näherungslösung u( x, t) Residuum setze u in die Differentialgleichung ein Abweichung vom exakten Wert R( x, t) = Lux [ (, t) ] f( x, t) Ansatz für die Nährungslösung Residuum = 0 für die exakte Lösung Bestimme die Entwicklungskoeffizienten so, dass das Residuum klein ist.

Ansatz für die Näherungslösung Finite Differenzen ersetze x und t durch diskrete Variable Finite Elemente x und t bleiben kontinuierliche Variable Entwicklung der Näherungslösung u( x, t) a i ()ϕ t i ( x) i ϕ i ( x) a i (t)... Basisfunktionen... Entwicklungskoeffizenten

Bestimmung der Entwicklungskoeffizenten Kollokationsverfahren Die Differentialgleichung ist an N Stützstellen erfüllt R(x i ) = 0, (i=1,...,n) Methode der kleinsten Fehlerquadrate bestimme a i so, dass das Integral minimal wird ------ R( x) a [ ] 2 d x = 0 i Galerkin Verfahren multipliziere das Residuum mit N Gewichtsfunktionen w i R ( x )w d x = 0 i

Galerkin Verfahren Galerkin-Verfahren Für beliebige Gewichtsfunktionen folgt nicht, dass R klein ist spezielle Wahl Ansatzfunktionen = Gewichtsfunktionen w i ( x) u = ------- = ϕ a i ( x) i Methode der gewichteten Residuen ist unabhängig von finiten Elementen beliebige Basisfunktionen möglich analytische und numerische Methoden Beispiel: Fourier-Galerkin Verfahren der Konvektion in der Atmosphere führte auf die Lorenz Gleichungen

Finite Elemente Spezielle Wahl der Basisfunktionen Basisfunktion ϕ i 0 nur in kleinen Bereichen Entwicklungskoeffizenten a i entsprechen den Feldwerten an den Punkten x i Basisfunktionen interpolieren die Lösung zwischen den Punkten hat-function u u i-1 u i+1 u i 1 ϕ i-1 ϕ i ϕ i+1 x i-1 x i x i+1

1D stückweise lineare Funktionen u u i-1 u i+1 u i 1 ϕ i-1 ϕ i ϕ i+1 x i-1 x i x i+1 ϕ i = x x -------------------- i 1, x x i x i 1 x < x i i 1 x i + 1 x --------------------, x i x < x i + 1 x i + 1 x i 0, x< x i 1 or x x i + 1

Eigenschaften des FE-Ansatzes Entwicklungskoeffizenten entsprechen den Feldwerten es gilt: ϕ i ( x i ) = 1 ϕ j ( x i ) 0 und = für j i ux ( i ) = a i ϕ i ( x i ) = i a i Innerhalb eines Elementes linear interpoliert x i x < x i + 1 wird das Feld x i + 1 x x x ux ( ) u i -------------------- i = + u i + 1 -------------------- x i + 1 x i x i + 1 x i ϕ i ϕ i+1

Beispiel: 1D Transportgleichung u ----- = v u ----- t x für jede Basisfunktion ϕ i erhält ergibt sich eine Gleichung multipliziere Residuum mit Basisfunktion integriere über das Grundgebiet u ----- + v u ----- ϕj dx = 0 t x Einsetzen der Näherungslösung uxt, ( ) = a i ()ϕ t i ( x) i da iϕi dϕ i ϕ dt j + va i ϕ j dx = 0 dx i

Matrix-Schreibweise da iϕi dϕ i ϕ, i = 1,...,N dt j + va i ϕ j dx = 0 dx i A da + B a = 0 dt mit A ij und = ϕ i ϕ j dx dϕ i B ij = v ϕ j dx dx System von N gewöhnlichen Differentialgleichungen

2 ------- u t 2 v 2 u = -------- x 2 Beispiel: Wellengleichung im Intervall (0,1) Galerkin Verfahren: 1 2 da iϕi dt 2 ϕ j v 2 d ϕ i aiϕ j dx 2 dx = 0 i 0 partielle Integration 2 1 2 da i dt 2 ϕ i ϕ j dx a i v 2 dϕ i dϕ j dx v 2 dϕ i + ϕ dx dx j = 0 dx i 0 A ij i 1 0 B ij 1 0

Matrix-Schreibweise A d 2 d a t 2 + B a = 0 periodische Zeitabhängigkeit: Ansatz Eigenwertproblem ω 2 A a 0 = B a 0 a() t = a 0 e iωt

Finite element resources http://www.engr.usask.ca/%7emacphed/finite/fe_resources/fe_resources.html Internet Finite Element Resources Conferences and Courses Conferences Courses Newsgroups Particular Newsgroup Items Listservers and Bulletin Board Systems FE Homepages and FTP Sites Societies FE Links CFD Links Mesh Generation Links Public Domain FE Programs FREEWARE with source code Programs for the PC SHAREWARE or no source code Mesh Generators FE Visualisation Public Domain Mathematical Software and Miscellaneous Items. Books FE Code from Books

Gittererzeugung, Visualisierung: GID http://gid.cimne.upc.es/intro/index.html

User-Interface

Beispiel: Strömungsprobleme Aerodynamic analysis of the new large telescope facility in La Palma island, Islas Canarias, Spain.

Zusammenfassung Finite Element Links Software, Bücher, Vergleiche,... http://www.engr.usask.ca/%7emacphed/finite/ fe_resources/fe_resources.html GETDP / GMSH is a scientific computation software for the numerical solution of integro-differential equations, using finite element and integral type methods. http://www.geuz.org/getdp/ OFELI (Object Finite Element LIbrary) is a library of finite element C++ classes for multipurpose developement of finite element software. I http://sourceforge.net/projects/ofeli/ GID Modellierung, Gittererzeugung, Visualisierung http://gid.cimne.upc.es/intro/index.html

Beispiel: Getdp

Gebiete GROUP Gebiete und Regionen Group { /* The numbers correspond to physical regions defined in 'thermal.geo' (the input to the GMSH meshing tool) */ Hot = Region[10] ; Cold = Region[20] ; Domain = Region[100];

Funktionen FUNCTION Funktionen auf den Gebieten (groups) Function { alpha = 0.1 ; // Waermeleitfahigkeit

CONSTRAINT Randbedingungen, Anfangsbedingungen Randbedingungen Constraint { { Name Boundary; Case { { Region Hot; Value 200. ; { Region Cold ; Value 0. ;

FUNCTIONSPACE Funktionenraum FunctionSpace { { Name H1 ; Type Form0 ; BasisFunction { { Name wi ; NameOfCoef fi ; Function BF_Node ; Support Domain ; Entity NodesOf[ All ] ; Constraint { { NameOfCoef fi; EntityType NodesOf ; NameOfConstraint Boundary;

JACOBIAN geometrische Transformationen Geometrische Transformationen Randbedingungen im Unendlichen Jacobian { { Name JVol ; Case { { Region All ; Jacobian Vol ;

Integrale INTEGRATION Integrationsmethoden für Integrale über finite Elemente Integration { { Name I1 ; Case { { Type Gauss ; Case { { GeoElement Triangle ; NumberOfPoints 12 ;

FORMULATION Aufstellung der Gleichungen Eingabe der Gleichung in symbolischer Form Formulation { { Name HeatEquation; Type FemEquation ; Quantity { { Name temperature; Type Local ; NameOfSpace H1; Equation { Galerkin { [ alpha*dof{d temperature, {d temperature ] ; In Domain ; Jacobian JVol ; Integration I1 ;

RESOLUTION Beschreibung der Lösungsschritte Zeitabhängigkeit Lösen der Gleichungen Resolution { { Name HeatEquation; System { { Name HeatSys ; NameOfFormulation HeatEquation ; Operation { Generate[HeatSys] ; Solve[HeatSys] ; SaveSolution[HeatSys] ;

Beispiel: Zeitabhängigkeit FORMULATION Equation { Galerkin { Dt [ Dof{temperature, {temperature ] ; In Domain ; Jacobian JVol ; Integration I1 ; Galerkin { [ alpha*dof{d temperature, {d temperature ] ; In Domain ; Jacobian JVol ; Integration I1 ; theta Methode y t + 1 y t RESOLUTION dy dt = fyt (, ) -------------------- = f( ( 1 θ)y t t θy t + 1, t) Operation { InitSolution[HeatSys]; SaveSolution[HeatSys] ; TimeLoopTheta[tmin,tmax,dt,theta] { Generate[HeatSys] ; Solve[HeatSys] ; SaveSolution[HeatSys] ;

POSTPROCESSING Definition der Ausgabedaten (Integrale, Ableitungen der FE-Lösung) Ausgabedaten PostProcessing { { Name TemperaturePost ; NameOfFormulation HeatEquation ; Quantity { { Name Temp ; Value { Local { [ {Grad temperature ] ; In Domain;

Exportieren der Ausgabedaten POSTOPERATION wie werden die Postprocessing daten in files geschrieben? PostOperation { { Name TemperaturePlot; NameOfPostProcessing TemperaturePost; Operation { Print[ Temp, OnElementsOf Domain, File "thermaltime.pos"] ;