计算地球圆盘负荷产生的位移

1.研究背景

        计算受表面载荷影响的弹性体变形问题有着悠久的历史,涉及到许多著名的数学家和物理学家(Boussinesq 1885;Lamb 1901;Love 1911,1929;Shida 1912;Terazawa 1916;Munk & MacDonald 1960;Longman 1962,1963;Farrell 1972)。对于均匀弹性半空间的易于编码的解(Becker & Bevis 2004)在许多工程上是有用的,但在地球不是半空间也不均匀的地质应用中通常不适用。对于分层弹性半空间的解(例如Pan等人2007),其适用范围更广,但即使对于具有大孔径的载荷或在离载荷远处计算变形时,这种形式主义也不合适,这里所说的“大”是指距离地球半径的显著部分或更大的距离。因此,大多数地球物理应用的首选框架是分层、弹性、自引力球体和液态核心(Farrell 1972)。这个问题的解通常都是在涉及到载荷Love数的展开式中进行的,通常用符号(hn,kn,ln)来表示。通常,这些处理采用点载荷或圆盘载荷

        点载荷形式主义在接近标称点载荷时通常无效,从这个意义上讲,圆盘载荷形式主义更灵活。我们注意到,Love数形式主义适用于球对称弹性地球模型;横向各向同性可以被纳入这类模型,只要它符合球对称性(Pan等人2015)。弹性各向异性的一般类(例如横向变化的各向异性和方位各向异性)打破了球对称性,在这种情况下,Love数方法就不再适用,需要使用不同的方法,例如自引力有限元模型(虽然它们也存在自己的计算困难)。使用本文中提出的方程和代码的人隐含地假设弹性各向异性要么不存在,要么是一种非常特殊的类型。

        在1975年至2000年期间,关于表面载荷问题的地球物理文献大多集中在冰盖地壳回弹调整方面,这涉及到粘弹性地球的载荷响应(Haskell 1935;Cathles 1975;Peltier & Andrews 1976;Wu & Peltier 1982)。几乎所有这些研究都使用了Love数形式主义。载荷问题的线性弹性和线性粘弹性版本通过弹性-粘弹性对应原理相互联系(Alfrey 1944;Read 1950;Lee 1950)。但是,从2001年开始,人们对纯弹性问题的兴趣出现了重大回升。这是因为(i)大地测量学家和地球物理学家开始意识到全球定位系统(GPS)接收机网络,以及更近期的全球导航卫星系统(GNSS)接收机,正在记录由水、雪和冰固体地球上的载荷所引起的季节性、持续性(即渐进性)和瞬时的弹性位移(Blewitt等人2001;Heki 2001;Mangiarotti等人2001;Dong等人2002;Bevis等人2005),而且在较小程度上还受大气压力变化的影响(Vandam等人1994),以及(ii)地球弹性变形的大地测量观测可以用于监测冰体质量(Khan等人2010;Bevis等人2012;Nielsen等人2012;Spada等人2012;Nielsen等人2013)和陆地水储存(Bevis等人2004,2005;Steckler等人2010;Fu等人2013;Borsa等人2014)的变化,从而提供了一种研究气候循环和气候变化的新方法。

        描述分布式负荷或负荷变化的最简单方法是将其表示为一个或多个圆盘载荷。如果将质量或质量变化场表示为网格,则每个网格单元都可以视为圆盘载荷,或者,如果它位于地球测量站的远场,也可以视为点载荷。点载荷的响应可以与小而远的圆盘载荷的相似,除了在接近点载荷的地方,点载荷概念本身几乎从来不是现实的。本文的目的是两方面的。首先,我们希望向非专家提供有关如何计算地球对圆盘载荷的弹性响应的简单但完整的讨论。我们提供了一个MATLAB函数(diskload)来实现这个算法。我们的代码,像几乎所有它的同类一样,都会求和由载荷Love数和勒让德多项式或勒让德多项式的导数组成的项,其阶数为n = 0、1、2,...nmax。选择适当的nmax值非常重要,特别是当我们寻求对载荷的近场响应的准确估计时。当建模地球测量站位于质量损失正在使用非常高分辨率数字高程模型(DEM)进行表征的冰盖边缘附近时,由于DEM的网格间距或单元大小,适当的nmax值会按比例增加,所以这个错误相当容易发生。我们希望这个讨论能帮助我们的读者更自信地使用更一般和更复杂的软件包,例如REAR(Melini等人2015a,b)。(翻译自Bevis et al., 2016)

2.基本原理

        所谓的“圆盘载荷”是一种特定类型的表面质量密度分布,其特点是(i)均匀施加的压力,即恒定的载荷“厚度”,和(ii)轴对称,这两个特点使得圆盘载荷在与冰盖等高度调整相关的应用中非常方便(例如Spada等人2012;Melini等人2015b),因为它可以直接展开成一系列球谐函数。此处略去基本理论公式的推导,感兴趣的参考原文。

引用至Bevis et al.(2016)

水平和垂直位移计算公式:

3.matlab程序

下载地址:

https://github.com/danielemelini/diskload

diskload MATLAB函数用于计算对施加在给定角半径和高度的圆盘上的均匀表面压力载荷的响应。弹性响应是在地球表面上的一个或多个点上计算的,这些点与圆盘载荷的中心指定的角距离相隔。弹性响应是使用用户提供的弹性负载Love数(h,k,l)计算的,这些Love数是使用特定的地球弹性结构模型生成的。下图是h l k负荷love数随阶数的变化示意图。

下载程序,在matlab中打开:

clear
%% Load LNs
load REF_6371_loading_love_numbers
%% Set some constants
alpha = 0.1;                           % Disk radius (degrees)
theta = linspace(0,alpha*5,100);       % Range of colatitudes wrt disk center
Tw    = 1;                             % Disk height (equivalent water height, m)
nmin  = 0;                             % Minimum degree
nmax  = [100:100:40000];               % Range of maximum degrees
imass = 1;                             % choose imass,0 or 1 (uncompensated/compensated load)
%%
if imass==1
    fprintf('invoking a globally compensated load (icomp=1)\n')
else
    fprintf('invoking an uncompensated load (icomp=0)\n')
end
% Compute the disc response for the maximum value of nmax
[U,V,G]= diskload(alpha,imass,theta,Tw,nmin,nmax(end),h_love,k_love,l_love);
%% FIG 1
figure(1); clf;
plot(  theta./alpha, U, 'b', 'LineWidth', 1.5 );  hold on;
plot(  theta./alpha, V, 'r', 'LineWidth', 1.5 );
plot(  theta./alpha, G, 'g', 'LineWidth', 1.5 );
xlabel( '\theta/\alpha', 'FontSize', 16 );
ylabel( 'mm', 'FontSize', 16 );
xlim([0 5]);
ylim([-2.5 1]);
grid on;
legend( 'U', 'V', 'G', 'Location', 'best' );
tit1=['Disk radius \alpha = ',sprintf('%4.2f^\\circ',alpha)];
tit2=['  Load = ',sprintf('%4.2f',alpha),' m w.e.'];
title([tit1 tit2])
%% Examine sensitivity to nmax
% Compute the disc response vs nmax at theta = K * alpha
K=1.5;
[U,V,G]= diskload(alpha,imass,K*alpha,Tw,nmin,nmax,h_love,k_love,l_love);
%% FIG 2
figure(2); clf;
semilogx(  nmax, U, 'b', 'LineWidth', 1.5 );  hold on;
semilogx(  nmax, V, 'r', 'LineWidth', 1.5 );
semilogx(  nmax, G, 'g', 'LineWidth', 1.5 );
xlabel( 'n_{max}', 'FontSize', 16 );
ylabel( 'mm', 'FontSize', 16 );
xlim( [100 40000] );
ylim( [-1 .3] );     % good for K=1.5 or 1.25
%ylim( [-0.6 .2] );   % good for K = 1.75 or 2
nROT = 360 / alpha;
semilogx( [nROT nROT], ylim, 'k-.', 'LineWidth', 1.1 ); 
semilogx( 2*[nROT nROT], ylim, 'k-.', 'LineWidth', 1.1 ); 
semilogx(xlim,U(end)*[1 1],'b--', 'LineWidth', 0.5 );
semilogx(xlim,V(end)*[1 1],'r--', 'LineWidth', 0.5 );
semilogx(xlim,G(end)*[1 1],'g--', 'LineWidth', 0.5 );
hold off;
grid on;
legend( 'U', 'V', 'G', 'Location', 'best' );
Ht=title(['Truncation error: Loading response computed',...
       ' at \theta = ',sprintf('%4.2f',K),...
       ' \alpha as a function of n_{max} ']);
set(Ht,'FontSize',14)

运行结果:

关于负荷变形计算的问题,可以参见以下的文章:

利用GRACE和地表质量模型计算地球弹性负荷变形实验笔记_如何通过grace数据导出位移-CSDN博客

LoadDef python工具包学习二:利用JPL mascon数据计算地表站点的三维变形(E、N、V)_loaddef grace]-CSDN博客

参考文献:

Bevis, M., Melini, D., & Spada, G. (2016). On computing the geoelastic response to a disk load. Geophysical Journal International205(3), 1804-1812.

♥欢迎点赞收藏♥

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/763788.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Springboot笔记-02

1.使用Spring Initializer快速创建Spring Boot项目 IDE都支持使用Spring的项目创建向导快速创建一个Spring Boot项目; 选择我们需要的模块; 向导会联网创建Spring Boot项目; 默认生成的Spring Boot项目; 主程序已经生成好了&#…

Linux:http协议初步认识

文章目录 OSI七层模型http协议域名路径信息请求和响应 编写一个httpserver OSI七层模型 在结束了前面对于序列化反序列化等内容的学习后,重新回到对于OSI模型的部分 如上所示的是对于OSI接口的示意图,在这当中可以看到会话层的概念,会话层的…

QPainter绘图和QChart图表和QCustomplot绘曲线图

一,QPainter绘图 Qt里的所有绘图,比如一个按钮和一个Label的显示,都有绘图系统来执行。绘图系统基于 QPainter、QPaintDevice和QPainEngine类。QPainter是可以直接用来操作绘图的类,而 QPaintDevice和QPainEngine都比QPainter更底…

FinalShell下载安装 + FinallShell远程连接Linux 教程

图形化,命令行 : windows下的图形化和命令行 : Linux下的图形化和命令行 : 在Linux下,普遍使用命令行操作 ; FinalShell 下载安装: 下载安装地址 : 这里给出windows的下载地址 : install 双击打开就好 : 然后点击下一步 : 选择下载地址&…

vue学习日记14:工程化开发脚手架Vue CLI

一、概念 二、安装 1.全局安装&查看版本 注意启动cmd输入命令 要以管理员运行哦 安装了一次就行以后不用再创建了 yarn global addvue/cli vue --version 显示了版本号即可 2.创建项目架子 创建项目的路径在哪 项目就在哪 项目名字不能用中文 vue create project-n…

搜索测试题题解(3月19号总结)

目录 1.Dungeon Master 2.Oil Deposits 3.Find a way 1.Dungeon Master Sample InputcopyOutputcopy 3 4 5 S.... .###. .##.. ###.###### ##### ##.## ##...##### ##### #.### ####E1 3 3 S## #E# ###0 0 0Escaped in 11 minute(s). Trapped! 这道题与普通的bfs模板题就是…

线性顺序表算法库

list.cpp 具体函数实现 #include <stdio.h> #include "list.h" #include <malloc.h>/************************************************** ①函数名: CreateList 功 能: 用数组构建顺序表 参 数: ①SqList *&L:传入的线性表 ②ElemType a[]:使用…

Iterator对象功能学习

package config;import java.util.Iterator; import java.util.Properties; import java.util.Set;/*** 这个类演示了如何使用Properties类来存储和访问键值对。* Properties类继承自Hashtable&#xff0c;因此它可以用来存储键值对数据&#xff0c;且支持同步。*/ public clas…

代码随想录算法训练营第十七天|110、平衡二叉树

110. 平衡二叉树 已解答 简单 相关标签 相关企业 给定一个二叉树&#xff0c;判断它是否是 平衡二叉树 示例 1&#xff1a; 输入&#xff1a;root [3,9,20,null,null,15,7] 输出&#xff1a;true示例 2&#xff1a; 输入&#xff1a;root [1,2,2,3,3,null,null,4,4] 输出&am…

计算机二级Python基础操作题

题目来源&#xff1a;计算机二级Python半个月抱佛脚大法&#xff08;内呈上真题版&#xff09; - 知乎 第4&#xff0c;5&#xff0c;6&#xff0c;7&#xff0c;9&#xff0c;10&#xff0c;11套 1. 基础题1 sinput() print("{:\"^30x}".format(eval(s))) b …

安科瑞ADL逆流检测多功能电表

ADL400 导轨式多功能电能表&#xff0c;是主要针对电力系统&#xff0c;工矿企业&#xff0c;公用设施的电能统计、管理需求而设计的一款智能仪表&#xff0c;产品具有精度高、体积小、安装方便等优点。集成常见电力参数测量及电能计量及考核管理&#xff0c;提供上48月的各类电…

2024年企业上云补贴,免费领取腾讯云优惠代金券

腾讯云代金券领取渠道有哪些&#xff1f;腾讯云官网可以领取、官方媒体账号可以领取代金券、完成任务可以领取代金券&#xff0c;大家也可以在腾讯云百科蹲守代金券&#xff0c;因为腾讯云代金券领取渠道比较分散&#xff0c;腾讯云百科txybk.com专注汇总优惠代金券领取页面&am…

springboot实现文件上传

SpringBoot默认静态资源访问方式 首先想到的就是可以通过SpringBoot通常访问静态资源的方式&#xff0c;当访问&#xff1a;项目根路径 / 静态文件名时&#xff0c;SpringBoot会依次去类路径下的四个静态资源目录下查找&#xff08;默认配置&#xff09;。 在资源文件resour…

3·15日,上海飞北京,东航全球首架C919亲测初体验

引言&#xff1a;“望闻问切”亲测 感受C919机型的航班 【阿明观察 &#xff5c; 科技热点关注】 赶巧了&#xff01;2024年3月15日消费者权益日这天&#xff0c;上海飞北京&#xff0c;我选择了采用C919的东方航空公司航班。 真赶巧了&#xff01;上了飞机后我才知道&…

Spring Boot项目中使用MyBatis连接达梦数据库6

在开发中&#xff0c;使用Spring Boot框架结合MyBatis来操作数据库是一种常见的做法。本篇博客将介绍如何在Spring Boot项目中配置MyBatis来连接达梦数据库6&#xff0c;并提供一个简单的示例供参考。(达梦六不仅分表还分模式.) 我拿SYSTEM表的LPS模式下面Student表做案例。 1.…

GPT中使用的Gaussian Error Linear Unit (GELU)

GPT中使用的Gaussian Error Linear Unit (GELU) flyfish 高斯误差线性单元&#xff08;GELU&#xff09; Gaussian Error Linear Unit g e l u ( x ) x P ( X ≤ x ) X ∼ N ( 0 , 1 ) gelu(x) x P(X \leq x) \qquad X \sim \mathcal{N}(0, 1) gelu(x)xP(X≤x)X∼N(0,1) G…

2024年【山东省安全员C证】免费试题及山东省安全员C证作业考试题库

题库来源&#xff1a;安全生产模拟考试一点通公众号小程序 山东省安全员C证免费试题根据新山东省安全员C证考试大纲要求&#xff0c;安全生产模拟考试一点通将山东省安全员C证模拟考试试题进行汇编&#xff0c;组成一套山东省安全员C证全真模拟考试试题&#xff0c;学员可通过…

【Qt】使用Qt实现Web服务器(六):QtWebApp用户名密码登录

1、示例 1)演示 2)登录 3)显示 2、源码 示例源码Demo1->LoginController void LoginController::service(HttpRequest& request, HttpResponse& response) {

【Linux】/proc文件系统

&#x1f525;博客主页&#xff1a;PannLZ &#x1f618;欢迎关注&#xff1a;&#x1f44d;点赞&#x1f64c;收藏✍️留言 文章目录 /proc文件系统1.获取与进程相关的信息:/proc/ID2./proc 目录下的系统信息3. 访问/proc 文件4.动态创建/proc文件系统4.1创建目录4.2创建proc…

双点双向路由引入实验

双点双向路由引入实验 1、OSPF和ISIS路由协议的优先级分别是什么&#xff1a;OSPF&#xff1a;10&#xff0c;150&#xff0c;ISIS&#xff1a;15&#xff0c;15 2、加表原则&#xff1a;当不同的路由协议学习到相同的目的网络&#xff0c;比较优先级&#xff0c;优先级数值小…