Norwegian University of Science and Technology NTNU Norwegian University of Science and Technology NTNU

Faculty of Engineering Science and Technology Faculty of Engineering Science and Technology
Department of Structural Engineering Page 1 of 5 Department of Structural Engineering Page 2 of 5
Responsible contact during exam: ENGLISH Problem 1 (weight ca. 50 %)
Professor Odd Sture Hopperstad, 73594703/73594700
y, v v3 K
EXAM IN COURSE TKT4191 4 7 3 u3 K 1
4 7 3
[
6
FINITE ELEMENT METHODS 1 8 6 8
5
1 2 1 5 2
K 1
Thursday 16 December 2004 [ 1 [ 1
09.00 – 13.00 x, u
Figure 1a Figure 1b
Authorized equipment: D – No printed or hand-written papers allowed. Approved, sim-
ple calculator allowed. Figure 1a shows a quadrilateral element with arbitrary shape (Q8). The element has 8 nodes
Two sheets of equations appended. and 2u 8 = 16 degrees of freedom. The element is mapped on to a square in the reference
coordinates [ and K that vary between 1 and 1, see Figure 1b. The relation between
”physical” coordinates x and y and reference coordinates [ and K is expressed by the map-
ping
This exam consists of 7 pages; 1 front page, 4 pages with problems, and 2 pages with
equations ­x ½
® ¾ > N @^c`
¯y¿
Note that most of the questions are independent of each other. The displacements are given by
­u ½
® ¾ > N @^d`
¯v¿
The grades will be announced 10 January 2005 or earlier (unless something unforeseen
occurs).
a) The element in Figure 1 is a C0 element. What does that mean?
b) What general criteria must the shape functions satisfy in order to guarantee convergence
towards exact solution when the element size is gradually reduced? What particular re-
quirements do these criteria entail for the shape functions of the C0 element in Figure 1?
Norwegian University of Science and Technology NTNU Norwegian University of Science and Technology NTNU
Faculty of Engineering Science and Technology Faculty of Engineering Science and Technology
Department of Structural Engineering Page 3 of 5 Department of Structural Engineering Page 4 of 5
Problem 2 (weight ca. 25 %)
c) Establish the shape functions N1 and N5 of the element in Figure 1.
d) The strain energy of the element may be expressed as
1 1
^` ^` dV ^d` >k @^d`
2³
T T
U
2
Use this expression for the strain energy to establish a relation for the element stiffness
matrix >k @ expressed in terms of matrix symbols.
e) Assume that the element has uniform thickness h and write the volume integral for >k @
C
L
in d) as a surface integral in reference coordinates. Establish the expression for numeri-
cal integration of >k @ and explain briefly the quantities that are involved.
f) What do you understand by full, reduced and selectively reduced integration? How
many integration points must be used in full integration of the element in Figure 1?
Substantiate the answer. (Note: A polynomial of degree 2n  1 is integrated exactly by
n-point Gauss quadrature.)
g) The stresses are calculated based on the nodal displacements of the element. Establish Figure 2
the matrix expression for calculation of the stresses and comment upon the accuracy.
What can you say about equilibrium? How can the nodal stresses of the element in Fig-
ure 1 be calculated based on the stresses in the Gauss points? a) A two-dimensional structure, see Figure 2, is completely fixed along the foundation,
and is subjected to wind loading that gives a uniform distributed pressure on one side
h) Consider an element of the type shown in Figure 1, but the shape is rectangular and the and an equally large suction on the opposite side. The structure is symmetric about a
edges are parallel to the x and y axes. The edges parallel to the x axis have length 2a, vertical axis through the apex. A static analysis is to be conducted for the given struc-
while the other edges have length 2b. The element is loaded by a constant body force ture, and you have available an isoparametric Q4 element having 4 nodes and 8 d.o.f.
where the component in the y direction is Fy per unit volume, while the component in Make use of the symmetry properties of the structure and sketch a proper mesh (be-
tween 20 and 30 elements) where you clearly specify the boundary conditions.
the x direction is zero. The thickness is uniform and equal to h. Determine the consis-
tent nodal loads in nodes 1 and 5 caused by the body force Fy . Show the directions of b) Explain briefly the main ingredients of the Patch test. What can be said about an ele-
the nodal forces in a sketch and comment upon the result. ment that fulfils the Patch test?
c) Give an example on how the Patch test may be conducted for the Q4 element. Make a
sketch of the mesh, and how the load distribution is for the different cases that you
would consider.
Norwegian University of Science and Technology NTNU Norwegian University of Science and Technology NTNU
Faculty of Engineering Science and Technology Faculty of Engineering Science and Technology
Department of Structural Engineering Page 5 of 5 Department of Structural Engineering Page 1 of 7
Problem 3 (weight ca. 25 %) Responsible contact during exam: ENGLISH
Professor Odd Sture Hopperstad, 93037174/73594703/73594700
C
EXAM IN COURSE TKT4191
b
y p0 FINITE ELEMENT METHODS 1
A B
x
a/2
Tuesday 13 December 2005
a/2
09.00 – 13.00
Figure 3
Figure 3 shows a triangular two-dimensional structure with uniform thickness h. The structure
is fixed in A and B as shown in the figure. The structure is made of an elastic material with Authorized equipment: D – No printed or hand-written papers allowed. Approved, sim-
Young’s modulus E, while Poisson’s ratio Q is equal to zero. Assume the displacement field ple calculator allowed.
Two sheets of equations appended.
u a1x  a 2 y, v a 3 y
where a , a and a are constants.
1 2 3 This exam consists of 9 pages; 1 front page, 6 pages with problems, and 2 pages with
equations
a) Start from the given generalized displacement field, and show that the displacements
may be expressed by direct interpolation in the form
­u B ½ Note that all of the questions are independent of each other.
­u ½ ª[  K / 2 K 0 º ° °
^u` ® ¾ « ®u C ¾ > N @^d`
¯v¿ ¬ 0 0 K¼» ° °
¯ vC ¿
The grades will be announced 4 January 2006 or earlier (unless something unforeseen
where u B is the displacement of B in the x direction, u C and vC are respectively the x occurs).
and y displacement of C, and [ x / a and K y / b .
b) Find the stiffness matrix >k @ for the triangular element expressed by Young’s modulus
E, area A and thickness h.
c) Calculate the consistent nodal load vector for a uniform line load p 0 along the edge BC.
Norwegian University of Science and Technology NTNU Norwegian University of Science and Technology NTNU
Faculty of Engineering Science and Technology Faculty of Engineering Science and Technology
Department of Structural Engineering Page 2 of 7 Department of Structural Engineering Page 3 of 7
Problem 1 (weight ca. 10 %) Problem 4 (weight ca. 10 %)
y, v v3 K y
4 7 3 u3 K 1
4 7 3
[ 3
8 6 8 6
9
9
5 b
1 2 1 5 2
K 1 v2
[ 1 [ 1 1 2
x, u u2 x
a
Figure 1a Figure 1b
Figure 2
Figure 1a shows a rectangular plane stress element with arbitrary shape (Q9). The element has
9 nodes and 2u 9 = 18 d.o.f. The element is mapped onto a square in the reference coordi-
nates [ and K that vary between 1 and 1, see Figure 1b. The displacement field for the element in Figure 2 is given by
Use the method of zero-lines to establish the shape functions for nodes 1, 5 and 9. § x y· x y § x y· x y
u ( x, y ) ¨1   ¸ u1  u2  u3 , v( x, y ) ¨ 1   ¸ v1  v2  v3
© a b¹ a b © a b¹ a b
Problem 2 (weight ca. 10 %) Assume that the element is given the nodal displacements
Give an account of u1 0, v1 0, u2 0, v2 0.001a, u3 0.001b, v3 0
x why it is necessary to use numerical integration to determine the stiffness matrix,
x what quadrature rule is mostly used to determine the stiffness matrix and why it is The element is in plane stress and the constants in Hooke’s law are E 70000 MPa and
used, and
Q 0.3 . The yield strength is f y 220 MPa .
x what is meant by full integration of an element.
Calculate the stresses V x , V y and W xy in the element, and use von Mises’ yield criterion to
Problem 3 (weight ca. 10 %) check if plastic yielding occurs for the actual stress state. According to von Mises’ yield crite-
rion, yielding occurs when
Give an account of what is meant by zero-energy modes (hourglass modes), when they may
occur, and to what extent they create problems in stress analysis.
Ve V x2  V y2  V xV y  3W xy2 t f y
Sketch the zero-energy modes of the Q4 and Q8 elements.
Norwegian University of Science and Technology NTNU Norwegian University of Science and Technology NTNU
Faculty of Engineering Science and Technology Faculty of Engineering Science and Technology
Department of Structural Engineering Page 4 of 7 Department of Structural Engineering Page 5 of 7
Problem 5 (weight ca. 10 %) Problem 7 (weight ca. 10 %)
What is meant by consistent loads in the finite element method? vi
1
­ Px ½ 8 ui
The element in Figure 2 is loaded by a concentrated load ^P` ® ¾ at the point
¯ Py ¿
­ x ½ ­a / 2 ½ 2
¾ . Calculate the consistent load vector ^re ` for this load case.
7
® ¾ ® y
¯ y ¿ ¯b / 2 ¿
3
6
Problem 6 (weight ca. 10 %) 4 5
x
y
H Figure 4
I
G
Figure 4 shows a rectangular plane stress element with 8 nodes. Each node has two d.o.f. ui
and vi as shown in the figure.
D F
E Using the method of generalized d.o.f., the displacements are interpolated within the element
as
u ¬« X ¼» ^au ` , v ¬« X ¼» ^av `
A B C
x where the column vectors ^au ` and ^av ` contain the generalized d.o.f., while the row vector
¬« X ¼» contains polynomials.
Figure 3
Determine the polynomials in ¬« X ¼» for the element in Figure 4 so that the compatibility
Convergence towards correct solution in the finite element method requires that the com- criterion is fulfilled across inter-element boundaries between elements of this type.
pleteness and compatibility criteria are fulfilled. Give an account of these two criteria for C0 Substansiate the answer. Sketch without any calculations the shape functions N1 and N 2 in
elements. 8 8
Figure 3 shows a patch of elements of the type CST, Q4, LST and Q8. Is the compatibility
the relations u ¦Nu
i 1
i i and v ¦N v
i 1
i i .
criterion generally fulfilled along all the inter-element boundaries in Figure 3? Substantiate
the answer.
Norwegian University of Science and Technology NTNU Norwegian University of Science and Technology NTNU
Faculty of Engineering Science and Technology Faculty of Engineering Science and Technology
Department of Structural Engineering Page 6 of 7 Department of Structural Engineering Page 7 of 7
Problem 8 (weight ca. 10 %) Problem 10 (weight ca. 10 %)
A D
3
5
y
B C 2
x
Figure 5 y 6 vi
Figure 5 shows a test coupon that is used for tensile testing of metals. The coupon is fixed
4 ui
along the boundary A-B and loaded by a uniform traction force in the positive x-direction x
along the boundary C-D. The coupon is modelled in FEMplate using plane stress elements of
the type Q8. Symmetry is invoked, and only ¼ of the coupon is modelled.
1
Sketch a mesh with about 60 Q8 elements for the part of the coupon that is modelled and
show the boundary conditions that you would apply for the actual problem. Try to make a
Figure 6
mesh that is as good as possible given the available number of elements. You do not need to
show the mid-side nodes of the Q8 element in the sketch.
Figure 6 shows an isoparametric plane stress element with two d.o.f. in each node. The thick-
ness of the element is constant and equal to t , and the area is A . The element is subjected to
a uniform body force
Problem 9 (weight ca. 10 %)
­ Fx ½
^F` ® ¾
The strain energy of an element is expressed in the form ¯ Fy ¿
1
³ 2 ^` ^` dV
T
Ui Determine the shape functions for node 1 and 4 using area coordinates, and calculate the con-
sistent nodal forces in the two nodes for the actual load.
The displacement field of the element is given by ^u` > N @^d`i , where ^d`i is the nodal Use that for area coordinates we have:
d.o.f. of element i, while the strain field is expressed by ^` > B @^d`i . Determine the strain A A A
³[ ³[ ³ [ [ dA
2
dA , dA ,
energy U of a patch of elements and identify the system stiffness matrix > K @ . The relation- A
i
3 A
i
6 A
i j
12
ship between local and global d.o.f. is given as ^d`i >L@i ^D` .
Norges teknisk-naturvitenskapelige universitet NTNU Norges teknisk-naturvitenskapelige universitet NTNU
Fakultet for ingeniørvitenskap og teknologi Fakultet for ingeniørvitenskap og teknologi
Institutt for konstruksjonsteknikk Page 1 of 9 Institutt for konstruksjonsteknikk Page 2 of 9
Responsible contact during exam: ENGLISH
Problem 1 (weight ca. 10 %)
Professor Odd Sture Hopperstad, 73594703/73594700
a) Discuss briefly the main computational ingredients in the finite element method.
b) The element stiffness relation can be expressed as >k @^d`  ^r` and >k c@^dc`  ^rc` in
two different coordinate systems. Assume that the relation between ^dc` and ^d` is given
EXAM IN COURSE TKT4191 by the transformation
^dc` > T@^d`
FINITE ELEMENT METHODS 1 where > T@ is a transformation matrix, and determine the relation between the nodal loads
^r` and ^rc` . Establish the relationship between the element stiffness matrices >k @ and
>k c@ in the two coordinate systems.
c) Show in detail how boundary conditions can be introduced into the global stiffness rela-
Monday 4. December 2006 tion
Kl 09.00 – 13.00 >K @^D` ^R`
Problem 2 (weight ca. 10 %)
Authorized equipment: D – No printed or hand-written papers allowed. Approved, sim-
ple calculator allowed.
Two sheets of equations appended.
3
b
This exam consists of 9 pages; 1 front page, 6 pages with problems, and 2 pages with y
equations 6 5
b
Note that most of the questions are independent of each other. 1 4 2
x
a a
The grades will be announced 4 January 2007 or earlier (unless something unforeseen
Figure 1
occurs).
Figure 1 shows an isosceles triangular plane stress element, i.e. two of the edges have the
same length. Geometry and nodes are illustrated in the figure. Origin in the (x,y)-coordinate
system is localized in node 4. Let the x-axis be coincident with the edge running through
nodes 1, 4 and 2, and let the y-axis pass through node 3.
Use the “method of zero lines” to determine the shape functions for nodes 3, 5 and 2.
Norges teknisk-naturvitenskapelige universitet NTNU Norges teknisk-naturvitenskapelige universitet NTNU
Fakultet for ingeniørvitenskap og teknologi Fakultet for ingeniørvitenskap og teknologi
Institutt for konstruksjonsteknikk Page 3 of 9 Institutt for konstruksjonsteknikk Page 4 of 9
Problem 3 (weight ca. 10 %) Problem 5 (weight ca. 10 %)
a) Discuss whether static equilibrium and kinematic compatibility are fulfilled within the y
element, on the boundary between elements and in the nodes in the finite element method. P
b) State the criteria for convergence towards the exact solution in C0 elements. q
4 7 3
c) If the highest degree polynomial in the shape functions for a C0 element is of order p , h/2
x
which order has the error in displacement and strain and what is the convergence rate of 8 6
the same quantities?
h/2
1 5 2
Problem 4 (weight ca. 10 %)
a a
x Figure 3
Figure 3 shows a rectangular Q8 element. The element is subjected to a uniform distributed
load q in positive y-direction along edge 4-7-3 and a concentrated load P in the positive y-
45
direction in the point  x, y   a / 2, h / 2  . Calculate the consistent nodal loads for this load
a case.
-t-
Problem 6 (weight ca. 10 %)
a y y y
V0 W0
a a C P
Figure 2
Figure 2 shows a triangular plane structure. Along the edge x = 0 a normal stress V 0 y / a and b
a shear stress W 0 y / a are acting as shown in Figure 2. The structure is in a state of plane
y
stress, and the stress field is given by
x y x y x y
V x W0 V0 , V y (3W 0  2V 0 )  (2W 0  V 0 ) , W xy (2W 0  V 0 )  W 0 A B
a a a a a a
x
Establish the static boundary conditions for the structure along the vertical and inclined edges,
a/2 a/2
and check whether the structure is in static equilibrium. It is assumed that body forces can be
neglected.
Figure 4
Norges teknisk-naturvitenskapelige universitet NTNU Norges teknisk-naturvitenskapelige universitet NTNU
Fakultet for ingeniørvitenskap og teknologi Fakultet for ingeniørvitenskap og teknologi
Institutt for konstruksjonsteknikk Page 5 of 9 Institutt for konstruksjonsteknikk Page 6 of 9
Figure 4 shows a triangular plane structure with constant thickness t. The structure is fixed in ­ w ½ ­w½
points A and B as shown in the figure. It is made of an elastic material with Young’s modulus °° w[ °° ° wx °
E, while the Poisson number  is zero. ® ¾ > J @ ®° w °¾
° w ° ° °
Assume the displacement field ¯° wK ¿° ¯° wy ¿°
a1 x  a2 y, v b) Let the inverse Jacobian matrix >  @ >J@
1
u a3 y be given by
where a1 , a2 and a3 are constants. Determine the expression for the strain energy U and the
ª *11 *12 º
>@ >J @
1
load potential : for the structure when it is loaded by a concentrated force P in positive x- « »
direction in point C. Determine the displacement of point C in terms of a , b , E and P . ¬* 21 * 22 ¼
so that
Problem 7 (weight ca. 10 %)
­w½ ­ w ½
°° wx °° ° w[ °
For incompatible elements the Patch test is used to determine whether the element converges ® ¾ > @ ®° ¾°
towards the exact solution with mesh refinement. Discuss how the Patch test can be done for a °w° ° w °
plane stress element, and explain why an element that fulfils the test also satisfies the conver- ¯° wy ¿° °¯ wK °¿
gence requirements.
Determine the expressions for the normal strains H x and H y as functions of the reference co-
Problem 8 (weight ca. 10 %)
ordinates and the nodal displacements.
For isoparametric plane stress elements the relation between ”physical” coordinates x and y
and reference coordinates [ and K is expressed by the mapping
n n
Problem 9 (weight ca. 10 %)
x ¦ N [ ,K  x ,
i i y ¦ N [ ,K  y
i i
i 1 i 1 Assume that the Q8 element in Figure 3 has thickness with quadratic variation over the ele-
ment in accordance with the shape functions, i.e.
where  xi , yi  are nodal coordinates, N i N i [ ,K  are shape functions and n is the number
8
of nodes in the element. The displacements u and v are interpolated by use of the same
shape functions, i.e.
t ¦ N [ ,K  t
i 1
i i
n n
where ti is the nodal thickness. Determine the quadrature rule (Gauss) so that the stiffness
u ¦ N [ ,K  u ,
i 1
i i v ¦ N [ ,K  v
i 1
i i
coefficients kij (i.e. the elements of the stiffness matrix >k @ ) are integrated exactly. Substan-
tiate the answer.
where  ui , vi  are nodal displacements.
a) Determine the Jacobian matrix > J @ for an arbitrary isoparametric plane stress element,
when it is given that
Norges teknisk-naturvitenskapelige universitet NTNU Norges teknisk-naturvitenskapelige universitet NTNU
Fakultet for ingeniørvitenskap og teknologi Fakultet for ingeniørvitenskap og teknologi
Institutt for konstruksjonsteknikk Page 7 of 9 Institutt for konstruksjonsteknikk Page 8 of 9
Problem 10 (weight ca. 10 %)
Formler/Formulas
|
Elastisitetsteori/Theory of elasticity
y
Likevekt/Equilibrium:
wV x wW xy wW zx wW xy wV y wW yz wW zx wW yz wV z
M    Fx 0,    Fy 0,    Fz 0
M wx wy wz wx wy wz wx wy wz
x
>w @ ^`  ^F` ^0`
T
Kinematisk kompatibilitet/Kinematic compatibility:
Figure 5 wu wv ww wu wv wv ww ww wu
Hx , Hy , Hz , J xy  , J yz  , J zx 
wx wy wz wy wx wz wy wx wz
Figure 5 shows a plane specimen with a circular hole subjected to a moment M. The specimen
is thin, and thus plane stress may be assumed. Owing to the symmetry of the specimen only
the shaded quarter part in the figure shall be analyzed by the finite element method. ^` >w @^u`
Sketch an element model with only 3 isoparametric Q8 elements for the reduced structure,
and show the boundary conditions in the sketch. What quadrature rule would you choose for
the element? Substantiate the answer. Material lov/Material law:
Assume that you have the choice between Q8 and Q9 elements for the actual problem. Which 1 1 1
Hx (V x Q (V y  V z )), H y (V y Q (V z  V x )), H z (V z Q (V x  V y ))
element type would you choose? Substantiate the answer. E E E
1 1 1 E
J xy W xy , J yz W yz , J zx W zx ; G
G G G 2 1  Q 
^` > E@ ^` ^` >E@^`
1
or
Statiske randkrav/Static boundary conditions:
)x V x l  W xy m  W zx n, ) y W xy l  V y m  W yz n, ) z W zxl  W yz m  V z n
^` > @^n`
Norges teknisk-naturvitenskapelige universitet NTNU Norges teknisk-naturvitenskapelige universitet Studentnummer:
Fakultet for ingeniørvitenskap og teknologi Institutt for konstruksjonsteknikk
Student number:
Institutt for konstruksjonsteknikk Page 9 of 9
Faglig kontakt under eksamen:
Energiprinsipper/Energy principles: Institutt for konstruksjonsteknikk, Gløshaugen
Professor Kolbein Bell, tlf 94679 / mobil 45 666 123
Virtuelle forskyvningers prinsipp/Principle of virtual displacements:
³ ^G ` ^` dV ³ ^G u` ^F` dV  ³ ^G u` ^` dS
T T T
Eksamen i emne
Prinsippet om stasjonær potensiell energi/Principle of stationary potential energy: TKT4191 Elementmetoden 1
G3 p ³ ^G ` >E@^`  ^G ` >E@^ ` dV
T T
0
Mandag 17. desember, 2007
Tid (time): kl 0900 - 1300 (4 timer - 4 hours)
 ³ ^G u` ^F` dV  ³ ^G u` ^` dS  ^G D` ^P`
T T T
0
Elementmetoden/Finite element method: Hjelpemidler: Ingen trykte eller håndskrevne hjelpemidler tillatt.
Bestemt enkel kalkulator tillatt.
Forskyvninger/Displacements: Aid: No printed or handwritten notes are permitted.
A specific, simple calculator is permitted.
^u` > N @^d`
Sensur: Om intet uforutsatt inntreffer vil sensuren foreligge innen 11.01.2007.
Censoring: Examination results should be available not later than 11.01.2007.
Tøyninger/Strains:
^` > B@^d` Før svarene direkte på oppgavearket - bruk baksiden om du får for lite plass. På side
17 finner du en del formler til bruk hvor det måtte passe; dette arket kan rives fra
(behøver ikke leveres).
Merk at for noen av spørsmålene skal du bare gi svaret (i en boks som er avsatt for
Stivhetsmatrise/Stiffness matrix:
dette formål). Eventuelle utregninger gjøres på kladd, men de skal ikke tas med i
besvarelsen.
>k @ ³ > B@ >E@> B@ dV
T
Betydningen (vekten) av hvert problem er angitt i %.
Answer on these pages - use the back if you need more space. On page 17 you will
find some formulas for use at your discretion; this sheet need not be handed in.
Konsistent lastvektor/Consistent nodal load vector:
Note that for some questions you are only asked to give the final answer (in a box pro-
vided for this purpose). Computations, if any, are carried out on auxiliary paper not
^re ` ³ >B @ > E@^ 0 ` dV  ³ > N @ ^F` dV  ³ > N @ ^` dS
T T T
to be handed in.
The importance (the weight) of each problem is indicated in %.
Stivhetsrelasjonen/Stiffness relation:
>k @^d` ^re ` Husk å skrive studentnummeret ditt i boksen øverst til høyre på hvert ark.
Remember to write your student number in the box at the top of each page.
Side 1 av 18
TKT4191 Elementmetoden 1 Studentnummer: TKT4191 Elementmetoden 1 Studentnummer:
Eksamen - desember 2007 Student number: Eksamen - desember 2007 Student number:
Problem 1 (10 %) b) Still opp systemets totale potensielle energi og etabler stivhetsrelasjonen, P = KG, ved
hjelp av prinsippet om minimum potensiell energi (PMPE)
P Establish the total potential energy of the system and determine the stiffness relation,
A B P = KG, by the principle of minimum potential energy (PMPE)
k k G
L/3 L/3 L/3
Figur 1
Bjelke A-B i figur 1 antas å være fullstendig stiv. Systemet har dermed 1 frihetsgrad, f.eks G.
Beam A-B in figure 1 is completely rigid. Hence the system has 1 degree og freedom, e.g. G.
a) Etabler stivhetsrelasjonen for systemet, P = KG, ved en likevektsbetraktning og bestem G.
Use equilibrium to establish the stiffness relation, P = KG, for the system and determine G.
Problem 2 (10 %)
Samme system som i figur 1, men nå er bjelken fleksibel med bøyestivhet EI.
Same system as in figure 1, but the beam is now flexible with bending stiffness EI.
a) Problemet skal løses ved hjelp av metoden til Rayleigh-Ritz. Anta en (1) fri parameter
(generalisert koordinat) og velg den enkleste tillatte form for tverrforskyvningen v som
også tar hensyn til bøyestivheten. Tøyningsenergien ved bøyning er gitt i formelarket bak-
erst. Sett EI = 0,1975 ˜ kL og bestem nedbøyningen G.
3
The problem is to be solved by the Rayleigh-Ritz method. Assume one (1) free parameter
(generalized coordinate) and choose the simplest permissible shape for the transverse dis-
placement v that also take into consideration the bending stiffness of the beam. You find
the strain energy due to bending amongst the formulas at the back.
Let EI = 0,1975 ˜ kL and determine the displacement G.
3
Side 2 av 18 Side 3 av 18
TKT4191 Elementmetoden 1 Studentnummer: TKT4191 Elementmetoden 1 Studentnummer:
Eksamen - desember 2007 Student number: Eksamen - desember 2007 Student number:
Problem 3 (30 %)
y v v3 K
5 u3
3
4 3
4 4
avbildning 1
[
Q4
mapping
1
1 2 [2 = 1
1 2 x 1 1 K2 = 1
u
3 6
Formfunksjoner (shape functions): Ni = 1 (1 + [i[)(1 + KiK)
4
Figur 2
Figur 2 viser det enkle firkantelementet Q4 (med 8 frihetsgrader) i både fysiske og naturlige
koordinater. Formfunksjonene for dette C0 elementet, i naturlige koordinater, er også vist.
Figure 2 shows the simple quadrilateral C0 element Q4 (with 8 degrees of freedom) in both
physical and natural coordinates. The shape functions, in natural coordinates, are also shown.
a) Hvordan uttrykkes elementets geometri ved en isoparametrisk formulering, dvs hvordan
avbildes elementet fra fysiske til naturlige koordinater? Sett opp avbildningen på matrise-
form og forklar symbolene.
How is the geometry of the element expressed in an isoparametric formulation, or in other
words, how is the element mapped from physical to natural coordinates? Write the map-
ping in matrix notation and explain the symbols.
b) Hvilken forskyvningsantakelse vil du velge for v med 2 frie parametre? Gi bare svaret.
Which displacement assumption will you choose for v with 2 free parameters? Give the
answer only.
v=
Side 4 av 18 Side 5 av 18
TKT4191 Elementmetoden 1 Studentnummer: TKT4191 Elementmetoden 1 Studentnummer:
Eksamen - desember 2007 Student number: Eksamen - desember 2007 Student number:
b) Tøyningene er definert som - the strains are defined as
­ H ½
° x ° ­ u ½ N 1 0 ..... 0
® Hy ¾ = >w @® ¾ = > w @> N @^ d ` = > B @^ d` where >N@ = and
° ° ¯ v ¿ 0 N 1 ..... N 4
¯ J xy ¿
T
^d` = u1 v1 . . . . u4 v4 er forskyvningsvektoren (is the displacement vector).
Etabler operatormatrisen > w @ - se formelarket bakerst.
Establish the operator matrix > w @ - see formulas at the back.
d) Bestem Jacobi-determinanten i knutepunktene 2 og 4 for elementet i figur 2. Angi bare
svarene.
Determine the Jacobi determinant at nodes 2 and 4 of the element in figure 2. Answers
c) For å kunne bestemme matrisen [B] trenger vi Jacobi-matrisen [J] definert som only.
­ ½ ­ ½
° w-
----- ° ° w
----- ° J2 = J4 =
° w[ ° ° wx °
® ¾ = >J@® ¾
° w ° ° w °
° ------ ° ° ----- °
wK wy T T
¯ ¿ ¯ ¿ e) Etabler via vektorene u ,x u ,y v ,x v ,y og u ,[ u ,K v ,[ v ,K et uttrykk , på
matriseform, for matrisen [B] for et generelt isoparametrisk C0 element.
Bruk kjerneregelen (se formlene bakerst) til å etablere generelle uttrykk for elementene i
[J], og bestem Jacobimatrisen for elementet i figur 2. Establish, via the two vectors above, an expression, in matrix form, for the matrix [B] of a
In order to determine matrix [B] we need the Jacobi matrix [J] as defined above. general isoparametric C0 element.
Use the chain rule (see formulas at the back) to establish general expressions for the
elements of [J] and determine the Jacobi matrix for the element in figure 2.
Side 6 av 18 Side 7 av 18
TKT4191 Elementmetoden 1 Studentnummer: TKT4191 Elementmetoden 1 Studentnummer:
Eksamen - desember 2007 Student number: Eksamen - desember 2007 Student number:
g) Q4 elementet har følgende uavhengige forskyvningsmønstre:
The Q4 element has the following independent displacement modes:
1 2 3 4
f) Hvorfor benytter vi numerisk integrasjon ved beregning av stivhetsmatrisen [k] for iso-
parametriske elementer? 5 6 7 8
Why do we use numerical integration for the computation of the stiffness matrix [k] for
isoparametric elements?
Hva representerer disse ved full integrasjon?
What do these modes represent in case of full integration?
Hva skjer ved redusert integrasjon og hvorfor?
Hvor mange integrasjonspunkter må vi ha for å integrere et 6. grads polynom i 1 dimen- What happens in case of reduced integration and why?
sjon eksakt med Gauss kvadratur?
How many integration points do we need in order to integrate a polynomial of degree 6 in
1 dimension exactly by use of Gauss quadrature?
n=
Hva mener vi med full og redusert integrasjon av [k], og hvor mange integrasjonspunkter
krever disse integrasjonene for Q4?
What does full and reduced integration of [k] mean, and how many integration points do
these integrations require for the Q4 element?
Side 8 av 18 Side 9 av 18
TKT4191 Elementmetoden 1 Studentnummer: TKT4191 Elementmetoden 1 Studentnummer:
Eksamen - desember 2007 Student number: Eksamen - desember 2007 Student number:
Problem 4 (20 %)
Vi betrakter en feltstørrelse I (f.eks forskyvningskomponenten i x- eller y-retning) som må til-
fredsstille C0-kontinuitet mellom naboelementer.
We consider a field variable I (e.g. displacement component in x- or y-direction) that is requi-
red to satisfy C0-continuity between neighbouring elements. N1 =
a) Hva mener vi med C0-kontinuitet?
What do we mean by C0-continuity?
N3 =
b)
K
K=1 N5 =
4 3
8 7 [ Figur 3
1 5 6 2
K= 1
[= 1 1
[=
[= 1 3 [=1
3
Figur 3 viser et spesielt firkantelement avbildet i naturlige koordinater [,K. Elementet har 8
frihetsgrader, I1, I2, ...., I8, definert i de viste knutepunktene. N7 =
Figure 3 shows a special quadrilateral element mapped onto a square in natural coordinates
[,K. The element has 8 degrees of freedom, I1, I2, ...., I8 , at the nodes shown.
Bestem formfunksjonene N1, N3, N5 og N7 slik at kravene til kontinuitet mellom naboelementer
med felles knutepunkter er tilfredsstilt. Angi kun funksjonene i boksene nedenfor og skisser
hver funksjon som en “flate” på elementet.
Determine the shape functions N1, N3, N5 og N7 in such a way that interelement continuity is
maintained between elements with common nodal points. Give the functions only, and sketch
each function as a “surface” on the element.
Side 10 av 18 Side 11 av 18
TKT4191 Elementmetoden 1 Studentnummer: TKT4191 Elementmetoden 1 Studentnummer:
Eksamen - desember 2007 Student number: Eksamen - desember 2007 Student number:
c) Problem 5 (8 %)
Trekantelementet i figur 4 har 6 frihetsgrader fordelt på tre For triangulære C0-elementer er det sterkt ønskelig å interpolere den feltvariable størrelsen,
3 hjørnenoder og tre midtside noder. eventuelt størrelsene, med fullstendige polynomer i x og y. Dette gir knutepunktsplasseringer
5 Bestem alle 6 formfunksjonene uttrykt ved arealkoordinater. for elementer basert på lineære, kvadratiske og kubiske polynomer som vist i figur 5. Alle
6 disse elementene tilfredsstiller også C0-kontinuitet mellom elementene.
Angi bare svarene.
4 2 The triangular element in figure 4 has 6 degrees of freedom, For triangular C0-elements it is very desirable to interpolate the field variable(s) by complete
1 polynomials in x and y. This gives the nodal point arrangements for elements based on linear,
one at each corner node and one at each mid-side node.
Determine all 6 shape functions in terms of area coordinates. quadratic and cubic polynomials shown in figure 5. All three elements satisfy C0- continuity
Figur 4 between elements.
Answers only.
N1 = N4 =
N2 = N5 =
lineær - linear kvadratisk - quadratic kubisk - cubic
N3 = N6 =
Figur 5
Hvor mange knutepunkter vil du måtte ha for det neste fullstendige polynom, dvs et fullstendig
d) Et element for spenningsanalyse må kunne beskrive stivlegemebevegelsene korrekt. 4. grads polynom, og hvor vil du plassere knutepunktene når du også ønsker full C0-kontinuitet
Hvilke krav stiller dette til formfunksjonene for et C0-element? mellom naboelementer? Vis plasseringe på figuren.
Kontroller at formfunksjonene i spørsmål c) tilfredsstiller kravet i trekantens flatesenter. How many nodal points will you need for the next complete polynomial, that is a complete 4th
An element for stress analysis must represent rigid body movements correct. order polynomial, and where will you place them in order to also maintain inter-element C0-
Which requirement does this impose on the shape functions for a C0-element? continuity? Indicate your proposal on the figure.
Check that the shape functions in question c) satisfy this requirement at the element center.
Side 12 av 18 Side 13 av 18
TKT4191 Elementmetoden 1 Studentnummer: TKT4191 Elementmetoden 1 Studentnummer:
Eksamen - desember 2007 Student number: Eksamen - desember 2007 Student number:
Problem 6 (10 %) Problem 7 (12 %)
For Q4 elementet i figur 2 er sidekanten 2-3 en ytre D2
y v rand. På den virker det en jevnt fordelt last (traksjon)
p
p (kraft per lengdeenhet), som vist på figur 6. D1
5
3 Bestem elementets lastvektor {re} for den viste
4 4 belastning. skive stav - bar
(a) ( Q4 )
Side 2-3 of the Q4 element in figure 2 is an external plate
Q4 D7
edge on which is applied a uniform load (traction) of D4 bjelke
intensity p (force per unit length), as shown in figure stav
6. bar D6 beam
1 x D3 D5 D8
2 Determine the element load vector {re} for the
u loading shown.
3 6
d8 d1
d6
Figur 6 d2
d5 C
d7 (b)
d4
A d3
d2 d5
d2 d4 D
d1 d3 B d1 d4
d1 d3 d3
d2 d4 d6
Figur 7
Figur 7a viser en plan konstruksjon som består av 4 elementer, ett skiveelement (Q4), to stav-
elementer og ett bjelkeelement, med totalt 8 ukjente frihetsgrader, D1 til D8 som vist i figuren.
Lokal nummerering av elementfrihetsgradene er vist i figur 7b.
Figure 7a shows a plane structure consisting of 4 elements, one plate element (Q4), two bar
elements and one beam element, with a total of 8 degrees of freedom, D1 to D8, as shown in the
figure. Local numbering of the element degrees of freedom is shown in figure 7b.
a) Anta at stivhetsmatrisene for de fire elementene (med frihetsgrader referert til globale
akser) er angitt som:
[kA] = [ aij ] [kB] = [ bij ] [kC] = [ cij ] [kD] = [ dij ]
Er noen av elementene i disse elementmatrisene lik 0 (null), og i så fall hvilke?
Assume that the stiffness matrices of the four elements (with degrees of freedom referred to
global axes) are designated:
[kA] = [ aij ] [kB] = [ bij ] [kC] = [ cij ] [kD] = [ dij ]
Are any of the elements of these element matrices equal to 0 (zero), and if so which?
Side 14 av 18 Side 15 av 18
TKT4191 Elementmetoden 1 Studentnummer: TKT4191 Elementmetoden 1
Eksamen - desember 2007 Student number: Eksamen - desember 2007
Formler - Formulas
Likevekt - equilibrium:
wV x wW xy wW zx wW wV wW wW wW wV
-------- + --------- + --------- + F x = 0 , --------xy- + --------y + --------yz- + F y = 0 , --------zx- + --------yz- + --------z + F z = 0
wx wy wz wx wy wz wx wy wz
Kinematikk - kinematics:
wu wv ww wu wv wv ww ww wu
H x = ------ , H y = ----- , H z = ------- , J xy = ------ + ----- , J yz = ----- + ------- , J zx = ------- + ------
wx wy wz wy wx wz wy wx wz
Elastisk tøyningsenergi - elastic strain energy:
³
1 T
Generelt - in general: U = --- ^ H ` > E @ ^ H ` dV Pascals trekant
2
V Pascal’s triangle
b) Angi, i rutene nedenfor, elementene i systemstivhetsmatrisen [K], uttrykt ved element-
1
³
stivhetene aij, bij, cij og dij (erstatt i og j med tall). Det er tilstrekkelig å angi alt på og over 1 2
Bøyning av bjelke - bending of beam: U = --- EIv ,xx dx
hoveddiagonalen. 2 x y
L
Indicate, in the boxes below, the elements of the system stiffness matrix [K], expressed in x2 xy y2
T
terms of the element stiffnesses aij, bij, cij og dij (replace i and j by numbers). It is sufficient Potensiell energi av ytre belastning: : = –^ R ` ^ D `
to show the content of the boxes on and above the main diagonal. Potential of external loading x3 x2y xy2 y3
.
Total potensiell energi: 3p = U + : .
1 2 3 4 5 6 7 8 .
Total potential
1 Prinsippet om minimum potensiell energi: G3 = 0
Principle of minimum potetial energy
2 w- w wx w wy w w wx w wy
Kjerneregelen - chain rule : ----- = ----- ------ + ----- ------ ------ = ----- ------ + ----- ------
w[ wx w[ wy w[ wK wx wK wy wK
Element matriser - element matrices:
3
³ > B @ > E @ > B @ dV
T
Stivhet - stiffness: > k @ =
4 V
K =
³ ³ ³
T T T
Lastvektor - load vector: ^ r ` = > B @ > E @ ^ H ` dV + > N @ ^ F ` dV + > N @ ^ ) ` dS
e 0
5
V V S
Arealkoordinater - area coordinates:
6
7
8
Side 16 av 18 Side 17 av 18
TKT4191 Elementmetoden 1 Studentnummer:
Eksamen - desember 2007 Student number:
EKSAMEN I EMNE TKT4191
ELEMENTMETODEN 1
3 Arealkoordinatene til punkt P:
Area coordinates of point P:
[1 = 0
Høsten 2004
A
[ i = ----i i = 1,2,3 og / and A = A1 + A2 + A3
[2 = 0 A1 A Løsningsforslag
A2
P Ÿ [1 + [2 + [3 = 1
A3 2
Oppgave 1
1 [3 = 0
a) Et C0-element må tilfredsstille såkalt C0-kontinuitet som betyr at vi mellom
naboelementer må ha kontinuitet i selve forskyvningen(e), men ikke i dens deriverte
(egentlig kontinuitet i til og med den 0. deriverte). Egentlig Cm-1 hvor m er
variasjonsorden (graden av derivert i funksjonalen) og 2m er orden på problemets
differensialligning.
b) Kontinuitetskriteriet:
Formfunksjonene må være kontinuerlige og m-deriverbare innen elementet, og de må
sikre Cm-1 kontinuitet mellom naboelementer (m er variasjonsordenen).
Formfunksjonene må være entydige og slik at for C0-elementer:
x Ni = 1 i knutepunkt i og lik null i alle andre knutepunkter.
x Ni må være identisk lik null langs alle sidekanter/flater som ikke omfatter
knutepunkt i.
Fullstendighetskriteriet:
Formfunksjonene må være slik at elementet kan representere vilkårlige konstante
tøyningstilstander, og det må kunne beskrive alle stivlegemebevegelser korrekt (uten
at det oppstår spenning, såkalt “self-straining”).
For det aktuelle problemet (gjelder for u og v:):
x u c0  c1x  c 2 y Ÿ ui c0  c1x i  c 2 yi
x u ¦N u i i c0 ¦ Ni  c1 ¦ Ni x i  c2 ¦ Ni yi for isoparametrisk formulering.
x De er identiske dersom ¦N i 1
1
Side 18 av 18
c) Formfunksjonene: La [i og Ki være koordinater til integrasjonspunkt nr i langs [ -aksen og nr j
 langs K -aksen og la
N1 c1 (K  1)([  1)([  K  1)
¬ª  [i , Ki  ¼º ¬ª B  [i , Ki  ¼º > E@ ¬ªB  [i , Ki  ¼º
 -1=0 T
ª¬) ij º¼
c (2)(2)(1) 1
1
da er
Ÿ c1  14

1 1 n n
>k @ h³ ³ >B@ >E@>B@ Jd[dK | h ¦¦ W W ) J  [ , K 
T
1 5 N5 c5 ([  1)([  1)(K  1) i j ij i j
1 1 i 1 j 1
 -1=0 c5 (1)(1)(2) 1
 +1=0  ++1=0
Ÿ c5 1
2
Her er Wi og Wj “vekttallene” i en integrasjonsregel med n punkter og J([i , K j ) er
verdien av Jakobi-determinanten i integrasjonspunktet ([i , K j ) .
N1 1
4 (1  [)(K  1)([  K  1) f) Full integrasjon innebærer at vi integrerer >k @ eksakt når elementet har basisform, det
Altså
N5 1
(1  K)(1  [ 2 ) vil si en form som gir konstant Jakobi-determinant – velger da selvsagt den “laveste”
2
Gauss-regel som sikrer dette. Redusert integrasjon betyr at vi bruker en
integrasjonsregel med færre punkt enn det full integrasjon krever (men vi bruker same
regel for alle bidrag til >k @ .
d) Elementstivhetsmatrisen >k @ :
Selektiv integrasjon betyr at vi bruker forskjellig integrasjons-regel for forskjellige
^d` >k @^d` ³ ^` ^` dV ³ ^` >E@^` dV
T T T
bidrag til >k @ (f.eks. full integrasjon av bøyningsbidrag og redusert integrasjon av
1 1 1
2 2 2
V V
skjærbidrag).
T ª º
2 ^d` « ³ > B @ > E @> B @ dV » ^d`
1 T
¬V ¼ Full integrasjon av vårt 8-noders element krever en 3 u 3 regel. I > B@ -matrisen har vi
nemlig ledd med både [ 2 og K2 , og dette gir 4.grads ledd i >@ . Detter betinger 3
>k @ ³ >B @ >E@> B@ dV
T
Ÿ integrasjonspunkter (som integrerer en 5.grads funksjon eksakt, mens 2 punkter bare
V integrerer en 3.grads funksjon eksakt.)
g) Spenning:
e) Stivhetsmatrisen:
^` >E@^` >E@> w @^u` >E@>w @> N@^d` >E@> B@^d`
>k @ ³ > B@ > E@> B @ dV ³ h > B@ > E@> B @ dA
T T
V A
Vi avleder altså spenningene ved å derivere (via > w @ ) et forskyvningsfelt
1 1
( ^u` > N@^d` ) som allerede er beheftet med feil, og dermed må vi regne med at feilen
h ³ ³ > B @ > E@> B @Jd[dK
T
1 1 i ^` blir større enn feilen i ^u` ( ^d` ).
hvor dA Jd[dK . Likevekten tilfredsstilles ikke punktvis (heller ikke inne i elementet), men vi har en
midlere (integrert) likevekt.
Spenningene beregnet etter ^` >E@>B@^d` er mest nøyaktige i integrasjonspunktene
i en integrasjonsregel som er en orden lavere enn det som tilsvarer full integrasjon
(henger sammen med integrasjonspunktenes egenskaper, og bekreftes av resultater).
2 3
 Beregner spenningene i 2 u 2 Gauss-punkter, og ekstrapolerer verdiene ut til
knutepunktene via et bilineært felt gjennom de 4 integrasjonspunktene. Fx1 Fx5 0
1 1 1 1
Fy1 hab ³ ³ N1Fy d[dK 1
4 Fy hab ³ ³ ([2  [ 2K  K2 [K2  [K  1)d[dK
1 1 1 1
u u 1
1
4 F hab ³ (2[
y
2
  [  2)d[
2
3
2
3  13 Fy hab
1
1 1 1 1
Fy5 hab ³ ³ N1Fy d[dK 1
2 Fy hab ³ ³ (1  [ 2  K [ 2 K)d[dK
u u 1 1 1 1
1
1
2 Fy hab ³ (2  2[ 2 )d[ 4
3 Fy hab
1
2×2 Gauss
Vi får samme resultat i 2, 3 og 4 som i 1, og i 6, 7 og 8 som i 5, dvs.
h)
y, 
P Q P
4 3
7 K y
b Ÿ dy bdK P 1
3 Fy hab
b [ x
Ÿ dx ad[
Fy 6 a Fy Q 4
3 Fy hab
8 Q Q
x,  dA dxdy abd[dK 4Q  4P 4Fy hab
b 0
1 5 2
P Q P
a a
Konsistent (kinematisk ekvivalent) lastvektor: Likevekten stemmer altså, men vi har fått det ”pussige” resultat at, selv med vanlige
1 1
(standard) frihetsgrader går kreften i hjørnepunktene den ”gale” veien!
^re ` ³ > N @ ^F` dV ³ h > N @ ^F` dA ³ > N @ ^F` d[dK
T T
hab ³
T
V A 1 1
ªN0 0º
>N@ > N0 @ > N1 N 2 " N8 @
T
« N 0 »¼
¬0
­0½
^F` ® ¾
¯Fy ¿
4 5
Oppgave 2 c) Elementnett:
a) Elementnettet:
Parallellogramform
ved oppleggene
Case 1: V x konst.
b) Patchtesten er en numerisk test som utføres for en gruppe av elementer utsatt for et Vx 1 Vx 1
lineært forskyvningsfelt, og som benyttes til å undersøke om et element er gyldig
med hensyn til konvergens mot korrekt løsning ved nettforfining. 3 6 9
Prosedyre:
x Lag en gruppe av elementer med irregulær form med minst en intern node.
x Innfør opplagerbetingelser som er akkurat tilstrekkelig til å unngå
stivlegemebevegelser. 5 8
x Foreskriv laster på en eller flere av rendene til elementgruppen som er 2
arbeidsekvivalente med en konstant spenningstilstand.
x Indre noder i elementgruppen er verken belastet eller fastholdt.
x Løs ligningssystemet og sammenlign resultatene med eksakt løsning. 1 4 7
x Hvis de beregnede resultatene er korrekte I all punkter, er patchtesten
oppfylt.
Et stabilt element, dvs. et element uten timeglassmoder, som tilfredsstiller 3 6 9 2.5
3
patchtesten kan beskrive:
x Stivlegemebevegelser uten tøyning
x Konstante tøyningstilstander
x Kompatibilitet med naboelementer for konstante tøyningstilstander.
5 8 5
2
5
1 4 7
6 7
Case 2: V y konst. Case 3: W xy konst.
W xy
Vy 1
3 6 9
3 6 9
5 8
5 8 2
2
1 4 7
1 4 7
W xy
Vy 1
2.5 5 2.5
3 6 9 2.5
3
5
2.5 2.5
3 6 9
5
5 8
5
2
5 8
2 1 4 7 2.5
5
1 4 7
5
8 9
 T
ª b 0 0 º ª1 0 0 º ª b 0 0º
>k @ ³ «« 0 a »» E ««0 1 0 »» «« 0 0 a »» dV
Oppgave 3 1 1
Ÿ 0
ab ab
a) Interpolasjonspolynom > N @ : ¬« a / 2 a 0 ¼» ¬«0 0 1/ 2 ¼» ¬« a / 2 a 0¼»
V
ª b 2  a 2 / 8 a 2 / 4 0 º ª1  8b 2 / a 2 2 0 º
­u B ½ ­ u(a, 0) ½ ª a 0 0 º ­ a1 ½ ­ a1 ½ EAh « » EAh « »
° ° ° ° « a / 2 b 0 » °a ° a 2 / 4 a2 / 2 0 » 2
> A @ °®a 2 °¾
4 0»
^d` ®u C ¾ ® u(a / 2, 0) ¾ a 2b2 « 8b 2 «
« »® 2¾ « 0 0 a 2 ¼» « 0 0 8 ¼»
° ° ° ° ° ° ° ° ¬ ¬
¯ vC ¿ ¯ v(a / 2, b) ¿ ¬« 0 0 b ¼» ¯ a 3 ¿ ¯a 3 ¿
­ a1 ½ ª b 0 0 º ­u B ½
° ° 1 « ° °
a 0 »» ®u C ¾
c) Lastvektor:
Ÿ ®a 2 ¾  a / 2
° ° ab «
¯a 3 ¿ ¬« 0 0 a ¼» ¯° v C ¿° L  x½
^r ` ³ ¬ª N ¼º ^` dS ³ L1 ®­ ­1½ p 0 L
L
T

e ¾ p 0 dx ®¾
S ¯ 0
x ¿ ¯1¿ 2
(Komponentene i ^re ` har samme retning som p 0 .)
­ a1 ½ ª b 0 0º ­u B ½
ªx y 0º ° ° ªx y 0º 1 « »° °
^u` « » ®a 2 ¾ « » « a / 2 a 0 » ® u C ¾
¬ 0 0 y¼ ° ° ¬ 0 0 y ¼ ab «
¯a 3 ¿ ¬ 0 0 a ¼» ¯° v C ¿°
p0L/2
­u B ½
1 ª xb  ya / 2 ya 0 º° °
®u C ¾
ab ¬« 0 0 ya ¼» ° °
¯ vC ¿ C
p0
[ x /a½ 1 ª xb  ya / 2 ya 0º ª[  K / 2 K 0 º
¾ Ÿ >N@ p0L/2
K y / b¿ ab ¬« 0 0 ya ¼» «
¬ 0 0 K¼» L

b) Stivhetsmatrisen >k @ :
B
>k @ ³ >B @ >E@>B @dV
T
V
a b
sin D cosD =
ªw 0º ª b 0 0º 2L L
« wx » 1 ª xb  ya / 2 ya 0º 1 «
>B @ > w @> N @ «0
w
wy » « 0 0 a »»
«w » ab ¬ 0 0 ya ¼» ab « ­ p0 L b ½
¬ wy
w
wx ¼ ¬« a / 2 a 0 ¼» ° 2 L°
­re2 cos D ½ ° ° ­2b ½
° ° ° p0 L b ° p0 ° °
ª1 0 0 º Ÿ ^re ` ® re1cosD ¾ ® ¾ ®2b ¾
>E@ E ««0 1 0 »» ° r sin D ° ° 2 L° 4 ° °
¯ e1 ¿ ° p0 L a ° ¯a¿
¬«0 0 1/ 2 ¼» °
¯ 2 2L ¿
°
10 11
 EKSAMEN I EMNE TKT4191 Oppgave 3
ELEMENTMETODEN 1
x Nullenergimoder er deformasjonsmønster som kan opptre uten at
elementet/elementene opptar energi.
Høsten 2005
x Ved å bruke redusert integrasjon kan man få nullenergimoder. Man skiller mellom
Løsningsforslag nullenergimoder som kan bre seg utover hele nettet og nullenergimoder som ikke gjør
det.
x I en spenningsanalyse der timeglassmønster oppstår vil deformasjonen bli altfor stive.
Oppgave 1
Nullenergimoder for Q4:
1
N1 [ ([  1)n(n  1)
4
1 2
N5 ([  1)n(n  1)
2
N9 ([ 2  1)( n 2  1)
Oppgave 2
Ved å bruke 1×1 Gauss får man 2 ˜ 4  3  3 2 nullenergimoder.
Ved bruk av geometrisk avbildning kan komponentene i > B @ > E@> B @ J inneholde
T
x
Nullenergimode for Q8:
både polynomer i telleren og nevneren. Disse kan være vanskelig å integrere analytisk
og derfor vil det være hensiktsmessig å bruke numerisk integrasjon.
x Når elementet innehar sin basis form / det  J  konstant , består komponentene i
> B@ > E@> B@ J av polynomer. Gauss kvadratur blir brukt da denne er den mest
T
optimale kvadraturregelen for polynomer i 1 dimensjon. Den fungerer også bra i 2D
og 3D.
x Ved full integrasjon brukes en kvadraturregel slik at > k @ integreres eksakt for et
element som innehar sin basis form.
Ved å bruke 2×2 Gauss får man 2 ˜ 8  3  4 ˜ 3 1 nullenergimode.
Oppgave 4 Oppgave 5
u1 0, v1 0, u2 0, v2 0.001a, u3 0.001b, v3 0 Konsistente laster er laster som er fremkommet ved bruk av VFP. De er arbeidsekvivalente
y med de jevnt fordelte laster.
Ÿ u  x, y  ˜ 0.001b=0.001y
b
­ Px ½ ­x ½ ­a 2 ½
v  x, y 
x
˜ 0.001a=0.001x ^P` ® ¾ i punktet ® ¾ ® ¾
a ¯ Py ¿ ¯ y¿ ¯b 2 ¿
^re ` > N x @ ^P`
T
Finner tøyningene:
wu
Hx 0 ­ a 2 b 2½
wx °1  a

b °
wv ­ r1x ½ ° ° ­0 ½
Hy 0 ° ° ° a 2 ° ° °
wy ® r2x ¾ ® ¾ Px ®1 2 ¾ Px
° ° ° a ° °1 2 °
wu wv ¯ r3x ¿ ¯ ¿
J xy  2 ˜ 0.001 0.002 ° b2 °
wy wx
¯° b ¿°
Finner spenningene:
­ a 2 b 2½
°1   °
­ r1y ½ a b ­0 ½
Vx
E
1  2
 H x  H y  0 ° °
°
° a 2
°
° ° °
® r2y ¾ ® ¾ Py ®1 2 ¾ Py
° ° ° a ° °1 2 °
Vy
E
1  2
 H y  H x  0 ¯ r3y ¿ ° b2 ° ¯ ¿
E 70000 ¯° b ¿°
W xy GJ xy J xy ˜ 0.002 53.85MPa
2 1    2 1  0.3
Flyt?
Ve 3 ˜  53.85MPa 
2
93.27MPa<220MPa=f y
Ÿ Det vil ikke oppstå plastisk flyt.
Oppgave 6 Oppgave 7
x Fullstendighetskriteriet: Langs rendene 1-2-3-4 og 5-6-7-8 kan vi bestemme et tredjeordens polynom av y :
u v q1  q 2 y  q 3 y 2  q 4 y 3
Forskyvningsfunksjonene må være slik at stivlegemebevegelser kan beskrives uten at
og tilsvarende for v .
det oppstår spenninger og forskyvningsfunksjonene må være slik at konstante
tøyninger/spenninger kan beskrives innen elementene.
Langs rendene 1-8 og 4-5 kan vi bestemme et førsteordens polynom av x :
u v q5  q6 x
For C 0 innebærer dette at forskyvningsfunksjonene u and v må minst inneholde et
fullstendig polynom av orden 1. og tilsvarende for v .
x Kompatibilitetskriteriet for C 0 : u ( x, y ) kan da bestemmes som:
u ( x, y ) q 1  q 2 y  q3 y 2  q 4 y 3   q5  q6 x 
Til og med den 0’te deriverte av forskyvningsfunksjonene må være kontinuerlig
mellom elementene, dvs. forskyvningsfunksjonene selv. =q1q5  q1q 6 x  q 2q5 y  q 2q 6 xy  q 3q5 y 2
+q 3q 6 xy 2  q 4q 5 y 3  q 4q 6 xy 3
=a u1  a u2 x  a u3 y  a u4 xy  a u5 y 2
Kompatibilitet:
 a u6 xy 2  a u7 y 3  a u8 xy 3
y
H (tilsvarende for v )
I
G
Ÿ ¬« x ¼» ¬«1 x y xy y xy y xy ¼»
2 2 3 3
Skisserer N1 og N 2 :
D F 1
E
N1
Ikke oppfylt Ikke oppfylt
A x
B C
Q4: Lineært forskyvningsfelt langs rander – 2 node
Q8: Kvadratisk forskyvningsfelt langs rander – 3 node
1
N2
Tilsvarende:
CST: Lineært forskyvningsfelt langs rander – 2 node
LST: Kvadratisk forskyvningsfelt langs rander – 3 node
Oppgave 8 Oppgave 9
Randbetingelser: Tøyningsenergien til en samling av elementer:
nel
U ¦U
i 1
i
³ 2 >B@^d`  >E@>B@^d`  dV
1 T
Ui i i
Veli
1
=³ ^d`i >B@ >E@> B@^d`i dV
T T
Veli
2
1
^d`i > k @i ^d`i
T
2
Elementnett, Q8:
n el
1
¦ 2 ^d` >k @ ^d`
T
U i i i
i 1
n el
¦ 2 >L@ ^D` >k @ >L@ ^D`
1 T
i i i
i 1
n el
1
¦ 2 ^D` >L@ >k @ >L@ ^D`
T T
i i i
i 1
T§ ·
n el
1
^D` ¨ ¦ >L@i > k @i >L@i ¸ ^D`
T
U=
2 ©i1 ¹
n el
> K @ ¦ > L@i T >k @i > L@i
i 1
> k @i ³ >B@ >E@>B@dV
T
Veli
 Norges teknisk-naturvitenskapelige universitet NTNU
Fakultet for ingeniørvitenskap og teknologi
Institutt for konstruksjonsteknikk Side 1 av 15
Oppgave 10
s
3
r [2 EKSAMEN I TKT4191 ELEMENTMETODEN 1
s [3
1  r  s [1 Mandag 4. desember 2006
Kl 09.00 – 13.00
5
6
Oppgave 1 (vekt ca. 10 %)
a)
1 4 2 r
Beregningsmessige hovedpunkter i elementmetoden:
§1 · § 1·
N1 C1 1  r  s  ¨  r  s ¸ C1[1 ¨ [1  ¸ , C1 2
©2 ¹ © 2¹ 1. Diskretisering: Et kontinuerlig legeme med uendelig mange frihetsgrader
ŸN 1 [  2[  1
1 1 erstattes med en tilnærmet diskret modell med et endelig antall frihetsgrader.
N4 C4 r 1  r  s  C4[ 2[1 , C4 4 Legemet deles inn i elementer og knutepunkter (noder) med definerte
frihetsgrader.
Ÿ N4 4[1[ 2 2. Elementanalyse: Forskyvningsfeltet ^u` innen elementet tilnærmes ved
Belastning: ^u` > N @^d` , der er > N @ interpolasjonsfunksjoner og ^d` er frihetsgrader i
nodene. Bruk av kinematisk kompatibilitet, lineært elastisk materiallov og
­ Fx ½ virtuelle forskyvningers prinsipp (eller prinsippet om stasjonær potensiell
^F` ® ¾ energi) gir elementstivhetsrelasjonen >k @^d` ^re ` , hvor >k @ er elementets
¯ Fy ¿
stivhetsmatrise og ^re ` er konsistent lastvektor.
­ rx1 ½ ­ Fx ½ 3. Systemanalyse: Bruk av kompatibilitet i nodene (relasjon mellom elementets
® ¾
¯ ry1 ¿
³ [  2[
1 1  1 ® ¾ tdA
¯ Fy ¿
og systemets frihetsgrader) og virtuelle forskyvningers prinsipp (eller
prinsippet om stasjonær potensiell energi) for et system av elementer gir
§ A A · ­ Fx ½ systemets stivhetsrelasjon > K @^D` ^R` , hvor > K @ , ^D` og ^R` er
=¨2  ¸® ¾t
© 6 3 ¹ ¯ Fy ¿ henholdsvis stivhetsmatrisen, frihetsgradene og lastvektoren til systemet.
­0 ½ 4. Innføring av randkrav: Systemets stivhetsmatrise er singulær dersom
=® ¾ konstruksjonen ikke er fastholdt mot stivlegemebevegelser. Innføring av
¯0 ¿ tilstrekkelige randkrav gjør stivhetsmatrisen inverterbar, og stivhetsrelasjonen
kan løses for de ukjente forskyvningene i nodene.
­ rx4 ½ ­ Fx ½ 5. Løsning av ligningssystemet: Direkte eller iterative løsningsmetoder brukes for
® ¾
¯ ry4 ¿
³ 4[1[2 ¯® Fy ¿¾ tdA å løse det resulterende ligningssystemet: ^D` > K @ ^R`
1
4At ­ Fx ½ 6. Beregning av spenninger: Elementets frihetsgrader ^d` er entydig gitt av
= ® ¾
12 ¯ Fy ¿ systemets frihetsgrader ^D` , og tøyningene og spenningene innen et element
At ­ Fx ½ kan da beregnes ved
= ® ¾
3 ¯ Fy ¿
Norges teknisk-naturvitenskapelige universitet NTNU Norges teknisk-naturvitenskapelige universitet NTNU
Fakultet for ingeniørvitenskap og teknologi Fakultet for ingeniørvitenskap og teknologi
Institutt for konstruksjonsteknikk Side 2 av 15 Institutt for konstruksjonsteknikk Side 3 av 15
c)
b)
Systemstivhetsrelasjonen > K @^D` ^R` kan skrives på formen
Virtuelle forskyvningers prinsipp uttrykt i de to koordinatsystemene:
ª K11 K12 º ­D x ½ ­ R c ½
^G d` ^r` ^G dc` ^rc` ^G d` > T@ ^rc`
T T T T
« »® ¾ ® ¾
¬K 21 K 22 ¼ ¯ Dc ¿ ¯R x ¿
hvor det er brukt at
hvor ^D x ` = frie forskyvninger, mens ^Dc ` = foreskrevne forskyvninger.
^dc` > T@^d`
Det følger at:
Det følger da at
>K11 @^D x `  > K12 @^Dc ` ^R c `
^G d` ^r`  >T@ ^rc`
T T
0
>K 21 @^D x `  >K 22 @^Dc ` ^R x `
Dette må holde for enhver ^G d` , og relasjonen mellom nodelastene ^r` og ^rc` blir
Hvis ^Dc ` forhindrer stivlegemebevegelser, er > K11 @ inverterbar, og de ukjente forskyvning-
^r` > T@ ^rc`
T
ene er bestemt ved
>K11 @^Dx ` ^R c `  >K12 @^Dc `
Gitt at
eller
>k @^d`  ^r` and >k c@^dc`  ^rc`
^D x ` >K11 @ ^R c `  > K12 @^Dc `
1
Det følger da at
Reaksjonskreftene ^R x ` bestemmes ved bruk av den andre ligningen.
>k @^d`  ^r`  > T@ ^rc`
T
>T@ >k c@^dc` >T@ >k c@>T@^d`
T T
dvs.
>k @ >T@ >k c@>T@
T
Norges teknisk-naturvitenskapelige universitet NTNU Norges teknisk-naturvitenskapelige universitet NTNU
Fakultet for ingeniørvitenskap og teknologi Fakultet for ingeniørvitenskap og teknologi
Institutt for konstruksjonsteknikk Side 4 av 15 Institutt for konstruksjonsteknikk Side 5 av 15
Oppgave 3 (vekt ca. 10 %)
Oppgave 2 (vekt ca. 10 %)
a)
Node 3:
Statisk likevekt er oppfylt i nodene, men generelt ikke innen et element eller på randen mel-
N3 C3 y  y  b  lom to elementer.
1
N 3  0, 2b  1 Ÿ C3 2b  2b  b  1 Ÿ C3 Kinematisk kompatibilitet er oppfylt i nodene og innen et element, og vanligvis også på ran
2b 2
1
N3 y  y  b
2b 2 b)
Konvergenskriterier for C0-elementer:
Node 5:
1) Innen hvert element må den antatte forskyvningsfunksjonen og dens deriverte (tøyninge-
§ 2b · ne) være kontinuerlige
N 5 C5 y ¨ x  y  2b ¸
© a ¹
2) Over randen mellom to elementer må den antatte forskyvningsfunksjonen være kontinuer-
§a · § 2b a · 1
N 5 ¨ , b ¸ 1 Ÿ C5b ¨  b  2b ¸ 1 Ÿ C5 lig
©2 ¹ © a 2 ¹ 2b 2
1 § 2b · 3) Innen hvert element må den antatte forskyvningsfunksjonen inneholde et komplett poly-
N 5 y ¨ x  y  2b ¸ nom av orden 1. Dette betyr at elementet eksakt kan beskrive vilkårlige stivlegeme-
2b 2 © a ¹
bevegelser og konstante tøyningstilstander
c)
Node 2:
Forutsetter et komplett polynom av grad p . La h være elementlengden.
§ 2bx · § 2b ·
N 2 C2 ¨  y ¸ ¨ x  y  2b ¸
© a ¹© a ¹ Feilen i forskyvning er O  h p 1  , mens feilen i tøyning er O  h p  .
§ 2ba · § 2b · 1
N 2  a , 0  1 Ÿ C2 ¨  0 ¸ ¨ a  0  2b ¸ 1 Ÿ C2
© a ¹© a ¹ 8b 2 Konvergenshastigheten er p  1 for forskyvningene og p for tøyningene.
§
1 2bx · § 2b ·
N2 ¨  y ¸ ¨ x  y  2b ¸
8b 2 © a ¹© a ¹
Norges teknisk-naturvitenskapelige universitet NTNU Norges teknisk-naturvitenskapelige universitet NTNU
Fakultet for ingeniørvitenskap og teknologi Fakultet for ingeniørvitenskap og teknologi
Institutt for konstruksjonsteknikk Side 6 av 15 Institutt for konstruksjonsteknikk Side 7 av 15
Oppgave 4 (vekt ca. 10 %) y-retning:
wW yx wV y 2W0  V 0 2W0  V 0
Statiske randkrav langs vilkårlig rand:   0
wx wy a a
V l W m ) ; W lV m )
x xy x xy y y
Konklusjon: Likevekt er oppfylt innen skiven.
For vertikal rand, x = 0:
y y
l 1, m 0 :  V x )x V 0 ;  W xy )y W0
a a For vertikal rand, x = 0:
For skrå rand, x = y: y y y
Vx V 0 , Vy (2W 0  V 0 ) , W xy W 0
a a a
2 2
l ,m  :
2 2
For skrå rand, x = y:
2 2
Vx  W xy )x 0 Ÿ Vx W xy
2 2 x x x
2 2 Vx W 0  V 0  , V y (W 0  V 0 ) , W xy (W 0  V 0 )
W xy V y )y 0 Ÿ W xy Vy a a a
2 2 Ÿ Vx W xy Vy
Likevektsligningene for skiven:
Konklusjon: Likevekt er oppfylt på randen.
wV x wW yx wW xy wV y
 0,  0
wx wy wx wy
x-retning:
wV x wW yx W0 W 0
  0
wx wy a a
Norges teknisk-naturvitenskapelige universitet NTNU Norges teknisk-naturvitenskapelige universitet NTNU
Fakultet for ingeniørvitenskap og teknologi Fakultet for ingeniørvitenskap og teknologi
Institutt for konstruksjonsteknikk Side 8 av 15 Institutt for konstruksjonsteknikk Side 9 av 15
Oppgave 5 (vekt ca. 10 %)
Generelt uttrykk for konsistent lastvektor for konsentrert last:
Aktuelle formfunksjoner bestemmes ved null-linjemetoden: ­Fx ½
^re ` > N @
T
® ¾
1 1 ¯Fy ¿
N4 C 4 x  x  a  , C4 Ÿ N4 x  x  a
2a 2 2a 2
Konsistent lastvektor for aktuell konsentrert last:
N7 C7  a  x  a  x  , C7
1
a2
Ÿ N7
1 2
a2
 a  x2 
­ry4 ½ ­1½
1 1 ° ° P
ª x  x  a  2  a 2  x 2  x  x  a º
T P° °
N3 C3 x  a  x  , C3 Ÿ N3 x  x  a ®ry7 ¾ ®6¾
2a 2 2a 2 ° ° 2a 2 ¬ ¼x a/2 8° °
¯ ry3 ¿ ¯3¿
Generelt uttrykk for konsistent lastvektor for flatelast: Alle de andre elementene i ^re ` er lik null.
^re ` ³ > N @ ^` dS
T
Samlet konsistent lastvektor (kun ledd ulik null):
Konsistent lastvektor for aktuell flatelast:
­ry4 ½ ­1 ½ ­1½
­ry4 ½ ° ° qa ° ° P ° °
a ­1 ½ ®ry7 ¾ ®4 ¾  ® 6 ¾
° ° 1
ª x  x  a  2  a 2  x 2  x  x  a  º qdx
T qa ° ° ° ° 3 ° ° 8° °
®ry7 ¾ ³ 2a 2 ¬ ¼ ®4 ¾
3 ° ° ¯ ry3 ¿ ¯1 ¿ ¯3¿
° ° a
¯ ry3 ¿ ¯1 ¿
Alle de andre elementene i ^re ` er lik null.
Norges teknisk-naturvitenskapelige universitet NTNU Norges teknisk-naturvitenskapelige universitet NTNU
Fakultet for ingeniørvitenskap og teknologi Fakultet for ingeniørvitenskap og teknologi
Institutt for konstruksjonsteknikk Side 10 av 15 Institutt for konstruksjonsteknikk Side 11 av 15
Oppgave 6 (vekt ca. 10 %) T ­ Fx ½ ­P½ § a ·
:  ^u ` ® ¾  ^uC vC ` ® ¾ uC P  ¨ a1  a2b ¸ P
¯ Fy ¿ ¯0 ¿ © 2 ¹
Gitt forskyvningsfeltet
u a1 x  a2 y, v a3 y Potensiell energi:
abtE § 2 1 2 2· § a ·
hvor a1 , a2 og a3 er konstanter. 3 U : ¨ a1  a2  a3 ¸  ¨ a1  a2b ¸ P
4 © 2 ¹ © 2 ¹
Tøyningene er:
Stasjonær verdi av 3 :
wu wv wu wv
H a, H a, J  a w3 abtE a P
x
wx
1 y
wy
3 xy
wy wx
2
0 Ÿ a1  P 0 Ÿ a1
wa1 2 2 btE
w3 abtE 4P
0 Ÿ a2  bP 0 Ÿ a2
Tøyningsenergien: wa2 4 atE
w3 abtE
1 0 Ÿ a3  0 0 Ÿ a3 0
³ 2 ^` >E@^` dV wa3
T
U 2
Siden tverrkontraksjonen Q er lik null, har vi at Forskyvningsfeltet for skiven:
ªE 0 0 º P P
« x
0 »»
u y, v 0
> E@ «0 E btE atE
¬« 0 0 E / 2 ¼»
Tøyningsenergien for skiven: I punkt C, x a / 2, y b:
ªE 0 0 º ­ a1 ½ Pa Pb
uC  , vC 0
° ° abt
a ` «« 0 0 »» ® a3 ¾
1 1
³ 2 ^` >E@^` dV ^a1
U
T
a E 2btE atE
3 2
2 2
¬« 0 0 E / 2 ¼» ¯°a2 ¿°
abtE § 2 1 2 2·
¨ a1  a2  a3 ¸
4 © 2 ¹
Lastpotensialet:
 ³ ^u` ^` dS
T
:
For punktlasten P i C gir dette:
Norges teknisk-naturvitenskapelige universitet NTNU Norges teknisk-naturvitenskapelige universitet NTNU
Fakultet for ingeniørvitenskap og teknologi Fakultet for ingeniørvitenskap og teknologi
Institutt for konstruksjonsteknikk Side 12 av 15 Institutt for konstruksjonsteknikk Side 13 av 15
Oppgave 7 (vekt ca. 10 %) Oppgave 8 (vekt ca. 10 %)
a)
Patchtesten er en numerisk test som utføres for en gruppe av elementer utsatt for et lineært
forskyvningsfelt, og som benyttes til å undersøke om et element er gyldig med hensyn til kon- For en funksjon M M [ ,K  har vi at
vergens mot korrekt løsning ved nettforfining.
wM wM wx wM wy
Prosedyre: 
w[ wx w[ wy w[
1. Stivhetsmatrisen > K @ dannes for en samling av elementer med vilkårlig form som omslut- wM wM wx wM wy

ter en node i. wK wx wK wy wK
2. Et sett av nodeforskyvninger ^D` beregnes for de ytre nodene som tilsvarer en vilkårlig I
eller
konstant tøyningstilstand for elementene.
3. Tøyningene i elementene ^` og nodekreftene ^R` beregnes på vanlig måte.
4. Testen er tilfredsstilt hvis og bare hvis ­ wM ½ ­ wM ½
a. Tøyningene ^` i elementene er konstante ° °
° w[ ° ° wx ° ª x,[ y,[ º ª ¦ N i ,[ xi ¦ N [ y º»
® ¾ > J @ ®° wM ¾° , >J@ « y,K ¼»
«
i, i
b. Indre nodekrefter ^R ` er lik null, dvs. node i er i statisk likevekt
i ° wM ° ° ° ¬ x,K ¬ ¦ N i ,K xi ¦N Ky ¼ i, i
¯° wK ¿° ¯° wy ¿°
b)
Normaltøyningene:
Hx u, x *11u,[ *12u,K *11 ¦ Ni ,[ ui  *12u,K ¦ Ni ,K ui
Hy v, y * 21v,[ * 22v,K * 21 ¦ Ni ,[ vi  * 22 ¦ Ni ,K vi
Elementer som oppfyller Patch-testen tilfredsstiller konvergenskravene fordi:
x Punkt a) er en kontroll på at fullstendighetskriteriet er tilfredsstilt
x Punkt b) er en kontroll på at integrert likevekt er oppfylt på randen mellom elemente-
ne (dvs. at integralene i virtuelle forskyvningers prinsipp konvergerer mot korrekt
verdi) selv om kompatibilitetskriteriet ikke er oppfylt.
Norges teknisk-naturvitenskapelige universitet NTNU Norges teknisk-naturvitenskapelige universitet NTNU
Fakultet for ingeniørvitenskap og teknologi Fakultet for ingeniørvitenskap og teknologi
Institutt for konstruksjonsteknikk Side 14 av 15 Institutt for konstruksjonsteknikk Side 15 av 15
Oppgave 9 (vekt ca. 10 %) Oppgave 10 (vekt ca. 10 %)
For et Q8-element inneholder formfunksjonene følgende polynomledd (fra Pascals trekant):
y
1, [ , K , [ 2 , [K , K 2 , [ 2K , [K 2
Stivhetsmatrisen >k @ er definert som
M M
1 1
>k @ ³ >@ >E@> @ dV ³ ³ >@ > E@> @ tJd[ dK
T T
1 1
x
Produktet >  @ > E@>  @ inneholder polynomer i 4. grad i [ og K , mens J er konstant.
T
Tykkelsen inneholder polynomer i 2. grad i [ og K , slik at integranden >  @ > E@>  @ tJ inne-
T
holder polynomer av 6. grad i [ og K .
For Q8-elementet representerer en 3 ganger 3 Gauss-regel full integrasjon. Elementet har kun
ikke-kommuniserbare timeglassmoder for en 2 ganger 2 Gauss-regel, som ofte gir bedre re-
En n u n Gauss-regel integrere et polynom av grad n  1 eksakt. For det aktuelle elementet
sultater. Dette taler for å velge en 2 ganger 2 Gauss-regel for dette tilfellet.
trengs det derfor en 4 u 4 Gauss-regel for eksakt integrasjon av >k @ .
Q9-elementet kan beskrive ren bøyning selv når de ikke er rektangulære. Dette er ikke tilfelle
for Q8. Derfor bør Q9 velges for det aktuelle tilfellet.