基于三维6节点矩形单元的结构变形分解方法技术

本发明专利技术属于力学分析领域，公开了一种基于三维6节点矩形单元的结构变形分解方法，包括以下步骤：构造出6节点矩形单元的基本位移和基本变形基向量，得到三维6节点矩形单元的正交力学基矩阵；建立空间结构模型，得到6节点矩形单元各个节点的位移向量；将6节点矩形单元的位移向量投影到三维6节点矩形单元的正交力学基矩阵上，得到每个矩形单元的基本变形和基本位移的投影系数向量；根据投影系数向量中投影系数的大小，得出矩形单元发生的主要基本变形和次要基本变形，即可实现对结构模型的变形分解和变形识别。本发明专利技术可以精确识别厚度较宽度和长度来说很小的单相受力薄板类构件的面内变形情况，还能精确识别出面外弯曲变形等面外变形。外变形。外变形。

【技术实现步骤摘要】
基于三维6节点矩形单元的结构变形分解方法


[0001]本专利技术属于力学分析
，具体涉及一种基于三维6节点矩形单元的结构变形分解方法。
技术背景
[0002]三维4节点矩形单元正交分解出的变形基中主要是面内的基本变形，面外的变形只有翘曲变形，这对于面外基本变形的分解是不准确的。考虑到4节点矩形单元的缺点，以及异形柱、跨高比较大的连梁等构件单相受力的特点，构造出三维6节点矩形单元。
[0003]目前已经公布的变形分解方法均为对某一构件的分解，在对体型庞大的结构模型进行变形分解时，可算性尤为重要，三维6节点矩形单元作为一种宏观的矩形单元，可准确识别构件面内及面外变形情况，并且6节点矩形单元自由度较少、计算速度较快，可快速完成复杂结构的计算，并识别出各个部位的构件的面内变形与面外变形情况。

技术实现思路

[0004]本专利技术的目的在于提出一种基于三维6节点矩形单元的结构变形分解方法，不仅可以精确识别厚度较宽度和长度来说很小的单相受力薄板类构件的面内变形情况，还能精确识别出面外弯曲变形等面外变形。r/>[0005]为本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于三维6节点矩形单元的结构变形分解方法，其特征在于，包括以下步骤：步骤1：在三维直角坐标系中，依据受力平衡条件和正交分解理论，构造出6节点矩形单元的基本位移和基本变形基向量，从而得到三维6节点矩形单元的正交力学基矩阵P；步骤2：在空间直角坐标系下，建立空间结构模型，并采用四边形壳单元对空间结构模型进行划分，选择四边形壳单元的4个角点和两条对边上的中点组成三维6节点矩形单元，得到6节点矩形单元各个节点的位移向量s；步骤3：将6节点矩形单元各个节点的位移向量s投影到三维6节点矩形单元的正交力学基矩阵P上，得到每个矩形单元的基本变形和基本位移的投影系数向量r；步骤4：根据投影系数向量r中投影系数的大小，得出每个矩形单元发生的主要基本变形和次要基本变形，即可实现对空间结构模型的变形分解和变形识别。2.根据权利要求1所述的基于三维6节点矩形单元的结构变形分解方法，其特征在于，所述的6节点矩形单元的基本位移和基本变形包括：X轴向的刚体平动位移、Y轴向的刚体平动位移、Z轴向的刚体平动位移、绕X轴向刚体转动位移、绕Y轴向刚体转动位移、绕Z轴向刚体转动位移、XOY面X轴向拉压变形、XOY面Y轴向拉压变形、XOY面X轴向面内弯曲变形、XOY面Y轴向面内弯曲变形、XOY面剪切变形、XOY面翘曲变形、XOY面面外弯曲变形、XOY面胀缩变形、XOY面冲切变形、XOY面反向翘曲变形、XOY面反对称弯曲变形以及XOY面反对称拉伸
‑
压缩变形。3.根据权利要求2所述的基于三维6节点矩形单元的结构变形分解方法，其特征在于，所述步骤1具体包括以下步骤：步骤1.1：在空间直角坐标系下，依据受力平衡条件、力矩平衡条件以及正交理论，针对6节点矩形单元构造出单元的18种基本位移和基本变形基向量p1～p
18
，具体如下：p1表示X轴向刚体平动位移基向量：p2表示Y轴向刚体平动位移基向量：p3表示Z轴向刚体平动位移基向量：p4表示绕X轴向刚体转动位移基向量：p5表示绕Y轴向刚体转动位移基向量：p5＝(0,0,
‑
0.5,0,0,0.5,0,0,0.5,0,0,
‑
0.5,0,0,0,0,0,0)
T
，p6表示绕Z轴向刚体转动位移基向量：
p7表示XOY面X轴向拉压变形基向量：p7＝(
‑
0.5,0,0,0.5,0,0,0.5,0,0,
‑
0.5,0,0,0,0,0,0,0,0)
T
，p8表示XOY面Y轴向拉压变形基向量：p9表示XOY面X轴向面内弯曲变形基向量：p9＝(
‑
0.5,0,0,0.5,0,0,
‑
0.5,0,0,0.5,0,0,0,0,0,0,0,0)
T
，p
10
表示XOY面Y轴向面内弯曲变形基向量：p
10
＝(0,
‑
0.5,0,0,0.5,0,0,
‑
0.5,0,0,0.5,0,0,0,0,0,0,0)
T
，p
11
表示XOY面剪切变形基向量：p
12
表示XOY面翘曲变形基向量：p
12
＝(,0,0,
‑
0.5 0,0,0.5,0,0,
‑
0.5,0,0,0.5,0,0,0,0,,00)
T
，p
13
表示XOY面面外弯曲变形基向量：p
14
表示XOY面胀缩变形基向量：p
15
表示XOY面冲切变形基向量：p
16
表示XOY面反向翘曲变形基向量：p
17
表示XOY面反对称弯曲变形基向量：p
18
表示XOY面反对称拉伸
‑
压缩变形基向量：其中，m为矩形单元在Y轴向的长度，n为矩形单元在X轴向的长度；步骤1.2：将步骤1.1所示的6个基本位移基向量与12个基本变形基向量组合构造出6节点矩形单元的完备正交力学基矩阵P：P＝[p
1 p
2 p
3 p
4 p
5 p
6 p
7 p
8 p
9 p
10 p
11 p
12 p
13 p
14 p
15 p
16 p
17 p
18
]。4.根据权利要求1所述的基于三维6节点矩形单元的结构变形分解方法，其特征在于，
步骤2具体包括以下步骤：步骤2.1：利用有限元软件，在三维直角坐标系下建立空间结构模型，并采用四边形壳单元对空间结构模型进行划分，然后选择任一四边形壳单元的4个角点和X轴方向上的两个中点组成6节点矩形单元，则第i个6节点矩形单元的节点坐标d
i
：d
i
＝(x
i1
,y
i1
,z
i1
,x
i2
,y
i2
,z
i2
,x
i3
,y
i3
,z
i3
,x
i4
,y
i4
,z
i4
,x
i5
,y
i5
,z
i5
,x
i6
,y
i6
,z
i6
)，其中：d
i
表示空间模型中第i个矩形单元的节点坐标，i＝1,2,3...,k，k为矩形单元的数量；x
ij
表示第i个矩形单元第j个节点的X轴上的坐标值；y
ij
表示第i个矩形单元第j个节点的Y轴上的坐标值；z
ij
表示第i个矩形单元第j个节点的Z轴上的坐标值；j＝1,2,3,4,5,6；步骤2.2：空间模型在任意荷载作用下发生位移和变形，产生变形后的第i个6节点矩形单元节点坐标d
′
i
：d
′
i
＝(x
′
i1
,y
′
i1
,z
′
i1
,x
′
i2
,y
′
i2
,z
′
i2
,x
′
i3
,y
′
i3
,z
...

【专利技术属性】
技术研发人员：孙攀旭，赵瑞青，杨强，严亚丹，王东炜，
申请(专利权)人：郑州大学，
类型：发明
国别省市：