CN105808505B - 一种计算弧形地连墙侧位移的方法 - Google Patents

一种计算弧形地连墙侧位移的方法 Download PDF

Info

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
Application number
CN201610183960.5A
Other languages
English (en)
Other versions
CN105808505A (zh
Inventor
耿大新
余小强
童立红
闵世超
谢建华
郑焕
曾宇凡
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
East China Jiaotong University
Original Assignee
East China Jiaotong University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by East China Jiaotong University filed Critical East China Jiaotong University
Priority to CN201610183960.5A priority Critical patent/CN105808505B/zh
Publication of CN105808505A publication Critical patent/CN105808505A/zh
Application granted granted Critical
Publication of CN105808505B publication Critical patent/CN105808505B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/13Differential equations
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/02Measuring 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βα为在β面上,作用于单位宽度上的扭矩;F为在α面上,作用于中面单位宽度上的拉应力;F为在β面上,作用于中面单位宽度上的拉应力;FTαβ为在α面上,作用于单位宽度上的平错力;FTβα为在β面上,作用于中面单位宽度上的平错力。
(4)将式(1-1)及式(1-2)代入上式(1-3)内力方程中,即可得到下式(1-4):
其中:D为薄壳的弯曲刚度,
(5)取壳体的任意一微元,建立弧形地连墙的内力与其所受荷载之间的关系;
(6)由于在柱壳中,横向剪力F对环向平衡的影响较小,因此弧形地连墙的平衡微分方程为下式(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βα为在β面上,作用于单位宽度上的扭矩;F为在α面上,作用于中面单位宽度上的拉应力;F为在β面上,作用于中面单位宽度上的拉应力;FTαβ为在α面上,作用于单位宽度上的平错力;FTβα为在β面上,作用于中面单位宽度上的平错力。
步骤四:将式(1-1)及式(1-2)代入上式(1-3)内力方程中,即可得到下式(1-4):
其中:D为薄壳的弯曲刚度,
步骤五:取壳体的任意一微元,如图3和图4所示,为表示简明,图中仅显示这一微元的中面。将其薄膜内力、横向剪力与弯矩、扭矩分开表示。图3中为薄膜内力与横向剪力,用单箭头矩矢表示,图4中为弯矩及扭矩,采用双箭头矩矢表示。图中的q1、q2、q3为每单位中面面积范围内的荷载,包括体力和面力在内。
根据图3、4,建立弧形地连墙的内力与其所受荷载之间的关系。如下式(1-5):
其中,F为在α面上,作用于单位宽度上的横向剪力;、F为在β面上,作用于单位宽度上的横向剪力;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βα为在β面上,作用于单位宽度上的扭矩;F为在α面上,作用于中面单位宽度上的拉应力;F为在β面上,作用于中面单位宽度上的拉应力;FTαβ为在α面上,作用于单位宽度上的平错力;FTβα为在β面上,作用于中面单位宽度上的平错力;
(4)将式(1-1)及式(1-2)代入上式(1-3)内力方程中,即可得到下式(1-4):
其中:D为薄壳的弯曲刚度,
(5)取壳体的任意一微元,建立弧形地连墙的内力与其所受荷载之间的关系;
其中,F为在α面上,作用于单位宽度上的横向剪力;、F为在β面上,作用于单位宽度上的横向剪力;q1、q2、q3为壳体微元上每单位中面面积范围内的荷载,即分别为α、β、γ方向的外力;
(6)由于在柱壳中,横向剪力F对环向平衡的影响较小,因此弧形地连墙的平衡微分方程为下式:
再将(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βα为在β面上,作用于单位宽度上的扭矩;F为在α面上,作用于中面单位宽度上的拉应力;F为在β面上,作用于中面单位宽度上的拉应力;FTαβ为在α面上,作用于单位宽度上的平错力;FTβα为在β面上,作用于中面单位宽度上的平错力;F为在α面上,作用于单位宽度上的横向剪力;、F为在β面上,作用于单位宽度上的横向剪力;
(11)最后根据相应工况的边界条件,由微分方程(1-11)解出F,再将其带入(1-10)式,即可求得中面位移。
CN201610183960.5A 2016-03-29 2016-03-29 一种计算弧形地连墙侧位移的方法 Expired - Fee Related CN105808505B (zh)

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)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108717013B (zh) * 2018-05-28 2019-10-18 华南理工大学 通过圆柱体侧向压缩力位移响应确定材料参数的方法
CN110866362A (zh) * 2019-12-06 2020-03-06 广西交通科学研究院有限公司 一种地连墙设计厚度的计算方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
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