CN114707389A - 一种基于插值小波的多分辨梁单元的构造方法 - Google Patents

一种基于插值小波的多分辨梁单元的构造方法 Download PDF

Info

Publication number
CN114707389A
CN114707389A CN202210485933.9A CN202210485933A CN114707389A CN 114707389 A CN114707389 A CN 114707389A CN 202210485933 A CN202210485933 A CN 202210485933A CN 114707389 A CN114707389 A CN 114707389A
Authority
CN
China
Prior art keywords
unit
matrix
nodes
displacement
interpolation
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
CN202210485933.9A
Other languages
English (en)
Other versions
CN114707389B (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.)
Lanzhou University
Original Assignee
Lanzhou 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 Lanzhou University filed Critical Lanzhou University
Priority to CN202210485933.9A priority Critical patent/CN114707389B/zh
Publication of CN114707389A publication Critical patent/CN114707389A/zh
Application granted granted Critical
Publication of CN114707389B publication Critical patent/CN114707389B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • 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/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了一种基于插值小波的多分辨梁单元的构造方法，从而显式建立了一种具有埃尔米特插值性质与精确重构低阶多项式能力的小波多分辨逼近格式。继而，以此逼近格式为位移插值函数，基于拉格朗日运动学定理推导得到了描述梁单元运动的控制方程，从而提出了一种基于插值小波的多分辨梁单元的构造方法，用于相关结构的力学分析。可通过调节单元的分辨率水平调控整体或局部的求解精度，无需复杂的网格划分；轴向自由度数量与横向自由度数量可独立配置，以实现不同工况下的最优求解；本发明提出的小波多分辨梁单元的求解效率明显优于目前常用的基于三次多项式插值的梁单元，尤其是在自由振动分析与屈曲分析中优势非常明显。

Description

一种基于插值小波的多分辨梁单元的构造方法
技术领域
本发明属于有限元分析技术领域，具体涉及一种基于插值小波的多分辨梁单元的构造方法。
背景技术
杆件是土木工程、机械与航天航空等工程领域中常用的一种结构类型，在仿真模拟分析时，主要使用基于低阶多项式插值的梁单元对其进行力学分析。由于此类梁单元自身的分辨率水平无法调节，因此在求解中只能通过调整网格划分来调控结构分析精度，导致在实际应用中需反复调整网格划分才能获得满足要求的结果。而网格划分需消耗大量的时间，由此造成整体分析效率不高。特别是对于涉及裂纹或复杂载荷等存在局部大梯度的问题，生成所需的局部细化网格颇为耗时。同时，由于低阶多项式无法高效地逼近三角函数，导致上述以低价多项式为插值函数的梁单元在自由振动与屈曲分析中的效率不高。此外，在许多结构中杆件的轴向变形远比弯曲变形简单，且转角可由横向位移的一阶导数有效表征(欧拉伯努利梁理论)，显然对于此类问题只需使用较少的轴向自由度(相较于横向自由度)，并在杆件连接处与结构支撑处设置转角自由度即可满足要求。但在上述目前最常用的梁单元中，其轴向、横向与转角自由度数量总是相等，由此造成在这些结构分析中计算资源的浪费。
小波分析作为一种新型的数学工具，其最大的优势是具有局部多分辨分析特性。以小波级数为位移近似格式的有限元单元可以通过调控小波基函数的分辨率水平来便利地调节整体和局部分析，无需反复进行网格划分。特别是对于存在局部大梯度的问题，可以通过在相应的局部区域添加高分辨率水平的小波基函数来有效求解，无需进行局部细化有限元网格。这是一种优于传统网格加密与阶次升高的多分辨分析技术，具有优良的分析效率、求解精度与数值稳定性。然而，构造要求位移插值函数具有连续导数的小波梁单元并非易事。目前主要有两种构造途径，一是通过转换矩阵建立小波展开系数与节点位移之间的联系，如区间B样条小波梁单元，但这一方式需计算转换矩阵的逆，导致计算量大且存在数值失稳的风险(转换矩阵的可逆性依赖于单元内部节点分布)。另一方式是采用具有埃尔米特插值性质的三角小波来构造梁单元，但由于基于三角小波的位移插值函数无法精确重构低阶多项式，违背了有限元理论中所要求的位移插值函数需可精确表征刚体位移与常应变状态这一基本要求，因此存在解不收敛的风险，且其分析精度将会依赖于具体问题，如对于结构中经常出现的无横向载荷作用的梁，其分析精度从理论上来说不如采用低阶多项式插值的传统梁单元，因为此时后者可以获得数值意义上的精确解。
发明内容
本发明基于插值小波分析理论，提出了一种新型边界延拓方式，从而建立了一种具有埃尔米特插值性质与精确重构低阶多项式能力的小波多分辨逼近格式。继而，以此逼近格式为位移插值函数，基于拉格朗日运动学定理推导得到了描述梁单元运动与变形的控制方程，从而建立了一种基于插值小波的多分辨梁单元的构造方法，用于相关结构的力学分析。本发明的详细技术方案如下：
一种基于插值小波的多分辨梁单元的构造方法，计算过程按以下步骤进行：
1）生成描述单元拉伸变形的内部节点：
a.选定分析拉伸变形的基础分辨率水平
Figure DEST_PATH_IMAGE001
，生成描述拉伸的基础节点，
Figure 690991DEST_PATH_IMAGE002
；
b.在需要提高精度的区域添加局部节点
Figure DEST_PATH_IMAGE003
；
c.识别所有节点的二分点式表征参数p t (n)与k t (n)，其中函数p t (n)为不小于
Figure 171913DEST_PATH_IMAGE004
并使得
Figure DEST_PATH_IMAGE005
为整数的最小整数；
d.合并所有位于同一位置的节点，并将所有节点按坐标从小到大依次编号为
Figure 13967DEST_PATH_IMAGE006
；
2）生成描述单元弯曲变形的内部节点：
a.选定描述弯曲变形的基础分辨率水平
Figure DEST_PATH_IMAGE007
，生成描述弯曲的基础节点
Figure 471493DEST_PATH_IMAGE008
；
b.在需要提高精度的区域添加局部节点
Figure DEST_PATH_IMAGE009
；
c.识别所有节点的二分点式表征参数p t (n)与k b (n)，其中函数p t (n)为不小于
Figure 817024DEST_PATH_IMAGE010
并使得
Figure DEST_PATH_IMAGE011
为整数的最小整数；
d.合并所有位于同一位置的节点，并删除所有
Figure 436224DEST_PATH_IMAGE012
以及
Figure DEST_PATH_IMAGE013
或
Figure 808650DEST_PATH_IMAGE014
的局部节点
Figure DEST_PATH_IMAGE015
，最后将所有节点按坐标从小到大依次编号为
Figure 69867DEST_PATH_IMAGE016
；
3）生成单元形函数：
基于单元内部节点，分别生成描述拉伸变形的单元形函数：
Figure DEST_PATH_IMAGE017
式中：
Figure 535483DEST_PATH_IMAGE018
为二阶插值小波函数；
同时生成描述弯曲变形的单元形函数：
Figure DEST_PATH_IMAGE019
其中
Figure 60005DEST_PATH_IMAGE020
式中
Figure DEST_PATH_IMAGE021
为四阶插值小波函数；
4）计算单元刚度矩阵：
Figure 640766DEST_PATH_IMAGE022
其中子矩阵
Figure DEST_PATH_IMAGE023
各元素具体为：
Figure 174515DEST_PATH_IMAGE024
其中
Figure DEST_PATH_IMAGE025
与
Figure 494638DEST_PATH_IMAGE026
分别为梁单元的抗拉刚度与抗弯刚度，j ku 和j kw 为控制单元刚度矩阵积分精度的独立参数，l e 为单元长度；其中，当抗拉刚度与抗弯刚度分别为常数时，有：
Figure DEST_PATH_IMAGE027
和
Figure 721220DEST_PATH_IMAGE028
5）计算单元质量矩阵：
Figure DEST_PATH_IMAGE029
其中子矩阵
Figure 370637DEST_PATH_IMAGE030
各元素具体为：
Figure DEST_PATH_IMAGE031
其中
Figure 442498DEST_PATH_IMAGE032
为梁单元的线密度，j u 和j w 为控制单元质量矩阵积分精度的独立参数；其中，当线密度为常数时，有：
Figure DEST_PATH_IMAGE033
6）计算单元几何刚度矩阵：
Figure 617128DEST_PATH_IMAGE034
其中子矩阵
Figure DEST_PATH_IMAGE035
各元素具体为：
Figure 483452DEST_PATH_IMAGE036
式中
Figure DEST_PATH_IMAGE037
为梁单元的轴力，j g 为控制单元几何刚度矩阵积分精度的独立参数；其中，当轴力为常数时，有：
Figure 743532DEST_PATH_IMAGE038
7）计算单元内部载荷的广义等效节点载荷：
Figure DEST_PATH_IMAGE039
其中：
Figure 852041DEST_PATH_IMAGE040
式中
Figure DEST_PATH_IMAGE041
其中
Figure 881176DEST_PATH_IMAGE042
、
Figure DEST_PATH_IMAGE043
与
Figure 183982DEST_PATH_IMAGE044
分别为梁单元上作用的轴向分布力，横向分布力与分布弯矩，j p 、j q 和j θ 为控制载荷积分精度的独立参数，
Figure DEST_PATH_IMAGE045
、
Figure 665779DEST_PATH_IMAGE046
与
Figure DEST_PATH_IMAGE047
为梁单元内部的轴向集中力、横向集中力与集中弯矩，
Figure 79443DEST_PATH_IMAGE048
分别为其作用点；其中，当轴向分布力、横向分布力或分布弯矩均匀分布时，分别有：
Figure DEST_PATH_IMAGE049
，
Figure 963085DEST_PATH_IMAGE050
或
Figure DEST_PATH_IMAGE051
8）生成全局坐标系下的单元运动控制方程：
Figure 436792DEST_PATH_IMAGE052
式中整体坐标系下的质量矩阵、刚度矩阵、几何刚度矩阵与广义等效节点载荷分别为：
Figure DEST_PATH_IMAGE053
广义节点位移与端部集中载荷的等效节点载荷分别为：
Figure 405885DEST_PATH_IMAGE054
Figure DEST_PATH_IMAGE055
其中坐标转换矩阵为
Figure 859125DEST_PATH_IMAGE056
，
Figure DEST_PATH_IMAGE057
式中
Figure 862853DEST_PATH_IMAGE058
，
Figure DEST_PATH_IMAGE059
，这里
Figure 773040DEST_PATH_IMAGE060
与
Figure DEST_PATH_IMAGE061
分别为单元起点与终点坐标，
Figure 229429DEST_PATH_IMAGE062
分别为单元起点处沿全局坐标系x方向的位移、y方向的位移与转角，
Figure DEST_PATH_IMAGE063
分别为单元终点处沿全局坐标系x方向的位移、y方向的位移与转角，
Figure 984896DEST_PATH_IMAGE064
分别为作用在单元起点处沿全局坐标系x方向的集中力、y方向的集中力与集中弯矩，
Figure DEST_PATH_IMAGE065
分别为作用在单元终点处沿全局坐标系x方向的集中力、y方向的集中力与集中弯矩；
9）基于步骤8）获取的单元运动方程，依据标准有限元的分析流程直接组装获得结构的总体运动方程，求解便可获得整体广义节点位移，继而获得相应的单元节点位移，由此基于形函数重构获得单元位移场。
一、[0, 1]区间上C 1型小波多分辨格式的构造
基于插值小波多分辨理论，在[0, 1]区间上逼近连续函数f(x)的一个多分辨格式可写为：
Figure 577551DEST_PATH_IMAGE066
，
Figure DEST_PATH_IMAGE067
其中θ(x) 为四阶插值小波，奇数集
Figure 393060DEST_PATH_IMAGE068
可为集合
Figure DEST_PATH_IMAGE069
的任意子集，d j,k 为展开系数。
上述展开格式中需用到实际中并不存在的域外节点值
Figure 841140DEST_PATH_IMAGE070
，
Figure DEST_PATH_IMAGE071
，需结合边界延拓技术用域内节点值插值给出此四个虚拟节点值。但目前所有边界延拓技术只能保证改进后的插值格式具有C 0连续性，即保证在x=0和1处，插值格式恒等于原函数，而无法保证其一阶导数值相等。
在区间
Figure 400297DEST_PATH_IMAGE072
与区间
Figure DEST_PATH_IMAGE073
上分别基于节点值f(0)，f'(0)，
Figure 113038DEST_PATH_IMAGE074
构造三次艾米特插值多项式可得：
Figure DEST_PATH_IMAGE075
将上述值代入原小波多分辨近似格式可得：
Figure 833870DEST_PATH_IMAGE076
其中基函数
Figure DEST_PATH_IMAGE077
从上述推导可以看出，为了保证可以构造三次埃尔米特插值多项式，基础分辨率水平j 0必须大于等于1。定义节点为
Figure 264851DEST_PATH_IMAGE078
，将上述近似格式中所有用到的节点按坐标从小至大排序后，可将该近似格式改写为常用格式：
Figure DEST_PATH_IMAGE079
其中展开系数：
Figure 362120DEST_PATH_IMAGE080
基函数为
Figure DEST_PATH_IMAGE081
。
上式即为本发明所构建的小波多分辨插值格式，可以证明其可精确重构所有不高于三阶的多项式(相较于三角小波梁单元的优点)，且具有埃尔米特插值性质(相较于其他小波插值格式的区别)，即满足关系：
Figure 194947DEST_PATH_IMAGE082
，
Figure DEST_PATH_IMAGE083
，
Figure 352259DEST_PATH_IMAGE084
，
Figure DEST_PATH_IMAGE085
Figure 506422DEST_PATH_IMAGE086
同时从上述过程可以看出，本发明所提出的C 1型插值基函数的构造无需计算任何逆矩阵，其可显式写为标准小波基函数的线性组合(相较于区间B样条小波梁单元的优点)。
二、小波多分辨梁单元的构造
采用上述本发明所构造的埃尔米特型小波多分辨插值格式表征梁单元的横向位移
Figure DEST_PATH_IMAGE087
可得(独有的插值格式)：
Figure 672961DEST_PATH_IMAGE088
式中形函数
Figure DEST_PATH_IMAGE089
，展开系数
Figure 94715DEST_PATH_IMAGE090
、
Figure DEST_PATH_IMAGE091
、
Figure 422928DEST_PATH_IMAGE092
以及
Figure DEST_PATH_IMAGE093
，其中
Figure 562923DEST_PATH_IMAGE094
为梁截面的转角，
Figure DEST_PATH_IMAGE095
为沿梁轴向的归一化坐标，l e 为梁单元的长度。
同时，用基于二阶插值小波的多分辨近似格式表征梁单元的轴向位移
Figure 267573DEST_PATH_IMAGE096
可得：
Figure DEST_PATH_IMAGE097
式中形函数
Figure 809413DEST_PATH_IMAGE098
，展开系数
Figure DEST_PATH_IMAGE099
以及
Figure 308528DEST_PATH_IMAGE100
，其中
Figure DEST_PATH_IMAGE101
为二阶插值小波基函数。
用上述两位移插值函数分别表征梁单元的动能、弹性势能以及外力虚功，继而基于拉格朗日运动学定理可得梁单元无阻尼振动方程为：
Figure 699932DEST_PATH_IMAGE102
式中单元广义节点位移向量
Figure DEST_PATH_IMAGE103
单元质量矩阵：
Figure 942695DEST_PATH_IMAGE104
其中子矩阵
Figure DEST_PATH_IMAGE105
上式中各元素具体为：
Figure 339041DEST_PATH_IMAGE106
其中
Figure DEST_PATH_IMAGE107
为梁单元的线密度，j u 和j w 为控制单元质量矩阵积分精度的独立参数。在上述计算格式中，将线密度用位移插值函数所采用的小波尺度级数进行了展开。其中，当线密度为一常数时，有：
Figure 9057DEST_PATH_IMAGE108
单元刚度矩阵为：
Figure DEST_PATH_IMAGE109
其中子矩阵
Figure 123643DEST_PATH_IMAGE110
上式中各元素具体为：
Figure DEST_PATH_IMAGE111
其中
Figure 904518DEST_PATH_IMAGE112
与
Figure DEST_PATH_IMAGE113
分别为梁单元的抗拉刚度与抗弯刚度，j ku 和j kw 为控制单元刚度矩阵积分精度的独立参数。在上述计算格式中，将抗拉刚度与抗弯刚度用位移插值函数所采用的小波尺度级数进行了展开。其中，当抗拉刚度与抗弯刚度分别为常数时，有：
Figure 686529DEST_PATH_IMAGE114
和
Figure DEST_PATH_IMAGE115
单元几何刚度矩阵为：
Figure 28911DEST_PATH_IMAGE116
其中子矩阵
Figure DEST_PATH_IMAGE117
其中各元素具体为：
Figure 630793DEST_PATH_IMAGE118
式中
Figure DEST_PATH_IMAGE119
为梁单元的轴力，j g 为控制单元几何刚度矩阵积分精度的独立参数。在上述计算格式中，将轴力用位移插值函数所采用的小波尺度级数进行了展开。其中，当轴力为常数时，有：
Figure 215358DEST_PATH_IMAGE120
单元内部载荷(除端点集中载荷)的广义等效节点载荷为：
Figure DEST_PATH_IMAGE121
其中：
Figure 586297DEST_PATH_IMAGE122
式中
Figure DEST_PATH_IMAGE123
其中
Figure 598115DEST_PATH_IMAGE124
分别为梁单元上作用的轴向分布力，横向分布力与分布弯矩，
Figure DEST_PATH_IMAGE125
为控制载荷积分精度的独立参数，
Figure 421715DEST_PATH_IMAGE126
为梁单元内部的轴向集中力、横向集中力与集中弯矩，
Figure DEST_PATH_IMAGE127
分别为其作用点。在上述计算格式中，将分布载荷用位移插值函数所采用的小波尺度级数进行了展开。其中，当轴向分布力、横向分布力或分布弯矩均匀分布时，分别有：
Figure 809971DEST_PATH_IMAGE128
，
Figure DEST_PATH_IMAGE129
或
Figure 35416DEST_PATH_IMAGE130
单元端点处集中载荷的广义等效节点载荷为：
Figure DEST_PATH_IMAGE131
式中
Figure 710811DEST_PATH_IMAGE132
为作用在梁单元起点处的轴向集中力、横向集中力与集中弯矩，
Figure DEST_PATH_IMAGE133
为作用在梁单元终点处的轴向集中力、横向集中力与集中弯矩。
上述计算过程中所有用到的函数值、导数值与积分值均可由基础数据库获得，无需进行数值积分，具体计算过程与基础数据库的生成方法已由发明人给出(Xiaojing Liu,Guirong Liu, Jizeng Wang, Youhe Zhou, A wavelet multiresolution interpolationGalerkin method for targeted local solution enrichment, Computational Mechanics, 2019, 64: 989–1016)。
设梁单元的起点与终点在全局坐标系下的坐标分别为
Figure 287286DEST_PATH_IMAGE134
，定义坐标转换矩阵
Figure DEST_PATH_IMAGE135
，
Figure 213654DEST_PATH_IMAGE136
其中
Figure DEST_PATH_IMAGE137
，
Figure 293605DEST_PATH_IMAGE138
。由单元坐标系下梁的运动控制方程可得全局坐标系下的运动控制方程为：
Figure DEST_PATH_IMAGE139
式中整体坐标系下的质量矩阵、刚度矩阵、几何刚度矩阵与广义等效节点载荷分别为：
Figure 912805DEST_PATH_IMAGE140
，
Figure DEST_PATH_IMAGE141
，
Figure 976576DEST_PATH_IMAGE142
，
Figure DEST_PATH_IMAGE143
广义节点位移与端部集中载荷的等效节点载荷分别为：
Figure 706635DEST_PATH_IMAGE144
Figure DEST_PATH_IMAGE145
其中
Figure 142558DEST_PATH_IMAGE146
分别为单元起点处沿全局坐标系x方向的位移、y方向的位移与转角，
Figure DEST_PATH_IMAGE147
分别为单元终点处沿全局坐标系x方向的位移、y方向的位移与转角，
Figure 932659DEST_PATH_IMAGE148
分别为作用在单元起点处沿全局坐标系x方向的集中力、y方向的集中力与集中弯矩，
Figure DEST_PATH_IMAGE149
分别为作用在单元终点处沿全局坐标系x方向的集中力、y方向的集中力与集中弯矩。
本发明的有益效果在于：
1.可通过调节单元的分辨率水平调控整体或局部的求解精度，无需复杂的网格划分；
2.轴向自由度数量与横向自由度数量可独立配置，以实现不同工况下的最优求解；
3.只在单元的端点处设置转角自由度以施加连续性条件或位移约束，相较于传统梁单元大幅减少了总自由度数量；
4.所使用的位移插值函数可精确重构所有不高于三阶的多项式，满足有限元理论的一致性要求，且具有优良的逼近三角函数的能力；
5.位移插值函数的构造无需计算任何逆矩阵，保证了方法的效率与稳定性；
6.本发明构建的小波多分辨梁单元全面支持变材料密度、变材料弹性模量与变杆件截面，以及常见的载荷与边界条件；
7.本发明提出的小波多分辨梁单元的求解效率明显优于目前常用的基于三次多项式插值的梁单元，尤其是在自由振动分析与屈曲分析中优势非常明显。
附图说明
图1实施例1等截面简支梁受横向均布载荷作用示意图；
图2实施例1本发明方法所使用的5组单元弯曲分辨率水平对应的单元内部弯曲节点分布图；
图3是实施例1本发明方法在不同单元分辨率水平下的绝对误差；
图4是实施例2变截面悬臂梁受横向集中载荷作用示意图；
图5是实施例3等截面简支梁示意图；
图6是实施例4阶梯梁受端部集中载荷作用示意图。
具体实施方式
下面结合具体实施例对本发明作进一步说明：
本部分提供4个实施例来直观展示本发明方法的有效性与性能，其中前三个算例的对照结果直接来自商业有限元软件ANSYS (2022 R1)的测试手册，最后一个算例的对照结果同样基于该测试手册上所提供的相关命令流修改而得(ANSYS, Inc., AnsysMechanical APDL Verification Manual, Release 2022 R1, Canonsburg: ANSYS,Inc., 2022)。
实施例1: 横向均布载荷作用下等截面简支梁的弯曲分析
如图1所示，梁长L=200 in，梁截面尺寸b=h=2.5 in，材料弹性模量E=30 Mpsi，横向分布载荷p=1.79253 lb/in。
表1. 实施例1本发明方法与商业有限元软件求解结果的比较
Figure 218147DEST_PATH_IMAGE150
在计算中将整根梁划分为1个单元，其中单元的拉伸分辨率始终设置为
Figure DEST_PATH_IMAGE151
(即采用两个拉伸自由度)，分别使用5组弯曲分辨率(如图2所示)进行求解，所得结果的绝对误差如图3所示。从图3中可以看出，当单元的整体或局部分辨率水平逐步升高时，解的误差在相应的区间内逐步减小，说明本发明方法可以通过调整单元的整体或局部分辨率水平有效调控整体或局部的分析精度，而无需重新划分网格。表1给出了本发明方法与目前常用商业有限元软件(使用基于三次多项式插值的梁单元)求解结果的比较，从中可以看出本发明方法具有优良的求解效率。如二者均使用均匀节点分布时，本发明方法使用21个自由度的相对误差仅为商业有限元软件使用27个自由度(划分为8个单元)的43%。
实施例2：横向集中载荷作用下连续变截面悬臂梁的弯曲分析
如图4所示，梁长L=20 in，高h=0.5 in，宽为一底边长b=3 in的等腰三角形，材料弹性模量E=30 Mpsi，横向集中载荷P=10 lbs。
表2. 本发明方法与商业有限元软件求解结果的比较
Figure 751897DEST_PATH_IMAGE152
在计算中将整根梁划分为1个单元，其中单元的拉伸分辨率设置为0，弯曲分辨率设置为1，无局部节点，此时单元共7个自由度(所允许的最少自由度数量)。表2分别给出了本发明方法与商业有限元软件的求解结果，其中后者将梁划分为10个等长单元。从中可以看出，本发明方法获得了数值意义上的精确解(即只存在数值计算的舍入误差)，证明了本发明梁单元所使用的位移插值函数可精确重构所有不高于3阶的多项式(本问题的精确解为二次多项式)。此外，表2的结果也说明本发明方法在变截面梁分析中相较于目前常用的商业有限元更为高效。
实施例3：等截面简支梁的自由振动分析
如图5所示，梁长L=80 in，梁截面积A=4 in2，惯性矩I=1.3333 in4，材料密度ρ=0.000728 lb-sec2/in4，弹性模量E=30 Mpsi。
表3. 本发明方法与商业有限元软件求解结果的比较
Figure DEST_PATH_IMAGE153
在计算中将整根梁划分为1个单元，其中单元的拉伸分辨率设置为0，弯曲分辨率设置为4，无局部节点，此时单元共21个自由度。表3给出了基于本发明方法与商业有限元软件所得的前5阶固有频率，其中后者使用8个等长单元。从表3可以看出，本发明方法在自由振动分析中相较于目前常用的商业有限元具有更高的分析效率，尤其是对于高阶自由振动模态。如对于第5阶模态，前者使用21个自由度时的相对误差仅为后者使用27个自由度时的7.76%。
实施例4：阶梯梁屈曲分析
如图6所示，梁长L=6 m，左半部分与右半部分的横截面尺寸分别为b 0=2h 0=0.2 m与b 1=h 1=0.1 m，材料弹性模量E=200 GPa。
表4. 本发明方法与商业有限元软件求解结果的比较
Figure DEST_PATH_IMAGE155
在计算中将整根梁划分为2个等长的单元，其中两个单元的拉伸分辨率均设置为0，弯曲分辨率均设置为4，无局部节点，此时每个单元各有21个自由度，整体共39个自由度。表4给出了基于本发明方法与商业有限元软件所得的临界屈曲载荷，其中后者使用12个等长单元，整体也共39个自由度。从表4可以看出，本发明方法在屈曲分析中相较于目前常用的商业有限元具有更高的分析效率。如对于两端简支阶梯梁，二者使用相同数量的自由度时，本发明方法的相对误差仅为商业有限元的3.36%。

Claims (1)

1.一种基于插值小波的多分辨梁单元的构造方法，其特征在于，计算过程按以下步骤进行：
1）生成描述单元拉伸变形的内部节点：
a.选定分析拉伸变形的基础分辨率水平
Figure DEST_PATH_IMAGE002
，生成描述拉伸的基础节点，
Figure DEST_PATH_IMAGE004
；
b.在需要提高精度的区域添加局部节点
Figure DEST_PATH_IMAGE006
；
c.识别所有节点的二分点式表征参数p t (n)与k t (n)，其中函数p t (n)为不小于
Figure DEST_PATH_IMAGE008
并使得
Figure DEST_PATH_IMAGE010
为整数的最小整数；
d.合并所有位于同一位置的节点，并将所有节点按坐标从小到大依次编号为
Figure DEST_PATH_IMAGE012
；
2）生成描述单元弯曲变形的内部节点：
a.选定描述弯曲变形的基础分辨率水平
Figure DEST_PATH_IMAGE014
，生成描述弯曲的基础节点
Figure DEST_PATH_IMAGE016
；
b.在需要提高精度的区域添加局部节点
Figure DEST_PATH_IMAGE018
；
c.识别所有节点的二分点式表征参数p t (n)与k b (n)，其中函数p t (n)为不小于
Figure DEST_PATH_IMAGE020
并使得
Figure DEST_PATH_IMAGE022
为整数的最小整数；
d.合并所有位于同一位置的节点，并删除所有
Figure DEST_PATH_IMAGE024
以及
Figure DEST_PATH_IMAGE026
或
Figure DEST_PATH_IMAGE028
的局部节点
Figure DEST_PATH_IMAGE030
，最后将所有节点按坐标从小到大依次编号为
Figure DEST_PATH_IMAGE032
；
3）生成单元形函数：
基于单元内部节点，分别生成描述拉伸变形的单元形函数：
Figure DEST_PATH_IMAGE034
式中：
Figure DEST_PATH_IMAGE036
为二阶插值小波函数；
同时生成描述弯曲变形的单元形函数：
Figure DEST_PATH_IMAGE038
其中
Figure DEST_PATH_IMAGE040
式中
Figure DEST_PATH_IMAGE042
为四阶插值小波函数；
4）计算单元刚度矩阵：
Figure DEST_PATH_IMAGE044
其中子矩阵
Figure DEST_PATH_IMAGE046
各元素具体为：
Figure DEST_PATH_IMAGE048
其中
Figure DEST_PATH_IMAGE050
与
Figure DEST_PATH_IMAGE052
分别为梁单元的抗拉刚度与抗弯刚度，j ku 和j kw 为控制单元刚度矩阵积分精度的独立参数，l e 为单元长度；其中，当抗拉刚度与抗弯刚度分别为常数时，有：
Figure DEST_PATH_IMAGE054
和
Figure DEST_PATH_IMAGE056
5）计算单元质量矩阵：
Figure DEST_PATH_IMAGE058
其中子矩阵
Figure DEST_PATH_IMAGE060
各元素具体为：
Figure DEST_PATH_IMAGE062
其中
Figure DEST_PATH_IMAGE064
为梁单元的线密度，j u 和j w 为控制单元质量矩阵积分精度的独立参数；其中，当线密度为常数时，有：
Figure DEST_PATH_IMAGE066
6）计算单元几何刚度矩阵：
Figure DEST_PATH_IMAGE068
其中子矩阵
Figure DEST_PATH_IMAGE070
各元素具体为：
Figure DEST_PATH_IMAGE072
式中
Figure DEST_PATH_IMAGE074
为梁单元的轴力，j g 为控制单元几何刚度矩阵积分精度的独立参数；其中，当轴力为常数时，有：
Figure DEST_PATH_IMAGE076
7）计算单元内部载荷的广义等效节点载荷：
Figure DEST_PATH_IMAGE078
其中：
Figure DEST_PATH_IMAGE080
式中
Figure DEST_PATH_IMAGE082
其中
Figure DEST_PATH_IMAGE084
、
Figure DEST_PATH_IMAGE086
与
Figure DEST_PATH_IMAGE088
分别为梁单元上作用的轴向分布力，横向分布力与分布弯矩，j p 、j q 和j θ 为控制载荷积分精度的独立参数，
Figure DEST_PATH_IMAGE090
、
Figure DEST_PATH_IMAGE092
与
Figure DEST_PATH_IMAGE094
为梁单元内部的轴向集中力、横向集中力与集中弯矩，
Figure DEST_PATH_IMAGE096
分别为其作用点；其中，当轴向分布力、横向分布力或分布弯矩均匀分布时，分别有：
Figure DEST_PATH_IMAGE098
，
Figure DEST_PATH_IMAGE100
或
Figure DEST_PATH_IMAGE102
8）生成全局坐标系下的单元运动控制方程：
Figure DEST_PATH_IMAGE104
式中整体坐标系下的质量矩阵、刚度矩阵、几何刚度矩阵与广义等效节点载荷分别为：
Figure DEST_PATH_IMAGE106
广义节点位移与端部集中载荷的等效节点载荷分别为：
Figure DEST_PATH_IMAGE108
Figure DEST_PATH_IMAGE110
其中坐标转换矩阵为
Figure DEST_PATH_IMAGE112
，
Figure DEST_PATH_IMAGE114
式中
Figure DEST_PATH_IMAGE116
，
Figure DEST_PATH_IMAGE118
，这里
Figure DEST_PATH_IMAGE120
与
Figure DEST_PATH_IMAGE122
分别为单元起点与终点坐标，
Figure DEST_PATH_IMAGE124
分别为单元起点处沿全局坐标系x方向的位移、y方向的位移与转角，
Figure DEST_PATH_IMAGE126
分别为单元终点处沿全局坐标系x方向的位移、y方向的位移与转角，
Figure DEST_PATH_IMAGE128
分别为作用在单元起点处沿全局坐标系x方向的集中力、y方向的集中力与集中弯矩，
Figure DEST_PATH_IMAGE130
分别为作用在单元终点处沿全局坐标系x方向的集中力、y方向的集中力与集中弯矩；
9）基于步骤8）获取的单元运动方程，依据标准有限元的分析流程直接组装获得结构的总体运动方程，求解便可获得整体广义节点位移，继而获得相应的单元节点位移，由此基于形函数重构获得单元位移场。
CN202210485933.9A 2022-05-06 2022-05-06 一种基于插值小波的多分辨梁单元的构造方法 Active CN114707389B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210485933.9A CN114707389B (zh) 2022-05-06 2022-05-06 一种基于插值小波的多分辨梁单元的构造方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210485933.9A CN114707389B (zh) 2022-05-06 2022-05-06 一种基于插值小波的多分辨梁单元的构造方法

Publications (2)

Publication Number Publication Date
CN114707389A true CN114707389A (zh) 2022-07-05
CN114707389B CN114707389B (zh) 2024-04-12

Family

ID=82176215

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210485933.9A Active CN114707389B (zh) 2022-05-06 2022-05-06 一种基于插值小波的多分辨梁单元的构造方法

Country Status (1)

Country Link
CN (1) CN114707389B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101114308A (zh) * 2007-08-24 2008-01-30 西安交通大学 一种三类变量区间b样条小波梁单元构造方法
CN110924904A (zh) * 2018-09-20 2020-03-27 中国石油化工股份有限公司 一种由井底泵功图调整抽油机电机转速的方法
CN112180592A (zh) * 2020-10-14 2021-01-05 兰州大学 一种用于光学逆向设计的快速优化方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101114308A (zh) * 2007-08-24 2008-01-30 西安交通大学 一种三类变量区间b样条小波梁单元构造方法
CN110924904A (zh) * 2018-09-20 2020-03-27 中国石油化工股份有限公司 一种由井底泵功图调整抽油机电机转速的方法
CN112180592A (zh) * 2020-10-14 2021-01-05 兰州大学 一种用于光学逆向设计的快速优化方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张兴武;陈雪峰;杨志勃;何正嘉;: "基于多变量小波有限元的一维结构分析", 工程力学, no. 08, 25 August 2012 (2012-08-25) *

Also Published As

Publication number Publication date
CN114707389B (zh) 2024-04-12

Similar Documents

Publication Publication Date Title
CN108345741B (zh) 基于无网格rkpm各向异性材料二维热变形和热应力分析方法
Xing et al. A differential quadrature finite element method
Qian Full analytical sensitivities in NURBS based isogeometric shape optimization
Gofuku et al. Point-tangent/point-normal B-spline curve interpolation by geometric algorithms
Liu et al. Isogeometric method based in-plane and out-of-plane free vibration analysis for Timoshenko curved beams
Li et al. Beam structure optimization for additive manufacturing based on principal stress lines
Gang et al. Mesh deformation on 3D complex configurations using multistep radial basis functions interpolation
Yan et al. Structure/material concurrent optimization of lattice materials based on extended multiscale finite element method
Wang et al. Voronoi polygonal hybrid finite elements with boundary integrals for plane isotropic elastic problems
CN111125963A (zh) 基于拉格朗日积分点有限元的数值仿真系统及方法
Petrolito Hybrid-Trefftz quadrilateral elements for thick plate analysis
Ding et al. Exact and efficient isogeometric reanalysis of accurate shape and boundary modifications
Roh et al. Integration of geometric design and mechanical analysis using B‐spline functions on surface
Kirsch et al. Accurate displacement derivatives for structural optimization using approximate reanalysis
Trebotich et al. A Projection Method for Incompressible Viscous Flow on Moving QuadrilateralGrids
CN114707389A (zh) 一种基于插值小波的多分辨梁单元的构造方法
Devarajan Vibration Analysis of Timoshenko Beams using Isogeometric Analysis
Zhang et al. A T-splines-oriented isogeometric topology optimization for plate and shell structures with arbitrary geometries using Bézier extraction
Chen et al. An edge center based strain-smoothing element with discrete shear gap for the analysis of Reissner–Mindlin shell
Wang et al. Kirchhoff–Love shell representation and analysis using triangle configuration b-splines
He et al. Shape optimization of a flexible beam with a local shape feature based on ANCF
CN117390778A (zh) 一种考虑薄板应变强化效应的塑性安定上限载荷计算方法
Xia et al. Analysis-aware modelling of spacial curve for isogeometric analysis of Timoshenko beam
CN115453873A (zh) 非线性结构动力学系统瞬态响应无条件稳定时间积分方法
Ding et al. An improved explicit-implicit precise integration method for nonlinear dynamic analysis of structures

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