CN107871049B - 曲面单层网壳梁单元欧拉角的计算方法 - Google Patents
曲面单层网壳梁单元欧拉角的计算方法 Download PDFInfo
- Publication number
- CN107871049B CN107871049B CN201711194465.5A CN201711194465A CN107871049B CN 107871049 B CN107871049 B CN 107871049B CN 201711194465 A CN201711194465 A CN 201711194465A CN 107871049 B CN107871049 B CN 107871049B
- Authority
- CN
- China
- Prior art keywords
- node
- vector
- axis
- local coordinate
- curved surface
- 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
Images
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 (3)
1.曲面单层网壳梁单元欧拉角的计算方法,其特征在于,包括如下步骤:
S1、将曲面单层网壳进行三角化,即增加适当的虚梁,使得网格均为三角形网格,保证虚梁处在最接近曲面的面上;
S2、建立节点与坐标文件node.txt及element.txt;
S3、找出与某一节点相连的节点集合,集合中的每一个节点与该节点组成了梁单元,对任意梁单元,p1、p2为梁单元两节点坐标矢量,求定向第三点坐标矢量p3,具体包括:
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;
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()为反余弦函数。
3.如权利要求1所述的曲面单层网壳梁单元欧拉角的计算方法,其特征在于,S4具体包括以下步骤:
(1)当梁处在铅直平面中,
梁的局部坐标轴对应的矢量xr0由下式求得:
xr0=xr/norm(xr),其中xr=p2-p1;
梁的局部坐标轴对应的矢量yr0由下式求得:
yr0=cross(Z,xr0);
梁的局部坐标轴对应的矢量zr0由下式求得:
zr0=cross(xr0,yr0);
(2)当梁不处在铅直平面中,
梁的局部坐标轴对应的矢量xr0由下式求得:
xr0=xr/norm(xr),其中xr=p2-p1;
梁的局部坐标轴对应的矢量yr由下式求得:
yrt=cross(zrt,xr);yr=yrt/norm(yrt);
其中zrt=p3-p1;
梁的局部坐标轴对应的矢量zr由下式求得:
zr=cross(xr,yr);
其中,上述公式中,norm()为求矢量范数的函数;cross()为矢量叉积运算函数。
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 CN107871049A (zh) | 2018-04-03 |
CN107871049B true 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) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109115695B (zh) * | 2018-08-07 | 2020-05-19 | 华中科技大学 | 一种各向异性体材料光学常数和欧拉角的提取方法 |
CN111611742B (zh) * | 2020-06-11 | 2022-03-18 | 上海卫星工程研究所 | 平面变形绝对-相对欧拉角计算方法及系统 |
Citations (6)
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 | 哈尔滨工程大学 | 一种面向大姿态机动的四元数全角域转换欧拉角的获取方法 |
-
2017
- 2017-11-24 CN CN201711194465.5A patent/CN107871049B/zh active Active
Patent Citations (6)
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)
Title |
---|
Unified kinematic framework for a non-nominal Euler axis/angle rotation;E.L.de Angelis;《Acta Astronautica》;20151231;第333-338页 * |
局部双层柱面网壳的几何非线性稳定性分析;陈东 等;《空间结构》;19990930;第5卷(第3期);第33-39页 * |
考虑杆件失稳的单层球面网壳稳定性分析与试验研究;顾磊 等;《建筑结构学报》;20170731;第38卷(第7期);第25-33页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107871049A (zh) | 2018-04-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Deng et al. | Synthesis of deployable/foldable single loop mechanisms with revolute joints | |
Chen et al. | Intelligent computational design of scalene-faceted flat-foldable tessellations | |
Morterolle et al. | Numerical form-finding of geotensoid tension truss for mesh reflector | |
Wei et al. | General approach for inverse kinematics of nR robots | |
Gan et al. | Mobility change in two types of metamorphic parallel mechanisms | |
Zhang et al. | A kirigami-inspired 8R linkage and its evolved overconstrained 6R linkages with the rotational symmetry of order two | |
Qin et al. | Multi-furcation in a derivative queer-square mechanism | |
CN107871049B (zh) | 曲面单层网壳梁单元欧拉角的计算方法 | |
Ghorbel et al. | Dynamic stiffness formulation for free orthotropic plates | |
Parigi et al. | Three-dimensionality in reciprocal structures: concepts and generative rules | |
Stephan et al. | Reticulated structures on free-form surfaces | |
Du et al. | Stiffness and dexterous performances optimization of large workspace cable-driven parallel manipulators | |
CN108875132B (zh) | 一种多工况下应力约束类桁架连续体结构的拓扑优化方法 | |
Xu et al. | Discrete schemes for Gaussian curvature and their convergence | |
Zhu et al. | Three-dimensional geometry of magnetic reconnection induced by ballooning instability in a generalized Harris sheet | |
CN211468825U (zh) | 姿态控制装置和平板式卫星 | |
CN106021930B (zh) | 扭转力矩作用下的最小质量张拉整体结构设计方法 | |
Abramczyk | Building Structures Roofed with Multi-Segment Corrugated Hyperbolic Paraboloid Steel Shells | |
Chen et al. | Bearing rigidity and formation stabilization for multiple rigid bodies in SE (3) | |
US10628533B2 (en) | Global optimization of networks of locally fitted objects | |
Yang et al. | A class of quasi-quartic trigonometric Bézier curves and surfaces | |
Senba et al. | Vibration reduction by natural frequency optimization for manipulation of a variable geometry truss | |
Liu et al. | An effective self-collision detection algorithm for multi-degree-of-freedom manipulator | |
Jeong et al. | Virtual joint method for kinematic modeling of wheeled mobile manipulators | |
Lu et al. | Dynamics analysis of a novel 5-DoF 3SPU+ 2SPRR type parallel manipulator |
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 |