水动力计算软件AQWA 培训ppt课件.ppt

《水动力计算软件AQWA 培训ppt课件.ppt》由会员分享，可在线阅读，更多相关《水动力计算软件AQWA 培训ppt课件.ppt（129页珍藏版）》请在三一办公上搜索。

1、25th 29th september 2006,AQWA TRAINING COURSE,目录,培训内容:AQWA程序分析能力理论背景模拟技术modelling techniques分析过程分析数据的准备结果的输出内容和解释实例,客户,Regulatory BodiesABS(US),Lloyds Register(UK)OperatorsBP,US Navy,NavionEngineering CompaniesJohn Brown Hydrocarbons,Kellog Brown and RootArup Energy,Friede and GoldmanMitsubishi Heav

2、y Industries,Kawasaki Heavy IndustriesPI Rauma,IntecWorley Engineering,IHC Gusto,SBM,客户,Consultants:Armstrong TechnologyJP KennyJJ McMullenMoffatt&NicholNoble DentonHan PadronAtkinsAMOGD2MAtlantiaMoss MaritimeInocean,AQWA 程序介绍,AQWA的程序结构和分析能力AQWA LINE-3-D diffraction eigen mode and dynamic stability

3、analysisAQWA FER-Spectral analysis of structure motion(wave frequency or/and drift frequency)and mooring tension in irregular waves,AQWA 程序介绍,AQWA NAUT-Time domain program for wave frequency structure motion and mooring tension analyses in regular wavesAQWA DRIFT-Time domain program for drift freque

4、ncy and wave frequency structure motion and mooring tension analysis in irregular wavesAQWA Graphical Supervisor(AGS)-AQWA pre and post processor;on-line analysisAQWA WAVE-Interface between hydrodynamic program and FE program to transfer wave loads from AQWA LINE to a FE model for structural strengt

5、h analysis,AQWA 程序介绍,LINE,LIBRIUM,FER,NAUT,DRIFT,WAVE,ASAS(FE model),General Relations between Programs,Typical AQWA Model(1)-Moored Tanker,Typical AQWA Model(2)-Semi Sub,AQWA Result(1)-Force&Response Curves,AQWA Result(2)-Pressure Contour,AQWA Result(3)-Shear Force&Bending Moment,AQWA LINE 介绍,AQWA-

6、LINE 是3D衍射和辐射分析程序通过大量的板格来描述结构在每个板格中心有源点（a source），程序会求解受到如下边界条件时源点的强度:没有水流穿过船体没有水流穿过海底板格为自由面这样就可以在任意点计算出压力和速度,AQWA LINE理论,AQWA LINE的3D线性波理论假设流体无旋，不可压小波（small wave elevation）速度势能控制方程体的边界条件,AQWA LINE理论,自由表面运动条件 线性化之后自由表面动力条件,AQWA LINE理论,线性化之后组合自由表面条件 海底边界条件,AQWA LINE理论,辐射条件辐射条件是一种物理状态，使波不会在错误的方向上传递，如水

7、穿透船体，从而避免数上的不确定性。数值方法势能项的线性叠加 式中：下标I是入射波，d是衍射波,j=1,2,6 是6个自由度的辐射波，x是单位波幅下结构的运动.,AQWA LINE理论,对于有限水深d，入射波的势能为式中k是波数，定义为:使用脉冲源分布（pulsating source distribution）求解衍射和辐射波势能,AQWA LINE理论,式中：是辐射体强度（source strength）;S 是入水结构面;(x,y,z)为流域中域点的坐标;(,)是S上源点坐标;G 为格林函数，是求解Laplace方程的基础，并且满足所有边界条件（除了体边界条件）.G 可以表达为:,AQWA

8、 LINE理论,式中pv 表示积分的principal value;J0 是Bessel函数的第一项,AQWA LINE理论,在每个板格结构表面上的源强度（source strength）假设为常数，通过体边界条件求解积分方程计算:对于衍射势能，在结构表面由于入射势能产生的法向速度会变为零，而所减少的法向速度会转化为结构运动.,AQWA LINE理论,压力和第一阶波浪力的计算 一旦计算了源点强度和势能，从线性化的Bernoulli方程中可以计算每个板格的水动力压力:第一阶波浪力是通过在体表面上进行积分得到的,二阶平均漂浮力计算远场求解（Far field solution）(动量守恒方法，只用

9、于水平力）):式中 SR 是在流场结构周围的垂直圆柱边界其半径为R，是SR的流体体积，The mean force is the time average of the above expression and the first order term becomes zero.The pressure includes the non-linear term in the Bernoulli equation and therefore will not disappear.,AQWA LINE理论,AQWA LINE理论,近场求解(压力对于六个自由度进行积分方法):式中 WL 为水线位置;

10、r 是相关波表面（relative wave surface elevation）;S0 水下结构表面;X 结构表面的运动;MS 结构的质量;R是结构的旋转矩阵;是结构中心加速度.,AQWA LINE理论,Quadratic Transfer Function(QTF)Components at both difference and sum frequenciesEach with real and imaginary parts,AQWA LINE理论,In which,AQWA LINE理论,运动方程 结构在波浪中的响应X是通过计算下面方程得到:式中 Ms 是结构质量,Ma 是附加质量,

11、C 是阻尼,K 流体刚度,F 是波浪力（包括衍射和辐射力）.,AQWA波概述,波的类型:1)小波即线性波（Airy Waves）a=A cos(-t+kx)(:为频率，弧度/秒;k:波数)在以下模块都有：AQWA LINE,LIBRIUM,FER,DRIFT,NAUT 2)2阶波（Stokes 2nd Order Waves）a=A cos(-t+kx)+0.5 k A cos2(-t+kx)AQWA NAUT模块缺省值用2阶波,AQWA波概述,波的类型(续):3)流量函数波(只用于深水和Morison单元)a)载波表面压力为零；或者b)在波表面不运动。,AQWA波概述,波浪形式:1)规则波用

12、于AQWA LINE,NAUT(缺省)2)不规则波通过定义波谱确定，用于AQWA LIBRIUM,FER,DRIFT.不规则波也可以在NAUT中定义(with an IRRE card in the JOB deck).,AQWA波概述,波浪形式(续):在AQWA中的波谱类型:a.P-M spectrumb.JONSWAP spectrumc.User defined spectrumd.Gaussian spectrum for Cross Swell不规则波可以是一下形式:a.长波;或b.短波,也就是 a spread sea(只用于 AQWA LIBRIUM and FER),结构上的波

13、力,对于衍射结构(通过板单元来模拟)-入射波力(Froude-Krylov force):静止波压力产生（from the pressure in the undisturbed waves）.-衍射力:由于固定结构干扰入射波产生.-辐射力：由于结构摆动产生的波产生.-平均漂浮力(net force due to high order effect),结构上的波力,对于Morison结构(用Morison单元模拟,eg TUBEs)-Morison 力(包括阻力，drag)是用Morison方程计算的。.,Morison 力,Morison力的计算方程对于细长圆柱单元，例如管道单元，单元上单位

14、长度的水动力可以用Morison方程计算:式中Ca 和 Cd是附加质量和单元的阻力（drag）系数；;是每单位长度单元的体积，D是单元的直径.,结构上的波力,AQWA 的总体坐标系,定向特性:原点在静水线上Z轴的正方向是垂直向上X轴的正方向是右舷Y轴的正方向是根据右手法则,z,0,z,y,x,W.L.,AQWA 的总体坐标系,AQWA的总体坐标系是固定的参考轴，而与东、南、西、北等方向无关,AQWA的方向,AQWA中定义的波、风和水流方向是作为船体的行进方向.在定义的波浪(或水流、风)之间有角度时，角度是与X正方向逆时针转为正.中的方向是安度输入输出.,X axis,Wave directio

15、n(or current,wind),positive angle,AQWA坐标系的方向,对于非对称结构，输入的方向的角度必须在范围内-180 度到+180度内；对于X方向对称的结构，角度范围在0度到180度之间；对于X、Y方向对称结构，角度范围在0度到90度之间.,y,x,v,q,相位角,In AQWA LINE,the phase angle()of a parameter defines the time difference(dt)from the time when wave crest is at the CoG of the structure to the time when

16、the parameter reaches its peak value.(dt=*T/360,where T is the wave period).正的相位角表示滞后于波.,AQWA 的文件名,每个AQWA文件有三个部分:文件前缀（两个字母）-表示程序类型alLINEabLIBRIUMafFERadDRIFTanNAUT运行的文件(最多六个字母),AQWA 的文件名,文件扩展名(三个字母)文件扩展名说明文件类型(eg,.dat)例如:altest.dat,输入数据(.dat)文件,This is an ASCII text file containing all the input dat

17、a necessary for the Stages of Analysis about to be executed.The data is in fixed format and must be entered using a text editor.An editor which allows 80 columns of data to be displayed,and which indicates the column number of the current cursor position,is highly recommended.Note:A graphical user i

18、nterface,which allows interactive data input in a user-friendly environment,is currently under development.,列表文件(.lis),This is an ASCII text file containing all the output data(in text form)from the Stages of Analysis which have just been executed.It can be examined using a text editor.,重启动文件(.res),

19、This is a binary file,which is written by all the AQWA programs,and which contains a summary of the data associated with all the Stages of Analysis which have so far been executed.Examples:If Stages 1 to 4 have been executed,the restart file will contain:-model definition-hydrodynamic database-main

20、analysis parameters,水动力文件(.hyd),This is a binary file which is created by AQWA-LINE at the end of the diffraction/radiation analysis(Stage 3).It contains the hydrodynamic database produced by the AQWA-LINE run.At the end of Stage 3,the contents of the restart and hydrodynamics files are:restart file

21、(RES)hydrodynamics file(HYD)model definition hydrodynamic databasehydrodynamic database,水动力文件(.hyd),The hydrodynamics file is therefore a sub-set of the restart file.The other programs do not create a hydrodynamics file,but they can read and use this file(only when ALDB option is on,or a FILE card w

22、ith*.HYD file name in deck 6).,定位文件Position(.pos and.eqp),Both are binary files.eqp file:Created by AQWA LIBRIUM which stores the equilibrium positions of a system of structures.This file can be read in by FER,DRIFT and NAUT as the start position for the analysis(an option RDEP is required).pos file

23、:Created during a time domain analysis by DRIFT or NAUT in which the positions,velocities,etc of a system of structures are stored for every time step.,AGS 显示文件(.plt),This is a binary file,which is created during a main analysis(Stage 5)and contains either:the time history of forces and motions(NAUT

24、,DRIFT)orforces and responses as a function of frequency(FER,LINE),AQWA 文件类型,每次运行AQWA都会包括多个文件。文件主要是后缀不同，一般的扩展名如下:.dat输入数据文件(ASCII).lis输出列表文件(ASCII).res重启动文件(binary),输入、输出.hyd水动力文件(binary),输入、输出.pos输出定位（position）文件(binary).plt输出图形文件(binary).eqp平衡位置文件(binary),输入、输出.uss源强度文件(binary),输入、输出,AQWA 文件类型,.me

25、s output message file(ASCII).potpotential file(binary),input and output.pac pressures at centroids(binary),output.vacvelocities at centroids(binary),output.msd input mass distribution for BM/SF(ASCII).sfm input mass distribution for splitting forces(ASCII).lininput file for AGS mesh generator(ASCII)

26、.wht input file specifying a wave height time history(ASCII).wvt input file specifying a wind velocity time history(ASCII).xft input file specifying an external force time history acting on a structure(ASCII),AQWA Restart Stages,Stage 1定义模型 Decks 1 to 5Stage 2水动力数据 Decks 6 to 8Stage 3衍射/辐射分*Stage 4主

27、要分析参数 Decks 9 to 20Stage 5主要分析*只是计算步骤,Stage 1,Deck 1COOR节点坐标Deck 2ELM*单元定义Deck 3MATE材料特性Deck 4GEOM几何特性Deck 5GLOB整体约束*结构编号,Stage 2,Deck 6FDR*规则波定义(频率和方向)Deck 7WFS*水动力特性静水特性Analysis PositionDeck 8DRC*漂流力系数*结构编号注:在第二步输入的水动力特性是用修改或替代第三步AQWA-LINE 的计算结果,Stage 3,该步是AQWA-LINE的主要分析，但只是一个程序计算过程而已.注:在第二步输入的水动力

28、特性是用修改或替代第三步AQWA-LINE 的计算结果,Stage 4,Deck 9DRM*漂流运动参数Deck 10HLD*船体阻力（Drag）系数Deck 11ENVR环境参数Deck 12CONS约束Deck 13SPEC谱参数WAVE规则波参数(N)Deck 14MOOR停泊线定义*Structure Number,Stage 4(cont),Deck 15STRT初始条件Deck 16TINT 时间积分参数(D,N)LMTS迭代参数(B)Deck 17HYDC管道的水动力参数Deck 18PROP打印输出选项Deck 19NONEReserved for future useDeck

29、 20NONEReserved for future use,Stage 5,该步骤只是程序求解计算.,分析方法,一般的分析方法是在AQWA-LINE中运行Stages 1 to 3，然后在其它程序运行 Stages 4 to 5，例如 AQWA-NAUT.实例:Step 1:AQWA-LINE(Restart 1 to 3)runInput FilesOutput Filesaltest.dataltest.lisaltest.resaltest.hydaltest.plt,分析方法(续),Step 2:文件拷贝copy altest.resabtest.resStep 3:AQWA-LIB

30、RIUM(Restart 4 to 5)runInput FilesOutput Filesabtest.databtest.listabtest.resabtest.resabtest.eqp,AQWA 单元类型,The following types of elements are accepted in AQWA:QPPL:Quadrilateral panelTPPL:Triangular panelTUBE:Tube element(circular cross section)STUB:Slender tube element(non-circular cross section

31、allowed)PMAS:Point mass,AQWA 单元类型,PBOY:Point buoyancyFPNT:Field point(for wave surface calculation),should only be defined after QPPL and TPPLDISC:Circular disc with no thickness.An additional parameter DIFF is needed for diffracting QPPL and TPPL elements.,AQWA LINE警告检查,General warningsNo.of elemen

32、ts 6000:solution timeNormals point out:modelling conventionNo gaps:force balanceFacets cannot cut surface:solution requirementDimensions KR:good practice,AQWA LINE警告检查,Stage#1 checksArea ratio of adjacent elements Element centres at least one facet radius apart Shape factor 0.2 warningShape factor 0

33、.02 fatal error,Warnings in AQWA LINE,Stage#2 checksLongest side 0.000273 D for reasonable no.of steps to convergence.Depth below WL 0.0000132.dittoDepth below WL 0.5.rf.nz Distance above sea bed must be 0.5.rf Warning if nodes not connected to another element.,AQWA-LINE 的实例输入数据,AQWA-LINE 实例输入数据(续),

34、定义Morison单元,AQWA LINE的常用选项(1),DATA用于检查输入数据(LBDFN)GOON用于忽略一般警告(L)REST定义开始过程(LBDFN)LDOP载荷输出-输出.POT 和.USS 文件 用于压力计算(例如压力显示,SF/BM)(L)PRCE为 Decks 1 5(LBDFN)设置打印项PPEL为每个单元设置打印属性(LBDFN),AQWA LINE的常用选项(2),NPPP无压力后后处理阻止节点连接警告ALDBAqwa-Line 数据库-reads in AQWA-LINE database TAKE CARENQTF使用近场求解，用于漂浮力系数CQTF计算 QTF

35、矩阵PFIXPartially FIXed,AQWALINE 后处理,波表面云图显示,AQWALINE 后处理,Air Gap,AQWALINE后处理,压力云图,AQWALINE后处理,衍射波表面,AQWA 数据处理,处理原因:为了添加附加节点、单元、阻尼等，或者把多个数据库合并成一个，而不重新完全运行AQWA-LINE分析.处理过程:1)拷贝已有数据文件AL1.DAT，改名为 AL2.DAT;2)添加新节点，（非衍射）单元，附加阻尼 等(CRNM 选项需要用 RAO 重新计算)到AL2.DAT文件。,AQWA 数据处理(续.1),3)删除deck 6中所有波频率和方向并改为:06 FDR10

36、6FILE AL1.HYD06CSTR N(N is the structure number in AL1)END06CPDB4)在 AQWA-LINE中运行AL2.DAT,AQWA 数据处理(续.2),For对于复杂结构数据和合并，新的AL*.DAT文件应该包括所有的结构定义，并且每个结构都有一个deck 6，包括文件卡片中对于的文件名.以上过程也可以应用于其它AQWA程序.,AQWA 数据处理(续.3),JOB MESH NAUTTITLE A SHIP+TUBES AND DISCSOPTIONS REST ENDRESTART 1 5 01 COOR 010001 10.0 0.00

37、0 2.000 010002 10.0 0.000 0.000.END010999 0.00 0.000-0.500 02 ELM1 02QPPL DIFF(1)(1)(2)(12)(11)02QPPL DIFF(1)(21)(22)(32)(31).02TUBE(1)(101)(102)(1)(1)02DISC(1)(1)(2)(3)END02PMAS(1)(999)(2)(2)03 MATE 03 1 7.8E3 END03 2 1.0E7 04 GEOM 04TUBE 1 1.00 0.05 0.00 04CONT 0.75 1.00 04DISC 3 1.20 04CONT 1.14

38、1.00 END04PMAS 2 1.0E9 0.00 0.00 9.0E9 0.0 9.0E9,AQWA 数据处理(续.4),相互作用的结构(1),PFIX方法把漂浮模型和固定模型合并到一个结构中利用特殊的组编号999使用PFIX选项设置固定结构的基础,相互作用的结构(2),相互作用水动力AQWA能够计算水动力系数，可以完全的考虑水动力之间的相互作用.可以包括达到10浮体结构相互作用.,结构间的水动力响应作用,对称性,SYMX 表示AQWA 可以假设结构是X轴（FRA）对称的。利用对称性可以减少计算时间.RMXS不用对称（remove symmtry），创建完整模型（即使模型是对称的）。该方

39、法只用于T/QPPL单元.MSTR把结构移动到新定义的位置。只用于单元和与单元关联的节点，而不是所有在STRC卡片中列出的节点.实际上在FRA中移动节点和单元，而与Deck 15卡片中的位置不同.SYMY 和RMYS 对于Y轴有相同的影响,AQWA LIBRIUM 介绍,AQWA-LIBRIUM 用于确定平衡位置，静稳定性和动稳定性结构是以很小的步骤移动直到力平衡，从而找到平衡位置然后，程序计算线性化刚度矩阵的特征值，得到静稳定性通过阻抗矩阵的特征值得到动稳定性,AQWA LIBRIUM理论,确定静平衡位置的方程.结构的X位置通过下式计算:式中K是系统的刚度矩阵，F是力.程序迭代直到X小于所定

40、义的容差,分析方法,Step 2:File Copyingcopy altest.resabtest.resStep 3:AQWA-LIBRIUM(Restart 4 to 5)runInput FilesOutput Filesabtest.databtest.listabtest.resabtest.resabtest.eqp,分析方法(续),在AQWA-LIBRIUM主分析开始时(Stage 5)，需要提供如下数据:主模型定义restart file水动力数据 restart file主分析参数 input data fileNote:As the data from the resta

41、rt file corresponds to decks 1 to 8,only decks 9 to 20 are required in the input data file(decks 1 to 8 must be omitted).,AQWA-LIBRIUM 实例数据(1),AQWA-LIBRIUM实例数据(1)(cont.),AQWA-LIBRIUM实例数据(2),AQWA-LIBRIUM实例数据(2)(cont.),AQWA LIBRIUM常用选项,STAT只用于静稳定性分析DYNA只用于动力稳定性分析LSTF线性刚度（Linear STiFness），用从AQWA-LINE得到

42、的水静力刚度.PBIS打印迭代步,Drag in AQWA(deck 10),1)水流和风力系数CUFX Dir1 DirN C1.CNWIFX Dir1 DirN C1.CN etc(Dir:方向编号;C1:阻力（Drag）系数在方向的相对水流速度force in X direction=CUFX.V2force in Y direction=CUFY.V2yaw moment=CURZ.V2,Drag in AQWA(deck 10)(cont 1),2)Morison阻力（drag）系数 MDIN Nrow Ncol C1.C6Nrow:阻力矩阵的行数Ncol:阻力矩阵的列数C1 C6:

43、阻力（Drag）系数,AQWA约束(deck 12),1)忽略用户指定的运动自由度DACF Ns Ndof(Ns:结构数;Ndof:自由度数)2)结构间的机械连接 不允许关联的平移运动，但允许关联的选择运动.DCON Nt Ns1 Nd1(Nd3)Ns2 Nd2(Nd4)Nt:number of D.O.F.being locked by this constraint.Nt=0:Ball and Socket,rotation in 3 D.O.F.Nt=1:Universal joint,rotation in 2 D.O.F.Nt=2:Hinge,rotation in 1 D.O.F.

44、Nt=3:Rigid connection,no rotation.,AQWA约束(deck 12)(cont.1),AQWA约束(deck 12)(cont.2),AQWA约束(deck 12)(cont.3),AQWA停泊索(deck 14),AQWA 程序除了AQWA LINE 外，都可以定义停泊索。下面是常用的停泊类型(更详细的内容参考AQWA手册):LINE:线弹性线（无重量）LINE Ns1 Nd1 Ns2 Nd2 K L(Ns1,Ns2:structure numbers;Nd1,Nd2:node numbers;K:stiffness;L:unstretched length)

45、POLY:多项（Polynomial）弹性线（无重量）.POLY K1 K2 K3 K4 K5NLIN Ns1 Nd1 Ns2 Nd2(Ts)L(Fw)(Fp)(K1,.,K5:stiffness;Ts:winch tension;Fw:winch winding in friction factor;Fp:winch paying out friction factor;winding in tension Tw=Ts*(1-Fw);paying out tension Tp=Ts*(1+Fp);Ts,Fw and Fp are only needed when the POLY line

46、is used as a winch),AQWA停泊索(deck 14)(cont.1),COMP/ECAT:合成弹性悬链线，有重量.14COMP Nz Nx Ne Zmin Zmax Slope 14ECAT M1 A1 EA1 Tmax1 L1 14ECAT M2 A2 EA2 Tmax2 L2 14ECAT M3 A3 EA3 Tmax3 L3 14NLIN Ns1 Nd1 Ns2 Nd2式中 Nz,Nx-number of database points within z and x ranges.Ne-number of ECAT in this COMP line.Zmin,Zma

47、x-Z range(measured from the anchor)for the attachment node.Slope-sea bed slope(in degrees;positive for slope going up from anchor towards attachment point).M1,M2,M3-mass per unit length for ECAT 1,2,3.A1,A2,A3-equivalent cross section area.EA1,EA2,EA3-Youngs modulus x area.Tmax1-3-maximum tension.L1

48、,L2,L3-length of ECAT 1,2,3.Ns,Nd-structure number and node number.,停泊数据,AQWA停泊索(deck 14)(cont.2),AQWA 览索动力学(只用于COMP/ECAT):14COMP 20 30 3 770.0 800.0 3.00 14ECAT 170.00 A 6.0000E8 7.500E6 250.0 14ECAH Ca Cd De Cx 14ECAT 200.00 A 9.0000E8 7.500E6 350.0 14ECAT 170.00 A 6.0000E8 7.500E6 600.0 14NLID 1

49、9910 0 9911 式中 ECAH:card name for catenary hydrodynamic coefficient definition.A:equivalent cross section area(must be non-zero).Ca:added mass coefficient.Cd:transverse drag coefficient.De:effective diameter for drag.Cx:inline drag coefficient.NLID:same as NLIN but with cable dynamics.,AQWA停泊索(deck

50、14)(cont.3),防护装置（Fenders）14POLY K1 K2 K3 K4 K5 14FEND Size(Kf)(Kc)14FLIN Type Ns1 Nd1 Nd2 Ns2 Nd3 Nd4式中 K1 K5-non-linear stiffness coefficientsSize-uncompressed size of fenderKf-tangential friction coefficientKc-normal damping coefficientType-1=fixed fender,2=floating fenderNs1-Structure to which fe