基于显-隐式格式的三维时域土-结相互作用分析的异步并行算法

来源 :中国科学:技术科学 | 被引量 : 1次 | 上传用户:dfqq209
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
提出了一种分析三维时域土-结构动力相互作用的并行计算方法.该方法利用显式和隐式积分格式各自的优点,采用集中质量显式有限元和黏弹性人工边界模拟无限域地基,采用隐式积分格式有限元法计算上部结构的响应,两者可采用不同的时间步距.在此基础上,采用MPI通信协议,提出了土-结构系统中土体区域的并行计算方法,以及土体区域与结构间的异步并行计算方法,并通过自编的Fortran源程序实现了地震作用下土-结相互作用分析的并行计算.以某一核电结构模型为例,分别采用串行算法和并行算法分析了土-核电结构体系在SV波输入下的
其他文献
本文提出了一种基于分层复杂性度量的形状描述方法 (简称HCMD),以应用于一般的形状图像检索任务.该方法属于基于区域的形状描述方法,不需要专门提取形状的边界完成特征抽取任务.本文的主要贡献是提出了一种分层的形状描述框架,基于该框架,形状沿各个方向被迭代地切割成小的区块,每一层级的形状区块被施以各种几何特性的度量,以刻画其形状复杂性.这种分层抽取形状复杂性特征的描述机制使得HCMD具有由粗到细的形状
期刊
本文讨论非PH(Pythagorean hodograph)曲线的一类五次OR曲线(curves with rational offsets)的构造方法.OR曲线是具有有理形式的等距线的一类参数曲线,在CAD中有着广泛的应用.本文采用参数曲线的复数表示形式,根据导矢曲线的因式分解的不同,将非PH曲线的五次OR曲线分为两种类型,并分别给出这两类曲线的构造方法.在给定C1连续的初始条件下,通过指定一个
期刊
~~
期刊
本栏目新闻资讯从2016年10期起以全新模式呈现给读者。以国内、国际两个视点出发,多领域体会行业发展的魅力所在,进而思考行业激烈竞争态势的发展方向。丰富内容敬请关注和浏览,欢迎扫码参与!(温馨提示:此码扫描仅适用于安卓系统手机,苹果系统版即将推出,敬请期待!)
期刊
信息传播算法求解可满足问题时有良好的有效性,使得难解区域变窄.然而,信息传播算法不总有效,常表现为不收敛.对于这种现象,至今缺少系统的理论解释.调查传播(survey propagation,SP)算法是最为有效的信息传播算法,对SP算法的收敛性研究是设计其他信息传播算法的重要基础,并为信息传播算法的广泛应用提供理论依据.为了分析SP算法的收敛性,通过对消息更新方程取双曲正切,将消息取值从[0,1
期刊
本文采用Legendre Galerkin-Chebyshev配置方法求解广义空间分数阶Burgers方程.该方法基于Legendre Galerkin变分形式,但是非线性项与右端源项采用Chebyshev-Gauss插值逼近.首先,通过在空间方向采用Legendre Galerkin-Chebyshev配置方法离散,时间方向采用leap-frog/Crank-Nicolson格式离散,得到了方程
期刊
研究了L(0,1)模态导波在管道弯头处的反射特性,提出了一种新的基于导波模态置信准则计算的L(0,1)模态导波弯头反射分析方法.使用数值模拟方法研究了L(0,1)模态导波在弯头处的模态转换,并分析了导波激励频率、弯曲半径和弯曲角度3种因素对L(0,1)模态导波反射的影响;利用导波模态置信准则表示导波弯头反射程度,反射程度随频率和弯曲半径变化趋势与数值模拟结果进行了对比,最后进行了实验验证.结果表明
期刊
二次正交Arnoldi方法的存储量小并能保持与标准Arnoldi方法类似的数值稳定性和收敛性,因此成为求解二次特征值问题的重要方法.算法的运行过程中计算量和存储量会不断增加,将算法进行重新启动是算法在实际使用中的必然需求,该方法的特殊分解形式对算法的重启提出了新要求.本文分析了该方法所形成子空间的性质和重启时子空间应具有的形式和性质,提出了一种能够保持算法特殊子空间结构且简便易实现的重启方法.在此
期刊
本文研究了Stokes方程基于混合有限元离散的预处理方法.基于三角形剖分上速度场的连续P_2元以及压力场的P_0元,我们采用局部子区域逆矩阵建立了预处理子.通过分析预处理矩阵的特征值分布,我们证明了该预处理子大大改善原问题系数矩阵的条件数.本文给出了相应的数值结果.
期刊
通过有限元空间和有限体积元空间的一种双射投影得到了不可压缩流问题低次等阶稳定有限体积元方法.该方法采用低次等阶元P1-P1(或Q1-Q1)对Navier-Stokes(N-S)方程进行数值求解,利用局部压力投影技术进行稳定化处理.通过有限元和有限体积元方法的等价性进行有限体积元方法的理论分析.发现不可压缩流N-S问题在f∈H~1时,稳定有限体积元方法与稳定有限元方法之间具有O(|logh|~(1/
期刊