- +1-530-264-8006
- info@tutorsglobe.com

18,76,764

Questions

Asked

21,311

Experts

9,67,568

Questions

Answered

Start Excelling in your courses, Ask an Expert and get answers for your homework and assignments!!

Submit Assignment2015 © Tutors Globe. All rights reserved.

## Finite Elements and Differential Equations

Finite Elements::Triangulating a RegionA drawback of finite difference methods is that they require a very regular grid as well as thus a very regular region either rectangular or a regular part of a rectangle. Finite elements is a manner that works for any shape region because it isn’t built on a grid but on triangulation of the region that is cutting the region up into triangles as we did in a previous lecture.

The following figure depicts a triangularization of a region.

An annular region with a triangulation note down that the nodes and triangles are very evenly spaced.

This figure was formed by the script program mywasher.m. Notice that the nodes are consistently distributed. This is better for the finite element process where we will use it.

Open the program mywasher.m. This program describes a triangulation by defining the vertices in a matrix V in which every row contains the x and y coordinates of a vertex. Notice that we list the interior nodes first afterwards the boundary nodes.

Triangles are described in the matrix T. Every row of T has three integer numbers indicating the indices of the nodes that form a triangle. For example the first row is 43 42 25, so T1 is the triangle.

The finite elements for the triangulation of a one dimensional object.

with vertices v

_{43}, v_{42}and v_{25}. The matrix T in this case was produced by the Mat lab commands Delaunay. The command produced more triangles than want and the unwanted ones were deleted.We a three dimensional plot of the region as well as triangles is produced by having the last line be trimesh (T,x,y,z).

What is a finite element?The restricted element method is a mathematically complicated process. Nevertheless a finite element is actually a very simple object. To every node vj we associate a function Φ

_{j}that has the properties Φ_{j}(v_{j}) = 1 and Φ_{j}(v_{i}) = 0 for i /= j. Between nodes, Φ_{j }is a linear function. This function is the finite element. There are fancier kinds of finite elements that we will not discuss.If we consider one dimension subsequently a triangulation is just a subdivision into subintervals. We illustrate an uneven subdivision of an interval and the finite elements corresponding to every node in the above figure.

In two dimensions Φ

_{j}is calm of triangular piece of planes. Therefore Φ_{j}is a function whose graph is a pyramid with its peak over node v_{j}.What is a finite element solution?A finite element (approximate) explanation is a linear combination of the elements:

Therefore finding a finite element solution amounts to finding the best values for the constants {C

_{j}}^{n}_{j}=1. In the program mywasher.m the vector z contains the node values C_{j}. These values give the height at every node in the graph. For illustration if we set all equal to 0 except one equal to 1then the function is a finite element. Do this for one boundary node afterwards for one interior node.Notice that a total of linear functions is a linear function. Therefore the solution using linear elements is a piecewise linear function.

As well notice that if we denote the j-th vertex by v

_{j}, thenU(v

_{j}) = C_{j}Therefore we see that the constants C

_{j}are just the values at the nodes.Take the one-dimensional case since we know that the solution is linear on every subinterval knowing the values at the endpoints of the subintervals that are the nodes gives us complete knowledge of the solution. Could be a restricted element solution since it is piecewise linear.

A probable finite element solution for a one dimensional object. Values are assigned at every node and a linear interpolant is used in between.

In the two-dimensional case the solution is linear on every triangle and so again if we know the values {C

_{j}}^{n}_{j}=1 at the nodes afterwards we know everything.:Experiment with finite elementsBy changing the values in c in the program we are able to produce different three dimensional shapes based on the triangles. The point subsequently of a finite element solution is to find the values at the nodes that best approximate the true solution. This task is able to be subdivided into two parts- (1) assigning the values at the boundary nodes and (2) assigning the values at the interior nodes.

:Values at boundary nodesOnce a triangulation as well as a set of finite elements is set up the next step is to incorporate the boundary conditions of the problem. Presume that we have fixed boundary conditions that is of the form

u(x) = g(x), for x ∈ ∂D,

Where D is the object (domain) as well as ∂D is its boundary. After that the boundary condition directly determines the values on the boundary nodes.

In particular presume that vlis a boundary node. Since Φl(vl) = 1 and all the other elements are zero at node vl, then to create

We should choose:

Cl= g(vl).

Therefore the constants Cjfor the boundary nodes are set at exactly the value of the boundary condition at the nodes.

Therefore if there are m interior nodes then,

C

_{j}= g(vj), for all m + 1 ≤ j ≤ n.In the program mywasher the initial 32 vertices correspond to interior nodes and the last 32 correspond to boundary nodes. By setting the last 32 values of z we get the boundary conditions. We could do this by adding up the following commands to the program

z(33:64) = .5;

Or more elaborately we might utilize functions:

z(33:48) = x(33:48).^2 - .5*y(33:48).^2

z(49:64) = .2*cos(y(49:64))

Latest technology based Matlab Programming Online Tutoring AssistanceTutors, at the

www.tutorsglobe.com, take pledge to provide full satisfaction and assurance inMatlab Programming helpviaonline tutoring. Students are getting 100% satisfaction byonline tutorsacross the globe. Here you can get homework help for Matlab Programming, project ideas and tutorials. We provide email basedMatlab Programming help. You can join us to ask queries 24x7 with live, experienced and qualified online tutors specialized in Matlab Programming. ThroughOnline Tutoring, you would be able to complete your homework or assignments at your home. Tutors at theTutorsGlobeare committed to provide the best qualityonline tutoringassistance forMatlab Programming Homework helpandassignment helpservices. They use their experience, as they have solved thousands of the Matlab Programming assignments, which may help you to solve your complex issues of Matlab Programming.TutorsGlobeassure for the best quality compliance to your homework. Compromise with quality is not in our dictionary. If we feel that we are not able to provide thehomework helpas per the deadline or given instruction by the student, we refund the money of the student without any delay.