Dynamic analysis of moving mass by Newmark beta method
14 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
Shabnam Sulthana Mohamed Isaque
el 4 de Sept. de 2019
Respondida: Samsun Anadolu
el 4 de Dic. de 2020
Hi,
I am new to matlab and have couple of doubts regarding my coding for the dynamic analysis of a beam of moving mass.
I am modelling a beam of 11 nodes and moving mass with a velocity of 1m/s. i wouldlike to plot the displacement of node5 when the load is in each nodes . for example displacement of node5 when load is in node 1,displacement of node5 load is in node 2,displacement of node5 when load is in node3..........displacement of node5 when load is in ith node.
For each time step, i should get the output of displacement of one node. i have attached the code which i have done,i am getting displacement as incremental for each nodes and the displacement is not changing for the before nodes.
could anyone help me in this issue.
Thanks
3 comentarios
Jim Riggs
el 5 de Sept. de 2019
Editada: Jim Riggs
el 5 de Sept. de 2019
A few things need clarification:
1) Number of elements is defined 2 different times;
ne1 = 10
numberElements = size(elementNodes,1)
2) elementNodes is defined more than once:
elementNodes = [...] % hard coded
elementNodes(:,1) = ii % computed
elementNodes(:,2)=ii+1
3) P is defined at the top as a scalar, then as an array, then re-defined as a scalar:
P = 10000
....
P = zeros(n,1)
P = 10000
Is P supposed to be a vector or a scalar?
4) You define a set of six constants (a1, a2,...a6) but then don't use them where you could to simplify the equations; e.g.
kbar = ks + C*a1 + M*a2;
A = M*a3 + C*a4;
B = M*a5 + C^2*a6;
...and inside the for loop:
veldv = Dv*a1 - vel(i)*a4 - accl(i)*a6;
accldv = DC*a2 - vel(i)*a3 - accl(i)*a5;
5) The calculation of A overwrites the initial definition of A
A = 0.00767; % area in m^2
...
A = M*a3 + C*a4;
So, appart from this confusion in the code, I still do not understand the question. Perhaps a drawing would help describe the relationship of the moving mass and the beam.
Respuestas (1)
Samsun Anadolu
el 4 de Dic. de 2020
Good Job, thanks for sharing.
I am a beginner numerical analysis problem.
I have checked similar codes; I couldn't find any clue about the boundary conditions of the system.
Do these systems have a boundary condition? Could you please explain it?
0 comentarios
Ver también
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!