CN105005662B - 分析结果精确的输电铁塔杆件应力计算方法 - Google Patents
分析结果精确的输电铁塔杆件应力计算方法 Download PDFInfo
- Publication number
- CN105005662B CN105005662B CN201510423172.4A CN201510423172A CN105005662B CN 105005662 B CN105005662 B CN 105005662B CN 201510423172 A CN201510423172 A CN 201510423172A CN 105005662 B CN105005662 B CN 105005662B
- Authority
- CN
- China
- Prior art keywords
- mrow
- mover
- msub
- node
- overbar
- 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.)
- Active
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 20
- 238000004364 calculation method Methods 0.000 title claims abstract description 9
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims abstract description 90
- 239000011159 matrix material Substances 0.000 claims abstract description 72
- 238000006073 displacement reaction Methods 0.000 claims abstract description 71
- 229910000831 Steel Inorganic materials 0.000 claims abstract description 53
- 239000010959 steel Substances 0.000 claims abstract description 53
- 229910052742 iron Inorganic materials 0.000 claims abstract description 45
- 239000000463 material Substances 0.000 claims description 40
- 238000000034 method Methods 0.000 claims description 15
- 238000010008 shearing Methods 0.000 claims description 15
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 230000000694 effects Effects 0.000 claims description 8
- 230000007246 mechanism Effects 0.000 claims description 6
- 238000012546 transfer Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 5
- 238000011156 evaluation Methods 0.000 abstract description 4
- 239000000203 mixture Substances 0.000 abstract 2
- 238000012067 mathematical method Methods 0.000 abstract 1
- 230000005540 biological transmission Effects 0.000 description 10
- 230000005611 electricity Effects 0.000 description 9
- 238000013461 design Methods 0.000 description 4
- 238000009826 distribution Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000005489 elastic deformation Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 210000002683 foot Anatomy 0.000 description 3
- 238000003491 array Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000000205 computational method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 230000007704 transition 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Abstract
本发明公开了一种分析结果精确的输电铁塔杆件应力计算方法，首先利用有限单元法建立铁塔结构的有限元力学模型，并根据铁塔结构组成进行整体结构的离散化；再对每一单元生成单元刚度矩阵，结合铁塔结构中各杆件之间的空间角度关系和连接关系，叠加生成整体刚度矩阵；然后，根据铁塔所受载荷生成载荷阵列，并以节点位移阵列作为未知量，与整体刚度矩阵，载荷阵列组成矩阵方程；最后，通过求解矩阵方程，获取节点应变，最终得到铁塔结构中每根杆件的应力。本发明通过数学方法直接对输电铁塔杆件应力进行精确求解，可为铁塔结构安全评价提供科学依据。
Description
本申请是申请号201410423136.3、申请日2014.8.26、名称“一种基于有限元分析的输电铁塔杆件应力计算方法”的分案申请。
技术领域
本发明涉及一种对输电铁塔进行结构力学分析并精确计算每根杆件应力的方法，属于输电线路运行状态安全评价领域。
背景技术
电能生产与传输是国家经济社会发展的命脉，输电线路更是关系国计民生的“生命线”。输电线路作为重要的生命线工程，为电力能源的高效供给和合理分配提供了有力保证。经济社会的持续健康发展，对输电线路运行的安全和可靠性提出了更高要求。然而，由于我国幅员辽阔，输电线路路径长、塔架高、大量穿越深山大河，甚至是无人区，自然地理气象条件复杂多变，输电线路作为大电网的主干，其自身是一个巨大的环境灾害承载体。长期运行资料表明线路安全事故大多由机械力学原因造成。虽然在线路规划设计中，已基于铁塔塔材理论计算强度值，按特定工况条件计算允许应力并设置安全系数，在一定程度保证线路的安全运行。但实际运行工况复杂且时变，且无法完全由设计规范所预计，部件实际运行应力可能超出安全阈值；因此实际输电线路虽已严格按规范设计、施工、运行，但断线、倒塔等安全事故仍时有发生。
针对08年特大冰雪灾害，业界和学界专家普遍认为电网缺乏科学有效的防灾减灾事故预警综合体系等是停电事故发生的重要原因之一。建立这一体系的关键是对输电线路的安全状态给予科学评价。当前相关研究主要通过应用拉力传感器，设计在线应力监测装置。由于硬件成本、装置本体易受外力破坏等影响因素，其进一步推广应用受到了制约。开发一种零硬件成本、无环境制约、可复用在所有运行输电线路的结构安全评价系统成为了电网安全运行的当务之急，其中对铁塔杆件进行精确的应力分布计算是其中的关键。
针对输电铁塔的应力分布分析计算，见诸公开报导的技术路线均是对铁塔结构进行多次简化等效，将铁塔所受载荷进行简单的分解或叠加，并忽略铁塔结构的弹性形变。此种分析方法给出的结果，无法为安全评价提供准确地应力数据。从力学分析的角度考虑，不精确。故本发明提出，对铁塔结构进行有限元建模，考虑到铁塔结构的弹性形变，对每一根杆件进行力学分析，并根据铁塔结构对单根杆件的力学性质进行合理叠加，给出整体铁塔结构的精确力学性质，以得到铁塔结构在复杂载荷下的精确应力分布。
发明内容
为了解决上述问题，本发明提出一种基于有限元分析的输电铁塔杆件应力计算方法，用于分析计算铁塔每根杆件的应力，包括如下步骤：
步骤1：根据有限单元法建立铁塔结构的有限元力学模型，根据得到的力学模型获取铁塔节点；并将铁塔结构离散化，得到与铁塔两节点间杆件数目相同的单元，单元彼此之间仅靠节点(铁塔中主材与主材、主材与斜材、斜材与斜材的交汇点视为节点)相连；
步骤2：针对每个单元和节点位移，分别生成单元刚度矩阵[k](e)和节点位移阵列并根据杆件之间的空间角度关系、杆件之间的连接关系，转换叠加出铁塔整体刚度矩阵
步骤3：将铁塔所受均布载荷和非节点载荷等效移置到节点上，形成节点载荷阵列(其中)，并考虑到矩阵方程中的奇异性，引入铁塔塔腿四节点的位移约束条件，求解矩阵方程，得到节点位移阵列
步骤4：计算铁塔节点的应变、应力，获取铁塔结构中薄弱的杆件。
其中，步骤1中的铁塔结构离散化方法是将铁塔中主材与主材、主材与斜材、斜材与斜材的交汇点归结为节点，每两个节点之间的杆件视为一个杆单元。
步骤2中的节点位移阵列和单元刚度矩阵[k](e)分别为：
1)点位移阵列的表达式为：
其中，为整体坐标系中的节点位移阵列；为第1个节点的位移阵列；为第2个节点的位移阵列；以此类推为第n个节点的位移阵列；ui，vi，wi为第i节点在局部坐标系中三个方向的线位移；θxi，θyi，θzi为第i节点处截面绕三个坐标轴的转动，θxi代表截面的扭转，θyi，θzi分别代表截面在xz及xy坐标面内的转动。
2)单元刚度矩阵[k](e)的表达式为：
其中，[k](e)为杆单元在单元局部坐标系内的刚度矩阵；A为杆单元横截面面积；Iy为在xz面内截面惯性矩；Iz为在xy面内的截面惯性矩；Ip为单元的扭转惯性矩；l为长度；E和G分别为材料的弹性模量和剪切模量。
根据所获得的单元刚度矩阵[k](e)并依据杆件之间的空间角度关系、杆件之间的连接关系，即可得到铁塔整体刚度矩阵
步骤3中的节点载荷矩阵为：
其中，为整体坐标中所有节点载荷阵列；为整体坐标中第i个节点的载荷列阵；Nxi为第i个节点的轴向力，Nyi、Nzi分别为第i个节点在xy及xz面内的剪力；Mxi为第i个节点的扭矩，Myi、Mzi为第i个节点在xz及xy面内的弯矩。
考虑到矩阵方程中为奇异矩阵，为使得方程组有解，本发明中引入铁塔4个塔脚与基础连接部分固定，限制铁塔结构的刚性位移这一约束条件，从而保证整体刚度方程具有唯一解。
步骤4中计算铁塔节点的应变和应力公式为：
其中，σx，σy，σz为坐标轴x，y，z方向上的3个正应力分量；τxy，τyz，τzx为在xy平面，yz平面，zx平面内的3个切应力分量；分别为铁塔空间结构在发生形变时，节点处在坐标轴x，y，z三个方向上的产生的线应变分量； 为铁塔空间结构在形变时，节点处在xy平面，yz平面，zx平面内产生的3个剪应变分量。；u＝u(x,y,z)，v＝v(x,y,z)，w＝w(x,y,z)分别为铁塔空间结构在形变时，节点在坐标轴x，y，z方向的位移；E为杆件的弹性模量；μ为泊松比。
本发明的技术效果：
1)充分考虑铁塔结构的弹性形变，对每一根杆件进行力学分析，根据铁塔结构对单根杆件的力学性质进行叠加，使得分析结果更加精确。
2)不仅可分析单一工况下铁塔结构受力，对分析复杂工况下的铁塔结构受力分析也具有优越的性能。
3)计算中引入铁塔4个塔脚与基础连接部分固定，限制铁塔结构的刚性位移这一约束条件，从而保证整体刚度方程具有唯一解，巧妙地解决了矩阵方程中的奇异性矩阵的奇异性问题。
附图说明
图1为220kV干字型铁塔示意图。
图2为铁塔模型中塔脚施加约束条件示意图。
具体实施方式
一种基于有限元分析的输电铁塔杆件应力计算方法，可概括为四个阶段：前期处理、有限元力学分析、工况荷载处理和后期处理。前期处理包括建立铁塔结构的有限元力学模型和对整体结构进行离散化处理；有限元力学分析即对铁塔结构的有限元力学模型进行分析，结合铁塔结构中各杆件之间的空间角度关系和连接关系，叠加生成整体刚度矩阵；工况荷载处理即将铁塔所受均布载荷、非节点载荷等效移置到节点上，形成节点载荷阵列；后期处理即以节点位移阵列作为未知量，与整体刚度矩阵，载荷阵列组成矩阵方程，并求解矩阵方程，给出节点应变，最终找出铁塔结构中薄弱的杆件。该方法主要包括如下步骤：
步骤1：根据有限单元法建立铁塔结构的有限元力学模型，根据得到的力学模型获取铁塔节点；并将铁塔结构离散化，得到与铁塔两节点间杆件数目相同的单元，单元彼此之间仅靠节点(铁塔中主材与主材、主材与斜材、斜材与斜材的交汇点视为节点)相连；
步骤2：针对每个单元和节点位移，分别生成单元刚度矩阵[k](e)和节点位移阵列并根据杆件之间的空间角度关系、杆件之间的连接关系，转换叠加出铁塔整体刚度矩阵
步骤3：将铁塔所受均布载荷和非节点载荷等效移置到节点上，形成节点载荷阵列(其中)，并考虑到矩阵方程中的奇异性，引入铁塔塔腿四节点的位移约束条件，求解矩阵方程，得到节点位移阵列
步骤4：计算铁塔节点的应变、应力，获取铁塔结构中薄弱的杆件。
下面对每个步骤作进一步详细说明：
步骤1中：根据有限单元法建立铁塔结构的有限元力学模型，根据得到的力学模型获取铁塔节点；并将铁塔结构离散化，得到与铁塔两节点间杆件数目相同的单元，单元彼此之间仅靠节点(铁塔中主材与主材、主材与斜材、斜材与斜材的交汇点视为节点)相连，其实施过程为：
首先建立铁塔结构的有限元力学模型，以铁塔的横担方向作为整体坐标系的x轴，线路方向作为y轴，竖直方向作为z轴，并满足右手定则；以杆单元所在直线作为单元局部坐标系的x轴，杆件与局部坐标系下的x轴方向重合，其正方向与整体坐标系x轴正方向一致。
输电铁塔作为空间杆件系统，用杆-梁单元(以后简称杆单元)进行离散，本发明中将铁塔中主材与主材、主材与斜材、斜材与斜材的交汇点化为节点，每两个节点之间的杆件视为一个杆单元。由于铁塔的结构比较复杂，杆件数目比较多，若在结构离散化即划分单元时，单元分的越小，单元数目就越多，计算时间就越长，因此考虑按照铁塔本身的自然结构划分，既提高了计算精度，又减少了计算工作量。
步骤2中：针对每个单元和节点位移，分别生成单元刚度矩阵[k](e)和节点位移阵列并根据杆件之间的空间角度关系、杆件之间的连接关系，转换叠加出铁塔整体刚度矩阵
对铁塔结构离散化后，要对单元进行力学特性分析，即确定单元节点力和节点位移之间的关系。为了分析和确定这一关系，需要选择位移模式，位移函数是单元上点的位移对点的坐标的函数，本方法用单元内部点的坐标的多项式来表示，空间中的杆件，每个节点具有6个自由度，即杆件除了承受一维轴力、两维剪力和两维弯矩的作用外，还可能承受一维扭矩的作用。并且，空间杆单元承受一维轴力、两维剪力、两维弯矩、一维扭矩，即对应着节点的6个自由度，分别为3个方向上的线位移和在节点处截面绕3个坐标轴的转动，因此单元内部点的坐标的多项式可表示为δ＝k1u+k2v+k3w+k4θx+k5θy+k6θz，据此，可形成所有节点的位移阵列
其中，为整体坐标系中的节点位移阵列；为第1个节点的位移阵列；为第2个节点的位移阵列；以此类推为第n个节点的位移阵列；ui，vi，wi为第i节点在局部坐标系中三个方向的线位移；θxi，θyi，θzi为第i节点处截面绕三个坐标轴的转动，θxi代表截面的扭转，θyi，θzi分别代表截面在xz及xy坐标面内的转动。
建立单元刚度方程的基本步骤为：在假定单元位移函数的基础上，根据弹性力学理论，来建立应变、应力与节点位移之间的关系式。然后根据虚位移原理，求得单元节点力与节点位移之间的关系，从而得出如下单元刚度矩阵[k](e)：
其中，[k](e)为杆单元在单元局部坐标系内的刚度矩阵；A为杆单元横截面面积；Iy为在xz面内截面惯性矩；Iz为在xy面内的截面惯性矩；Ip为单元的扭转惯性矩；l为长度；E和G分别为材料的弹性模量和剪切模量。
单元单元刚度矩阵[k](e)得到后，根据杆件之间的空间角度关系、杆件之间的连接关系，转换叠加出铁塔整体刚度矩阵其具体实施过程如下：
首先，假设局部坐标为x，y，z；整体坐标为以从局部坐标系的x轴正方向顺时针转到整体坐标系的x轴正方向为正，则x轴的方向余弦为：
y轴的方向余弦为
z轴的方向余弦为
令
则整体坐标系内的单元刚度矩阵表示为：
其中，[λ0]为节点转换阵，[λ]为坐标转换矩阵，[k](e)为杆单元在单元局部坐标系内的刚度矩阵。
整体刚度矩阵的集成规则如下：
1)先求出每个单元的刚度矩阵[k](e)；
2)将的每个子块进行换号，换成对应的整体编号；
3)将换成以后的字块送到整体刚度矩阵中的对应位置上；
4)若在同一位置上有几个单元的相应子块送到时，则应进行叠加。
经过上述步骤后，就可以得到整体刚度矩阵中的每个子块，从而形成了整体刚度矩阵
其中，为第j个节点产生单位位移时，在第i个节点上引起的节点力，称为刚度子矩阵。值得注意的是：中每个子块都是6×6阶矩阵，如果整体结构具有n个节点，那么整体刚度矩阵的阶数6n×6n。
步骤3将铁塔所受均布载荷和非节点载荷等效移置到节点上，形成节点载荷阵列(其中)，并考虑到矩阵方程中的奇异性，引入铁塔塔腿四节点的位移约束条件，求解矩阵方程，得到节点位移阵列其具体实施过程为：
首先将铁塔所受均布载荷、非节点载荷等效移置到节点上，形成节点载荷阵列空间中的杆件，每个节点具有6个自由度，即杆件除了承受轴力、剪力和弯矩的作用外，还可能承受扭矩的作用。并且，空间杆单元承受一维轴力、两维剪力、两维弯矩、一维扭矩，即对应着节点的6个自由度。输电铁塔的杆单元正是空间杆单元。
其中，为整体坐标中所有节点载荷阵列；为整体坐标中第i个节点的载荷列阵；Nxi为第i个节点的轴向力，Nyi、Nzi分别为第i个节点在xy及xz面内的剪力；Mxi为第i个节点的扭矩，Myi、Mzi为第i个节点在xz及xy面内的弯矩。
由于矩阵方程中为奇异矩阵，方程组无解，若要求解该方程，必须引入约束条件，限制铁塔结构的刚性位移，保证整体刚度方程有唯一解。位移约束条件的作用是使结构上的节点的位移分量为常数值，即δi＝δ0。引入位移约束条件，就是要将δi＝δ0引入到结构总体刚度方程中。
本发明中，采用对角元素置1法，将δi＝δ0引入整体刚度矩阵
针对输电铁塔的4个塔腿中，与基础连接的部分是固定端约束，因此δ0＝0；将K的第i行的主对角线元素Kii置1，其余元素清零，且将第i行的载荷项Ri用δ0代替，上式变为
通过置1法将位移约束条件δi＝δ0引入到整体刚度方程中，并没有改变矩阵K和R中的各元素储存顺序，而且矩阵K仍然为对称矩阵。
本发明计算方法中，对与基础相连的4个塔腿节点的子矩阵采用置1法，即可代入24个位移边界条件，消除整体刚度矩阵的奇异性，从而采用高斯消元法进行矩阵方程求解。
步骤4：计算节点应变、应力等，找出铁塔结构中最薄弱的杆件，其实施过程如下：
由铁塔结构的整体刚度矩阵方程求解出各节点的位移δ后，就可以得到各单元的节点位移δe。铁塔结构受力后，其内部各点将沿x，y，z三个坐标轴方向发生位移。如果各点沿x，y，z三个坐标轴方向的位移以u，v，w表示，它们是点的坐标函数，即u＝u(x,y,z)，v＝v(x,y,z)，w＝w(x,y,z)。
铁塔结构在形变时，内部任意一点处有3个线应变分量εx，εy，εz以及3对剪应变分量γxy＝γyx，γyz＝γzy，γzx＝γxz。由弹性力学可知，应变与位移之间的关系即几何方程为：
式中εx，εy，εz为分别为铁塔空间结构在发生形变时，节点处在坐标轴x，y，z三个方向上的产生的应变分量；γxy，γyz，γzx为铁塔空间结构在形变时，节点处在xy平面，yz平面，zx平面内产生的的3个应变分量。
铁塔结构受到作用时，内部任意一点处的应力状态也是三维的，有3个正应力分量σx，σy，σz以及三对切应力分量τxy＝τyx，τyz＝τzy，τzx＝τxz。
在线弹性范围内，应力和应变用如下方程表示：
其中，σx，σy，σz为坐标轴x，y，z方向上的3个正应力分量；τxy，τyz，τzx为在xy平面，yz平面，zx平面内的3个切应力分量；分别为铁塔空间结构在发生形变时，节点处在坐标轴x，y，z三个方向上的产生的线应变分量； 为铁塔空间结构在形变时，节点处在xy平面，yz平面，zx平面3个平面内产生的剪应变分量；E为杆件的弹性模量；μ为泊松比。
在本说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (1)
1.一种分析结果精确的输电铁塔杆件应力计算方法，其特征在于：包括如下步骤：
步骤1：根据有限单元法建立铁塔结构的有限元力学模型，根据得到的力学模型获取铁塔节点；并将铁塔结构离散化，得到与铁塔两节点间杆件数目相同的单元，单元彼此之间仅靠节点相连；铁塔中主材与主材、主材与斜材、斜材与斜材的交汇点视为节点；
步骤2：针对每个单元和节点位移，分别生成单元刚度矩阵[k](e)和节点位移阵列并根据杆件之间的空间角度关系、杆件之间的连接关系，转换叠加出铁塔整体刚度矩阵
步骤3：将铁塔所受均布载荷和非节点载荷等效移置到节点上，形成节点载荷阵列其中并考虑到矩阵方程中的奇异性，引入铁塔塔腿四节点的位移约束条件，求解矩阵方程，得到节点位移阵列
步骤4：计算铁塔节点的应变、应力，获取铁塔结构每根杆件应力；
步骤2中：针对每个单元和节点位移，分别生成单元刚度矩阵[k](e)和节点位移阵列并根据杆件之间的空间角度关系、杆件之间的连接关系，转换叠加出铁塔整体刚度矩阵
对铁塔结构离散化后，要对单元进行力学特性分析，即确定单元节点力和节点位移之间的关系；为了分析和确定这一关系，需要选择位移模式，位移函数是单元上点的位移对点的坐标的函数，本方法用单元内部点的坐标的多项式来表示，空间中的杆件，每个节点具有6个自由度，即杆件除了承受一维轴力、两维剪力和两维弯矩的作用外，还可能承受一维扭矩的作用；并且，空间杆单元承受一维轴力、两维剪力、两维弯矩、一维扭矩，即对应着节点的6个自由度，分别为3个方向上的线位移和在节点处截面绕3个坐标轴的转动，因此单元内部点的坐标的多项式可表示为δ＝k1u+k2v+k3w+k4θx+k5θy+k6θz，其中，u，v，w分别为3个方向上的线位移；θx，θy，θz分别为节点处截面绕3个坐标轴的转动；k1,k2，…，k6为6个自由度的系数；据此，可形成所有节点的位移阵列
其中，为整体坐标系中的节点位移阵列；为第1个节点的位移阵列；为第2个节点的位移阵列；以此类推为第n个节点的位移阵列；ui，vi，wi为第i节点在局部坐标系中三个方向的线位移；θxi，θyi，θzi为第i节点处截面绕三个坐标轴的转动，θxi代表截面在yz坐标面内的扭转，θyi，θzi分别代表截面在xz及xy坐标面内的转动；
建立单元刚度方程的基本步骤为：在假定单元位移函数的基础上，根据弹性力学理论，来建立应变、应力与节点位移之间的关系式；然后根据虚位移原理，求得单元节点力与节点位移之间的关系，从而得出如下单元刚度矩阵[k](e)：
其中，[k](e)为杆单元在单元局部坐标系内的刚度矩阵；A为杆单元横截面面积；Iy为在xz面内截面惯性矩；Iz为在xy面内的截面惯性矩；Ip为单元的扭转惯性矩；l为长度；E和G分别为材料的弹性模量和剪切模量；
单元单元刚度矩阵[k](e)得到后，根据杆件之间的空间角度关系、杆件之间的连接关系，转换叠加出铁塔整体刚度矩阵其具体实施过程如下：
首先，假设局部坐标为x，y，z；整体坐标为以从局部坐标系的x轴正方向顺时针转到整体坐标系的轴正方向为正，则x轴的方向余弦为：
<mrow>
<msub>
<mi>l</mi>
<mrow>
<mi>x</mi>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</msub>
<mo>=</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>x</mi>
<mover>
<mi>y</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</msub>
<mo>=</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mover>
<mi>y</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>x</mi>
<mover>
<mi>z</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</msub>
<mo>=</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mover>
<mi>z</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
</mrow>
y轴的方向余弦为
<mrow>
<msub>
<mi>l</mi>
<mrow>
<mi>y</mi>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</msub>
<mo>=</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>,</mo>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>y</mi>
<mover>
<mi>y</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</msub>
<mo>=</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>,</mo>
<mover>
<mi>y</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>y</mi>
<mover>
<mi>z</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</msub>
<mo>=</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>,</mo>
<mover>
<mi>z</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
</mrow>
z轴的方向余弦为
<mrow>
<msub>
<mi>l</mi>
<mrow>
<mi>z</mi>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</msub>
<mo>=</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>,</mo>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>z</mi>
<mover>
<mi>y</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</msub>
<mo>=</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>,</mo>
<mover>
<mi>y</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>z</mi>
<mover>
<mi>z</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</msub>
<mo>=</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>,</mo>
<mover>
<mi>z</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
</mrow>
令
<mrow>
<msub>
<mrow>
<mo>&lsqb;</mo>
<mi>&lambda;</mi>
<mo>&rsqb;</mo>
</mrow>
<mn>01</mn>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>l</mi>
<mrow>
<mi>x</mi>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>l</mi>
<mrow>
<mi>x</mi>
<mover>
<mi>y</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>l</mi>
<mrow>
<mi>x</mi>
<mover>
<mi>z</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>l</mi>
<mrow>
<mi>y</mi>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>l</mi>
<mrow>
<mi>y</mi>
<mover>
<mi>y</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>l</mi>
<mrow>
<mi>y</mi>
<mover>
<mi>z</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>l</mi>
<mrow>
<mi>z</mi>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>l</mi>
<mrow>
<mi>z</mi>
<mover>
<mi>y</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>l</mi>
<mrow>
<mi>z</mi>
<mover>
<mi>z</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
则整体坐标系内的单元刚度矩阵表示为：
<mrow>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mover>
<mi>k</mi>
<mo>&OverBar;</mo>
</mover>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>e</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mi>&lambda;</mi>
<mo>&rsqb;</mo>
</mrow>
<mi>T</mi>
</msup>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mi>k</mi>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>e</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>&lsqb;</mo>
<mi>&lambda;</mi>
<mo>&rsqb;</mo>
</mrow>
其中，[λ0]为节点转换阵，[λ]为坐标转换矩阵，[k](e)为杆单元在单元局部坐标系内的刚度矩阵；
整体刚度矩阵的集成规则如下：
1)先求出每个单元的刚度矩阵[k](e)；
2)将的每个子块进行换号，换成对应的整体编号；
3)将换成以后的字块送到整体刚度矩阵中的对应位置上；
4)若在同一位置上有几个单元的相应子块送到时，则应进行叠加；
经过上述步骤后，就可以得到整体刚度矩阵中的每个子块，从而形成了整体刚度矩阵
其中，为第j个节点产生单位位移时，在第i个节点上引起的节点力，称为刚度子矩阵；值得注意的是：中每个子块都是6×6阶矩阵，如果整体结构具有n个节点，那么整体刚度矩阵的阶数6n×6n。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510423172.4A CN105005662B (zh) | 2014-08-26 | 2014-08-26 | 分析结果精确的输电铁塔杆件应力计算方法 |
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510423172.4A CN105005662B (zh) | 2014-08-26 | 2014-08-26 | 分析结果精确的输电铁塔杆件应力计算方法 |
| CN201410423136.3A CN104281739B (zh) | 2014-08-26 | 2014-08-26 | 一种基于有限元分析的输电铁塔杆件应力计算方法 |
Related Parent Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410423136.3A Division CN104281739B (zh) | 2014-08-26 | 2014-08-26 | 一种基于有限元分析的输电铁塔杆件应力计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105005662A CN105005662A (zh) | 2015-10-28 |
| CN105005662B true CN105005662B (zh) | 2017-12-01 |
Family
ID=52256609
Family Applications (13)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710092257.8A Pending CN107220396A (zh) | 2014-08-26 | 2014-08-26 | 解决奇异性矩阵奇异性问题的输电铁塔杆件应力计算方法 |
| CN201510422952.7A Active CN105069191B (zh) | 2014-08-26 | 2014-08-26 | 充分考虑铁塔结构弹性形变的输电铁塔杆件应力计算方法 |
| CN201410423136.3A Active CN104281739B (zh) | 2014-08-26 | 2014-08-26 | 一种基于有限元分析的输电铁塔杆件应力计算方法 |
| CN201710092256.3A Pending CN107423457A (zh) | 2014-08-26 | 2014-08-26 | 分析精确的基于有限元分析的输电铁塔杆件应力计算方法 |
| CN201710092218.8A Pending CN107391780A (zh) | 2014-08-26 | 2014-08-26 | 考虑铁塔结构的弹性形变的输电铁塔杆件应力计算方法 |
| CN201510423536.9A Active CN105095576B (zh) | 2014-08-26 | 2014-08-26 | 一种输电铁塔杆件应力计算方法 |
| CN201710092245.5A Pending CN107391781A (zh) | 2014-08-26 | 2014-08-26 | 基于有限元分析的输电铁塔杆件应力计算方法 |
| CN201710092243.6A Pending CN107423456A (zh) | 2014-08-26 | 2014-08-26 | 一种输电铁塔杆件应力计算方法 |
| CN201710092231.3A Pending CN107229765A (zh) | 2014-08-26 | 2014-08-26 | 分析计算铁塔每根杆件应力的输电铁塔杆件应力计算方法 |
| CN201710092220.5A Pending CN107220395A (zh) | 2014-08-26 | 2014-08-26 | 充分考虑铁塔结构弹性形变的输电铁塔杆件应力计算方法 |
| CN201510423172.4A Active CN105005662B (zh) | 2014-08-26 | 2014-08-26 | 分析结果精确的输电铁塔杆件应力计算方法 |
| CN201710092219.2A Pending CN107357945A (zh) | 2014-08-26 | 2014-08-26 | 性能优越的基于有限元分析的输电铁塔杆件应力计算方法 |
| CN201510423520.8A Active CN105045998B (zh) | 2014-08-26 | 2014-08-26 | 可分析复杂工况的输电铁塔杆件应力计算方法 |
Family Applications Before (10)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710092257.8A Pending CN107220396A (zh) | 2014-08-26 | 2014-08-26 | 解决奇异性矩阵奇异性问题的输电铁塔杆件应力计算方法 |
| CN201510422952.7A Active CN105069191B (zh) | 2014-08-26 | 2014-08-26 | 充分考虑铁塔结构弹性形变的输电铁塔杆件应力计算方法 |
| CN201410423136.3A Active CN104281739B (zh) | 2014-08-26 | 2014-08-26 | 一种基于有限元分析的输电铁塔杆件应力计算方法 |
| CN201710092256.3A Pending CN107423457A (zh) | 2014-08-26 | 2014-08-26 | 分析精确的基于有限元分析的输电铁塔杆件应力计算方法 |
| CN201710092218.8A Pending CN107391780A (zh) | 2014-08-26 | 2014-08-26 | 考虑铁塔结构的弹性形变的输电铁塔杆件应力计算方法 |
| CN201510423536.9A Active CN105095576B (zh) | 2014-08-26 | 2014-08-26 | 一种输电铁塔杆件应力计算方法 |
| CN201710092245.5A Pending CN107391781A (zh) | 2014-08-26 | 2014-08-26 | 基于有限元分析的输电铁塔杆件应力计算方法 |
| CN201710092243.6A Pending CN107423456A (zh) | 2014-08-26 | 2014-08-26 | 一种输电铁塔杆件应力计算方法 |
| CN201710092231.3A Pending CN107229765A (zh) | 2014-08-26 | 2014-08-26 | 分析计算铁塔每根杆件应力的输电铁塔杆件应力计算方法 |
| CN201710092220.5A Pending CN107220395A (zh) | 2014-08-26 | 2014-08-26 | 充分考虑铁塔结构弹性形变的输电铁塔杆件应力计算方法 |
Family Applications After (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710092219.2A Pending CN107357945A (zh) | 2014-08-26 | 2014-08-26 | 性能优越的基于有限元分析的输电铁塔杆件应力计算方法 |
| CN201510423520.8A Active CN105045998B (zh) | 2014-08-26 | 2014-08-26 | 可分析复杂工况的输电铁塔杆件应力计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (13) | CN107220396A (zh) |
Families Citing this family (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104636603B (zh) * | 2015-01-20 | 2017-08-29 | 华北电力大学（保定） | 一种特高压单柱拉线塔扭振频率的两自由度计算方法 |
| CN106326504B (zh) * | 2015-06-16 | 2020-01-31 | 中国电力科学研究院 | 一种输电杆塔空间模型塔材定位方法 |
| CN106326505B (zh) * | 2015-06-16 | 2020-01-24 | 中国电力科学研究院 | 一种平面节点处理方法 |
| CN106777815A (zh) * | 2017-01-18 | 2017-05-31 | 江苏龙胜机床制造有限公司 | 基于有限元的数控折弯机机架分析方法 |
| CN107451371B (zh) * | 2017-08-11 | 2021-01-26 | 四川大学 | 提高三角单元计算精度的方法 |
| CN107944192A (zh) * | 2017-12-14 | 2018-04-20 | 青岛海西重机有限责任公司 | 一种组合连杆应力自动分析方法及装置 |
| CN109583025B (zh) * | 2018-10-30 | 2023-06-16 | 国网江西省电力有限公司电力科学研究院 | 一种覆冰腐蚀地线承载力分析方法及系统 |
| CN110750925B (zh) * | 2019-09-27 | 2024-02-06 | 哈尔滨工业大学（深圳） | 柔性机械手力变形模型建模方法及装置 |
| CN111339704B (zh) * | 2020-02-28 | 2023-07-18 | 四川电力设计咨询有限责任公司 | 输电铁塔错心节点的强度设计方法 |
| CN111753404B (zh) * | 2020-06-02 | 2021-02-09 | 中煤科工开采研究院有限公司 | 薄松散层采煤沉陷区双层介质地基附加应力计算方法 |
| CN112084683B (zh) * | 2020-07-22 | 2022-07-08 | 国网江西省电力有限公司电力科学研究院 | 一种输电线路用铝合金材料金具的塑性变形分析方法 |
| CN112182689B (zh) * | 2020-08-27 | 2023-03-31 | 山东电力工程咨询院有限公司 | 一种输电塔主材内力的监测方法及系统 |
| CN112363372B (zh) * | 2020-11-19 | 2023-03-10 | 东方晶源微电子科技（北京）有限公司深圳分公司 | 一种负显影光刻工艺的仿真方法、负显影光刻胶模型、opc模型及电子设备 |
| CN114638046B (zh) * | 2022-05-12 | 2022-08-09 | 中国铁路设计集团有限公司 | 一种铁路桥墩数字孪生变截面仿真计算方法 |
| CN115062432B (zh) * | 2022-06-29 | 2023-09-08 | 国网甘肃省电力公司电力科学研究院 | 一种基于数据驱动的输电杆塔倾斜预警方法及系统 |
| CN115758751B (zh) * | 2022-11-22 | 2023-09-26 | 中铁工程设计咨询集团有限公司 | 一种无缝线路设计方法、装置、设备及可读存储介质 |
| CN117057161B (zh) * | 2023-10-10 | 2024-01-30 | 清华大学深圳国际研究生院 | 一种基于数字孪生的风机安装力学分析方法及系统 |
| CN117609739B (zh) * | 2024-01-19 | 2024-04-05 | 北京云摩科技股份有限公司 | 基于多点形变数据联合分析的结构在线监测方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101446517A (zh) * | 2008-12-17 | 2009-06-03 | 中国电力科学研究院 | 一种输电线路高塔结构振动测试方法 |
| CN101979974A (zh) * | 2010-09-07 | 2011-02-23 | 东南大学 | 一种输电塔全场应力监测方法 |
| CN102254061A (zh) * | 2011-06-09 | 2011-11-23 | 山东电力研究院 | 输电塔线塔耦合体系的有限元建模及受力方法 |
| CN102968542A (zh) * | 2012-12-12 | 2013-03-13 | 河北省电力勘测设计研究院 | 应用ansys软件进行输电铁塔结构分析的方法 |
Family Cites Families (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2002041195A2 (en) * | 2000-11-17 | 2002-05-23 | Battelle Memorial Institute | Structural stress analysis |
| WO2007012200A1 (en) * | 2005-07-25 | 2007-02-01 | The University Of Manitoba | Design system for composite wind towers |
| WO2008105221A1 (ja) * | 2007-02-28 | 2008-09-04 | Keio University | 構造解析数値計算装置 |
| JP4253822B2 (ja) * | 2007-07-18 | 2009-04-15 | 株式会社テクニカルスルー | 粒子法の解析プログラム |
| US8209134B2 (en) * | 2008-12-04 | 2012-06-26 | Laura P. Solliday | Methods for modeling the structural health of a civil structure based on electronic distance measurements |
| CN101634619A (zh) * | 2009-08-20 | 2010-01-27 | 中国电力科学研究院 | 一种特高压钢管塔内力计算方法及装置 |
| CN101847855B (zh) * | 2010-05-25 | 2013-08-21 | 中国电力工程顾问集团中南电力设计院 | 架空输电线路铁塔抗导线舞动的方法 |
| CN102542125A (zh) * | 2012-03-15 | 2012-07-04 | 中国电力科学研究院 | 一种输电塔组装过程模拟方法 |
| CN103105257A (zh) * | 2013-01-29 | 2013-05-15 | 华北电力大学 | 一种基于动力特性的覆冰致输电塔线结构失效的预警方法 |
| CN103310072B (zh) * | 2013-06-28 | 2015-12-23 | 哈尔滨理工大学 | 基于力反馈的股骨生物力学有限元分析系统 |
| CN103455686A (zh) * | 2013-09-17 | 2013-12-18 | 重庆大学 | 架空输电塔-线耦合体系有限元模型的建模方法 |
| CN103838913B (zh) * | 2013-12-25 | 2017-03-08 | 广西科技大学 | 曲线箱梁弯桥的有限单元法 |
| CN103914590B (zh) * | 2014-03-12 | 2017-02-01 | 国家电网公司 | 一种电力铁塔三维实体模型生成方法 |
-
2014
- 2014-08-26 CN CN201710092257.8A patent/CN107220396A/zh active Pending
- 2014-08-26 CN CN201510422952.7A patent/CN105069191B/zh active Active
- 2014-08-26 CN CN201410423136.3A patent/CN104281739B/zh active Active
- 2014-08-26 CN CN201710092256.3A patent/CN107423457A/zh active Pending
- 2014-08-26 CN CN201710092218.8A patent/CN107391780A/zh active Pending
- 2014-08-26 CN CN201510423536.9A patent/CN105095576B/zh active Active
- 2014-08-26 CN CN201710092245.5A patent/CN107391781A/zh active Pending
- 2014-08-26 CN CN201710092243.6A patent/CN107423456A/zh active Pending
- 2014-08-26 CN CN201710092231.3A patent/CN107229765A/zh active Pending
- 2014-08-26 CN CN201710092220.5A patent/CN107220395A/zh active Pending
- 2014-08-26 CN CN201510423172.4A patent/CN105005662B/zh active Active
- 2014-08-26 CN CN201710092219.2A patent/CN107357945A/zh active Pending
- 2014-08-26 CN CN201510423520.8A patent/CN105045998B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101446517A (zh) * | 2008-12-17 | 2009-06-03 | 中国电力科学研究院 | 一种输电线路高塔结构振动测试方法 |
| CN101979974A (zh) * | 2010-09-07 | 2011-02-23 | 东南大学 | 一种输电塔全场应力监测方法 |
| CN102254061A (zh) * | 2011-06-09 | 2011-11-23 | 山东电力研究院 | 输电塔线塔耦合体系的有限元建模及受力方法 |
| CN102968542A (zh) * | 2012-12-12 | 2013-03-13 | 河北省电力勘测设计研究院 | 应用ansys软件进行输电铁塔结构分析的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105069191B (zh) | 2018-03-16 |
| CN105069191A (zh) | 2015-11-18 |
| CN107423457A (zh) | 2017-12-01 |
| CN107220395A (zh) | 2017-09-29 |
| CN104281739A (zh) | 2015-01-14 |
| CN105095576A (zh) | 2015-11-25 |
| CN107229765A (zh) | 2017-10-03 |
| CN105045998A (zh) | 2015-11-11 |
| CN107357945A (zh) | 2017-11-17 |
| CN104281739B (zh) | 2017-05-10 |
| CN105005662A (zh) | 2015-10-28 |
| CN107423456A (zh) | 2017-12-01 |
| CN107220396A (zh) | 2017-09-29 |
| CN107391781A (zh) | 2017-11-24 |
| CN107391780A (zh) | 2017-11-24 |
| CN105045998B (zh) | 2017-12-26 |
| CN105095576B (zh) | 2017-11-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105005662B (zh) | 分析结果精确的输电铁塔杆件应力计算方法 | |
| CN104298808B (zh) | 输电铁塔非线性柔性构件的应力计算方法 | |
| CN104091033B (zh) | 基于超单元结合虚拟变形法的桥梁静力有限元模型修正方法 | |
| CN104408286B (zh) | 基于正交分解理论的平面构件变形分解与振型识别方法 | |
| CN106021709A (zh) | 一种混凝土早期开裂风险的评估和控制方法 | |
| Li et al. | Experimental and numerical seismic investigations of the Three Gorges dam | |
| Gao et al. | Rock creep modeling based on discontinuous deformation analysis | |
| Gizejowski et al. | Beam-column in-plane resistance based on the concept of equivalent geometric imperfections | |
| Thinh et al. | Bending and vibration analysis of multi-folding laminate composite plate using finite element method | |
| Liao et al. | Second-order Effect Finite Element Analysis of Stainless Steel Frame Structure | |
| CN106339545A (zh) | 一种基于非连续变形分析的岩体流变模拟方法 | |
| Li et al. | Study on Algorithm and Parameters of Seismic Design Code for Medium and Long Span Bridges | |
| Jin et al. | Seismic response analysis of long span cable-stayed bridge by response spectrum method | |
| Wang et al. | Hysteresis curve simulation of RC frame based on finite element software ABAQUS | |
| Su et al. | Bond behavior of recycled coarse aggregate concrete with rebar after freeze–thaw cycles: Finite element nonlinear analysis | |
| CN105117586A (zh) | 一种基于直接微分法的输电导线响应敏感性分析方法 | |
| Xue et al. | The structure statics analysis and optimal design of hydraulic support | |
| Zhang et al. | Simulation of Multiscale Time and Space in Fracture Process of Three-Point-Bending Concrete Beam with Offset Notch | |
| Song et al. | Mechanical characteristics analysis and design optimization of the derrick | |
| Alves et al. | Modeling of the Thermo-Elastoplastic Behavior of Space Trusses | |
| Zhang et al. | Universal Format of Shape Function for Numerical Analysis using Multiple Element Forms | |
| Zhou et al. | Simulation Method on Progressive Collapse of Shell Model under Severe Earthquake | |
| Li et al. | FIBER BEAM-COLUMN ELEMENT FOR THE SEISMIC DAMAGE ANALYSIS OF RC FRAMES | |
| Yan | Overview of the status of research on the behavior of beam-column end plate connection joints of steel frame | |
| Manouchehri | Nonlinear seismic behaviour of concrete arch bridges |
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 |