DeepC 实用教程（三）环境数据

一、前言

  SESAM （Super Element Structure Analysis Module）是由挪威船级社（DNV-GL）开发的一款有限元分析（FEA）系统，主要用于海洋工程结构的极限强度评估、波浪荷载计算和系泊系统分析等。它由一系列模块/程序组成，以 GeniE、HydroD 和 DeepC 等模块为核心，是海工结构分析领域的标准软件。GeniE、HydroD 和 DeepC 分别是 SESAM 系统的建模及前处理模块、水动力计算模块和深水耦合分析模块。

SESAM 2017 Overview

  DeepC 是 Sesam 系统中的深水耦合分析工具，由 DNV（挪威船级社）和 Marintek（挪威海洋技术研究所）联合开发，主要用于 SPARs、FPSOs、TLPs 或类似类型的含有系泊/立管系统的浮式结构的分析。DeepC 包含三个独立的程序模块：DeepC Concept Modeller、Riflex 和 Simo，它们形成了一个建模、分析与后处理的集成环境。DeepC Concept Modeller 用于创建模型输入、结果后处理等。Riflex 和 Simo 由 Marintek 拥有、开发和维护，它们是 DeepC 的两个分析引擎（求解器），用于执行非线性时域有限元分析。

二、风谱/风剖

定义风谱与风剖面

// -----------------------------------------------------------
// Wind Spectrums
NPD1 = NPD(0.002);
// -----------------------------------------------------------
//Wind Profiles
//                       Href      Exp    V
windProf1 = WindProfile(37.49 m/s, 10 m, 0.11);
windProf2 = WindProfile(38.1 m/s, 10 m, 0.11);
windProf3 = WindProfile(23.29 m/s, 10 m, 0.11);
// -----------------------------------------------------------

三、洋流

定义洋流

// -----------------------------------------------------------
// Current Profiles
CurProf1_elevations = Array(-639.13 m,-80.5 m,-80.46 m,0 m);
CurProf1_angles = Array(0 deg,0 deg,0 deg,0 deg);
CurProf1_velocities = Array(0.122 m/s,0.122 m/s,0.945 m/s,0.945 m/s);
CurProf1 = CurrentProfile(CurProf1_elevations, CurProf1_angles, CurProf1_velocities);
CurProf2_elevations = Array(-639.13 m,-79.9 m,-79.85 m,0 m);
CurProf2_angles = Array(0 deg,0 deg,0 deg,0 deg);
CurProf2_velocities = Array(0.122 m/s,0.122 m/s,0.945 m/s,0.945 m/s);
CurProf2 = CurrentProfile(CurProf2_elevations, CurProf2_angles, CurProf2_velocities);
CurProf3_elevations = Array(-638.13 m,-624.8 m,-500 m,-400 m,-300 m,-200 m,-100 m,-50 m,0 m);
CurProf3_angles = Array(0 deg,0 deg,0 deg,0 deg,0 deg,0 deg,0 deg,0 deg,342 deg);
CurProf3_velocities = Array(0.232 m/s,0.232 m/s,0.241 m/s,0.25 m/s,0.329 m/s,0.658 m/s,1.161 m/s,1.439 m/s,1.713 m/s);
CurProf3 = CurrentProfile(CurProf3_elevations, CurProf3_angles, CurProf3_velocities);
// -----------------------------------------------------------

四、波浪

4.1 规则波浪

定义规则波

// -----------------------------------------------------------
// Regular Wave
RE_7x5_7x5 = RegularWaveHT(7.5 m,WavePeriod(7.5 s),rwHeight);
RE_7x5_15x5 = RegularWaveHT(7.5 m,WavePeriod(15.5 s),rwHeight);
RE_10x5_7 = RegularWaveHT(10.5 m,WavePeriod(7 s),rwHeight);
RE_10x5_15 = RegularWaveHT(10.5 m,WavePeriod(15 s),rwHeight);
RE_10x5_7x5 = RegularWaveHT(10.5 m,WavePeriod(7.5 s),rwHeight);
RE_10x5_15x5 = RegularWaveHT(10.5 m,WavePeriod(15.5 s),rwHeight);
RE_14x5_9 = RegularWaveHT(14.5 m,WavePeriod(9 s),rwHeight);
RE_14x5_15x5 = RegularWaveHT(14.5 m,WavePeriod(15.5 s),rwHeight);
RE_16x5_9 = RegularWaveHT(16.5 m,WavePeriod(9 s),rwHeight);
RE_17x5_9x5 = RegularWaveHT(17.5 m,WavePeriod(9.5 s),rwHeight);
RE_17x5_15x5 = RegularWaveHT(17.5 m,WavePeriod(15.5 s),rwHeight);
RE_17x5_15 = RegularWaveHT(17.5 m,WavePeriod(15 s),rwHeight);
// -----------------------------------------------------------

4.2 随机波浪谱

定义波浪谱

// -----------------------------------------------------------
// Wave Spectrums
//              Hs       Tp      Gamma
Jons1 = Jonswap(11.89 m, 13.8 s, 2.2);
Jons2 = Jonswap(11.8 m, 13.5 s, 2);
Jons3 = Jonswap(7.1 m, 10.8 s, 1.8);
Jonswap_7x5_11x5 = Jonswap(7.5 m,11.5 s,2.511581877);
Jonswap_7x5_19x5 = Jonswap(7.5 m,19.5 s,1);
Jonswap_17x5_11x5 = Jonswap(17.5 m,11.5 s,5);
Jonswap_17x5_19x5 = Jonswap(17.5 m,19.5 s,1.476094409);
// -----------------------------------------------------------

五、方向

‘

定义风浪流的主方向

// -----------------------------------------------------------
// Directions
Dircur1 = Direction(190 deg);
Dircur2 = Direction(219 deg);
Dircur3 = Direction(90 deg);
Dirwave1 = Direction(210 deg);
Dirwave2 = Direction(240 deg);
Dirwave3 = Direction(61 deg);
Dirwind1 = Direction(227 deg);
Dirwind2 = Direction(260 deg);
Dirwind3 = Direction(61 deg);
// -----------------------------------------------------------

六、海床属性

定义海床属性

// -----------------------------------------------------------
// Seabed
//                      Kn       Kl    Kt   Nyl Nyt
floor = SeaBedProperty(3550 Pa, 0 Pa, 0 Pa,  0,  0);
// -----------------------------------------------------------

七、位置

7.1 创建位置

Create New Loaction

// -----------------------------------------------------------
// Seabed
//                      Kn       Kl    Kt   Nyl Nyt
floor = SeaBedProperty(3550 Pa, 0 Pa, 0 Pa,  0,  0);
// -----------------------------------------------------------

7.2 规则波时域条件

7.3 随机波时域条件

(Single) Irregular Time Condition

7.4 波浪散布图

北大西洋波浪散布图（Sactter Diagram）

// -----------------------------------------------------------
// Scatter Diagrams
// A scatter diagram gives the number of occurences of seastates wihtin given wave height and wave period ranges
start_hs_array = Array(0.5 m,1.5 m,2.5 m,3.5 m,4.5 m,5.5 m,6.5 m,7.5 m,8.5 m,9.5 m,10.5 m,11.5 m,12.5 m,13.5 m,14.5 m,15.5 m,16.5 m,17.5 m);
start_tp_array = Array(1.5 s,2.5 s,3.5 s,4.5 s,5.5 s,6.5 s,7.5 s,8.5 s,9.5 s,10.5 s,11.5 s,12.5 s,13.5 s,14.5 s,15.5 s,16.5 s,17.5 s,18.5 s,19.5 s);
Block_A = ScatterHsTp(EDDY_location, start_hs_array, start_tp_array);
Block_A.setNumberOfOccurences(0.5 m, 1.5 s, 0);
Block_A.setNumberOfOccurences(0.5 m, 2.5 s, 0);
Block_A.setNumberOfOccurences(0.5 m, 3.5 s, 13);
Block_A.setNumberOfOccurences(0.5 m, 4.5 s, 1337);
Block_A.setNumberOfOccurences(0.5 m, 5.5 s, 8656);
// 代码未全部列出/以上代码仅作示例
// -----------------------------------------------------------

7.4.1 散布图分块

Hs/Tp Diagram

// -----------------------------------------------------------
// Scatter Discretizations
// The purpose of a scatter discretization is to subdivide the scatter diagram into entities that 
// can be analysed with a single time domain analysis.
ScatterDiscretization1 = ScatterDiscretizationHsTp(Block_A);
ScatterDiscretization1.addBlock("B7x5_11x5",7.5 m,11.5 s);
ScatterDiscretization1.addBin("B7x5_11x5","Bin1",1);
ScatterDiscretization1.addBlockRange("B7x5_11x5",0.5 m,3.5 s,1.5 m,4.5 s);
ScatterDiscretization1.addBlockRange("B7x5_11x5",0.5 m,4.5 s,4.5 m,5.5 s);
ScatterDiscretization1.addBlockRange("B7x5_11x5",0.5 m,5.5 s,7.5 m,11.5 s);
ScatterDiscretization1.addBlock("B7x5_19x5",7.5 m,19.5 s);
// 代码未全部列出/以上代码仅作示例
// -----------------------------------------------------------

7.4.2 时域条件

Multiple Irregular Time Conditions

Multiple Regular Time Conditions

八、参考文献

[1]. Sesam User Manual - DeepC (Deep water coupled floater motion analysis).

[2]. Sesam Theory Manual - DeepC (Deep water coupled floater motion analysis).

[3]. DeepC – Deep Water Coupled Analysis Tool, A White Paper.

[4]. 关于风机 叶片/荷载/控制 方面的介绍请访问：https://www.zhihu.com/column/c_1485646874003058688

[5]. 关于风机 有限元分析 方面的介绍请访问：https://blog.csdn.net/shengyutou

[6]. 联系作者 ，Email： [email protected]，WeChat/Weixin： 761358045