橋梁抗風(fēng)性能仿真解讀

1、圖1.1位重型支援橋某工作狀態(tài)模型利用CFD-ACE1行橋梁抗風(fēng)性能仿真一、計(jì)算實(shí)例簡(jiǎn)述本計(jì)算實(shí)例為重型支援橋假設(shè)時(shí)三個(gè)階段的風(fēng)橋耦合仿真。其模型結(jié)構(gòu)如圖所示：eyi前期使用GEO瞰件進(jìn)行模型建立，網(wǎng)格劃分等工作，并把此計(jì)算案例分為三個(gè)域進(jìn)行計(jì)算：1號(hào)域?yàn)闃蝮w本身；2號(hào)域?yàn)闃蝮w附近的空氣流域；3號(hào)域?yàn)闃蝮w遠(yuǎn)方的空氣流域。為了提高計(jì)算效率以及計(jì)算精確度，則需要在1號(hào)和2號(hào)域中網(wǎng)格相比3圖1.2域1和域2示意圖號(hào)域中更為精密。下圖為三個(gè)域的示意圖。其中如圖1.2所示，實(shí)體橋結(jié)構(gòu)內(nèi)部為1號(hào)域，橋體周圍的立方六面體（紫色六面體）到橋表面之間的空間為2號(hào)域。圖1.33號(hào)域示意圖3號(hào)域?yàn)樽畲蟮牧⒎搅骟w到

2、中間的立方六面體之間的空間、操作界面介紹打開CFD-ACE其工作界面如下圖所示:尸Gun口7豆|AAKpCUMri I Ujlw4 m IEVE3 nF - N Hs4 * 串 * * IM F . 0“ + M + 4 叮 g | g | gu | DC | 父 | m |MwiuiwwP n* Hwi rnmlvr iIImhI設(shè)JH面扳F rfcfWU! rruiM-r ChMiFA!的貨 IOwiiF Umk 51Ml iScilri r XMcniRMj r r IAmhi 3infchiMi F rrM SbrhrCH r FM-nfald nj-dli iN SuhsP 41 r

3、id De Icfwmm infita rrr I r fz i Elvaric iElKirl r 3gmk 隍qu r iO*ris r sEDwt.ffWT-niil Bevern4.T! Saw9DcL2加皿2川I則-IrivortHKkGumpChiq1iiZJFt-jU-hJAr#=1)mhHsgHNiapi,_di上J圖2CFD-ACE軟件工作界面其主要分為菜單區(qū)、實(shí)體（域）顯示區(qū)域、設(shè)置面板、操作對(duì)象選擇區(qū)等四個(gè)部分1 .菜單區(qū)File、Edit、View、Units、Tool、Windows分別為基本操作的操作菜單。度q圖mI為文件打開、存儲(chǔ)、撤銷快捷按鈕。力：為控制模型顯示

4、方式的選項(xiàng)按鈕，其中可選擇模型的可視角度、透明度、陰影等。|2q0日|”四對(duì)為cfd-ACE計(jì)算過(guò)程中對(duì)計(jì)算仿真過(guò)程或者結(jié)果操作的按鈕。其功能有時(shí)時(shí)觀察收斂曲線、打開與其關(guān)圖2.1.1菜單區(qū)工作界面2 .實(shí)體（域）顯示區(qū)域圖2.2.1實(shí)體（域）顯示區(qū)域其操作方式主要為鼠標(biāo)操作，其中鼠標(biāo)左鍵功能為拖動(dòng)，鼠標(biāo)右鍵功能為旋轉(zhuǎn)，滾輪功能為放大。同時(shí)可以通過(guò)此區(qū)域進(jìn)行操作對(duì)象的選擇3 .操作對(duì)象區(qū)域此區(qū)域?yàn)椴僮鲗?duì)象選擇區(qū)域，方便對(duì)其進(jìn)行選擇，選擇后，被選部分會(huì)在實(shí)體顯示區(qū)域顯示。然后進(jìn)行相關(guān)操作GenwtryTree冥Bjnda*yNameBCTyiwKSubType小串nkedCeetHZoneKey

5、口All口FavoritH0Block套Fluid9&didBC-c國(guó)，口GenJ.nterfareGe.l &e心1gliNo-SlipGeM &ct_3胃MlNo-SlipGen_Jginnle選才Opopen-打開相應(yīng)文件，或者直接點(diǎn)擊料UHJ-ACLULC;PLAbndge.4gongkuang-12UUIiuti-5.4.DIb-11 j|e Edit ViEiw Units To&li Window圖3.1.1 打開open圖3.1.2選擇打開文件 e :f d舉亞母父寶 。, Li R3, 蟲 司內(nèi)當(dāng)2、選擇設(shè)置面板中的PT選項(xiàng)，其名果如圖3.2.1O圖3.2.1 PT設(shè)置版面P

6、T(problemtype)為問(wèn)題類型設(shè)置面板，跟據(jù)特定的仿真實(shí)例選擇相應(yīng)的計(jì)算模塊，進(jìn)而求解。其中，F(xiàn)low為流體動(dòng)力學(xué)模塊；HeatTransfer為傳熱輻射模塊；Turbulence為湍流模塊；Chemistry/Mixing(chem)為化學(xué)模塊；Radiation(Rad)為輻射模塊；Spray為噴霧/顆粒模塊；FreeSurfaces(VOF為自由液面模塊；Two-Fluid(Fluid2)為兩相流模塊；Stress為結(jié)構(gòu)動(dòng)力學(xué)模塊；GridDeformation為網(wǎng)格形變模塊；Plasma為等離子體模塊；Electric為電模塊；Magnetic為電磁模根據(jù)不同的問(wèn)題，需要選擇相

7、應(yīng)的計(jì)算模塊進(jìn)行求解，針對(duì)本算例我們用到了Flow流體動(dòng)力學(xué)模塊、Turbulence湍流模塊、Stress為結(jié)構(gòu)動(dòng)力學(xué)模塊，所以需要的操作時(shí)在此三個(gè)模塊前面打勾選中，如圖3.2.1所示。3、選擇設(shè)置面板中的M餓項(xiàng)，其2果如圖3.3.1圖3.3.1MO設(shè)置面板MO(modeloption)為模型選項(xiàng)，主要設(shè)置相關(guān)模型的參數(shù)，根據(jù)前面PT設(shè)置中所選的模塊，此算例的MO會(huì)出現(xiàn)shared、flow、turb和stress四個(gè)設(shè)置分類，如圖3.3.1所示。3.1shared設(shè)置shared為公共參數(shù)設(shè)置,可設(shè)置描述收斂的標(biāo)題名稱，時(shí)間參數(shù)等。TransientConditionsTimeDepend

8、enceSteadyTransient口crlurCfkQC：圖3.3.1.1時(shí)間參數(shù)設(shè)置如圖3.3.1.1所示，此為時(shí)間參數(shù)的設(shè)置，若算例為瞬態(tài)，則需要加上時(shí)間參數(shù)，貝需要在transientconditions中的timedependence選擇transient。然后出現(xiàn)下圖3.3.1.2所示菜單。TransientConditionsTimeDependenceTransientI-Fas1TimesteppingTransientTimeStep-*StandardNo.ofSteps11000.三TimeStep|1E-005sTimeAccuracy1Euler圖3.3.1.2時(shí)

9、間參數(shù)設(shè)置版面其中在No.ofsteps中加入迭代計(jì)算步數(shù)，在TimeStep中輸入時(shí)間單位的時(shí)間長(zhǎng)度。3.2Flow設(shè)置Flow為模型的流體設(shè)置，其設(shè)置面板如圖3.3.2.1所示。PTMO|VC|BC|1C|SG|Out|Run|ShsfdPressureIFlowI|BRejerericePressureIlOOOOOZ/m人2TuroStressFamModelI-FanModelKflomientumResistanceModelI-IMomentumResistanceModelHemolysisModel(A(ShearABJCTimeACIl-Herriiolysis圖3.3.2

10、.1Flow設(shè)置面板針對(duì)此計(jì)算案例，只需在Flow設(shè)置面板中的ReferencePressure中輸入算例模型中實(shí)際環(huán)境的表面大氣壓力，其為一個(gè)atm,故輸入100000。其余則不選擇。3.3Turb設(shè)置PT MO VCBCCOut RunS hairedFlowTurbTurbulenceModelKEpsilonWallFunctionStandardWallStress圖3.3.3.1為Turb設(shè)置面板Turb設(shè)置中主要選擇湍流模型，和Wall的功能。其中TurbulenceModel中有如下模型選擇，按照算例的計(jì)算理論，我們需選擇K-Epsilon湍流模型計(jì)算。K EpsilonTur

11、bulenceModelWdlFunctionRNGKEPsilonKateLaunderKEpsilonTwoLayerKEpsilonRealizableKEpsilonKOmegaSST-K-OmegaLowRe(Chien)Spalart-AllmarasV2-FDetachedEddySimulation(DES)LargeEddySimulationLES)ConstEddyViscosityUserDefinedViscosity圖3.3.3.2TurbulenceModel菜單WallFunction中，我們選擇StandardWall,如圖3.3.3.3所示。Turbulen

12、ceModelKEpsilonWallFunctionStandardWallTwo-layerModelNon-equilibriumModelSharedFlowTurbStressPTMOIVCI3CIrcIscIOutIRun|Analysis-iLinearElementOrderFirstF-ShellElementsMassMatrixType-ConsistentDampingUStructuralDampingNumericalDampingModalAnalysislModalAnalysis圖3.3.4.1Stress面板根據(jù)本算例計(jì)算理論，Analysis應(yīng)選擇Line

13、ar,ElementOrder選才FFirst。如下圖所示。AnalysisLinearElementOrderNonlinear圖3.3.4.2Analysis菜單ElementOrderFirstSecond圖3.3.4.3ElementOrder菜單4、選擇設(shè)置面板中的VC選項(xiàng)PT|MOVC|BCIC|&C|Out|Run|VCSettingModePropertiesFlowPorousMediaTurbulencfSpeciesStress圖3.4.1 VC 設(shè)置面板VC(volumeconditions)為體/域設(shè)置，主要針對(duì)不同域中流體、湍流、應(yīng)力以及自身屬性進(jìn)行設(shè)置。VCSet

14、tingMode中菜單如圖3.4.1所示，針對(duì)此算例，我們只需要對(duì)算例中三個(gè)域的Properties和Stress進(jìn)行設(shè)置。4.1 Properties設(shè)置首先選擇Properties,然后點(diǎn)擊三個(gè)域進(jìn)行設(shè)置點(diǎn)擊選中其中的一個(gè)域。VoluinMNamt|0Zo-nt舂Yclumc-2Fluid1421.1Nc-NrsnieFt而d事;ZL“蠢wlumilM421.G口呼Llngrouip口圓)送然后右鍵單擊，則出現(xiàn)以下菜單,備volum#2Fluid1421.旦lankSelectedCtrl+BRotateAboutNoNameGroupingAddintoFavorites.CopyVol

15、umeconditionProperties.通過(guò)此菜單可進(jìn)行相關(guān)基本操作，如給不同的域進(jìn)行命名等操作，點(diǎn)擊Proterties ,則又會(huì)出現(xiàn)以下窗口，點(diǎn)擊 OK則完成操作同時(shí)，VC設(shè)置版面則會(huì)出現(xiàn)你選擇域的Properties設(shè)置:VC Setting Modei PropertiesProperties * FluidFluid Subtype hysDensityAirDefine,. D efine.GasMaterialNameAirGasMaterialName從數(shù)據(jù)庫(kù)選擇Air,如下圖所示。MaterialPropertySourcesjImportFromDatabaseGas

16、MaterialNameFluidAmmoniakg/m 1 30Argon岡CarbonDioxideQCarbcnMonoxide區(qū)|Chlorine同Ethane圈Flourine巨HeliumOHydrogen1網(wǎng)Methanm當(dāng)設(shè)置好流域?qū)傩院螅瑒t會(huì)出現(xiàn)流體的物理和流體參數(shù)，如果不符合算例要求，用戶可進(jìn)行修改，如下圖所示。viscosity為粘性PhysDensityFluidViscosity1 Constant(Dynamic)ConstantRho1.1614Mu1.846E-005對(duì)于1號(hào)域，其設(shè)置同理如下圖所示。VCSettingMod#PropertiesProperti

17、esSolidMaterialPropertySources-ImportFromDatabaseDefine,.SolidMaterialNameSteeLAISL1020 PhysStructDensityConstantRho 7900kg/mMaterialType-jLinearStructMaterialIsotropicMaterialYoungsModulus-JConstantN/mAE1200000000000Pois&OfisRatio-1ConstantThermalExpansionCoefficient，ConstantAlpha|LI癥0054.2 Stress設(shè)

18、置由于應(yīng)力模塊在本算例中只在1號(hào)域的橋梁實(shí)體中體現(xiàn)，所以只需進(jìn)行設(shè)置1號(hào)域中的Stress部分。選才VVCSettingMode中的Stress,然后選中一號(hào)域，則出現(xiàn)以下設(shè)置面板。activatestress激活壓力Stress“ActivateStressStressReference TemperatureConstantTref300ElementOptions(Quad,Brick,orPrism)StandardElements5、選擇設(shè)置面板中的BC選項(xiàng)BC(boundaryconditions)為邊界條件設(shè)置，主要設(shè)置計(jì)算邊界的參數(shù)。針對(duì)本算例需要設(shè)置以下五個(gè)部分：入口、出口、

19、橋體、邊界、固定點(diǎn)。針對(duì)對(duì)象選擇區(qū)域元素較多，所以我們首先應(yīng)該將其進(jìn)行分類，分為入口、出口、橋體、邊界、固定點(diǎn)五個(gè)部分。打包后為:5.1入口設(shè)置模型中入口示意圖如下圖所示:在對(duì)象區(qū)選擇入口對(duì)象，然后轉(zhuǎn)入到BC設(shè)置面板，如下圖所示PT|MO|VCBC|IC|&C|Out|Run|BCGettingMode1后也HtsralBCTypeE取(ExternalFaceonFluidVolume)一wSubTyp#1可比3njTurb-Pruir.iCanstantP01N/m2ReferencePre-sureLOOCXX)=N/mA2Temperature1CJIr*|T0CKMh口ire117d

20、0mVFixedPressureTurbPressure-ConstantP0N/mA2ReferencePressure100000jN/mA2BackflowTemperature-JConstantT300K5.3邊界設(shè)置邊界為3號(hào)域最外延除去入口和出口的部分，如下圖紅色部分所示:選擇邊界，然后將BCType設(shè)置為Wall,由于靠近邊界處的流體可忽略不計(jì),所以不進(jìn)行相應(yīng)的參數(shù)設(shè)置。其設(shè)置都保持默認(rèn)。5.4.2 Stress 設(shè)置Flow|，furblI “ i ” “ i j-a 11 ” 11StressTurlbuleniceRoughnessHeight3ConstantRHomF

21、lowTurbnimminaixijiii|:Stress:5.4橋?qū)嶓w設(shè)置橋?qū)嶓w如上圖所示部分，選中后其BC設(shè)置面板如圖所示，BCSettingMode設(shè)置為General,BCType設(shè)置為IntefaceBCSettingMode-1GeneralBCType|,imiwumi!*minAInterface(Fluid-SolidInterface)5.4.1Flow和Turb設(shè)置FlowSubType No-SlipTurbStressX-Direction VelocityConstantV-DirectionVelocityConstantm/sZ-DirectionVelocit

22、yConstantJFTl/5Flow Turb !Roughness Height* ConstantRH 0Stress由于是橋?qū)嶓w，所以其Flow以及Turb菜單無(wú)需設(shè)置。FlewSubType LoadTurbStress審 Applied PressureApplied Pressure-Constant“ImiplicrtPressureImplicitShearStressSpringForcePointLoad在SubType中選擇Load,將AppliedPressure選中，同時(shí)選中Implicit,如上圖所示。Pressure和ImplicitShearStress5.5

23、 Interface設(shè)置0(Fluid-Fluid Interface)StressInterface是指2號(hào)域和3號(hào)域之間的接觸面，如圖所示為紅色立方六面體表面。由于其為兩種大小網(wǎng)格的過(guò)渡面，不是實(shí)際的計(jì)算邊界，所以不需要做任何設(shè)置。只需將其選中，將其BCType設(shè)置為Interface5.6 橋梁固定部分的設(shè)置根據(jù)本算例實(shí)際情況，橋梁中有幾個(gè)部分為固定點(diǎn)，即在風(fēng)場(chǎng)中具形變、位移、撓度等都不會(huì)隨著風(fēng)流體的加載作用而改變，所以需要對(duì)這部分進(jìn)行特殊設(shè)置。選中固定部分，將其BCType選擇為Inteface。其Flow和Turb部分設(shè)置原理和操作與橋?qū)嶓w設(shè)置相同:FlowTurbStressSub

24、Type No-SlipX-Diir*ctiori Velocity ConstantU |om m/5V-Directicin Velocityi ConstantV 5m m 片Z-Direction VelocityConstantW flm m/sFlowi Turb iStressRoughness Heighti Consta ntRH mmStress設(shè)置如下圖所示，Sub Type 選擇 Prescribed Displacement ,根據(jù)本算例要求，需要固定三個(gè)維度上位移， 所以將Fixed X-Displacement、Fixed Y-Displacement 、Fixe

25、d Z-Displacement選中，同時(shí)用戶可根據(jù)實(shí)際情況設(shè)置 特殊大小的位移，填入 Delta X Y Z中的數(shù)值中。SubTypeiPrescribedDisplacement歹FixedX-Diaplacement17FixedYDisplacernentPFixedZDisplacementNodalPrescribedDisplacementX-Displacement1ConstantDeltaX|0mY-uispiacemefn:-1Cw式antDekaY|o二m4-UIsplacementCoristanl:Delta2|omiStress;Flow8Turb6、選擇設(shè)置面板

26、中的IC選項(xiàng)IC (initial conditions）為初始條件設(shè)置，主要設(shè)置計(jì)算初始參數(shù)。具體操作如下:1 、將 Initial Conditions設(shè)置為 User Specified ; IC Applied選擇為 Volume by VolumePT | M0 | VC | BC IC SC | Out | RunIC Option (For whole sinnulation, Apply button not applicable)Initial Condition -t User SpecifiedIC Applied ( Volume by Volume由于2號(hào)域和3號(hào)域都

27、為空氣流體域，所以兩者可以共同設(shè)置，同時(shí)選擇2號(hào)和3號(hào)域。Volume NameVC TypeBlanks 1 ICG roupZoneKey0Molume2Fluid1一volurnelFluid242 L.YolumelSolid3421.,需要對(duì)其進(jìn)行Shared、Flow、Turb三部分設(shè)置，其設(shè)置方法與邊界條件設(shè)置中的入口設(shè)置相同，如下圖所示:|Shedl: TemperatureFlowConstantTurbT 300Shared下嬴1XDirectioni VelodlyTurb-11c55 知 nt u gYDireetionVelocitySharedFlowI Turb

28、Ki ConstantmA2/s20.011 Constant口0.01-Kinetic EnergyDissipation RateConstantVwtn/sz-Uirectionveiocrty一匚m/sw.7.9JPressurejConstantN/mA2N/mA2pReferencePressure0zl1MOOO3、選擇1號(hào)域。VolumeN31TleVCTypeBlankeICGroupZoneKey工volume2Fluid142L二volumelFluid23“*volumelSolid3Shared設(shè)置面板為SharedTemperatureStressConstantT

29、300Stress設(shè)置面板為:SharedX-DirectionResidualStressf+ve=tensile)Stress!*:1ConstantSigm3xx2N/mA2Y-DirectionResidualStress(fve=tensile)-JConstantSigmayy|oJN/mA2Z-DirertionResidualStress(+-ve=tensile)iConstantSigmazzoN/mA2備注：如果計(jì)算需要從已有的計(jì)算結(jié)果開始計(jì)算則需要將InitialConditions設(shè)置為PreviousSolution，同時(shí)選擇上AppendResidualFile和

30、AppendOutputFile,然后在RestartFilename中添加初始計(jì)算文件，如下圖所示ICOption(Forwholesimulation,Applybuttonnotapplicable)InitialCondHion-iPreviousSolution7AppendResidualFileAppendOutputFileRestartFilenamegk2-liugudingcharrg-1200-7i9,DTFBrowse(ICAppliediForAllVolumes7、選擇設(shè)置面板中的SC選項(xiàng)SC（solvercontrol）為計(jì)算控制設(shè)置，其主要調(diào)節(jié)計(jì)算步數(shù)，收斂標(biāo)

31、準(zhǔn)，松弛因子等計(jì)算參數(shù)，通過(guò)設(shè)定，可以根據(jù)不同的算例得到準(zhǔn)確的模擬仿真結(jié)果。SC設(shè)置主要分為：IterSpatialSolversRelaxLimisAdv。因?yàn)楸舅憷裏o(wú)需設(shè)置過(guò)多下面只介紹Iter和Relax兩個(gè)設(shè)置方法，其余設(shè)置均保持默認(rèn)就能達(dá)到本算例的計(jì)算要求。Iter設(shè)置面板如下，主要設(shè)置最大計(jì)算迭代步數(shù)，和收斂標(biāo)準(zhǔn)，根據(jù)本算例實(shí)際情況我們?cè)O(shè)置Max.Iterations為10000步，Convergence（收斂）Crit為0.001,Min.Residual（殘余誤差）為1E-08。Relax設(shè)置面板如下，其主要是調(diào)節(jié)松弛因子，當(dāng)收斂曲線收斂不明顯時(shí)可IterSpatialSolv

32、ers進(jìn)行調(diào)節(jié)松弛因子，如下圖，如果將箭頭都向上調(diào)節(jié)，則收斂性就會(huì)改善I，mm、LimitsAdvTurbulencer Turbukntf Start Control,Relax!8、選擇設(shè)置面板中的OUT項(xiàng)OUT為輸出設(shè)置，主要控制仿真計(jì)算結(jié)果的輸出控制，其中菜單包含如下圖所示。GeneralRestsrtGraphicStiEEaF沁sMon征orPointMonitorPlaneUserAccessControl由于本算例輸出設(shè)置只涉及General、Graphic兩部分，所以針對(duì)此兩部分進(jìn)行講解。General設(shè)置版面如圖，OutputBy設(shè)置為Iteration;IntervalO

33、ptions設(shè)置為SpecfiedInterval（固定間隔）；由于本算例需要對(duì)沒(méi)間隔100步進(jìn)行結(jié)果分析，所以選中UniqueFile;設(shè)置Every為100。OutputControlsOutputByIterationIntervalOptionsSpecifiedIntervalDiagnosticsDiagnosticsGraphic為輸出結(jié)果控制，其輸出結(jié)果主要分為Shared、Flow、Turb、Stress四個(gè)部分，通過(guò)選擇使得仿真輸出相應(yīng)結(jié)果，具體操作為在相應(yīng)輸出結(jié)果前打鉤選中。針對(duì)本算例則需要輸出以下結(jié)果。-Shared“DensityStaticTemperaturerT

34、otalTemperatureFlowVelocityVector。Velocity1Magnitude“StaticPressure。TotalPressure歹LaminarViscosityrVorticityStrainRateI-X-DirectionVelocityResidualrY-DirectionVelocityResidualI號(hào)DirectionVelocityResidualI-PressureResidualTurbulenceTurbKineticEnergyyTurbDissipationRateTurbViscosityI-EffectiveVise.Y+I-

35、RGridParameter1 -TurbKineticEnergyResidualTurbDissipationRateResidualStress2 DisplacementPCartesianStressTensorFCartesianStrainTensor7 PrincipalStress審PrincipalStressDirectionQPrincipalStrainrReactionForces9、選擇設(shè)置面板中的RUN項(xiàng)其版面如下圖所示，在此設(shè)置版面中，控制仿真計(jì)算的運(yùn)行，有開始、暫停、終止、觀察其收斂曲線、輸出結(jié)果等功能。PT|MO|VC|BC|C|SC|OutI。|Paralle