File Exchange

image thumbnail

Flutter velocity of a suspension bridge

version 1.1 (67.6 KB) by E. Cheynet
The coupled flutter velocity of a single-span suspension bridge is computed in the frequency domain


Updated 03 Mar 2019

View License

The critical flutter velocity Vcr of a suspension bridge is estimated using a simple computational model accounting for the lateral, vertical and torsional motion of the bridge deck and using a multimodal approach. The computation is conducted in the frequency domain using the method proposed in [1]. The computed value of Vcr is compared to the famous analytical expressions from Selberg [2] and Rocard [3].
The present submission contains:

- The function fluterFd, which computes the critical flutter velocity following [1]

- The function VcrFlutter, which computes the critical flutter velocity following [2,3]

- An example file Example.m

- Two .mat file modalParameters_case1.mat and modalParameters_case2.mat that are used to load the eigen-frequencies and mode shapes of the two bridge models investigated

This is the first version of the submission. Some bugs may still exist. Any question, comments or suggestion is warmly welcomed.


[1] Jain, A., Jones, N. P., & Scanlan, R. H. (1996). Coupled aeroelastic and aerodynamic response analysis of long-span bridges. Journal of Wind Engineering and Industrial Aerodynamics, 60, 69-80.

[2] Selberg, A., & Hansen, E. H. (1966). Aerodynamic stability and related aspects of suspension bridges.

[3] Rocard, Y. (1963). Instabilite des ponts suspendus dans le vent-experiences sur modele reduit. Nat. Phys. Lab. Paper, 10.

Cite As

E. Cheynet (2019). Flutter velocity of a suspension bridge (, MATLAB Central File Exchange. Retrieved .

Comments and Ratings (1)


Perfect code, i have 3 questions .(1) Quasi-static assumption is used . so its no all the same in the references. every mode have three freedom, Every freedom need normalized?
(2)What does this code mean,why (jj,jj) is select.
dummyF = zeros(1,newN);
for jj=1:newN,
[dummyF(1,jj),~] = eigenFreqSyst(Mtot(jj,jj),Ktot(jj,jj),Ctot(jj,jj));
(3) In the reference ,is every mode is solved in one equation by generalized coordinates? IN this code every mode is soved in three degree of freedom ,That is a matrix 3(dof)*6(mode).
thank you very much



Added the project website

MATLAB Release Compatibility
Created with R2017b
Compatible with any release
Platform Compatibility
Windows macOS Linux

Discover Live Editor

Create scripts with code, output, and formatted text in a single executable document.

Learn About Live Editor