cfx和ansys的算例：水管的耦合问题算例2

1.问题描述：一根弯管，里面有流体入口流体速度10m/s ，开放出口压力（opening），管道两端固支。

现在想用ansys和cfx的MFX的流固耦合做个练习，观察在水流冲击下管道的变形情况。

2.模型描述：管道模型，网格，边界条件和接触面apdl/prep7!Selection tolerance!set element typeet,1,shell63 ! 3-D 20-Node StructuralR,1,0.01, , , , , ,!!set materialmp,ex,1,2.1E11 !Young modulusmp,prxy,1,0.3 !Poisson coefficientmp,dens,1,7800!simple pipe modelk,1,k,2,1k,3,0,1l,1,2l,1,3LFILLT,2,1,0.25, ,LPLOT WPSTYLE,,,,,,,,1 KWPA VE, 2 wpro,,,-90.000000 CSYS,4CYL4, , , ,0,0.1,90 CYL4, , , ,-90,0.1,0VDRAG,1,2 , ,, , ,1,3,2 vglue,allvdele,all,,0aplotFLST,2,6,5,ORDE,6 FITEM,2,4 FITEM,2,12 FITEM,2,20 FITEM,2,30 FITEM,2,33 FITEM,2,36 ADELE,P51X, , ,1FLST,2,7,5,ORDE,7 FITEM,2,5 FITEM,2,-6 FITEM,2,13 FITEM,2,-14 FITEM,2,21 FITEM,2,28 FITEM,2,31 ADELE,P51X, , ,1FLST,2,4,5,ORDE,4 FITEM,2,1 FITEM,2,22 FITEM,2,27 FITEM,2,34 ADELE,P51X, , ,1 aglue,alltype,1real,1esize,0.03MSHAPE,0,2DMSHKEY,1Amesh,all!boundary conditionnsel,s,loc,x,0D,all, , , , , ,ALL, , , , ,nsel,s,loc,z,0D,all, , , , , ,uz, , , , ,allsel,all!set fsi conditionsf,all,fsin,1allselsavecdwrite,db,example_shell,cdbfinish流体模型，网格，边界集合apdl/prep7/prep7et,2,fluid142,,,,1 !3D Fluid element with diplacement DOF optionet,3,shell63 !Mesh only element (3D quad 4 nodes) to mesh surfaces used in CFXpre !Fluid domain geometryk,1,k,2,1k,3,0,1l,1,2l,1,3LFILLT,2,1,0.25, ,LPLOTWPSTYLE,,,,,,,,1KWPA VE, 2wpro,,,-90.000000CYL4, , , ,0,0.1,90CYL4, , , ,-90,0.1,0VDRAG,1,2 , ,, , ,1,3,2vglue,allaplot!Fluid domain meshingallseltype,2mat,2esize,0.02vsweep,all!FSI interface surface meshasel,s,,,3asel,a,,,11asel,a,,,19asel,a,,,29asel,a,,,32asel,a,,,35ALLSEL,BELOW,AREAaplottype,3 !with mesh only elements amesh,allcm,fsi,elem !Create component named fsi allselASEL,S, , ,34ASEL,a, , ,22ALLSEL,BELOW,AREAtype,3 !with mesh only elements amesh,allcm,inlet,elem !Create component named inlet allsel,allASEL,S, , ,1ASEL,a, , ,27ALLSEL,BELOW,AREAtype,3 !with mesh only elementsamesh,allcm,outlet,elem !Create component named inletallsel,allasel,s,,,4asel,a,,,20asel,a,,,12asel,a,,,33asel,a,,,36asel,a,,,30ALLSEL,BELOW,AREAaplottype,3 !with mesh only elementsamesh,allcm,sym,elem !Create component named symallsel,allcdwrite,db,fluid,cdb !Create fluid.cdb file for CFXpre3.生成dat和defSet up the CFX Model and Create the CFX Definition FileSet up the example in the CFX preprocessor1.Start CFXpre from the CFX launcher.2.Create a new simulation and name it cfx_mfx3.Load the mesh from the ANSYS file named fluid.cdb. The mesh format is ANSYS.Accept the default unit of meters for the model.4.Define the simulation type:1.Set Option to Transient.2.Set Time duration - Total time to 0.5 s. Note: this value will be overridden byANSYS.3.Set Time steps - Timesteps to 0.005 s. Note: this value must be equal to the timestep set in ANSYS.4.Set Initial time - Option to Value, and accept the default of 0 s.5.Create the fluid domain and accept the default domain name. Use Assembly as thelocation.6.Edit the fluid domain using the Edit domain - Domain1 panel.1.Set Fluids list to Air at 25 C.2.Set Mesh deformation - Option to Regions of motion specified. Accept thedefault value of mesh stiffness.3.In the Fluid models tab, set Turbulence model - Option to None (laminar).4.Accept the remainder of the defaults.5.Initialize the model in the Initialisation tab. Click Domain Initialisation, and thenclick Initial Conditions. Select Automatic with value and set velocities and staticpressure to zero.7.Create the interface boundary condition. This is not a domain interface. Set Name toInterface1.1.In the Basic settings tab: - Set Boundary type to Wall. Set Location to FSI.2.In the Mesh motion tab: Set Mesh motion - Option to ANSYS Multifield.3.Accept the defaults for boundary details.8.Create the opening boundary condition. Set Name to Opening.1.In the Basic settings tab: Set Boundary type to Opening. Set Location to outle.2.In the Boundary details tab: Set Mass and momentum - Option to Static pres.(Entrain). Set Relative pressure to 0 Pa.3.In the Mesh motion tab: Accept the Mesh motion - Option default of Stationary. 9.Create the inlet boundary condition. Set Name to inlet. Edit the inlet boundary conditionusing Edit boundary: inlet in Domain: Domain1 panel.1.In the Basic settings tab: Set Boundary type to inlet. Set Location to inlet.2.In the Boundary details tab: Set Mass and momentum - Option to normal speed.Set normal speed value to 03.In the Mesh motion tab: Set Mesh motion to Stationary.10.Generate transient results to enable post processing through the simulation period.1.Click Output Control.2.Go to Trn Results tab.3.Create New. Accept Transient Results as the default name.4.Choose Time Interval and set to 0.005。