Question about Example "Dynamic Analysis of Clamped Beam"-Example from PDE Toolbox
2 views (last 30 days)
Show older comments
c = [2*G + mu; 0; G; 0; G; mu; 0; G; 0; 2*G + mu];
are being used in the linear case. I would like to understand how these coefficients correspond to the underlying equations, hoping to gain some insight in how the vector c must be specified. Because the problem is 2D (N=2) and there are 10 entries in vector c, i suppose the 2N(2N+1)/2-element form from specifyCoefficients is being used? But I don't understand where the entries in c are coming from. I would appreciate, if someone could explain (roughly) why the vector c looks this way.
0 Comments
Answers (1)
Alexis He
on 14 Apr 2023
The Lagrangian formulation of the equilibrium equations in index notation is :
For linear solution, we ignore the non-linear strain term in the Green-Lagrange strain tensor E and the first Piola-Kirchoff stress tensor P=F⋅S=F⋅C⋅E :
Thus, (1) will become :
For this 2-D problem, we have :
with using the constitutive relations in matrix form, we have :
replacing the with we have :
P.S. : means the c Coefficient of term in i-th equation
The 2N(2N+1)/2-Element Column Vector representation of c and the matrix representation of c coefficient ( with a “flattening” of the N-by-N-by-2-by-2 tensor to a 2N-by-2N matrix when N=2 ) are :
with using the material properties :
We obtained
0 Comments
See Also
Categories
Find more on General PDEs in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!