Numerical Investigation of the Dynamics for Low Tension Marine Cables

This paper presents a numerical investigation into the dynamics of marine cables which are extensively used in offshore industry. In this numerical study, the Euler-Bernoulli beam model is adopted to develop the governing equations of the cable. Bending stiffness is considered to cope with the low tension problem in local area of towing cable, and thus a more accurate solution with the consideration of the axial elongation can be given.The derived strongly-coupled and nonlinear governing equations are solved by a second-order accurate, implicit,and large time step stable central finite difference method. The quadratically convergent Newton-Raphson iteration method is applied to solving the discrete nonlinear algebraic equations. Then a towed array sonar system(TASS)problem is studied. The numerical solutions agree reasonably well with the experimental data and the simulated results of the references. The specified program of the present paper shows great robustness with high efficiency.     

Cable Sliding at Supports in Cable Structures

To develop an effective numerical method for the cable sliding problem in cable structures，two-node catenary cable element was built to model the cables based on analytical solution of elastic catenary．Cooperated with Newton method，continuation method Was used to solve the nonlinear equations．This approach is more efficient than using Newton method only and has a wider range to select initial values for the process to converge．The relationship between the tension on a cable segment and its unstrained length was derived and used to calculate the unbalanced cable tensions at the supports．An example is presented to show the correctness and efficiency of the proposed method．     

非线性方程组求解的新方法

提出了一种求解非线性方程组的方法.将非线性方程组的求解问题转化为最优化问题,对经典BFGS变尺度法进行改进:采用高精度通用数值算法求解函数梯度,采用新的数值解析法进行一维探索,进而有效地提高了BFGS变尺度法的效率和程序通用性;对遗传算法进行了改进.将改进的BFGS变尺度法与改进的遗传算法进行混合杂交,得到一种全局优化算法,数值测试表明该算法是可靠的.     

Computer simulation of thermochemical treatments: modelling diffusion and precipitation in metals

IntroductionThermochemicaltreatmentsarewidelyusedtoimProvethemechanicalprOPertiesofmetalliccomPonents.ForexamPle,inthecaseofsteel,hightemPeraturediffusionofcarbonand/ornitrogCnleadstoanincreaseofhardnessnearthesurface,aswellastolocalcompressiveresidualstfCsses.Bothcontributetoanincreaseofthelifetimeofthecomponent.Modellingdiffusioncanbeconsideredasrelativelystraightforward,sincethegoverningequationsarewellestablished.However,onemustaccountforthefactthatduringdiffusion,precipitationmayoccur…     

快速多极多域虚边界元法解不同材料组合结构

将快速多极算法和广义极小残值法（GMRES）的基本思想运用于虚边界元法的方程求解中,并构造了多域组合问题虚边界元法的快速多极展开的实施思路,且将此方法用于不同材料组合结构问题的求解.采用此方法能够使得原问题方程组求解的计算耗时量和储存量降至与所求问题的计算自由度数成线性比例.数值算例验证了方法的可行性、计算精度和计算效率.     

基于粒子群和蚁群算法的边界层微分方程求解

沿竖壁自然对流边界层微分方程组速度和温度均耦合,在打靶法中应用Newton求根的方法解对应的相似性微分方程组时对初值选择要求较高,在根值附近收敛变慢.将微分方程边值问题转化为初值问题求解的打靶过程可看作优化设计问题,用优化设计算法求解.将基于生物群信息传递规则和觅食规则提出的粒子群算法和蚁群算法应用到打靶法的求解过程中,并与其它优化算法计算结果进行了比较.结果表明,粒子群算法和蚁群算法用于沿变壁温竖壁自然对流层流边界层微分方程求解是可行的,计算过程稳定,对初值选择不敏感.     

二维薄板在超音速气流作用下的热弹性颤振

根据von-Karman几何非线性应变位移关系建立了二维薄板在超音速气流作用下的非线性动力学模型,其中气动力根据一阶准静态活塞理论计算.用假设模态法和伽辽金方法使方程离散化,然后用Runge-Kutta方法计算.分析结果用分叉图和相图描述.结果表明,温度效应降低了板的稳定性;随动压增大,系统经历稳定运动、周期运动和混沌运动状态.     

Dynamic response analysis of towed cable during deployment/retrieval

A numerical approach was developed to analyze the transient behavior of towed cable during ac- tively controlled deployment/retrieval (DR). The cable motion is described by the lumped parameter method, its corresponding boundary conditions are presented. In view of its varying length during DR, two auxiliary arguments are introduced to describe its continuous varying length and discrete number of nodes(equations), the length is determined by the pay out(or reel-in) rate, which is then used to determine the node number by a logic relation. For the discrete mathematical model of towed cable, an algorithm was developed to deal with the discrete governing equations. The simulation results indicate that the cable experiences more com- plex motions due to its varying length, and tension fluctuates seriously in the startup and ending stage of deployment/retrieval. The effect of towing ship's motion in waves on cable during deployment/retrieval is also considered via numerical simulation.     

基于Minimax目标函数的多宗量稳态热传导反演

基于Minimax目标函数,建立了求解稳态热传导反问题优化新模型,并应用BFGS算法进行求解.同时考虑了非均质的影响,对热物性参数和边界条件进行多宗量组合反演,并对信息误差和初值作了初步探讨,数值验证给出了令人满意的结果.     

基于Minimax目标函数的多宗量稳态热传导反演

基于Minimax目标函数,建立了求解稳态热传导反问题优化新模型,并应用BFGS算法进行求解.同时考虑了非均质的影响,对热物性参数和边界条件进行多宗量组合反演,并对信息误差和初值作了初步探讨,数值验证给出了令人满意的结果.     

改进弧长法求解屈曲问题

为了解决结构非线性有限元分析求解过程中结构失稳或材料出现软化时传统Newton-Raphson法无法通过极值点的问题,在传统弧长法的基础上,提出了求解过极值点问题的改进弧长法.该方法将非线性方程求解过程中出现的不平衡力向量分解为两个相互正交的向量,并建立其弧长的约束方程,求解得到非线性计算中的荷载因子.在求解过程中,给出了改进的确定弧长方法,通过弧长调整避免了求解荷载系数中出现复根的问题.通过两个拱型结构的屈曲分析算例,分别考虑几何非线性和几何材料双重非线性效应,对两个拱形结构进行了非线性屈曲分析,结果表明:结构在出现屈曲时发生急跳现象,验证了改进的弧长法在结构出现材料软化和失稳时能通过极值点.     

EFFICIENT STEADY-STATE ANALYSIS METHOD FOR CLOSED-LOOP PWM SWITCHING CONVERTERS

Since the early70s,many s mall-signal fre-quency-domain modeling technique have been de-veloped to study the dynamic perfor mances and sta-bilities ofPWMs witching converters around thestatic operating point.However,the accuracy of thederived models relies on accurate deter mination ofthe steady state operating point,and usually manyperfor mances of circuits are given based on steady-state,so it is very essential to analyze the steadystate of circuits.Analysis methods for steadystate ofnonlin…     

Dynamic Effect Research of Cable-Lead-In Rod on Towed System

This paper aims to research the cable-lead-in rod effect on a towed system through mathematical modeling and numerical simulations. The rod dynamics, as a key part of this study, is modeled using the combination of cable node governing equations and kinematic constraint conditions. As the first attempt to analyze such a problem, the rod is simply treated as an elastic cable segment so as to be incorporated into the dynamics of the cable, and a set of algorithm is then proposed based on the kinematic constraint conditions to fully describe its motions. Meanwhile, the cable and the underwater vehicle are modeled by the traditional lumped mass method and the 6 degree-of-freedom maneuverability equations for submarines respectively; the coupling boundary conditions besides the rod dynamics are also given to form the whole system's model. Several numerical cases are performed to investigate the rod effect on the system in different maneuver situations. Some meaningful conclusions are drawn through comparative analysis.     

基于MATLAB的公路桥梁车桥耦合数值计算方法

应用达朗贝尔原理推导了两自由度车辆和桥梁的振动平衡方程,提出用龙格-库塔法和NEWMARK法来求解车桥耦合振动问题。针对NEWMARK法提出了求解分离的车辆和桥梁运动方程组的分析策略:在每一个时间步长内进行迭代计算并将桥梁的振动平稳状态作为收敛条件。利用MATLAB结合两种数值计算原理分别编制了车桥耦合计算程序。算例分析表明:两种方法的计算精度都较高;采用NEWMARK法求解时,在每个时间步内迭代计算至桥梁振动平稳状态是有意义的。     

Direct algorithms for steady-state solution of long slender marine structures

The steady state solution of long slender marine structures simply indicates the steady motion response to the excitation at top of the structure.It is very crucial especially for deep towing systems to find out how the towed body and towing cable work under certain towing speed.This paper has presented a direct algorithm using Runge-Kutta method for steady-state solution of long slender cylindrical structures and compared to the time iteration calculation;the direct algorithm spends much less time than the time-iteration scheme.Therefore, the direct algorithm proposed in this paper is quite efficient in providing credible reference for marine engineering applications.     

移动最小二乘微分求积法解Helmholtz方程

用一种新型的数值方法--移动最小二乘微分求积法（MLSDQ）求解二维Helmholtz方程。MLSDQ方法是一种直接将微分方程离散的方法，它是将未知函数的各阶偏导数在离散点处的值用域内各配点的函数值加权组合来表示，权系数则直接用移动最小二乘Galerkin法中的形函数求导得到，通过MLSDQ技术将Helmholtz方程和相应的边界条件转化成为一组关于各配点位势的线性代数方程组，求解这组代数方程，便可得到各配点的位势，通过求解几个具有精确解的算例，讨论了方法的收敛性和数值精度，结果表明：该方法较适合于求解小波数的Helmholtz方程,对高波数的方程，需要设置大量的域内配点才能有较好的数值结果。     

一类分数阶非线性微分方程组的显式算法

讨论了一类分数阶微分方程组的一种数值算法,根据Caputo导数的性质,将分数阶微分方程组转化为Volterra积分方程组,再利用求解普通积分方程的Adams技巧,建立了分数阶微分方程组的一种显式数值算法,证明了该算法的收敛性与稳定性,并给出了数值仿真实例,证实了算法的有效性.     

直齿圆柱齿轮副动力学特性分析

利用拉格朗日方程建立了含间隙直齿圆柱齿轮副的动力学模型,通过齿轮轮齿弹性变形的原理数值计算建立了时变刚度的数学模型.利用4~5阶Runge-Kutta数值积分法对系统进行了数值求解.结合Poincaré映射图、相图、FFT频谱图、系统分岔图分析了系统随激励频率和阻尼变化时的动力学行为,发现了其稳定周期运动和倍周期运动及混沌运动.通过齿轮冲击模型数值计算,找出了不同初值情况下的冲击状态.     

Numerical Solution for Stiff Dynamic Equations of Flexible Multibody System

A nonlinear numerical integration method, based on forward and backward Euler integration methods, is proposed for solving the stiff dynamic equations of a flexible multibody system, which are transformed from the second order to the first order by adop- ring state variables. This method is of A0 stability and infinity stability. The numerical solutions violating the constraint equations are corrected by Blajer's modification approach. Simulation results of a slider-crank mechanism by the proposed method are in good a- greement with ones from other literature.     

Regularization Method to the Parameter Identification of Interfacial Heat Transfer Coefficient and Properties during Casting Solidification

The accurate material physical properties, initial and boundary conditions are indispensable to the numerical simulation in the casting process, and they are related to the simulation accuracy directly. The inverse heat conduction method can be used to identify the mentioned above parameters based on the temperature measurement data. This paper presented a new inverse method according to Tikhonov regularization theory. A regularization functional was established and the regularization parameter was deduced, the Newton-Raphson iteration method was used to solve the equations. One detailed case was solved to identify the thermal conductivity and specific heat of sand mold and interfacial heat transfer coefficient (IHTC) at the meantime. This indicates that the regularization method is very efficient in decreasing the sensitivity to the temperature measurement data, overcoming the ill-posedness of the inverse heat conduction problem (IHCP) and improving the stability and accuracy of the results. As a general inverse method, it can be used to identify not only the material physical properties but also the initial and boundary conditions' parameters.