Computational Solid Mechanics (CED501) – Mid-semester examination
Date and time: 16.09.2023 and 02:30 PM- 04:30 PM
Total points: 32
1. Consider the differential equations:
d dv
− a ( x) u + = f ( x)
dx dx
d du dv
− b( x ) + a ( x ) u + = g ( x )
dx dx dx
du du dv
with 0 x 1. Boundary conditions are: u(0) = u(1) = 0 , =0, =1, =0,
dx x =0 dx x =1 dx x =0
dv
=1.
dx x =1
(i) Considering two test functions p( x) and q( x) , write the weak form.
(ii) Considering a single element, approximate u , v , p , and q as: u = (1 − x)u1 + xu2 ,
v = (1 − x)v1 + xv2 , p = (1 − x) p1 + xp2 , and q = (1 − x)q1 + xq2 . Consider
a( x) = b( x) = f ( x) = g ( x) = 1 , determine the 4x4 stiffness matrix from the bilinear form.
(iii) Assume u1 = v1 = 0 , determine u 2 and v2 . [14]
2. Consider the functional
1 L
2
du 1 L
I (u ) =
2 0
a ( x )
dx
+ cu dx + [u ( L)]2 − f ( x)udx − u ( L)u − Q0u ( L)
2 0
c , , u , and Q0 are constants. Determine the Euler-Lagrange equations from the variational principle
I = 0 and the natural boundary conditions. [8]
P
3. Consider the equation v + v2 − h v = − . The symbol denotes SPH
t x h x
NP NP
averaging operator. Let v = mJ vJ ( x − xJ , h) , v 2 = mJ vJ2 ( x − xJ , h) , and
J =1 J =1
P NP
1 P
= mJ ( x − xJ , h) . Assume the smoothing length h is a function of time and
x J =1 J x J
substituting the formulas given above in the equation, using chain or product rule as applicable, deduce the
simplest form of the discretized version of the governing equation. [6]
4. Identity tensor I = ei ei . Consider a vector v = vk ek . Given e j e k = mjk e m . Show that, I v = v I
[4]
