CN107092730A - 适用于显式分析的三维无限元人工边界建立方法 - Google Patents
适用于显式分析的三维无限元人工边界建立方法 Download PDFInfo
- Publication number
- CN107092730A CN107092730A CN201710217997.XA CN201710217997A CN107092730A CN 107092730 A CN107092730 A CN 107092730A CN 201710217997 A CN201710217997 A CN 201710217997A CN 107092730 A CN107092730 A CN 107092730A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msubsup
- mover
- soil body
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 33
- 238000000034 method Methods 0.000 title claims abstract description 25
- 239000002689 soil Substances 0.000 claims abstract description 45
- 238000006073 displacement reaction Methods 0.000 claims description 15
- 230000008859 change Effects 0.000 claims description 4
- 230000009471 action Effects 0.000 claims description 2
- 238000004088 simulation Methods 0.000 description 12
- 230000001133 acceleration Effects 0.000 description 9
- 238000010586 diagram Methods 0.000 description 8
- 230000004044 response Effects 0.000 description 6
- 230000000694 effects Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000013016 damping Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 230000006641 stabilisation Effects 0.000 description 2
- 238000011105 stabilization Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种适用于显式分析的三维无限元人工边界建立方法,包括以下步骤:S1.创建三维土体几何模型,然后从三维土体几何模型中截取出部分的土体作为近场区域,剩余的土体作为远场区域;S2.建立三维无限元人工边界,使用无限元模拟远场区域,而采用有限元模拟近场区域;S3.将外源波动转化为等效的节点力,然后将其施加在构建的人工边界的节点上。
Description
技术领域
本发明涉及抗震分析领域,更具体地,涉及一种适用于显式分析的三维无限元人工边界建立方法。
背景技术
评估结构的安全性和完整性时,分析结构在地震作用下的响应是很重要的一项内容。结构的抗震分析受到结构、地震动、场地土特性等因素的影响,而与地基相连的结构体系,相较于地基往往具有刚度大的特点,此时刚性地基的假设并不能够反映结构真实的响应,所以研究中需要考虑土-结构相互作用的影响。
土-结构相互作用是由地基土和结构不同的材料性质在受力作用后引起形变,导致两者之间产生相互作用的问题,其动力分析不仅与结构动力学对结构响应的计算有关,还涉及到地基土的无限域模拟。土-结构相互作用动力分析研究的是波动经地基土的传播引起上部结构振动,此时上部结构产生的惯性力反作用于地基,引起地基新的振动又重新作用在结构上这样的过程。其分析方法主要分为两大类,即直接法和子结构法。子结构法基于叠加原理建立,只适用于分析线性模型,随着计算机技术的发展,建立结构复杂但精细的三维模型和考虑非线性等因素的影响,采用直接法模拟计算在SSI影响下的响应越来越受到关注。直接法是将土-结构作为整体进行动力分析,在有限元的数值模拟中,研究对象是有界的,而地基土是半无限性的,所以需要建立人工边界来模拟地基的辐射阻尼,吸收散射波在边界处产生的反射。
现阶段使用有限元法建立土-结构相互作用的模型时,通过设置人工边界模拟土的无限性,构建人工边界时主要采用粘弹性人工边界的方法。该方法基于无限空间中柱面和球面波理论推导得出,采用在截取的有限土体的边界处设置物理元件,即分别在边界的法向和切向施加离散的并联弹簧和阻尼器模拟。该法可以在ABAQUS的隐式动力分析中得以实现,但是对于大型复杂的、需要考虑非线性因素影响的三维结构模型,使用隐式动力分析计算时很难达到收敛,计算也很耗时,并且ABAQUS显式分析在保证计算稳定性的情况下可以提高计算效率,但是受到显式动力分析中无法设置接地弹簧-阻尼器的限制,影响了粘弹性人工边界在显式动力分析中的实现。
发明内容
本发明为了充分利用ABAQUS显式动力分析计算效率高的特点,克服上述现有技术在实现时存在的一些问题,提出了适用于显式分析的三维无限元人工边界建立方法。
为实现以上发明目的,采用的技术方案是:
适用于显式分析的三维无限元人工边界建立方法,包括以下步骤:
S1.创建三维土体几何模型,然后从三维土体几何模型中截取出部分的土体作为近场区域,剩余的土体作为远场区域;
S2.建立三维无限元人工边界,使用无限元模拟远场区域,而采用有限元模拟近场区域;
S3.将外源波动转化为等效的节点力,然后将其施加在构建的人工边界的节点上。
本发明在ABAQUS显式分析中采用无限元建立人工边界基于无穷远处位移为零的条件模拟地基土的无限性,截取感兴趣的三维土体,在有限土体的近场区域采用有限元模拟,而与土体的外侧即底面和四个侧面有一定距离的远场区域采用无限元模拟。该方法不仅可以保证计算的稳定性和精度,同时适用于显式动力分析,充分利用显式分析计算效率高的特点,解决了数值模拟计算复杂但结构精细的三维土-结构相互作用模型成本大的困难。
优选地,所述步骤S3中外源波动为位移场为up(t)的P波或S波,其垂直入射至人工边界处,此时人工边界上的节点l上施加的节点力fli(t)为:
其中表示输入波动的速度时程,可由up(t)导出,dt1、dt2表示行波延迟,λ为介质的拉梅常数,Al表示节点l的有效面积,Cln、Clτ分别为边界的法向、切向阻尼系数,并且Cln=ρcp,Clτ=ρcs,cp和cs分别为剪切波速和纵波波速,fli(t)表示节点l沿着i方向的节点力,i方向指的是坐标轴的x,y,z三个方向,fli(t)沿着坐标轴三个方向上的节点力分别表示为flx(t),fly(t),flz(t),节点力载荷公式中的下标lx,ly,lz表示载荷在节点l处的作用方向,上标±x,±y,±z表示节点l所在交界面的外法线方向,当与坐标轴方向相同时为正,反之为负。
与现有技术相比,本发明的有益效果是:
使用本发明提供的人工边界构建方法构建人工边界进行动力分析时不需要在边界处设置并联的弹簧-阻尼器,也可以精确的模拟地基的辐射阻尼,并且该边界可以在ABAQUS的显式动力分析中建立实施,充分的利用ABAQUS显式动力分析效率高的特点,对于大型的复杂结构在要求建立精细的三维模型和考虑非线性因素中,不仅保证计算的稳定性和精度还具有很好的计算效率。
附图说明
图1为应用本发明提供的边界建立方法的显式分析方法的流程图。
图2为本发明实例1和实例2所用的土体几何模型示意图。
图3为本发明中有限元土体四个侧面的三维无限元人工边界的几何体示意图。
图4为本发明中有限元土体底面的三维无限元人工边界的几何体示意图。
图5为本发明实例1为输入的入射波位移时程曲线的示意图。
图6为本发明实例1观察点的位移响应示意图。
图7为本发明实例1观察点的位移响应与远置边界数值解的比较示意图。
图8为本发明实例2输入的地震波加速度时程曲线示意图。
图9为本发明实例2观察点的加速度时程响应示意图。
具体实施方式
附图仅用于示例性说明,不能理解为对本专利的限制;
以下结合附图和实施例对本发明做进一步的阐述。
在ABAQUS的显式分析中建立三维的无限元人工边界进行土-结构相互作用的动力分析,如图1所示,具体几个关键实施过程如下:
(1)创建土体几何模型,截取有限的土体,将截取出的土体作为近场区域,剩余的土体作为远场区域;
(2)建立三维无限元人工边界,即在土体的外侧包括底面和四个侧面有一定距离的远场区域采用无限元模拟并赋予无限元单元特性,而近场区域的土体采用有限元模拟,网格划分时赋予有限单元特性。
(3)将外源波动输入转化为等效节点力载荷:
位移场为up(t)的P波垂直入射到有限元土体的边界处,此时在有限元与无限元交界面的节点l上,将波动转化为三维等效节点载荷的fli(t)计算公式如下:
由于考虑行波延迟,其中的l为交界面上的节点l到有限元体底面的距离,H为土体自由面到有限元体底面的距离,表示输入波动的速度时程,可由up(t)导出,dt1、dt2表示行波延迟,λ为介质的拉梅常数,Al表示节点l的有效面积,Cln、Clτ分别为边界的法向、切向阻尼系数,并且Cln=ρcp,Clτ=ρcs,cp和cs分别为剪切波速和纵波波速,fli(t)表示节点l沿着i方向的节点力,i方向指的是坐标轴的x,y,z三个方向,fli(t)沿着坐标轴三个方向上的节点力分别表示为flx(t),fly(t),flz(t),节点力载荷公式中的下标lx,ly,lz表示载荷在节点l处的作用方向,上标±x,±y,±z表示节点l所在交界面的外法线方向,当与坐标轴方向相同时为正,反之为负。同时关于S波的等效节点载荷转化以相同的原理可以推导得到。
(4)载荷施加,采用集中力的施加方法,加载在有限元与无限元交界处的每个节点上。
(5)赋予无限域的无限单元特性,还需调整无限元中的节点编号为逆时针规则编写且保证该处无限元编号后的第一个面是有限元与无限元的交界面,然后以inp文件形式进行分析计算。
(6)评估结果的合理性,由于模型的精确解很难求得,根据研究发现远置边界的直接有限元法数值解与精确解的误差在可接受的5%以内,所以最后需要将计算得到的土体响应与远置边界法得到的数值解进行比较。
实施例1:三维土的人工边界动力计算
如图2所示,采取三维土体模型进行分析截取土体的几何范围是400m×400m×300m,其中中间土体有限元近场区域的几何尺寸是:200m×200m×200m,该部分网格单元类型用C3D8R有限单元模拟,而有限单元土体的四个侧面和底面的外侧区域如图3、4所示采用无限元CIN3D8单元模拟。土体的材料参数特性为:密度ρ=2500kg/m3,弹性模量E=14.82GPa,泊松比ν=0.3,剪切波速Cs=1510m/s,压缩波速Cp=2825m/s。输入垂直入射的四分之一的正弦位移波,时长取为1s,其表达式为:当0≤t≤0.25s时u(t)=0.01sin(4πt),其余时刻位移为0,该入射波的位移时程如图5所示。取土的顶面即上表面的中心点作为观察点,经计算得到的位移时程如图6所示,并且与远置边界得到的数值解进行比较,比较的结果如图7呈现,从中可以观察到顶面中心点的位移时程不仅与远置边界计算得到的曲线变化趋势基本吻合,而且无限元边界得到的位移幅值0.022m与远置边界的位移幅值0.018m也很相近。根据波动理论自由面响应的幅值是输入位移波幅值的两倍,该方法得到的响应幅值与输入位移波幅值之间符合接近两倍的关系,充分体现了土体对波动的放大效应,也证实了在显式分析中使用三维无限元边界模拟的有效性。
实施例2:地震波载荷下三维土体的人工边界动力计算
采用实施例1中的土体模型,输入垂直入射的单向LBNS地震波,其加速度时程如图8所示,其中时长30s,加速度以重力加速度g为单位,并且加速度峰值是0.1g。将其转化为等效载荷施加在有限元与无限元交界面的节点上。取土的顶面即上表面的中心点作为观察点,计算得到其加速度时程响应如图9所示,该处只需要观察得到的土自由面中心点的加速度响应峰值是否与输入地震波峰值满足接近两倍的关系即可。结果发现得到的土自由面中心点的加速度峰值是0.239g是输入地震波峰值0.1g的2.39倍,该结果在误差允许的范围内。所以不管输入的是位移波还是加速度波,依然证实了在显式分析中使用三维无限元边界可以较准确的模拟地基土的无限性。
Claims (2)
1.适用于显式分析的三维无限元人工边界建立方法,其特征在于:包括以下步骤:
S1.创建三维土体几何模型,然后从三维土体几何模型中截取出部分的土体作为近场区域,剩余的土体作为远场区域;
S2.建立三维无限元人工边界,使用无限元模拟远场区域,而采用有限元模拟近场区域;
S3.将外源波动转化为等效的节点力,然后将其施加在构建的人工边界的节点上。
2.根据权利要求1所述的适用于显式分析的三维无限元人工边界建立方法,其特征在于:所述步骤S3中外源波动为位移场为up(t)的P波或S波,其垂直入射至人工边界处,此时人工边界上的节点l上施加的节点力fli(t)为:
<mrow>
<msubsup>
<mi>f</mi>
<mrow>
<mi>l</mi>
<mi>x</mi>
</mrow>
<mrow>
<mo>-</mo>
<mi>x</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>A</mi>
<mi>l</mi>
</msub>
<mfrac>
<mi>&lambda;</mi>
<msub>
<mi>c</mi>
<mi>p</mi>
</msub>
</mfrac>
<mo>&lsqb;</mo>
<msub>
<mover>
<mi>u</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<msub>
<mi>dt</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>u</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<msub>
<mi>dt</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<msubsup>
<mi>f</mi>
<mrow>
<mi>l</mi>
<mi>z</mi>
</mrow>
<mrow>
<mo>-</mo>
<mi>x</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>A</mi>
<mi>l</mi>
</msub>
<mo>&lsqb;</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>l</mi>
<mi>&tau;</mi>
</mrow>
</msub>
<msub>
<mover>
<mi>u</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<msub>
<mi>dt</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>l</mi>
<mi>&tau;</mi>
</mrow>
</msub>
<msub>
<mover>
<mi>u</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<msub>
<mi>dt</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<msubsup>
<mi>f</mi>
<mrow>
<mi>l</mi>
<mi>z</mi>
</mrow>
<mrow>
<mo>-</mo>
<mi>z</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>A</mi>
<mi>l</mi>
</msub>
<mo>&lsqb;</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>l</mi>
<mi>n</mi>
</mrow>
</msub>
<msub>
<mover>
<mi>u</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&rho;c</mi>
<mi>p</mi>
</msub>
<msub>
<mover>
<mi>u</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<msubsup>
<mi>f</mi>
<mrow>
<mi>l</mi>
<mi>y</mi>
</mrow>
<mrow>
<mo>-</mo>
<mi>y</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<msubsup>
<mi>f</mi>
<mrow>
<mi>l</mi>
<mi>y</mi>
</mrow>
<mrow>
<mo>+</mo>
<mi>y</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<msubsup>
<mi>f</mi>
<mrow>
<mi>l</mi>
<mi>x</mi>
</mrow>
<mrow>
<mo>+</mo>
<mi>x</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>f</mi>
<mrow>
<mi>l</mi>
<mi>x</mi>
</mrow>
<mrow>
<mo>-</mo>
<mi>x</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>f</mi>
<mrow>
<mi>l</mi>
<mi>z</mi>
</mrow>
<mrow>
<mo>-</mo>
<mi>y</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>f</mi>
<mrow>
<mi>l</mi>
<mi>z</mi>
</mrow>
<mrow>
<mo>+</mo>
<mi>y</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>f</mi>
<mrow>
<mi>l</mi>
<mi>z</mi>
</mrow>
<mrow>
<mo>+</mo>
<mi>x</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>f</mi>
<mrow>
<mi>l</mi>
<mi>z</mi>
</mrow>
<mrow>
<mo>-</mo>
<mi>x</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
其中表示输入波动的速度时程,可由up(t)导出,dt1、dt2表示行波延迟,λ为介质的拉梅常数,Al表示节点l的有效面积,Cln、Clτ分别为边界的法向、切向阻尼系数,并且Cln=ρcp,Clτ=ρcs,cp和cs分别为剪切波速和纵波波速,fli(t)表示节点l沿着i方向的节点力,i方向指的是坐标轴的x,y,z三个方向,fli(t)沿着坐标轴三个方向上的节点力分别表示为flx(t),fly(t),flz(t),节点力载荷公式中的下标lx,ly,lz表示载荷在节点l处的作用方向,上标±x,±y,±z表示节点l所在交界面的外法线方向,当与坐标轴方向相同时为正,反之为负。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710217997.XA CN107092730A (zh) | 2017-04-05 | 2017-04-05 | 适用于显式分析的三维无限元人工边界建立方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710217997.XA CN107092730A (zh) | 2017-04-05 | 2017-04-05 | 适用于显式分析的三维无限元人工边界建立方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107092730A true CN107092730A (zh) | 2017-08-25 |
Family
ID=59648632
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710217997.XA Pending CN107092730A (zh) | 2017-04-05 | 2017-04-05 | 适用于显式分析的三维无限元人工边界建立方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107092730A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110427680A (zh) * | 2019-07-26 | 2019-11-08 | 中国水利水电科学研究院 | 获取地震波斜入射作用下边坡地震动力放大效应的方法 |
CN114169043A (zh) * | 2021-11-15 | 2022-03-11 | 中国核电工程有限公司 | 时域分析模型的建模方法及装置、时域分析方法及系统 |
CN117473636A (zh) * | 2023-12-28 | 2024-01-30 | 天津风霖物联网科技有限公司 | 一种建筑物地下土体应力与形变关系测量方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130275101A1 (en) * | 2007-08-24 | 2013-10-17 | Sheng-Yuan Hsu | Method For Modeling Deformation In Subsurface Strata |
CN105224758A (zh) * | 2015-10-13 | 2016-01-06 | 河海大学 | 一种用于设计计算的有砟轨道高速铁路有限元建模分析方法 |
CN106339556A (zh) * | 2016-08-29 | 2017-01-18 | 东南大学 | 一种轨道交通桥梁和钢轨全空间噪声预测方法 |
CN106485012A (zh) * | 2016-10-20 | 2017-03-08 | 陈建永 | 一种软土基坑稳定性安全系数的计算方法 |
-
2017
- 2017-04-05 CN CN201710217997.XA patent/CN107092730A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130275101A1 (en) * | 2007-08-24 | 2013-10-17 | Sheng-Yuan Hsu | Method For Modeling Deformation In Subsurface Strata |
CN105224758A (zh) * | 2015-10-13 | 2016-01-06 | 河海大学 | 一种用于设计计算的有砟轨道高速铁路有限元建模分析方法 |
CN106339556A (zh) * | 2016-08-29 | 2017-01-18 | 东南大学 | 一种轨道交通桥梁和钢轨全空间噪声预测方法 |
CN106485012A (zh) * | 2016-10-20 | 2017-03-08 | 陈建永 | 一种软土基坑稳定性安全系数的计算方法 |
Non-Patent Citations (2)
Title |
---|
苏欣: "时域显式方法模拟二维复杂场地地震反应的实现及其应用", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
魏道攀: "考虑土-路面-结构相互作用的城市高架桥地震反应分析", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110427680A (zh) * | 2019-07-26 | 2019-11-08 | 中国水利水电科学研究院 | 获取地震波斜入射作用下边坡地震动力放大效应的方法 |
CN110427680B (zh) * | 2019-07-26 | 2021-07-27 | 中国水利水电科学研究院 | 获取地震波斜入射作用下边坡地震动力放大效应的方法 |
CN114169043A (zh) * | 2021-11-15 | 2022-03-11 | 中国核电工程有限公司 | 时域分析模型的建模方法及装置、时域分析方法及系统 |
CN117473636A (zh) * | 2023-12-28 | 2024-01-30 | 天津风霖物联网科技有限公司 | 一种建筑物地下土体应力与形变关系测量方法 |
CN117473636B (zh) * | 2023-12-28 | 2024-03-12 | 天津风霖物联网科技有限公司 | 一种建筑物地下土体应力与形变关系测量方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102866637B (zh) | 一种基于二次降阶的带操纵面机翼非定常气动力模拟方法 | |
CN104281730B (zh) | 一种大转动变形的板壳结构动响应的有限元分析方法 | |
JPWO2014002977A1 (ja) | 空気と水と土骨格の連成計算装置および連成計算方法並びに連成計算プログラム | |
CN111914446B (zh) | 一种有限元数值分析中超临界角斜入射地震动输入方法 | |
CN104317985B (zh) | 一种基于界带有限元和拉格朗日坐标的流体仿真方法 | |
CN107092730A (zh) | 适用于显式分析的三维无限元人工边界建立方法 | |
CN102521514A (zh) | 深水桥墩考虑水—墩耦合效应的地震响应评估方法 | |
CN112949065B (zh) | 模拟层状岩体力学行为的双尺度方法、装置、存储介质及设备 | |
CN110717216A (zh) | 不规则波下带柔性气囊直升机横摇响应预报方法 | |
Guo et al. | Simulation of dynamic 3D crack propagation within the material point method | |
Zhu et al. | Simulation of large-scale numerical substructure in real-time dynamic hybrid testing | |
CN103886138A (zh) | 一种考虑流固耦合的核电厂外层安全壳外挂水箱结构的抗震设计方法 | |
Sun et al. | Benchmarking the material point method for interaction problems between the free surface flow and elastic structure | |
CN109101752B (zh) | 一种复杂水工建筑物局部结构自振频率计算方法 | |
CN101261651A (zh) | 主梁断面气动自激力的全过程数值模拟方法 | |
CN104318005A (zh) | 一种基于Virtools的多元数据融合虚拟试验显示系统 | |
Arbain et al. | Vibration analysis of Kenyir dam power station structure using a real scale 3D model | |
CN116432498A (zh) | 一种适用于边坡动力分析的地震动输入实现方法 | |
Fu et al. | Dual reciprocity boundary element based block model for discontinuous deformation analysis | |
Zhang et al. | Rotation errors in numerical manifold method and a correction based on large deformation theory | |
Han et al. | Transient analysis of three‐dimensional dynamic interaction between multilayered soil and rigid foundation | |
Shirkhanghah et al. | The Effects of Height and Plan on the Along-Wind Response of Structures Considering Wind-Soil-Structure Interaction | |
Werkle | Modeling of soil-foundation-structure interaction for earthquake analysis of 3D BIM models | |
Lane et al. | Vehicle-track-underground modeling of rail induced wave propagation | |
Kavand et al. | 3-Dimensional Numerical Modelling of Pile Group Response to Liquefaction-induced Lateral Spreading |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170825 |