- With Finite Element Method 1 (125.8 KiB) Viewed 20 times
L/3 L/3 L/3 10h TH(x) P ht. Xu с H(x=L/6) H(x=3L/6) u2 313 H(x=5L/6) 4 u4 e 1 e2 e3 - - 1
with finite element method
-
answerhappygod
- Site Admin
- Posts: 899604
- Joined: Mon Aug 02, 2021 8:13 am
with finite element method
with finite element method
Problem 3: Write a finite element program in Matlab that can be used to compute the displacement of the beam at point B in Figure 1. The program should be capable of using an arbitrary number of elements (m elements and n nodes), see Figure 2. The change in cross-section area (height) for each element must be included (element stiffness matrix, k). Make a convergence study with respect to the displacement in point B versus number of nodes/elements. Present the result so that it is easy for the reader to see how it converges. How many elements are required for a converged result? Supply the Matlab code in the report but also upload a separate m-file on the Canvas system. Hint: a for-loop over the elements is probably the easiest way (quick-and-dirty) to solve the updating of the cross-section area and assembly of the global stiffness matrix (although it's not a matlabish way). 1 2 3 m O 1 2 3 4 n-1 n L Figure 2. Illustration of nodes and elements.
L/3 L/3 L/3 10h TH(x) P ht. Xu с H(x=L/6) H(x=3L/6) u2 313 H(x=5L/6) 4 u4 e 1 e2 e3 - - 1
Problem 3: Write a finite element program in Matlab that can be used to compute the displacement of the beam at point B in Figure 1. The program should be capable of using an arbitrary number of elements (m elements and n nodes), see Figure 2. The change in cross-section area (height) for each element must be included (element stiffness matrix, k). Make a convergence study with respect to the displacement in point B versus number of nodes/elements. Present the result so that it is easy for the reader to see how it converges. How many elements are required for a converged result? Supply the Matlab code in the report but also upload a separate m-file on the Canvas system. Hint: a for-loop over the elements is probably the easiest way (quick-and-dirty) to solve the updating of the cross-section area and assembly of the global stiffness matrix (although it's not a matlabish way). 1 2 3 m O 1 2 3 4 n-1 n L Figure 2. Illustration of nodes and elements.
L/3 L/3 L/3 10h TH(x) P ht. Xu с H(x=L/6) H(x=3L/6) u2 313 H(x=5L/6) 4 u4 e 1 e2 e3 - - 1
Join a community of subject matter experts. Register for FREE to view solutions, replies, and use search function. Request answer by replying!