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

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

1、 2006 ANSYS,Inc.All rights reserved.1ANSYS,Inc.Proprietary2525th th 2929th th septemberseptember 20020066 2006 ANSYS,Inc.All rights reserved.2ANSYS,Inc.Proprietary目录培训内容:AQWA程序分析能力理论背景模拟技术modelling techniques分析过程分析数据的准备结果的输出内容和解释实例 2006 ANSYS,Inc.All rights reserved.3ANSYS,Inc.Proprietary客户客户Regulat

2、ory BodiesABS(US),Lloyds Register(UK)OperatorsBP,US Navy,NavionEngineering CompaniesJohn Brown Hydrocarbons,Kellog Brown and RootArup Energy,Friede and GoldmanMitsubishi Heavy Industries,Kawasaki Heavy IndustriesPI Rauma,IntecWorley Engineering,IHC Gusto,SBM 2006 ANSYS,Inc.All rights reserved.4ANSYS

3、,Inc.Proprietary客户客户Consultants:Armstrong TechnologyJP KennyJJ McMullenMoffatt&NicholNoble DentonHan PadronAtkinsAMOGD2MAtlantiaMoss MaritimeInocean 2006 ANSYS,Inc.All rights reserved.5ANSYS,Inc.ProprietaryAQWA AQWA 程序介绍程序介绍AQWA的程序结构和分析能力的程序结构和分析能力AQWA LINE-3-D diffraction&radiation analysis program

4、 for wave force and structure response calculations;hydrostatic analysisAQWA LIBRIUM-Structure equilibrium position and force balance calculations;eigen mode and dynamic stability analysisAQWA FER-Spectral analysis of structure motion(wave frequency or/and drift frequency)and mooring tension in irre

5、gular waves 2006 ANSYS,Inc.All rights reserved.6ANSYS,Inc.ProprietaryAQWA AQWA 程序介绍程序介绍AQWA NAUT-Time domain program for wave frequency structure motion and mooring tension analyses in regular wavesAQWA DRIFT-Time domain program for drift frequency and wave frequency structure motion and mooring ten

6、sion 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 strength analysis 2006 ANSYS,Inc.All rights reserved.7ANSYS,In

7、c.ProprietaryAQWA 程序介绍LINELIBRIUMFERNAUTDRIFTWAVEASAS(FE model)General Relations between Programs 2006 ANSYS,Inc.All rights reserved.8ANSYS,Inc.ProprietaryTypical AQWA Model(1)-Moored Tanker 2006 ANSYS,Inc.All rights reserved.9ANSYS,Inc.ProprietaryTypical AQWA Model(2)-Semi Sub 2006 ANSYS,Inc.All ri

8、ghts reserved.10ANSYS,Inc.ProprietaryAQWA Result(1)-Force&Response Curves 2006 ANSYS,Inc.All rights reserved.11ANSYS,Inc.ProprietaryAQWA Result(2)-Pressure Contour 2006 ANSYS,Inc.All rights reserved.12ANSYS,Inc.ProprietaryAQWA Result(3)-Shear Force&Bending Moment 2006 ANSYS,Inc.All rights reserved.1

9、3ANSYS,Inc.ProprietaryAQWA LINE 介绍AQWA-LINE 是3D衍射和辐射分析程序通过大量的板格来描述结构在每个板格中心有源点（a source），程序会求解受到如下边界条件时源点的强度:没有水流穿过船体没有水流穿过海底板格为自由面这样就可以在任意点计算出压力和速度 2006 ANSYS,Inc.All rights reserved.14ANSYS,Inc.ProprietaryAQWA LINE理论AQWA LINE的3D线性波理论假设流体无旋，不可压小波（small wave elevation）速度势能控制方程体的边界条件 2006 ANSYS,Inc.A

10、ll rights reserved.15ANSYS,Inc.ProprietaryAQWA LINE理论自由表面运动条件 线性化之后自由表面动力条件 2006 ANSYS,Inc.All rights reserved.16ANSYS,Inc.ProprietaryAQWA LINE理论线性化之后组合自由表面条件 海底边界条件 2006 ANSYS,Inc.All rights reserved.17ANSYS,Inc.ProprietaryAQWA LINE理论辐射条件辐射条件是一种物理状态，使波不会在错误的方向上传递，如水穿透船体，从而避免数上的不确定性。数值方法势能项的线性叠加 式中：

11、下标I是入射波，d是衍射波,j=1,2,6 是6个自由度的辐射波，x是单位波幅下结构的运动.2006 ANSYS,Inc.All rights reserved.18ANSYS,Inc.ProprietaryAQWA LINE理论对于有限水深d，入射波的势能为式中k是波数，定义为:使用脉冲源分布（pulsating source distribution）求解衍射和辐射波势能 2006 ANSYS,Inc.All rights reserved.19ANSYS,Inc.ProprietaryAQWA LINE理论式中：是辐射体强度（source strength）;S 是入水结构面;(x,y,

12、z)为流域中域点的坐标;(,)是S上源点坐标;G 为格林函数，是求解Laplace方程的基础，并且满足所有边界条件（除了体边界条件）.G 可以表达为:2006 ANSYS,Inc.All rights reserved.20ANSYS,Inc.ProprietaryAQWA LINEAQWA LINE理论理论 式中 pv 表示积分的principal value;J0 是Bessel函数的第一项 2006 ANSYS,Inc.All rights reserved.21ANSYS,Inc.ProprietaryAQWA LINE理论 在每个板格结构表面上的源强度（source strength

13、）假设为常数，通过体边界条件求解积分方程计算:对于衍射势能，在结构表面由于入射势能产生的法向速度会变为零，而所减少的法向速度会转化为结构运动.2006 ANSYS,Inc.All rights reserved.22ANSYS,Inc.ProprietaryAQWA LINE理论压力和第一阶波浪力的计算 一旦计算了源点强度和势能，从线性化的Bernoulli方程中可以计算每个板格的水动力压力:第一阶波浪力是通过在体表面上进行积分得到的 2006 ANSYS,Inc.All rights reserved.23ANSYS,Inc.Proprietary二阶平均漂浮力计算远场求解（Far fiel

14、d solution）(动量守恒方法，只用于水平力）):式中 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理论 2006 ANSYS,Inc.All rights

15、reserved.24ANSYS,Inc.ProprietaryAQWA LINE理论近场求解(压力对于六个自由度进行积分方法):式中 WL 为水线位置;r 是相关波表面是相关波表面（relative wave surface elevation）;S0 水下结构表面;X 结构表面的运动;MS 结构的质量;R是结构的旋转矩阵;是结构中心加速度.2006 ANSYS,Inc.All rights reserved.25ANSYS,Inc.ProprietaryAQWA LINE理论Quadratic Transfer Function(QTF)Components at both differe

16、nce and sum frequenciesEach with real and imaginary parts 2006 ANSYS,Inc.All rights reserved.26ANSYS,Inc.ProprietaryAQWA LINE理论In which 2006 ANSYS,Inc.All rights reserved.27ANSYS,Inc.ProprietaryAQWA LINE理论运动方程 结构在波浪中的响应X是通过计算下面方程得到:式中 Ms 是结构质量,Ma 是附加质量,C 是阻尼,K 流体刚度,F 是波浪力（包括衍射和辐射力）.2006 ANSYS,Inc.Al

17、l rights reserved.28ANSYS,Inc.ProprietaryAQWA波概述波的类型: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阶波 2006 ANSYS,Inc.All rights reserved.29ANSYS,Inc.ProprietaryAQWA波概述波的类型(续):3

18、)流量函数波(只用于深水和Morison单元)a)载波表面压力为零；或者b)在波表面不运动。2006 ANSYS,Inc.All rights reserved.30ANSYS,Inc.ProprietaryAQWA波概述波浪形式:1)规则波用于AQWA LINE,NAUT(缺省)2)不规则波通过定义波谱确定，用于AQWA LIBRIUM,FER,DRIFT.不规则波也可以在NAUT中定义(with an IRRE card in the JOB deck).2006 ANSYS,Inc.All rights reserved.31ANSYS,Inc.ProprietaryAQWA波概述波浪形

19、式(续):在AQWA中的波谱类型:a.P-M spectrumb.JONSWAP spectrumc.User defined spectrumd.Gaussian spectrum for Cross Swell不规则波可以是一下形式:a.长波;或b.短波,也就是 a spread sea(只用于 AQWA LIBRIUM and FER)2006 ANSYS,Inc.All rights reserved.32ANSYS,Inc.Proprietary结构上的波力对于衍射结构(通过板单元来模拟)-入射波力(Froude-Krylov force):静止波压力产生（from the pres

20、sure in the undisturbed waves）.-衍射力:由于固定结构干扰入射波产生.-辐射力：由于结构摆动产生的波产生.-平均漂浮力(net force due to high order effect)2006 ANSYS,Inc.All rights reserved.33ANSYS,Inc.Proprietary结构上的波力对于Morison结构(用Morison单元模拟,eg TUBEs)-Morison 力(包括阻力，drag)是用Morison方程计算的。.2006 ANSYS,Inc.All rights reserved.34ANSYS,Inc.Propriet

21、aryMorison 力 Morison力的计算方程对于细长圆柱单元，例如管道单元，单元上单位长度的水动力可以用Morison方程计算:式中Ca 和 Cd是附加质量和单元的阻力（drag）系数；;是每单位长度单元的体积，D是单元的直径.2006 ANSYS,Inc.All rights reserved.35ANSYS,Inc.Proprietary结构上的波力流体力水动力水静力波浪中的固定结构周围压力（入射波或Froude-Krylov 力)结构对波的影响（衍射）静水中的运动结构（辐射力）In-phase(附加质量)Out-of-phase(辐射阻尼)F()K.xMa().xC().x 20

22、06 ANSYS,Inc.All rights reserved.36ANSYS,Inc.ProprietaryAQWA 的总体坐标系定向特性:原点在静水线上Z轴的正方向是垂直向上X轴的正方向是右舷Y轴的正方向是根据右手法则z0zyxW.L.2006 ANSYS,Inc.All rights reserved.37ANSYS,Inc.ProprietaryAQWA 的总体坐标系AQWA的总体坐标系是固定的参考轴，而与东、南、西、北等方向无关 2006 ANSYS,Inc.All rights reserved.38ANSYS,Inc.ProprietaryAQWA的方向AQWA中定义的波、风和

23、水流方向是作为船体的行进方向.在定义的波浪(或水流、风)之间有角度时，角度是与X正方向逆时针转为正.中的方向是安度输入输出.X axisWave direction(or current,wind)positive angle 2006 ANSYS,Inc.All rights reserved.39ANSYS,Inc.ProprietaryAQWA坐标系的方向对于非对称结构，输入的方向的角度必须在范围内-180 度到+180度内；对于X方向对称的结构，角度范围在0度到180度之间；对于X、Y方向对称结构，角度范围在0度到90度之间.yxvq 2006 ANSYS,Inc.All rights

24、 reserved.40ANSYS,Inc.Proprietary相位角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 the parameter reaches its peak value.(dt=*T/360,where T is the wave period).正的相位角表示滞后于波.2006 ANSYS,Inc.All rig

25、hts reserved.41ANSYS,Inc.ProprietaryAQWA 的文件名每个AQWA文件有三个部分:文件前缀（两个字母）-表示程序类型alLINEabLIBRIUMafFERadDRIFTanNAUT运行的文件(最多六个字母)2006 ANSYS,Inc.All rights reserved.42ANSYS,Inc.ProprietaryAQWA AQWA 的文件名的文件名文件扩展名(三个字母)文件扩展名说明文件类型(eg,.dat)例如:altest.dat 2006 ANSYS,Inc.All rights reserved.43ANSYS,Inc.Proprietar

26、y输入数据输入数据(.(.datdat)文件文件This is an ASCII text file containing all the input data 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 numb

27、er of the current cursor position,is highly recommended.Note:A graphical user interface,which allows interactive data input in a user-friendly environment,is currently under development.2006 ANSYS,Inc.All rights reserved.44ANSYS,Inc.Proprietary列表文件列表文件 (.(.lislis)This is an ASCII text file containin

28、g 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.2006 ANSYS,Inc.All rights reserved.45ANSYS,Inc.Proprietary重启动文件重启动文件(.(.resres)This is a binary file,which is written by all the AQWA programs,and which contains a summ

29、ary 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 analysis parameters 2006 ANSYS,Inc.All rights reserved.46ANSYS,Inc.Proprietary水动力文件水动力文件

30、(.(.hydhyd)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(RES)hydrodynamics fi

31、le(HYD)model definition hydrodynamic databasehydrodynamic database 2006 ANSYS,Inc.All rights reserved.47ANSYS,Inc.Proprietary水动力文件水动力文件 (.(.hydhyd)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 fi

32、le(only when ALDB option is on,or a FILE card with*.HYD file name in deck 6).2006 ANSYS,Inc.All rights reserved.48ANSYS,Inc.Proprietary定位文件定位文件Position(.pos and.Position(.pos and.eqpeqp)Both are binary files.eqp file:Created by AQWA LIBRIUM which stores the equilibrium positions of a system of struc

33、tures.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: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.2006 ANSYS,Inc.All ri

34、ghts reserved.49ANSYS,Inc.ProprietaryAGS 显示文件(.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,DRIFT)orforces and responses as a function of frequency(FER,LINE)2006 ANSYS,Inc.All rights reserved.50ANSYS,Inc.Pro

35、prietaryAQWA AQWA 文件类型文件类型每次运行AQWA都会包括多个文件。文件主要是后缀不同，一般的扩展名如下:.dat输入数据文件(ASCII).lis输出列表文件(ASCII).res重启动文件(binary),输入、输出.hyd水动力文件(binary),输入、输出.pos输出定位（position）文件(binary).plt输出图形文件(binary).eqp平衡位置文件(binary),输入、输出.uss源强度文件(binary),输入、输出 2006 ANSYS,Inc.All rights reserved.51ANSYS,Inc.ProprietaryAQWA A

36、QWA 文件类型文件类型.mes 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 g

37、enerator(ASCII).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)2006 ANSYS,Inc.All rights reserved.52ANSYS,Inc.ProprietaryAQWA Restart Stage

38、sStage 1定义模型 Decks 1 to 5Stage 2水动力数据 Decks 6 to 8Stage 3衍射/辐射分*Stage 4主要分析参数 Decks 9 to 20Stage 5主要分析*只是计算步骤 2006 ANSYS,Inc.All rights reserved.53ANSYS,Inc.ProprietaryStage 1Deck 1COOR节点坐标Deck 2ELM*单元定义Deck 3MATE材料特性Deck 4GEOM几何特性Deck 5GLOB整体约束*结构编号 2006 ANSYS,Inc.All rights reserved.54ANSYS,Inc.Pr

39、oprietaryStage 2Deck 6 FDR*规则波定义(频率和方向)Deck 7 WFS*水动力特性静水特性Analysis PositionDeck 8 DRC*漂流力系数*结构编号注:在第二步输入的水动力特性是用修改或替代第三步AQWA-LINE 的计算结果 2006 ANSYS,Inc.All rights reserved.55ANSYS,Inc.ProprietaryStage 3该步是AQWA-LINE的主要分析，但只是一个程序计算过程而已.注:在第二步输入的水动力特性是用修改或替代第三步AQWA-LINE 的计算结果 2006 ANSYS,Inc.All rights

40、reserved.56ANSYS,Inc.ProprietaryStage 4Deck 9DRM*漂流运动参数Deck 10HLD*船体阻力（Drag）系数Deck 11ENVR环境参数Deck 12CONS约束Deck 13SPEC谱参数WAVE规则波参数(N)Deck 14MOOR停泊线定义*Structure Number 2006 ANSYS,Inc.All rights reserved.57ANSYS,Inc.ProprietaryStage 4(cont)Deck 15 STRT初始条件Deck 16 TINT 时间积分参数(D,N)LMTS迭代参数(B)Deck 17 HYDC

41、管道的水动力参数Deck 18 PROP打印输出选项Deck 19 NONEReserved for future useDeck 20 NONEReserved for future use 2006 ANSYS,Inc.All rights reserved.58ANSYS,Inc.ProprietaryStage 5该步骤只是程序求解计算.2006 ANSYS,Inc.All rights reserved.59ANSYS,Inc.Proprietary分析方法一般的分析方法是在AQWA-LINE中运行Stages 1 to 3，然后在其它程序运行 Stages 4 to 5，例如 AQ

42、WA-NAUT.实例:Step 1:AQWA-LINE(Restart 1 to 3)runInput FilesOutput Filesaltest.dataltest.lisaltest.resaltest.hydaltest.plt 2006 ANSYS,Inc.All rights reserved.60ANSYS,Inc.Proprietary分析方法(续)Step 2:文件拷贝copy altest.resabtest.resStep 3:AQWA-LIBRIUM(Restart 4 to 5)runInput FilesOutput Filesabtest.databtest.l

43、istabtest.resabtest.resabtest.eqp 2006 ANSYS,Inc.All rights reserved.61ANSYS,Inc.ProprietaryAQWA 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 sectio

44、n allowed)PMAS:Point mass 2006 ANSYS,Inc.All rights reserved.62ANSYS,Inc.ProprietaryAQWA 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 diffractin

45、g QPPL and TPPL elements.2006 ANSYS,Inc.All rights reserved.63ANSYS,Inc.Proprietary AQWA LINE警告检查General warningsNo.of elements 6000:solution timeNormals point out:modelling conventionNo gaps:force balanceFacets cannot cut surface:solution requirementDimensions KR:good practice 2006 ANSYS,Inc.All ri

46、ghts reserved.64ANSYS,Inc.ProprietaryAQWA LINE警告检查Stage#1 checksArea ratio of adjacent elements Element centres at least one facet radius apart Shape factor 0.2 warningShape factor 0.02 fatal error 2006 ANSYS,Inc.All rights reserved.65ANSYS,Inc.ProprietaryWarnings in AQWA LINEStage#2 checksLongest s

47、ide 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.2006 ANSYS,Inc.All rights reserved.66ANSYS,Inc.ProprietaryAQWA-LINE 的实例输入数据 2006 ANSYS,Inc.All right

48、s reserved.67ANSYS,Inc.ProprietaryAQWA-LINE 实例输入数据(续)2006 ANSYS,Inc.All rights reserved.68ANSYS,Inc.Proprietary定义Morison单元 2006 ANSYS,Inc.All rights reserved.69ANSYS,Inc.ProprietaryAQWA LINE的常用选项(1)DATA用于检查输入数据(LBDFN)GOON用于忽略一般警告(L)REST定义开始过程(LBDFN)LDOP载荷输出-输出.POT 和.USS 文件 用于压力计算(例如压力显示,SF/BM)(L)PRC

49、E为 Decks 1 5(LBDFN)设置打印项PPEL为每个单元设置打印属性(LBDFN)2006 ANSYS,Inc.All rights reserved.70ANSYS,Inc.ProprietaryAQWA LINE的常用选项(2)NPPP无压力后后处理阻止节点连接警告ALDBAqwa-Line 数据库-reads in AQWA-LINE database TAKE CARENQTF使用近场求解，用于漂浮力系数CQTF计算 QTF 矩阵PFIXPartially FIXed 2006 ANSYS,Inc.All rights reserved.71ANSYS,Inc.Proprie

50、taryAQWALINE AQWALINE 后处理后处理波表面云图显示波表面云图显示 2006 ANSYS,Inc.All rights reserved.72ANSYS,Inc.ProprietaryAQWALINE AQWALINE 后处理后处理Air Gap 2006 ANSYS,Inc.All rights reserved.73ANSYS,Inc.ProprietaryAQWALINEAQWALINE后处理后处理压力云图压力云图 2006 ANSYS,Inc.All rights reserved.74ANSYS,Inc.ProprietaryAQWALINEAQWALINE后处理后

51、处理衍射波表面衍射波表面 2006 ANSYS,Inc.All rights reserved.75ANSYS,Inc.ProprietaryAQWA 数据处理处理原因处理原因:为了添加附加节点、单元、阻尼等，或者把多个数据库合并成一个，而不重新完全运行AQWA-LINE分析.处理过程处理过程:1)拷贝已有数据文件AL1.DAT，改名为 AL2.DAT;2)添加新节点，（非衍射）单元，附加阻尼 等(CRNM 选项需要用 RAO 重新计算)到AL2.DAT文件。2006 ANSYS,Inc.All rights reserved.76ANSYS,Inc.ProprietaryAQWA 数据处理(

52、续.1)3)删除deck 6中所有波频率和方向并改为:06 FDR106FILE AL1.HYD06CSTR N (N is the structure number in AL1)END06CPDB 4)在 AQWA-LINE中运行AL2.DAT 2006 ANSYS,Inc.All rights reserved.77ANSYS,Inc.ProprietaryAQWA 数据处理(续.2)For对于复杂结构数据和合并，新的AL*.DAT文件应该包括所有的结构定义，并且每个结构都有一个deck 6，包括文件卡片中对于的文件名.以上过程也可以应用于其它AQWA程序.2006 ANSYS,Inc.

53、All rights reserved.78ANSYS,Inc.ProprietaryAQWA 数据处理(续.3)JOB MESH NAUTTITLE A SHIP+TUBES AND DISCSOPTIONS REST ENDRESTART 1 5 01 COOR 010001 10.0 0.000 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)

54、(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 1.00 END04PMAS 2 1.0E9 0.00 0.00 9.0E9 0.0 9.0E9 2006 ANSYS,Inc.All rights reserved.79ANSYS,Inc.ProprietaryAQWA 数据处理(续.4)2006 ANSYS,Inc.All rig

55、hts reserved.80ANSYS,Inc.Proprietary相互作用的结构相互作用的结构(1)(1)PFIX方法把漂浮模型和固定模型合并到一个结构中利用特殊的组编号999使用PFIX选项设置固定结构的基础 2006 ANSYS,Inc.All rights reserved.81ANSYS,Inc.Proprietary相互作用的结构相互作用的结构(2)(2)相互作用水动力AQWA能够计算水动力系数，可以完全的考虑水动力之间的相互作用.可以包括达到10浮体结构相互作用.2006 ANSYS,Inc.All rights reserved.82ANSYS,Inc.Proprietar

56、y结构间的水动力响应作用 2006 ANSYS,Inc.All rights reserved.83ANSYS,Inc.Proprietary对称性对称性SYMX 表示AQWA 可以假设结构是X轴（FRA）对称的。利用对称性可以减少计算时间.RMXS不用对称（remove symmtry），创建完整模型（即使模型是对称的）。该方法只用于T/QPPL单元.MSTR把结构移动到新定义的位置。只用于单元和与单元关联的节点，而不是所有在STRC卡片中列出的节点.实际上在FRA中移动节点和单元，而与Deck 15卡片中的位置不同.SYMY 和RMYS 对于Y轴有相同的影响 2006 ANSYS,Inc.

57、All rights reserved.84ANSYS,Inc.ProprietaryAQWA LIBRIUM 介绍AQWA-LIBRIUM 用于确定平衡位置，静稳定性和动稳定性结构是以很小的步骤移动直到力平衡，从而找到平衡位置然后，程序计算线性化刚度矩阵的特征值，得到静稳定性通过阻抗矩阵的特征值得到动稳定性 2006 ANSYS,Inc.All rights reserved.85ANSYS,Inc.ProprietaryAQWA LIBRIUM理论确定静平衡位置的方程.结构的X位置通过下式计算:式中K是系统的刚度矩阵，F是力.程序迭代直到X小于所定义的容差 2006 ANSYS,Inc.A

58、ll rights reserved.86ANSYS,Inc.Proprietary分析方法Step 2:File Copyingcopy altest.resabtest.resStep 3:AQWA-LIBRIUM(Restart 4 to 5)runInput FilesOutput Filesabtest.databtest.listabtest.resabtest.resabtest.eqp 2006 ANSYS,Inc.All rights reserved.87ANSYS,Inc.Proprietary分析方法(续)在AQWA-LIBRIUM主分析开始时(Stage 5)，需要提

59、供如下数据:主模型定义restart file水动力数据 restart file主分析参数 input data fileNote:As the data from the restart 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).2006 ANSYS,Inc.All rights reserved.88ANSYS,Inc.ProprietaryAQWA-LIBRIUM 实例数据(1)2006 ANS

60、YS,Inc.All rights reserved.89ANSYS,Inc.ProprietaryAQWA-LIBRIUM实例数据(1)(cont.)2006 ANSYS,Inc.All rights reserved.90ANSYS,Inc.ProprietaryAQWA-LIBRIUM实例数据(2)2006 ANSYS,Inc.All rights reserved.91ANSYS,Inc.ProprietaryAQWA-LIBRIUM实例数据(2)(cont.)2006 ANSYS,Inc.All rights reserved.92ANSYS,Inc.ProprietaryAQWA

61、LIBRIUM常用选项STAT只用于静稳定性分析DYNA只用于动力稳定性分析LSTF线性刚度（Linear STiFness），用从AQWA-LINE得到的水静力刚度.PBIS打印迭代步 2006 ANSYS,Inc.All rights reserved.93ANSYS,Inc.ProprietaryDrag in AQWA(deck 10)Drag in AQWA(deck 10)1)水流和风力系数水流和风力系数CUFX Dir1 DirN C1.CNWIFX Dir1 DirN C1.CN etc(Dir:方向编号方向编号;C1 :阻力（阻力（Drag）系数）系数在在 方向的相对水流速度

62、方向的相对水流速度 force in X direction=CUFX.V2force in Y direction=CUFY.V2yaw moment =CURZ.V2 2006 ANSYS,Inc.All rights reserved.94ANSYS,Inc.ProprietaryDrag in AQWA(deck 10)Drag in AQWA(deck 10)(cont 1)(cont 1)2)Morison阻力（阻力（drag）系数）系数 MDIN Nrow Ncol C1.C6Nrow:阻力矩阵的行数阻力矩阵的行数Ncol:阻力矩阵的列数阻力矩阵的列数C1 C6:阻力（阻力（Dr

63、ag）系数）系数 2006 ANSYS,Inc.All rights reserved.95ANSYS,Inc.ProprietaryAQWAAQWA约束约束 (deck 12)(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 con

64、straint.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.Nt=3:Rigid connection,no rotation.2006 ANSYS,Inc.All rights reserved.96ANSYS,Inc.ProprietaryAQWA约束(deck 12)(cont.1)2006 ANSYS,Inc.All rights reserved.97ANSYS,Inc.ProprietaryAQWA约束

65、(deck 12)(cont.2)2006 ANSYS,Inc.All rights reserved.98ANSYS,Inc.ProprietaryAQWA约束(deck 12)(cont.3)2006 ANSYS,Inc.All rights reserved.99ANSYS,Inc.ProprietaryAQWAAQWA停泊索停泊索(deck 14)(deck 14)AQWA 程序除了程序除了AQWA LINE 外，都可以定义停泊索。下面是常用的停泊外，都可以定义停泊索。下面是常用的停泊类型类型(更详细的内容参考更详细的内容参考AQWA手册手册):LINE:线弹性线（无重量）线弹性线（无

66、重量）LINE Ns1 Nd1 Ns2 Nd2 K L(Ns1,Ns2:structure numbers;Nd1,Nd2:node numbers;K:stiffness;L:unstretched length)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 is used as a winch)2006 ANSYS,Inc.All rights reserved.100ANSYS,Inc.ProprietaryA