Scribd
Upload
6 views1 page

Midsem

Uploaded by

degeca7075
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
6 views1 page

Midsem

Uploaded by

degeca7075
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 1

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]

You might also like