File Exchange

image thumbnail

Dynamic response of a line-like structure to a random load

version 3.1 (404 KB) by E. Cheynet
The dynamic displacement response of a line-like structure to an uncorrelated Gaussian white noise input is computed in the time domain.


Updated 29 Jul 2020

View Version History

GitHub view license on GitHub

The submission contains:

-the function dynaResp_TD that computes the time history of the displacement response of a line-like structure to a given load.
- The function eigenModes.m used here to compute the mode shapes and eigenfrequency of a cantilever beam
- The file bridgeModalProperties.mat that loads the mode-shapes and eigenfrequency of a single span suspension bridge.
- 2 example files Example1.mlx and Example2.mlx

To keep the analysis as simple as possible, the structure has only one type of motion. No modal coupling is introduced and no added mass, stiffness or damping is included.

Any comments, suggestion or question is welcomed.

Cite As

Cheynet, E. Time-Domain Dynamic Response of a Line-like Structure to a Random Load. Zenodo, 2020, doi:10.5281/ZENODO.3965470.

View more styles

Comments and Ratings (2)

E. Cheynet

Dear Maede,
Thank you for the feedback! The corresponding example file has been updated so that the model properties of the bridge are directly imported, without needing to call for the functions LysefjordBridge and eigenBridge. For your other comment, Ncoeff is the number of the coefficient for the from which are computed the mode shapes with the Galerkin method. These are neither the number of mode shapes nor the number of nodes, and this is why I use the variable name "Ncoeff". Ideally, Ncoeff should be much larger than the number of the highest mode. if Nyy is equal to Ncoeff in other submissions, it may be a coincidence. The absence of "t" in the function dYdT is indeed unnecessary (for the present case) and I have thus simplified the function. Note that in the updated submission, zetaStruct (the structural modal damping ratios) is also defined in a simpler way for the sake of clarity.

Maede Zolanvari

I found some bugs. First, I guess you forgot to include the LysefjordBridge and eigenBridge. However, they can be found easily through your other submissions, so no big deal :) Second, I guess the input of the eigenBridge(geometry,Ncoeff ), must be Nnodes (which is equal to Nyy in other submissions), rather than Ncoef. I don't think we need any Ncoef parameter here.
Another thing I wanted to mention is that in "dYdT" funstion, input "t" is not being used at all, not even in "Fun" or calculating the "k_i's," then why adding it as an input?
In Lysefjord.m, I think it's better to add "zetaStruct" as "Bridge.zetaStruct."
That's all I found for now. Keep the good work :)

MATLAB Release Compatibility
Created with R2019b
Compatible with R2012b and later releases
Platform Compatibility
Windows macOS Linux

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!