二维海底电缆铺设的微分方程模型与数值算法.docVIP

７ ９ ０ 　 数 学 物 理 学 报　 Ｖ ｏ １ ． ３ ５ Ａ 　 Ｇ ｉ （ ｕ） 　 ｎ ＋ ＋ 　 ｌ 州 蝌 ｎ ＋ 　 ｌ 一 　 Ｈ） 一 　 阻 　 （ 　¨ 一 ｉ 州１ 、 ） 　 － （ 　 ） 　 ｘ ｎ ＋ ｌ 　 ￣ ｎ － － １ 　 ＋ 　 （ 　 ） 　 ｘ ｎ ￣ ｌ 　 ￣ ｎ ＋ ｌ 　 一 （ ｐ —ｐ ｏ ） ｇ —ｆ ｌ ［ ￣ ／ ｎ ＋ ｌ （ ｙ ｎ ＋ 　 一 　） 一Ｏ ｎ ＋ １ （ 　一 　一 　 ） ］ 　 （ ４ ． ３ ） 　 ＝ ０ ， 　 ｉ （ ｕ ） ＝（ 　＿ １ １ 一 　 ｎ ｉ ＋ 　 ） 　 ＋（ 　Ｉ ｌ —ｙ ｎ ＋ 　 ） 。 一（ △ ｓ 件 １ ） 　 ＝０ ， 　 其 中， 　 ２ 　 — 　 ２ 　 （ ４ ． ４ ） 　 Ｏ Ｌ ｉ 　 Ａｓ — ｉ ＋ ］ （ ＡＳ ｉ — ＋ ｌ ￣Ａｓ ｉ ） ’ 　 ２ 　 — Ａｓ ｉ （ Ａｓ ｉ ＋ — ｌ ￣ Ａｓ ｉ ） ’ 　 ２ 　 ＋ 　 瓦　 丽 ’Ｏ ｎ ＋ １ 　— ／ ｋ ｔ （ ／ ｋ ｔ ｎ ＋ — ｌ ￣／ ｋ ｔ ｎ ） ， 　 ｎ Ｕ＝（ 　－ － ｌ 　 ｘ 　 …， 　 揣 ， 　 ¨， …， 　 ｙ 　 ｎ ＋ 　 ｌ 　 ， Ｔ ｉ ｎ ¨， … 　 Ｔ ｎ ＋ ｌ ｎ 　 ＋ 接下来，我们在初始节点 （ ８ ｉ ， ｔ ｏ ） 处，将 （ ２ ． １ ５ ） 一 （ ２ ． １ ８ ） 式离散化如下　 ＝ ０ ， Ｙ 　＝ ８ ｉ ， 　 —　 ＝ ０， Ｙ 　— Ｙ 　 ＝ ０， ＴＯ ＋ 　＝ （ ｐ—Ｐ ｏ ） ｇ ｓ ｉ ＋ 　 ， ｉ ＝０ ， １ ， ． 一， ／ ｏ ． 　 。 （ ４ ． ５ ） 　 设 ｒ 表 示 已经沉 到海底 的 节点 总数 ，则在 边界 节点处 将 （ ２ ． １ ９ ） 一 （ ２ ． ２ １ ） 式离散 化如 下　 ＋　 ＝ 　 ， ＋　 ｌ 　 Ｔ　 ｎ ＋ ｆ Ｊ ０ 。 ， 　 ￣ ｎ 　 ＋ ＋ ｌ 　一 　ｕ ｐ （ ｔ ｎ ＋ １ ） 　 ｙ ｎ ＋ 　 ｌ 　＝０ 　 Ｔ ｎ ＋ ＋ ｌ 卅　一 　 ｎ ＋ １ ） ， 札一０ ， １ ， …， Ⅳ 一１ ． 　 （ ４ ． ６ ） 　 注 ４ ． １ 对 于任 意 时刻 ｔ ， 怎样 获取 ｒ的值 ？首先 ，当 ｔ＝ ０时，令 ＿ 厂＝ ０ ， 并 设 定任 意　 小 的常量　 ＞ ０ ．假 设在 ｔ＝ ｔ 　 时刻 有 ｒ＝ ｒ 　 ．则 在 ｔ＝ ｔ 卅 １时刻 ，我 们 比较　 ＋ 　 （ 　＝ 　 ｒ ＋１ ， ｒ ＋２ ， …， 　 ＋几 ） 与 ￡的大小，并令 ｆ ＝｛ 　ｌ ｙ Ｌ ＋ 　＜Ｅ ， ｉ 　＝ｒ ＋１ ， ｒ ＋２ ， …， 　 ＋佗 ） ． 在　 ｔ ＝ｔ ｎ ＋ ｌ 时刻 ｒ的值 为 ｒ—ｒ 　 ＋ｆ ． 其 中 ，对任 意集 合　 ， 　 表 示　 中元 素的个 数 ． 　 由（ ４ ． ２ ）（ ４ ． ４ ） 式可知 ，在　 时１ 时刻 ，当 ｉ ＝ 　 ＋ｎ时 ， 　 ｎ 。 　 － ＋ １ 　 ， 　， ｎ 。 － ＋ １ 　 ， 　如 ｎ － ＋ １ 　 是 未知 的 （ 因　 为 ｎ—ｌ时刻 ，第　 ＋几个节 点还未 沉放 入海底 ） ． 因此 ，我们需 要给这 几个 变量 取定 一个 合　 适 的近 似值 ，根 据几 何关 系有　 ／ ￣８ 　 Ｉ 一 　 ＋　 ｏ ＋ｎ ． （ 　 ｎ … － １ 一 　＿ｘ … ｎ － １ 一 ） ， 　 一 （ ４ ． ７ ） 　 赫 一 　 ＋ 　 Ａ 　 ｓ Ｉ ￣ 　 ￣ ｎ 　 （ 　 ｎ － … １ 　一 　 ｎ ¨ － １ 一） ．　 （ ４ ． ８ ） 　 这 种近 似处理 以后 ， 　 ｎ ｆ ｎ － 　 １， 　ｎ ｆ ｎ － 　 １就 可 以参 与计 算 了． 　 由（ ２ ． １ ７ ） 和 （ ２ ． １ ８ ） 式 可知 ， 在每个 时刻 ｔ ， 方 程组 （ ２ ． １ ４ ） 式 的弧长参 数取值 范 围为 ａ （ ｔ ） 　 ８ 　 ６ （ 　 ） ， 其中口 （ 　 ）随着沉入 海底 的部 分越来 越 多而逐 渐增 大 ， 　６ （ 　 ） 随着探 测船 行走过 程 中　 释放缆 绳也 逐渐增 大 ，这就 涉及 到一个 自由边界 计算 的 问题 ［ ２ ］ 　 ． 自由边界 是个 很复杂 的 问题 ．在这 篇论 文 中，上边界 ｂ （ ｔ ）由 （ ２ ． １ ８ ） 式 定义 ，下边界 ０ （ 　 ） 　 是 未知 的．怎样取 ｎ （ ｔ ） 的值呢 ？本文采 用如下方 法 ： 为 了让 已经坠入 海底 的