CN105808505B - 一种计算弧形地连墙侧位移的方法 - Google Patents
一种计算弧形地连墙侧位移的方法 Download PDFInfo
- Publication number
- CN105808505B CN105808505B CN201610183960.5A CN201610183960A CN105808505B CN 105808505 B CN105808505 B CN 105808505B CN 201610183960 A CN201610183960 A CN 201610183960A CN 105808505 B CN105808505 B CN 105808505B
- Authority
- CN
- China
- Prior art keywords
- plane
- unit width
- arc
- per unit
- displacement
- 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.)
- Expired - Fee Related
Links
- 238000006073 displacement reaction Methods 0.000 title claims abstract description 48
- 238000000034 method Methods 0.000 title claims abstract description 16
- 230000005483 Hooke's law Effects 0.000 claims abstract description 6
- 230000014509 gene expression Effects 0.000 claims description 19
- 238000005452 bending Methods 0.000 claims description 14
- RZVAJINKPMORJF-UHFFFAOYSA-N Acetaminophen Chemical group CC(=O)NC1=CC=C(O)C=C1 RZVAJINKPMORJF-UHFFFAOYSA-N 0.000 claims description 8
- 238000010008 shearing Methods 0.000 claims description 5
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 1
- 238000009412 basement excavation Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
- G01B21/02—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Operations Research (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Buildings Adapted To Withstand Abnormal External Influences (AREA)
Abstract
一种计算弧形地连墙侧位移的方法,所述方法通过将弧形地连墙简化为圆柱壳,建立曲线坐标系,根据经典壳体理论,建立弧形地连墙中面应变与位移之间的关系、根据胡克定律确定弧形地连墙任意一点应力与位移之间的关系,以及弧形地连墙体内任意一微元的平衡微分方程,推导出弧形地连墙的平衡微分方程,根据相应工况的边界条件,从而确定弧形地连墙的侧位移。
Description
技术领域
本发明涉及一种计算弧形地连墙侧位移的方法,属岩土工程应用技术领域。
背景技术
在基坑工程中,由于地下连续墙具有施工振动小、墙体刚度大、整体性好及施工速度快等优点,因而被广泛应用于深基坑开挖围护结构工程中。但是,目前对弧形地连墙的侧位移计算还没有非常好的方法,都只是通过数值模拟及现场监测来确定的。
发明内容
本发明的目的是,针对现有弧形地连墙的侧位移计算问题尚无好的方法,本发明提出一种计算弧形地连墙侧位移的方法。
本发明的技术方案是,本发明通过将弧形地连墙简化为圆柱壳,建立曲线坐标系,根据经典壳体理论,建立弧形地连墙中面应变与位移之间的关系、根据胡克定律确定弧形地连墙任意一点应力与位移之间的关系,以及弧形地连墙体内任意一微元的平衡微分方程,推导出弧形地连墙的平衡微分方程,根据相应工况的边界条件,从而确定弧形地连墙的侧位移。
本发明一种计算弧形地连墙侧位移的方法包括以下步骤:
(1)首先将弧形地连墙简化为圆柱壳,然后采用曲线坐标系来规定其坐标系,设定沿柱壳的母线方向为α轴,沿柱壳的周向为β轴,沿柱壳的径向为γ轴。
(2)根据经典壳体理论,得出弧形地连墙中面应变与位移之间的关系;
式中:εα为中面内各点沿α方向的线应变;εβ为中面内各点沿β方向的线应变;εαβ为中面内各点沿α及β方向的切应变;α、β、γ为各点沿α、β、γ方向上的坐标;δ为柱壳的厚度;R为中面主曲率半径;u、ν、ω分别为α、β、γ方向的扰度;q1、q2、q3分别为α、β、γ方向的外力;E为弹性模量;μ为泊松比;kα、kβ主曲率kα=1/Rα、kβ=1/Rβ;χα、χβ为中面内各点的主曲率kα及kβ的改变;χαβ为中面内各点沿α及β方向的扭率的改变。
(3)应用弹性力学中的胡克定律,计算出弧形地连墙任意一点应力和位移的关系;表达式如下式(1-2):
其中,σα为在α面上,作用于单位宽度上的应力分量;σβ为在β面上,作用于单位宽度上的应力分量;
再对壳体沿径向进行积分,可以得到弧形地连墙中面内力;
弧形地连墙中面内力为下式(1-3):
其中,薄壳所以δ为壳体的厚度;τβα为在β面上,作用于单位宽度上的切应力;Mα为在α面上,作用于单位宽度上的弯矩;Mβ为在β面上,作用于单位宽度上的弯矩;Mαβ为在α面上,作用于单位宽度上的扭矩;Mβα为在β面上,作用于单位宽度上的扭矩;FTα为在α面上,作用于中面单位宽度上的拉应力;FTβ为在β面上,作用于中面单位宽度上的拉应力;FTαβ为在α面上,作用于单位宽度上的平错力;FTβα为在β面上,作用于中面单位宽度上的平错力。
(4)将式(1-1)及式(1-2)代入上式(1-3)内力方程中,即可得到下式(1-4):
其中:D为薄壳的弯曲刚度,
(5)取壳体的任意一微元,建立弧形地连墙的内力与其所受荷载之间的关系;
(6)由于在柱壳中,横向剪力FSβ对环向平衡的影响较小,因此弧形地连墙的平衡微分方程为下式(1-6):
再将(1-4)式的后三式代入(1-6)后二式,得到:
其中
(7)经过公式代换,从而就可以得到,用中面位移表示的弧形地连墙平衡微分方程,其表达式如下式(1-8):
其中,δ为柱壳的厚度;R为中面主曲率半径;u、ν、ω分别为α、β、γ方向的扰度;E为弹性模量;μ为泊松比。
(8)由于弧形地连墙一般仅受法向载荷作用,q1=q2=0
因此将其代入公式(1-8)中,得到基本微分方程如下式(1-9)
(9)引入位移函数,即可得到中面位移表达式如下(1-10):
由于,中面位移要满足基本微分方程式,即(1-10)式要满足(1-9)式;将其带入发现,(1-9)式前两个方程总能满足,而第三个方程要求:
(10)再将(1-10)式代入(1-4)式及(1-7)式,即可将内力用位移函数F表示如下:
(11)最后根据相应工况的边界条件,由微分方程(1-11)解出F,再将其带入(1-10)式,即可求得中面位移。
本发明的有益效果是,本发明通过建立弧形地连墙中面应变与位移之间的关系,确定弧形地连墙任意一点应力与位移之间的关系,以及弧形地连墙体内任意一微元的平衡微分方程,推导出弧形地连墙的平衡微分方程,从而确定弧形地连墙的侧位移;所述方法计算弧形地连墙的侧位移准确。
附图说明
图1为本发明计算弧形地连墙侧位移方法的流程框图;
图2为定义坐标系;
图3为壳体微元薄膜内力与横向剪力,用单箭头矩矢表示;
图4为壳体微元薄膜的弯矩及扭矩,采用双箭头矩矢表示。
具体实施方式
本发明是一种计算弧形地连墙侧位移的方法具体实施如图1所示,步骤如下:
步骤一:首先将弧形地连墙简化为圆柱壳,然后采用曲线坐标系来规定其坐标系,设定沿柱壳的母线方向为α轴,沿柱壳的周向为β轴,沿柱壳的径向为γ轴,如图2所示。
步骤二:根据经典壳体理论,得出弧形地连墙中面应变与位移之间的关系如下式(1-1):
式中:εα为中面内各点沿α方向的线应变;εβ为中面内各点沿β方向的线应变;εαβ为中面内各点沿α及β方向的切应变;α、β、γ为各点沿α、β、γ方向上的坐标;δ为柱壳的厚度;R为中面主曲率半径;u、ν、ω分别为α、β、γ方向的扰度;q1、q2、q3分别为α、β、γ方向的外力;E为弹性模量;μ为泊松比;kα、kβ主曲率kα=1/Rα、kβ=1/Rβ;χα、χβ为中面内各点的主曲率kα及kβ的改变;χαβ为中面内各点沿α及β方向的扭率的改变。
步骤三:应用弹性力学中的胡克定律,计算出弧形地连墙任意一点应力和位移的关系。表达式如下式(1-2):
再对壳体沿径向进行积分,可以得到弧形地连墙中面内力。弧形地连墙中面内力为下式(1-3):
其中,薄壳所以为壳体的δ厚度;τβα为在β面上,作用于单位宽度上的切应力;Mα为在α面上,作用于单位宽度上的弯矩;Mβ为在β面上,作用于单位宽度上的弯矩;Mαβ为在α面上,作用于单位宽度上的扭矩;Mβα为在β面上,作用于单位宽度上的扭矩;FTα为在α面上,作用于中面单位宽度上的拉应力;FTβ为在β面上,作用于中面单位宽度上的拉应力;FTαβ为在α面上,作用于单位宽度上的平错力;FTβα为在β面上,作用于中面单位宽度上的平错力。
步骤四:将式(1-1)及式(1-2)代入上式(1-3)内力方程中,即可得到下式(1-4):
其中:D为薄壳的弯曲刚度,
步骤五:取壳体的任意一微元,如图3和图4所示,为表示简明,图中仅显示这一微元的中面。将其薄膜内力、横向剪力与弯矩、扭矩分开表示。图3中为薄膜内力与横向剪力,用单箭头矩矢表示,图4中为弯矩及扭矩,采用双箭头矩矢表示。图中的q1、q2、q3为每单位中面面积范围内的荷载,包括体力和面力在内。
根据图3、4,建立弧形地连墙的内力与其所受荷载之间的关系。如下式(1-5):
其中,FSα为在α面上,作用于单位宽度上的横向剪力;、FSβ为在β面上,作用于单位宽度上的横向剪力;q1、q2、q3为壳体微元上每单位中面面积范围内的荷载,即分别为α、β、γ方向的外力;
步骤六:在柱壳中,由于横向剪力FSB对于环向平衡的影响较小,可以忽略不计。因此,上列(1-5)第二式中,FSB/R≈0。即该弧形地连墙的平衡微分方程可以改写为下式(1-6):
再将(1-4)式的后三式代入(1-6)后二式,得到:
其中,
步骤七:将(1-4)式及(1-7)代入(1-6)中的前三式,其中D=Eδ3/[12(1-μ2)],即可得到中面位移表示的弧形地连墙平衡微分方程如下式(1-8):
步骤八:由于弧形地连墙一般仅受法向荷载作用,因此公式(1-8)中q1=q2=0,将其代入,得到基本微分方程如下式(1-9):
步骤九:引入位移函数F=F(α,β),将其带入,中面位移即可表达如下(1-10):
由于中面位移要满足基本微分方程式,即(1-10)式要满足(1-9)式。将其代入发现,(1-9)式前两个方程总能满足,而第三个方程要求:
步骤十:将(1-10)式代入(1-4)式及(1-7)式,即可将内力用位移函数F表示如下:
由此发现,各种工况的边界条件总可以用F进行表示。
步骤十一:在相应工况的边界条件下由微分方程(1-11)解出F,再将其代入(1-10)式,即可求得中面位移。
Claims (1)
1.一种计算弧形地连墙侧位移的方法,其特征在于,所述方法通过将弧形地连墙简化为圆柱壳,建立曲线坐标系,根据经典壳体理论,建立弧形地连墙中面应变与位移之间的关系、根据胡克定律确定弧形地连墙任意一点应力与位移之间的关系,以及弧形地连墙体内任意一微元的平衡微分方程,推导出弧形地连墙的平衡微分方程,根据相应工况的边界条件,从而确定弧形地连墙的侧位移;
所述方法步骤如下:
(1)首先将弧形地连墙简化为圆柱壳,然后采用曲线坐标系来规定其坐标系,设定沿柱壳的母线方向为α轴,沿柱壳的周向为β轴,沿柱壳的径向为γ轴;
(2)根据经典壳体理论,得出弧形地连墙中面应变与位移之间的关系;
式中:εα为中面内各点沿α方向的线应变;εβ为中面内各点沿β方向的线应变;εαβ为中面内各点沿α及β方向的切应变;α、β为各点沿α、β方向上的坐标;R为中面主曲率半径;u、ν、ω分别为α、β、γ方向的扰度;χα、χβ为中面内各点的主曲率kα及kβ的改变,kα、kβ主曲率kα=1/Rα、kβ=1/Rβ;χαβ为中面内各点沿α及β方向的扭率的改变;
(3)应用弹性力学中的胡克定律,计算出弧形地连墙任意一点应力和位移的关系,表达式如下式(1-2):
其中,σα为在α面上,作用于单位宽度上的应力分量;σβ为在β面上,作用于单位宽度上的应力分量;E为弹性模量;μ为泊松比;γ为各点沿γ方向上的坐标;
再对壳体沿径向进行积分,可以得到弧形地连墙中面内力;
弧形地连墙中面内力为下式:
其中,薄壳所以ταβ为在α面上,作用于单位宽度上的切应力;τβα为在β面上,作用于单位宽度上的切应力;Mα为在α面上,作用于单位宽度上的弯矩;Mβ为在β面上,作用于单位宽度上的弯矩;Mαβ为在α面上,作用于单位宽度上的扭矩;Mβα为在β面上,作用于单位宽度上的扭矩;FTα为在α面上,作用于中面单位宽度上的拉应力;FTβ为在β面上,作用于中面单位宽度上的拉应力;FTαβ为在α面上,作用于单位宽度上的平错力;FTβα为在β面上,作用于中面单位宽度上的平错力;
(4)将式(1-1)及式(1-2)代入上式(1-3)内力方程中,即可得到下式(1-4):
其中:D为薄壳的弯曲刚度,
(5)取壳体的任意一微元,建立弧形地连墙的内力与其所受荷载之间的关系;
其中,FSα为在α面上,作用于单位宽度上的横向剪力;、FSβ为在β面上,作用于单位宽度上的横向剪力;q1、q2、q3为壳体微元上每单位中面面积范围内的荷载,即分别为α、β、γ方向的外力;
(6)由于在柱壳中,横向剪力FSβ对环向平衡的影响较小,因此弧形地连墙的平衡微分方程为下式:
再将(1-4)式的后三式代入(1-6)后二式,得到:
其中
(7)经过公式代换,从而就可以得到,用中面位移表示的弧形地连墙平衡微分方程,其表达式如下式(1-8):
其中,R为中面主曲率半径;u、ν、ω分别为α、β、γ方向的扰度;E为弹性模量;μ为泊松比;
(8)由于弧形地连墙一般仅受法向载荷作用,q1=q2=0
因此将其代入公式(1-8)中,得到基本微分方程如下式(1-9):
(9)引入位移函数,即可得到中面位移表达式如下(1-10):
由于,中面位移要满足基本微分方程式,即(1-10)式要满足(1-9)式;将其带入发现,(1-9)式前两个方程总能满足,而第三个方程要求:
(10)再将(1-10)式代入(1-4)式及(1-7)式,即可将内力用位移函数F表示如下:
其中,Mα为在α面上,作用于单位宽度上的弯矩;Mβ为在β面上,作用于单位宽度上的弯矩;Mαβ为在α面上,作用于单位宽度上的扭矩;Mβα为在β面上,作用于单位宽度上的扭矩;FTα为在α面上,作用于中面单位宽度上的拉应力;FTβ为在β面上,作用于中面单位宽度上的拉应力;FTαβ为在α面上,作用于单位宽度上的平错力;FTβα为在β面上,作用于中面单位宽度上的平错力;FSα为在α面上,作用于单位宽度上的横向剪力;、FSβ为在β面上,作用于单位宽度上的横向剪力;
(11)最后根据相应工况的边界条件,由微分方程(1-11)解出F,再将其带入(1-10)式,即可求得中面位移。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610183960.5A CN105808505B (zh) | 2016-03-29 | 2016-03-29 | 一种计算弧形地连墙侧位移的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610183960.5A CN105808505B (zh) | 2016-03-29 | 2016-03-29 | 一种计算弧形地连墙侧位移的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105808505A CN105808505A (zh) | 2016-07-27 |
CN105808505B true CN105808505B (zh) | 2018-12-21 |
Family
ID=56454036
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610183960.5A Expired - Fee Related CN105808505B (zh) | 2016-03-29 | 2016-03-29 | 一种计算弧形地连墙侧位移的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105808505B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108717013B (zh) * | 2018-05-28 | 2019-10-18 | 华南理工大学 | 通过圆柱体侧向压缩力位移响应确定材料参数的方法 |
CN110866362A (zh) * | 2019-12-06 | 2020-03-06 | 广西交通科学研究院有限公司 | 一种地连墙设计厚度的计算方法 |
-
2016
- 2016-03-29 CN CN201610183960.5A patent/CN105808505B/zh not_active Expired - Fee Related
Non-Patent Citations (3)
Title |
---|
圆柱壳形地下连续墙构筑物的计算方法;王守忠 等;《港工技术》;19931231;第1993年卷(第4期);第2段,第1、4节 * |
地下连续墙测试中弯矩的计算方法探讨;焦志斌 等;《岩土工程学报》;20061130;第28卷;全文 * |
基于中厚圆柱壳理论的地铁隧道结构振动特性分析;王鑫 等;《岩土工程学报》;20110531;第33卷(第5期);第1、2节 * |
Also Published As
Publication number | Publication date |
---|---|
CN105808505A (zh) | 2016-07-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yoo et al. | Vibration analysis of rotating pre-twisted blades with a concentrated mass | |
Banerjee | Dynamic stiffness formulation and free vibration analysis of centrifugally stiffened Timoshenko beams | |
EP2103760B1 (en) | Method to control the vibrations in an articulated arm for pumping concrete, and relative device | |
Banerjee | Explicit frequency equation and mode shapes of a cantilever beam coupled in bending and torsion | |
CN105808505B (zh) | 一种计算弧形地连墙侧位移的方法 | |
CN107255922A (zh) | 一种基于自适应双层滑模的遥操作系统快速力估计方法 | |
CN111783198B (zh) | 一种基于双梁模型的悬索桥精细化动力分析方法 | |
CN109063305B (zh) | 随机振动环境下输流直管道减振设计方法 | |
US20180073945A1 (en) | External Force Measurement System for Work Machine, and Work Machine | |
CN108717503A (zh) | 一种设计工况下埋入式抗滑桩承担的滑坡推力计算方法 | |
Zhang et al. | A finite element method for modelling the vibration of initially tensioned thin-walled orthotropic cylindrical tubes conveying fluid | |
JP3882014B2 (ja) | 構造物の振動試験装置およびその振動試験方法 | |
CN114357661A (zh) | 柱塞-滑靴组件耦合动力学模型建立方法 | |
CN109635500A (zh) | 航空管道三维流固耦合参数共振响应特性预测方法和装置 | |
Doudkin et al. | Contact force calculation of the machine operational point | |
Steen | Elastic buckling and postbuckling of eccentrically stiffened plates | |
Shahba et al. | Derivation of an efficient non-prismatic thin curved beam element using basic displacement functions | |
Hien et al. | Free vibration analysis of joined composite conical-cylindrical-conical shells containing fluid | |
CN111753361A (zh) | 一种爆炸荷载作用下围岩动力稳定性复合判据计算方法 | |
Muhammad et al. | Structural Analysis of a Knuckle Joint using different materials | |
CN108319761A (zh) | 确定群桩位移的方法和装置 | |
KR102338844B1 (ko) | 진동저감시스템 및 진동저감시스템에서 유연체 바 구조물의 진동을 저감하는 방법 | |
CN114386144A (zh) | 基于联合强度的黄土区竖井井壁空间主动土压力测算方法 | |
Soldatos | Generalization of variationally consistent plate theories on the basis of a vectorial formulation | |
Fan et al. | Simulation analysis of critical contact load of buckling rod in tubing string |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20181221 Termination date: 20190329 |
|
CF01 | Termination of patent right due to non-payment of annual fee |