CN107871049A - 曲面单层网壳梁单元欧拉角的计算方法 - Google Patents

曲面单层网壳梁单元欧拉角的计算方法 Download PDF

Info

Publication number
CN107871049A
CN107871049A CN201711194465.5A CN201711194465A CN107871049A CN 107871049 A CN107871049 A CN 107871049A CN 201711194465 A CN201711194465 A CN 201711194465A CN 107871049 A CN107871049 A CN 107871049A
Authority
CN
China
Prior art keywords
node
vector
curved surface
layer lattice
local coordinate
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.)
Granted
Application number
CN201711194465.5A
Other languages
English (en)
Other versions
CN107871049B (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.)
Shenzhen University Architectural Design And Research Institute Co Ltd
Original Assignee
Shenzhen University Architectural Design And Research Institute Co Ltd
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 Shenzhen University Architectural Design And Research Institute Co Ltd filed Critical Shenzhen University Architectural Design And Research Institute Co Ltd
Priority to CN201711194465.5A priority Critical patent/CN107871049B/zh
Publication of CN107871049A publication Critical patent/CN107871049A/zh
Application granted granted Critical
Publication of CN107871049B publication Critical patent/CN107871049B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • 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
    • 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/30—Polynomial surface description

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Computer Hardware Design (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Civil Engineering (AREA)
  • Architecture (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Mathematical Physics (AREA)
  • Image Generation (AREA)

Abstract

本发明公开了一种曲面单层网壳梁单元欧拉角的计算方法，包括如下步骤：将曲面单层网壳进行三角化；建立节点与坐标文件node.txt及element.txt；找出与某一节点相连的节点集合，集合中的每一个节点与该节点组成了梁单元，对任意梁单元，求定向第三点坐标矢量；假定全局坐标系中Z轴对应的矢量为Z＝[0,0,1]，求局部坐标系各轴矢量；求得空间梁的欧拉角BETA为局部坐标轴z与主平面为铅直面的空间梁之局部坐标轴z0之间的角度。本发明，可适用计算机编程运算，计算精度和计算效率均很高，可满足大规模曲面单层网壳设计的需求，适用于多种复杂曲面的处理，可处理大规模曲面单层网壳梁单元欧拉角的计算问题。

Description

曲面单层网壳梁单元欧拉角的计算方法
技术领域
本发明涉及，具体涉及曲面单层网壳梁单元欧拉角的计算方法。
背景技术
曲面单层网壳具有适用跨度大、美观等优点，因此该类型结构广泛用于大跨公共建筑屋盖形式。但曲面单层网壳中梁单元截面形式采用箱型、工字型等时，存在空间梁单元定向问题，而目前此问题一般采用梁单元第三点定向或梁单元欧拉角定向来处理。
以往求梁单元第三点或梁单元欧拉角均是采用手工算法，费工费时，效率极低，且容易出错，不适用大规模曲面单层网壳设计的需求。对曲面单层网壳，若要处理空间梁单元定向问题，除上述手工计算梁单元第三点或梁单元欧拉角外，还有两种处理方案：
(1)若梁单元采用圆形截面，对圆形截面，任一方向其抗弯性能均是相同的，但圆形截面的梁单元的抗弯性能不佳，且难以满足建筑美观的要求，这种处理方法只适用于梁单元采用圆形截面的方式，局限性太大。
(2)对空间曲面网壳，先求其曲面函数方程，Z＝f(x,y)，然后，对曲面中梁单元的中点，求曲面的法向，再求相应的欧拉角。这种方法仅适用球面、柱面等较为规则的类类型，而对更为复杂的任意曲面，难以实现。
由此可见，目前的曲面单层网壳的梁单元的定向问题存在不适用任意曲面的问题。
发明内容
本发明所要解决的技术问题是目前的曲面单层网壳的梁单元的定向问题存在不适用任意曲面的问题。
为了解决上述技术问题，本发明所采用的技术方案是提供了一种曲面单层网壳梁单元欧拉角的计算方法，包括如下步骤：
S1、将曲面单层网壳进行三角化，即增加适当的虚梁，使得网格均为三角形网格，保证虚梁处在最接近曲面的面上；
S2、建立节点与坐标文件node.txt及element.txt；
S3、找出与某一节点相连的节点集合，集合中的每一个节点与该节点组成了梁单元，对任意梁单元，p1、p2为梁单元两节点坐标矢量，求定向第三点坐标矢量p3；
S4、假定全局坐标系中Z轴对应的矢量为Z＝[0,0,1]，求局部坐标系各轴矢量xr、xr0、yr、yr0、zr、zr0；
S5、空间梁的欧拉角BETA为局部坐标轴z与主平面为铅直面的空间梁之局部坐标轴z0之间的角度，也是将空间梁主平面绕梁自身轴转成铅直平面所需的转角，其数学表达如下：
若x-z平面法向(y轴)朝上，则BETA＝180/pi*acos(dot(zr0,zr))；
若x-z平面法向(y轴)朝下，则：
BETA＝360-180/pi*acos(dot(zr0,zr))；
式中dot()为矢量点积运算函数，acos()为反余弦函数。
在另一个优选的实施例中，S2具体包括以下步骤：
S21、NODE.txt为：
S22、ELEMENT.txt为：
在另一个优选的实施例中，S3中，其第一节点编号为NO1,其第二节点编号为NO2，并具体包括以下步骤：
S31、通过ELEMENT.txt，找出与节点NO1相连的节点集合S1、与节点NO2相连的节点集合S2，即S1中的每一个节点与节点NO1组成了梁单元，S2中的每一个节点与节点NO2组成了梁单元；
S32、对S1与S2做交集运算得其结果s12，即s12＝S1∩S2；
s12为：
一个节点(NOa)，此梁必是边梁，此时，梁的法向n取三角形NO1,NO2,NOa的法向n1，即n＝n1；
或两个节点(Noa,Nob)，此梁必是中梁，其法向n与两个三角形相关,且梁的法向n取两个三角形法向的平均值，即n＝(n1+n2)/2；
则定向第三点坐标矢量p3＝(p1+p2)/2+1*n，(p1、p2为梁单元两节点坐标矢量)。
在另一个优选的实施例中，S4具体包括以下步骤：
(1)当梁处在铅直平面中，
梁的局部坐标轴对应的矢量xr0由下式求得：
xr＝p2-p1；xr0＝xt/norm(xr)；
梁的局部坐标轴对应的矢量yr0由下式求得：
yr0＝cross(Z,xr0)；
梁的局部坐标轴对应的矢量zr0由下式求得：
zr0＝cross(xr0,yr0)；
(2)当梁不处在铅直平面中，
梁的局部坐标轴对应的矢量xr0由下式求得：
xr＝p2-p1；xr0＝xt/norm(xr)；norm为求矢量范数的函数；
zrt＝p3-p1；
梁的局部坐标轴对应的矢量yr由下式求得：
yr＝cross(zrt,xr)；yr＝yr/norm(yr)；
梁的局部坐标轴对应的矢量zr由下式求得：
zr＝cross(xr,yr)；
其中，上述公式中norm()为求矢量范数的函数；cross()为矢量叉积运算函数。
本发明，可适用计算机编程运算，计算精度和计算效率均很高，可满足大规模曲面单层网壳设计的需求，适用于多种复杂曲面的处理，可处理大规模曲面单层网壳梁单元欧拉角的计算问题。
附图说明
图1为本发明的梁单元的欧拉角的含义示意图；
图2为本发明的正确计算欧拉角后的曲面单层网壳效果示意图。
具体实施方式
本发明提供了一种曲面单层网壳梁单元欧拉角的计算方法，可适用计算机编程运算，计算精度和计算效率均很高，可满足大规模曲面单层网壳设计的需求，适用于多种复杂曲面的处理，可处理大规模曲面单层网壳梁单元欧拉角的计算问题。下面结合具体实施例和说明书附图对本发明予以详细说明。
本发明提供的曲面单层网壳梁单元欧拉角的计算方法包括如下步骤：
S1、如图2所示，将曲面单层网壳进行三角化，即增加适当的虚梁(刚度接近0的梁)，使得网格均为三角形网格，保证虚梁处在最接近曲面的面上，否则，应对梁单元进行分割并增加相应的虚梁；
S2、建立节点与坐标文件node.txt及element.txt，具体地，可采用如下方式建立：
S21、NODE.txt为：
S22、ELEMENT.txt为：
S3、找出与某一节点相连的节点集合，集合中的每一个节点与该节点组成了梁单元，对任意梁单元，p1、p2为梁单元两节点坐标矢量，求定向第三点坐标矢量p3，
具体地，采用如下计算步骤：
任意梁单元，其第一节点编号为NO1,其第二节点编号为NO2。
S31、通过ELEMENT.txt，找出与节点NO1相连的节点集合S1、与节点NO2相连的节点集合S2，即S1中的每一个节点与节点NO1组成了梁单元，S2中的每一个节点与节点NO2组成了梁单元，
S32、对S1与S2做交集运算得其结果s12，即s12＝S1∩S2，s12为：
一个节点(NOa)，此梁必是边梁，此时，梁的法向n取三角形NO1,NO2,NOa的法向n1，即n＝n1；
或两个节点(Noa,Nob)，此梁必是中梁，其法向n与两个三角形相关,且梁的法向n取两个三角形法向的平均值，即n＝(n1+n2)/2；
则定向第三点坐标矢量p3＝(p1+p2)/2+1*n(p1、p2为梁单元两节点坐标矢量)；
S4、假定全局坐标系中Z轴对应的矢量为Z＝[0,0,1]，求局部坐标系各轴矢量xr、xr0、yr、yr0、zr、zr0，
具体地，采用如下计算步骤：
(1)当梁处在铅直平面中，
梁的局部坐标轴对应的矢量xr0由下式求得：
xr＝p2-p1；xr0＝xt/norm(xr)；
梁的局部坐标轴对应的矢量yr0由下式求得：
yr0＝cross(Z,xr0)；
梁的局部坐标轴对应的矢量zr0由下式求得：
zr0＝cross(xr0,yr0)；
(2)当梁不处在铅直平面中，
梁的局部坐标轴对应的矢量xr0由下式求得：
xr＝p2-p1；xr0＝xt/norm(xr)；
zrt＝p3-p1；
梁的局部坐标轴对应的矢量yr由下式求得：
yr＝cross(zrt,xr)；yr＝yr/norm(yr)；
梁的局部坐标轴对应的矢量zr由下式求得：
zr＝cross(xr,yr)；
其中，上述公式中norm()为求矢量范数的函数；cross()为矢量叉积运算函数。
S5、如图1所示，空间梁的欧拉角BETA为局部坐标轴z与主平面为铅直面的空间梁之局部坐标轴z0之间的角度，也是将空间梁主平面绕梁自身轴转成铅直平面所需的转角，其数学表达如下：
若x-z平面法向(y轴)朝上，则BETA＝180/pi*acos(dot(zr0,zr))；
若x-z平面法向(y轴)朝下，则：
BETA＝360-180/pi*acos(dot(zr0,zr))；
式中dot()为矢量点积运算函数，acos()为反余弦函数。
本发明，可适用计算机编程运算，计算精度和计算效率均很高，可满足大规模曲面单层网壳设计的需求，适用于多种复杂曲面的处理，可处理大规模曲面单层网壳梁单元欧拉角的计算问题。
本发明不局限于上述最佳实施方式，任何人应该得知在本发明的启示下作出的结构变化，凡是与本发明具有相同或相近的技术方案，均落入本发明的保护范围之内。

Claims (4)

1.曲面单层网壳梁单元欧拉角的计算方法，其特征在于，包括如下步骤：
S1、将曲面单层网壳进行三角化，即增加适当的虚梁，使得网格均为三角形网格，保证虚梁处在最接近曲面的面上；
S2、建立节点与坐标文件node.txt及element.txt；
S3、找出与某一节点相连的节点集合，集合中的每一个节点与该节点组成了梁单元，对任意梁单元，p1、p2为梁单元两节点坐标矢量，求定向第三点坐标矢量p3；
S4、假定全局坐标系中Z轴对应的矢量为Z＝[0,0,1]，求局部坐标系各轴矢量xr、xr0、yr、yr0、zr、zr0；
S5、空间梁的欧拉角BETA为局部坐标轴z与主平面为铅直面的空间梁之局部坐标轴z0之间的角度，也是将空间梁主平面绕梁自身轴转成铅直平面所需的转角，其数学表达如下：
若x-z平面法向(y轴)朝上，则BETA＝180/pi*acos(dot(zr0,zr))；
若x-z平面法向(y轴)朝下，则：
BETA＝360-180/pi*acos(dot(zr0,zr))；
式中dot()为矢量点积运算函数，acos()为反余弦函数。
2.如权利要求1所述的曲面单层网壳梁单元欧拉角的计算方法，其特征在于，S2具体包括以下步骤：
S21、NODE.txt为：
S22、ELEMENT.txt为：
3.如权利要求1所述的曲面单层网壳梁单元欧拉角的计算方法，其特征在于，S3中，其第一节点编号为NO1,其第二节点编号为NO2，并具体包括以下步骤：
S31、通过ELEMENT.txt，找出与节点NO1相连的节点集合S1、与节点NO2相连的节点集合S2，即S1中的每一个节点与节点NO1组成了梁单元，S2中的每一个节点与节点NO2组成了梁单元；
S32、对S1与S2做交集运算得其结果s12，即s12＝S1∩S2，s12为：
一个节点(NOa)，此梁必是边梁，此时，梁的法向n取三角形NO1,NO2,NOa的法向n1，即n＝n1；
或两个节点(Noa,Nob)，此梁必是中梁，其法向n与两个三角形相关,且梁的法向n取两个三角形法向的平均值，即n＝(n1+n2)/2；
则定向第三点坐标矢量p3＝(p1+p2)/2+1*n，(p1、p2为梁单元两节点坐标矢量)。
4.如权利要求1所述的曲面单层网壳梁单元欧拉角的计算方法，其特征在于，S4具体包括以下步骤：
(1)当梁处在铅直平面中，
梁的局部坐标轴对应的矢量xr0由下式求得：
xr＝p2-p1；xr0＝xt/norm(xr)；
梁的局部坐标轴对应的矢量yr0由下式求得：
yr0＝cross(Z,xr0)；
梁的局部坐标轴对应的矢量zr0由下式求得：
zr0＝cross(xr0,yr0)；
(2)当梁不处在铅直平面中，
梁的局部坐标轴对应的矢量xr0由下式求得：
xr＝p2-p1；xr0＝xt/norm(xr)；
zrt＝p3-p1；
梁的局部坐标轴对应的矢量yr由下式求得：
yr＝cross(zrt,xr)；yr＝yr/norm(yr)；
梁的局部坐标轴对应的矢量zr由下式求得：
zr＝cross(xr,yr)；
其中，上述公式中norm()为求矢量范数的函数；cross()为矢量叉积运算函数。
CN201711194465.5A 2017-11-24 2017-11-24 曲面单层网壳梁单元欧拉角的计算方法 Active CN107871049B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711194465.5A CN107871049B (zh) 2017-11-24 2017-11-24 曲面单层网壳梁单元欧拉角的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711194465.5A CN107871049B (zh) 2017-11-24 2017-11-24 曲面单层网壳梁单元欧拉角的计算方法

Publications (2)

Publication Number Publication Date
CN107871049A true CN107871049A (zh) 2018-04-03
CN107871049B CN107871049B (zh) 2021-02-19

Family

ID=61754589

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711194465.5A Active CN107871049B (zh) 2017-11-24 2017-11-24 曲面单层网壳梁单元欧拉角的计算方法

Country Status (1)

Country Link
CN (1) CN107871049B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109115695A (zh) * 2018-08-07 2019-01-01 华中科技大学 一种各向异性体材料光学常数和欧拉角的提取方法
CN111611742A (zh) * 2020-06-11 2020-09-01 上海卫星工程研究所 平面变形绝对-相对欧拉角计算方法及系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011047205A (ja) * 2009-08-27 2011-03-10 Nippon Steel Engineering Co Ltd ラチスシェル構造の節点構造及びその施工方法
CN105260341A (zh) * 2011-11-30 2016-01-20 西北工业大学 基于角速度的欧拉角勒让德近似输出方法
CN205296398U (zh) * 2016-01-05 2016-06-08 江苏沪宁钢机股份有限公司 一种大跨度单层网壳钢结构屋面
CN105806572A (zh) * 2016-03-14 2016-07-27 浙江大学 一种单层网壳结构节点半刚性测定方法
JP2016191562A (ja) * 2015-03-30 2016-11-10 住友金属鉱山株式会社 結晶方位の解析方法および解析装置
CN107247157A (zh) * 2017-05-10 2017-10-13 哈尔滨工程大学 一种面向大姿态机动的四元数全角域转换欧拉角的获取方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011047205A (ja) * 2009-08-27 2011-03-10 Nippon Steel Engineering Co Ltd ラチスシェル構造の節点構造及びその施工方法
CN105260341A (zh) * 2011-11-30 2016-01-20 西北工业大学 基于角速度的欧拉角勒让德近似输出方法
JP2016191562A (ja) * 2015-03-30 2016-11-10 住友金属鉱山株式会社 結晶方位の解析方法および解析装置
CN205296398U (zh) * 2016-01-05 2016-06-08 江苏沪宁钢机股份有限公司 一种大跨度单层网壳钢结构屋面
CN105806572A (zh) * 2016-03-14 2016-07-27 浙江大学 一种单层网壳结构节点半刚性测定方法
CN107247157A (zh) * 2017-05-10 2017-10-13 哈尔滨工程大学 一种面向大姿态机动的四元数全角域转换欧拉角的获取方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
E.L.DE ANGELIS: "Unified kinematic framework for a non-nominal Euler axis/angle rotation", 《ACTA ASTRONAUTICA》 *
陈东 等: "局部双层柱面网壳的几何非线性稳定性分析", 《空间结构》 *
顾磊 等: "考虑杆件失稳的单层球面网壳稳定性分析与试验研究", 《建筑结构学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109115695A (zh) * 2018-08-07 2019-01-01 华中科技大学 一种各向异性体材料光学常数和欧拉角的提取方法
CN111611742A (zh) * 2020-06-11 2020-09-01 上海卫星工程研究所 平面变形绝对-相对欧拉角计算方法及系统

Also Published As

Publication number Publication date
CN107871049B (zh) 2021-02-19

Similar Documents

Publication Publication Date Title
Richter-Gebert et al. First steps in tropical geometry
Zhu et al. Flocking of multi-agent non-holonomic systems with proximity graphs
Chen Stability of nonlinear systems
CN107871049A (zh) 曲面单层网壳梁单元欧拉角的计算方法
CN104361625A (zh) 一种基于射线原理的带边界保留的云数据精简算法
CN104933251B (zh) 一种翼面载荷处理方法
García-Naranjo Generalisation of Chaplygin’s reducing multiplier theorem with an application to multi-dimensional nonholonomic dynamics
CN106530228A (zh) 一种矢量多边形直角化纠正方法
CN106021781B (zh) 一种通用飞机机身轴线式载荷处理方法
Gaiko Global bifurcation analysis of the Lorenz system.
CN101739723A (zh) 自由曲面网络智能布局方法
CN115243270B (zh) 5g网络规划方法、装置、计算设备和存储介质
CN104794747A (zh) 一种基于射线原理的三维点云数据精简算法
CN105488088A (zh) 基于树形结构的二维网络角度分配布局方法
US9703897B2 (en) Modular design of complex tensegrity structures
Verhoeff et al. Regular 3d polygonal circuits of constant torsion
Gudi et al. Convergence analysis of the lowest order weakly penalized adaptive discontinuous Galerkin methods
Zhao et al. The Algorithm of Fast Collision Detection Based on Hybrid Bounding Box
Nikolajevic et al. A method to compute efficient 3d helicopters flight trajectories based on a motion polymorph-primitives algorithm
Ziv Napoleon-like configurations and sequences of triangles
Xu et al. Design and workspace optimization of a 6/6 cable-suspended parallel robot
Ouyang et al. Formatting of Julia sets of complex dynamic systems
CN207282874U (zh) 微机消谐装置低压线路接头
Eelbode Solutions for the Hyperbolic Dirac Equation on R1, m
Ishii et al. A finite type invariant of order at most 4 for genus 2 handlebody-knots is a constant map

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
GR01 Patent grant
GR01 Patent grant