Bentley Communities
Bentley Communities
  • Site
  • User
  • Site
  • Search
  • User
AutoPIPE
  • Product Communities
AutoPIPE
AutoPIPE Wiki Vortex Shedding
    • Sign In
    • -Pipe and Vessel Stress Analysis - Wiki
      • +.Bentley AutoPIPE
      • -.Bentley AutoPIPE Vessel (powered by Microprotol)
        • +Features and Bugs
        • +General Information
        • -Support Solutions and FAQ's
          • +Air Coolers
          • Boilers
          • +Calculations
          • +Coils
          • +Cost Estimation
          • +Drawings
          • +Expansion Joints
          • +Import / Export
          • +Interface
          • +Lifting Lugs
          • +Loads and Load Combinations
          • +Material, Material Database, and Default Libraries
          • +Nozzles and Flanges
          • +Pressure Vessels
          • +Report
          • +Shell and Tube Heat Exchangers
          • +Shells
          • +Soil
          • +Tanks
          • +Units
          • -Vortex Shedding
            • 01. How to apply Kanti Mahajan theory of Vortex Shedding
            • 02. How to apply Henry Bednar theory of Vortex Shedding
          • +Weld Joints
        • +System Messages
        • +Video Gallery
      • +ADL PIPE
      • +Bentley AutoPIPE Nozzle
      • +Bentley PlantFLOW
      • +Bentley PULS
    • +AutoPIPE Brand Learning Resources
    • +Attend a live meeting with an AutoPIPE Analyst

    You are currently reviewing an older revision of this page.

    • History View current version

    Vortex Shedding

    Wind induced vibrations on tall vessels are a concern and AutoPIPE Vessel currentley has two methods to address this design consideration.

    Method 1 - Kanti Mahajan

    1. Obtain first the natural frequency ‘f’ of your vessel

     

    1. Check for vibration possibility

     

    If W/LDr2 <=20, issue message in report “Vibration analysis necessary to be performed as W/LDr2 <=20”

    vibration analysis must be performed and proceed to step 3

     

    where W= corroded wt. in lb, L=length in ft., Dr=average internal diameter of top half of the structure in ft.

     

    If W/LDr2 >20, issue message in report “Vibration analysis not necessary to be performed as W/LDr2 >20”

     

    1. Perform stability investigation,

     

    Calculate damping factor DF=Wδ/LDr2

     

    Default value of δ = 0.03 (Table 3)

     

    If DF <= 0.75, issue message in report “system is unstable and vortex shedding vibration analysis is required as DF <= 0.75”, proceed to step 4

     

    If DF > 0.75, issue message in report “Vortex shedding vibration analysis not necessary to be performed as DF > 0.75”

     

    1. Calculate critical wind velocity being where S = 0.226 (Strouhal number)

     

    Justification: In the book this formula appears as   in mph, but Dr is introduced in feet, f in Hz, so the number of Strouhal is affected by a conversion factor from feet to miles and hours to seconds. To reverse this change and work in consistent units we have to divide by 1.466 (1mph = 5280/3600 = 1.466 ft/s). And to place S in the denominator, we inverse the coefficient:

     

    1. Calculate maximum wind velocity at top of the structure (H = 10m or 30ft)

    Vw =Vb(L/H)0.143

    Where Vb = wind velocity at H

     

    1. Apply gust factor, 1.3 to obtain maximum gust velocity = 1.3Vw

     

    1. If Vc <=maximum gust velocity, perform amplitude calculations go to step 8.

     

    If Vc > maximum gust velocity, issue message in report “Vibration amplitude no need to be calculated as Vc > maximum gust velocity”.

     

    1. Calculate amplitude of maximum dynamic deflection, Z=L5VC2/W δ Dr (10)-6 (0.00243) in.

     

    1. Check if Z < Maximum Deflection input by user, issue message in report “structure is safe:amplitude of maximum dynamic deflection < maximum allowable deflection”

     

    If Z> Maximum Deflection input by user, issue message in report “structure is unsafe:amplitude of maximum dynamic deflection exceeds maximum allowable deflection, please modify design”

     

    Reference: Kanti Mahajan, PE. (1979). Design of Process Equipment. Tulsa, OK: Pressure Vessel Handbook Publishing.

     

    Method 2: Bednar

    1. Check vibration criteria. See points 1 and 2 in Kanti Mahajan.

     

    1. Calculate maximum wind velocity at top of the structure. See points 5 and 6 in Kanti Mahajan.

     

    1. Calculate first and second critical wind velocity, V1 and V2

     

     where S = 0.2 (Strouhal number)

     

    V2 = 6.25.

     

    If the design wind velocity V>V1, perform the detailed vibration analysis from step 3 below with V1

    If the design wind velocity V>V2, perform the detailed vibration analysis from step 3 below with V2

    If V1 or V2 > 60 mph/27 m/s, ignore the calculation for that V value as there is cutoff for maximum critical wind speed for design consideration (page 120).

     

    1. Consider default value of lift coefficient, CL=0.5, normally lift coefficient is 0.4-0.6 (page 111) for practical engineering. Allow users to input CL as optional to default.

     

    1. Create table 4.3 (page 116) in 2 dropdown selections for evaluation of logarithmic decrement, δ and Magnification Factor, M.F. Allow 2 dropdowns selections

     

    • Selection of soil types (Soft – I, Stiff – II, Rock, Very Stiff – III)
    • Equipment type (Tall process columns, unlined stacks, lined stacks)

     

    Determine δ and M.F. from the above selected combinations

     

    1. Determine equivalent force F=0.00086 (CL x M.F.)(d x H x V12) (replace V1 by V2 as needed)

    Where, d = outside projected dia., ft., H= total column height, ft.

     

    1. This force F acts on the top of the column so calculate the moments at each level in same way as static wind load analysis.

    Note in the example only the moment at the junction of column and skirt is considered but APV has different approach of calculating at each individual level

     

    1. Perform fatigue analysis at the bottom tangent line weld joint

     

    • Evaluate the cyclical stress range SR=2S, where S is the bending stress due to F acting at the top. We’ll be taking the maximum of |σM| , |σz| , |σeq|  .The calculation of SR in the example is confusing and we need to follow the standard practice of calculating the bending stress as implemented in the program S= σ = M/section modulus, section modulus = π (r24-r14)/4r2
    • calculate number of cycles to failure, N = (K/βSR)n (page 118)
    • where, n=5, K=780,000 (these values are for carbon steel, no other value specified for any other grade or type of steel) 

    stress intensification factor, β as follows

    β = 1.2, shell plate with smooth finish,

    β = 1.8, butt weld

    β = 3.0 fillet weld

    • Allow users to select value for β from dropdown, default β = 2.0
    • Allow user to specify safety factor, S.F. Default S.F.=20.
    • Safe vibration time or service life considering weld fatigue due to wind vortex induced vibration = (N/S.F.)x T/3600 x(1/24)x(1/365) years
    • T=time period (seconds per cycle) = 1/f

    Reference: Henry H. Bednar, PE. (1991). Pressure Vessel Design Handbook Second Edition. Malabar, FL: Krieger Publishing Company.

    Communities
    • Home
    • Getting Started
    • Community Central
    • Products
    • Support
    • Secure File Upload
    • Feedback
    Support and Services
    • Home
    • Product Support
    • Downloads
    • Subscription Services Portal
    Training and Learning
    • Home
    • About Bentley Institute
    • My Learning History
    • Reference Books
    Social Media
    •    LinkedIn
    •    Facebook
    •    Twitter
    •    YouTube
    •    RSS Feed
    •    Email

    © 2023 Bentley Systems, Incorporated  |  Contact Us  |  Privacy |  Terms of Use  |  Cookies