板壳结构非线性分析的数值方法

根据连续介质力学和非线性有限元理论，推导了进行板壳结构非线性分析的退化壳元的完全拉格朗日（简称Ｔ．Ｌ．）法列式，建立了材料在大变形、大应变情况的本构关系。数值算例证明了本文方法的正确性和有效性。     

二维浮式结构物与完全非线性波相互作用数值分析

基于粘流理论、流固耦合理论和自由液面追踪技术——VOF(Volume of Fluid)法,以不可压缩流体的连续方程和N-S方程及结构动力学方程为基本控制方程,利用商用软件ANSYS Workbench的二次开发功能,建立了粘性数值波浪水池,数值模拟了完全非线性波的生成,计算分析了二维浮式结构物遭遇波浪的完全非线性现象,得到了甲板上浪水位高度和结构物的纵摇、垂荡、横摇等运动响应。结果表明:文中的方法可以用于波浪与浮式结构物相互作用的完全非线性数值模拟.     

解非线性船波问题的一个数值计算方法

本文提出了一种计算完全非线性三维船渡问题的数值计算方法。速度势用考虑了波动的自由面影响的船体湿表面和局部真实的自由面上的简单源(1/r)分布表示,并采用精确的自由面条件。求解过程是一个迭代过程,计算结果与试验值作了比较,吻合较好.     

港湾内非线性波浪对船剖面作用的时域计算模型(英文)

文中用边界元数值求解复速度势模拟了波浪对船剖面的时域作用力,采用满足运动学和动力学条件的半拉朗格朗日法跟踪流体的自由表面.推导了适用本数值模型的辐射边界条件,此边界条件和海绵层边界条件组成的出口边界条件,可以使波浪很好地被吸收或传播出去.进行了相关物理模型实验,数值计算结果和实验结果的对比表明此模型能精确模拟波浪对船的作用力.     

完全非线性波的数值模拟

本文利用MAC差分格式，对欧拉动学方程进行离散求解，利用VOF法处理自由表面的演变。数值计算主要针对矩形围域内驻波随时间的变化，给出了不同时刻的自由表面波形图。计算结果同通常势流理论的线性解作了对比，发现本法可以用于揭示自由波面运动的非线性特征，且可以作为研究完全非线性波的一种方法。     

基于CFD-非线性有限元双向耦合的集装箱船波浪下结构崩溃数值仿真

[目的]充分考虑载荷非线性和物面非线性因素的影响,研究集装箱船在波浪下结构响应及动态结构崩溃模式。[方法]首先,基于CFD平台建立船舶水动力模型,采用重叠网格法实现船体水动力模型动边界网格与远场流体域的欧拉网格间的匹配,在流体全域内采用流体体积法模拟自由面非线性,在流场全域内求解三维N–S方程,实时求解非线性波浪载荷;然后,建立可模拟船舯崩溃行为的船舶非线性有限元模型,基于显式动力学非线性有限元法计算包含塑性和屈曲的时域崩溃响应;最后,实现水动力模型与结构有限元模型在湿表面上的流体压力和节点位移的传递,以此进行CFD求解器与非线性有限元求解器间的双向迭代耦合,并实时计算4 600 TEU集装箱船结构崩溃过程中的非线性波浪载荷和结构崩溃响应。[结果]结果显示,在极端波浪下上部至中部结构广泛进入塑性状态,主甲板、舷侧板、甲板纵骨和舷侧纵骨等构件在波浪下出现明显的屈服及失稳,甲板纵骨和舷侧纵骨等骨材发生严重的侧向失稳,船体结构丧失了承载能力。[结论]所提方法可较准确求解结构响应及动态崩溃模式,可作为研究船舶结构崩溃响应的一种新方法。     

非线性波与可渗潜堤的相互作用数值模拟

把波浪场与可渗潜堤内部的孔隙流场作为一个整体问题,从Navier-Stokes方程出发,采用k-ε模型来封闭雷诺方程,作为波浪场的控制方程,以Forchheimer方程作为孔隙流的控制方程,用VOF方法跟踪自由表面,应用施主与受主单元模型求解流体体积函数F的输运方程。以椭圆余弦波为入射波,数值模拟了椭圆余弦波在可渗潜堤前后的波面变化过程,并得到了潜堤内部孔隙流压力的变化结果。     

FPSO水动力特性的完全非线性数值模拟

  目的  为提高海洋结构物的安全性能,针对波浪与结构物相互作用的问题开展完全非线性数值模拟研究。  方法  基于三维完全非线性时域势流理论及高阶边界元法（HOBEM）,建立波浪与结构物相互作用的开敞水域模型。采用速度势分离技术将整个问题分解为入射部分和散射部分,入射势由理论解给定。采用混合欧拉—拉格朗日（MEL）方法追踪瞬时自由水面的流体质点,并采用四阶龙格—库塔法对瞬时自由水面进行更新。引进虚拟函数计算波浪载荷,而非直接求解速度势的时间导数。在自由水面的外侧设置人工阻尼层,防止波浪从远场边界反射。自由水面网格仅在初始时刻生成一次,并采用弹簧近似法在不改变网格节点顺序的情况下对瞬时水面进行网格重构,以避免数值不稳定。  结果  在验证所提出数值模型有效性和精确性的基础上,针对某浮式生产储卸油轮（FPSO）模型的水动力特性进行数值模拟,发现考虑非线性影响时FPSO的运动响应在共振区段明显增大,证明了传统线性方法的预报结果趋于危险。  结论  研究成果既可为海洋浮式结构物的设计提供更可靠的预报工具,也可为其实际应用提供理论依据。     

超大型浮体结构考虑非线性项的波浪力响应计算

为求解超大型浮体结构波浪响应问题,将其简化成弹性梁板模型。以有限元软件为基础来建立超大型浮体结构的弹性梁板模型,在势流理论的基础上,将速度势分解为入射势、绕射势和辐射势三部分,并利用格林函数法推导出一阶、二阶波浪速度势的表达式。以此为理论基础,编制了相应的计算程序来计算浮体结构所受到的二阶波浪力,最后建立了梁板模型的水弹性响应方程,对非线性波浪力作用下的浮体结构的水弹性响应进行求解分析,并分析了波向对响应振幅的影响。     

二阶全非线性Boussinesq方程的二维数值模拟研究

在非交错网格下采用有限差分法首次对一组非线性精确至()Ομ2阶的全非线性Boussinesq方程数学模型进行了二维数值模拟分析。首先通过在方程的非线性项中引入缓坡假定,考察了其对模型数值精度的影响;其次,在模型中对二阶非线性项采用不同精度,考察了其对模型数值结果的影响。数值模拟结果表明,所建立的二阶完全非线性Boussinesq方程二维数值模型具有良好的适用性,模型非线性项中引入缓坡假定以及在二阶非线性项选用不同的精度对数值模拟结果影响不明显。     

船体湿表面脉动压力预报的三维非线性方法

本文给出了一种预报船体湿表面脉动压力响应的三维非线性时历模拟方法.计及瞬时物面变化引起的各种非线性因素,在时域内建立船体非线性运动方程.在每一时刻,认为船体以其瞬时平均湿表面为平衡位置做简谐振荡运动,利用三维线性频域解来简化运算.在此基础上,进一步将绕射和辐射问题在静平衡位置求解,仅考虑占主导地位的入射波和静恢复力变化的非线性,作为一种兼顾工程需要的计算精度和计算效率的简化方法.以此方法为基础编制了计算机程序,对实船进行了计算,并与模型试验进行了对比.分析比较表明,该方法能够给出令人满意的船体湿表面脉动压力响应预报结果.     

基于 mN LS方程的周期波群非线性演化数值模拟研究

论文基于控制深水复波包络演化的修正的四阶非线性薛定谔方程（mNLS）建立模型,利用离散时间步长的虚拟频谱方法（split‐step pseudospectral method）进行求解,模拟了周期波群的非线性演化,对演化过程中复包络幅值、波面位移、谱成分能量以及模态值的变化进行分析。结果表明：该模型可以有效地模拟周期波群的非线性演化。     

基于ANSYS的船闸混凝土非线性仿真分析

以长洲水利枢纽1#船闸上闸首为例,采用有限元软件,对施工期闸首混凝土浇筑模拟计算,考虑水化热、外界气温、浇筑层厚度之间的相互影响,对闸首混凝土温度变化进行分析。并以温度场的计算成果为前提,采用顺序耦合法对闸首结构混凝土的浇筑进行热-结构耦合计算,分析施工期混凝土结构的应力变化情况。通过对计算成果合理性的验证表明,利用有限元法可较好地对船闸结构进行整体稳定、应力、应变等全方位计算,且具有精度高、形象直观等特点,更好的模拟施工条件对船闸裂缝的影响。     

近水面椭球体非线性兴波的数值模拟

本文开发了一种去奇异面元法用于计算以定常前进速度运动的近水面椭球体引起的三维非线性自由面势流.采用求解非线性深水波的Stokes理论求解该定常非线性问题;物面用NURBS(非均匀有理B样条)表达;孤立的Rankine源分布在物体内部和自由面上方.给出了波型、波阻以及收敛特性随一些参数的变化规律.波阻计算结果和他人的其它数值结果进行了比较,其吻合程度令人满意.     

基于三维势流理论的耐波性数值预报

准确预报船体运动响应对于砰击等波浪载荷的计算以及合理结构设计具有重要意义。船舶在大幅波浪中的运动呈明显的非线性,而现阶段耐波性预报多采用线性切片方法。三维水动力分析软件WASIM基于时域势流理论,采用Rankine面元法预报船舶在波浪中的运动响应,并考虑了多种非线性因素。本文以标模DTMB5512为对象,采用WASIM预报其在不同航速下的耐波性,并与基于线性切片理论的计算结果和模型试验结果进行对比。结果表明:利用WASIM计算得到的船体运动响应比其他方法更接近试验值,合理体现了船舶在风浪中的实际耐波性能。因此,利用WASIM能够较好地评估船舶在波浪中的非线性耐波特性。     

纵摇容器中液体晃荡的非线性数值模拟

在海洋工程领域,液体晃荡是一种普遍存在的物理现象。对于船舶而言,转动比平动有着更重要的影响。该文针对纵摇容器中的液体晃荡问题,采用高阶边界元方法建立自由水面满足完全非线性边界条件的时域数学模型。通过大地坐标系和随体坐标系之间的坐标变换,使得计算域仅控制在随体坐标系内。求解中采用半混合欧拉—拉格朗日方法追踪流体瞬时水面,运用四阶龙格库塔方法更新下一时间步的波面和速度势。通过与已发表试验和数值结果的对比,验证了建立模型的准确性。进而开展大量数值试验研究容器纵摇运动频率、纵摇转动中心和容器中布置一竖直隔板对晃动波面与荷载的影响。     

采用对数应变的附加混合应变元的非线性分析

将Ｓｉｍｏ等提出的混合应变元应用到采用对数应变张量的有限应变分析中，给出了基于对数应变张量的几何非线性附加混合应变元的弹塑性模型，得出了简便的计算方法。数值算例表明，该方法能得到较好的结果，有较高的精度。     

基于三维势流理论的耐波性数值预报

准确预报船体运动响应对于砰击等波浪载荷的计算以及合理结构设计具有重要意义。船舶在大幅波浪中的运动呈明显的非线性,而现阶段耐波性预报多采用线性切片方法。三维水动力分析软件 WASIM基于时域势流理论,采用 Rankine面元法预报船舶在波浪中的运动响应,并考虑了多种非线性因素。本文以标模 DTMB5512为对象,采用 WASIM预报其在不同航速下的耐波性,并与基于线性切片理论的计算结果和模型试验结果进行对比。结果表明：利用 WASIM计算得到的船体运动响应比其他方法更接近试验值,合理体现了船舶在风浪中的实际耐波性能。因此,利用 WASIM能够较好地评估船舶在波浪中的非线性耐波特性。     

基于非线性有限元法的船舶冰区破冰数值模拟

应用非线性有限元法进行了破冰船冰区破冰数值模拟。通过比较数值模拟结果和试验结果,对冰体材料模型进行了验证;采用该冰体材料模型,对破冰船以不同航速在不同厚度的层冰中破冰航行时的动态响应进行了数值研究,给出了破冰过程中层冰的变形、冰力的大小以及冰的变形能和动能变化,分析了船速、冰层厚度对破冰阻力的影响。该研究结果对分析破冰船在层冰中破冰时的动态响应特性具有一定的参考价值。