CN106354929B - 基于刚度变化原理的承载结构载荷传递路径可视化方法 - Google Patents
基于刚度变化原理的承载结构载荷传递路径可视化方法 Download PDFInfo
- Publication number
- CN106354929B CN106354929B CN201610756951.0A CN201610756951A CN106354929B CN 106354929 B CN106354929 B CN 106354929B CN 201610756951 A CN201610756951 A CN 201610756951A CN 106354929 B CN106354929 B CN 106354929B
- Authority
- CN
- China
- Prior art keywords
- bearing structure
- point
- load
- load path
- strain energy
- 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
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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Architecture (AREA)
- Mathematical Analysis (AREA)
- Structural Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Civil Engineering (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Shafts, Cranks, Connecting Bars, And Related Bearings (AREA)
Abstract
基于刚度变化原理的承载结构载荷传递路径可视化方法,先计算出承载结构中各点的U*及其相应的梯度矢量,利用四阶龙格库塔数值积分方法对承载结构内各点的U*梯度矢量进行数值积分,得到的曲线即为承载结构中的载荷路径;然后从U*大小为零处开始,提取U*场等值线,沿各等值线对等值线上各点做断面分析,提取其上的脊点,顺次连接各个脊点,即得到基于刚度变化的主载荷路径;本发明可以明确载荷路径在结构中的传递方向和分布形态,为准确把握承载结构中载荷传递和分配规律提供指导。
Description
技术领域
本发明属于承载结构设计技术领域,具体涉及一种基于刚度变化原理的承载结构载荷传递路径可视化方法。
技术背景
承载结构是工业装备的重要基础构件,它的主要作用是支撑其它零部件、传递载荷等。在结构设计阶段,载荷路径的可视化对于结构承载效能的评估及材料的优化利用亦有很大的帮助。主应力迹线法是目前应用最为广泛的一种提取载荷路径的方法,它广泛应用于混凝土结构布筋及结构的外形设计等方面,主应力迹线上各点处应力的性质相同,但其大小一般不相等,并且在一般情况下主应力迹线不是从载荷的作用点连接到约束点,主应力迹线的这些性质使其无法准确描述结构内载荷的传递路径。
发明内容
为了克服上述现有技术的缺点,本发明的目的在于提供一种基于刚度变化原理的承载结构载荷传递路径可视化方法,可以明确显示载荷路径在结构中的传递方向和分布形态,为准确把握承载结构中载荷传递和分配规律提供指导。
为达到上述目标,本发明采取的技术方案为:
基于刚度变化原理的承载结构载荷传递路径可视化方法,包括以下步骤:
1)求解承载结构的U*场:通过公式(1)计算出承载结构中各点的U*值并通过插值方法得到U*场,承载结构中某点的U*值的求解方法如下:
式中,U为承载结构受原始约束和载荷作用下的应变能;U'为承载结构在原始约束和载荷的基础上,另行添加对目标点全约束后的结构应变能;
2)可视化载荷路径:由U*场确定承载结构的载荷路径和主载荷路径,具体为:
2.1)生成载荷路径:利用步骤1)中得到的U*场求出承载结构内各点处的U*梯度,利用四阶龙格库塔数值积分方法对承载结构内各点U*的梯度矢量进行数值积分,得到的曲线即为刚度线,刚度线为承载结构中载荷路径;
2.2)生成主载荷路径:从U*大小为零处开始,在U*场中提取等值线;沿各条等值线对线上各点做断面分析,以U*梯度的垂直面作为断面,断面上曲线的极大值点即为脊点,提取各条等值线上的脊点;顺次连接各脊点得到U*场的脊线,脊线即基于刚度变化的主载荷路径。
所述的步骤1)中采用有限元方法求解承载结构的U*场,具体步骤如下:
1.1)有限元建模:利用有限元方法,对承载结构建模并划分网格,得到承载结构的有限元模型,进而得到承载结构的总刚度矩阵K;
1.2)求解承载结构受载时的应变能:根据实际情况对步骤1.1)得到的有限元模型施加约束和载荷,利用静力学分析得到承载结构的节点位移向量X,并求得节点位移向量的转置XT,按照公式(2)求解结构受载时的应变能:
1.3)依次全约束各节点并求解承载结构的应变能:在步骤1.2)基础上,依次全约束目标结构有限元模型的各节点,求解承载结构受载下的应变能,定义全约束第i个节点时的结构应变能为U'i;
1.4)求解承载结构的U*场:承载结构中第i个节点的U*i为:
依据公式(3),分别求解各节点处的U*值,并通过插值方法得到承载结构的U*场。
本发明的有益效果为:
由于使用的参数U*是基于承载结构全局应变能的变化,所以最终结果具有全局性;由于依据U*来确定载荷路径,所以载荷路径必然与加载点和约束点相连,符合载荷的传递性;由于利用刚度变化原理并以脊线作为载荷路径,所以使得到的载荷路径具有清晰的构型。
附图说明
图1为实施例中悬臂梁结构的尺寸及边界条件示意图。
图2为实施例中包含各点梯度矢量的U*场示意图。
图3为实施例中悬臂梁结构的刚度线及主载荷路径示意图。
图4为U*场等值线分布图及U*场三维分布曲面图,其中图4(a)为U*场等值线分布示意图,图4(b)为U*场三维分布曲面图。
图5为实施例中利用断面分析法提取脊线的示意图。
具体实施方式
下面结合320mm*200mm*5mm的悬臂梁结构受载情况下载荷路径可视化的实施例对本发明作进一步说明。
基于刚度变化原理的承载结构载荷传递路径可视化方法,包括以下步骤:
1)求解承载结构的U*场:利用有限元方法对该悬臂梁结构进行建模,施加约束和载荷,并计算出承载结构中各节点的U*值,然后通过插值得到U*场,承载结构中某点的U*值的求解方法如下:
式中,U为承载结构受原始约束和载荷作用下的应变能;U'为承载结构在原始约束和载荷的基础上,另行添加对目标点全约束后的结构应变能;
具体步骤如下:
1.1)有限元建模:利用有限元方法,对该悬臂梁结构划分网格,这里使用的是四节点壳单元,并定义材料为结构钢,得到承载结构的总刚度矩阵K;
1.2)求解承载结构受载时的应变能:本实施例中悬臂梁结构的尺寸及边界条件如图1所示,依据图1施加约束和载荷,即全约束模型左侧边的一列节点,对右侧边上中间的节点施加向下的载荷P,此实施例中P为100牛顿;利用静力学分析得到承载结构的节点位移向量X,并求得节点位移向量的转置XT,按照下式求解得到承载结构受载时的应变能:
1.3)依次全约束各节点并求解承载结构的应变能:在步骤1.2)基础上,依次全约束目标结构有限元模型的各节点,求解承载结构受载下的应变能,定义全约束第i个节点时的结构应变能为U'i;
1.4)求解承载结构的U*场:承载结构中第i节点的U*值为:
依据上式,分别求解各节点处的U*值,并通过插值方法得到承载结构的U*场,如图2所示;
2)可视化载荷路径:由步骤1)中得到的U*场确定承载结构的载荷路径和主载荷路径,具体为:
2.1)生成载荷路径:U*的分布云图反映了结构内部载荷传递特性,其梯度向量代表了结构内部载荷传递的方向,本实施例利用步骤1)中得到的U*场求出承载结构内各点处的U*梯度矢量,如图2所示;利用四阶龙格库塔数值积分方法对承载结构内各点U*的梯度矢量进行数值积分,得到的曲线即为刚度线,刚度线为承载结构中载荷路径,如图3中的细实线所示;
2.2)生成主载荷路径:图4(a)为某二维结构的U*等值线示意图,其对应的U*等高图如图4(b)所示;无量纲参数U*即代表承载结构中各点的刚度大小,在U*值大小随位置变化的U*场中,如图4所示,脊线即为基于刚度变化的主载荷路径;
脊线的位置由承载结构中各点U*的分布所决定,使用断面分析法来提取脊线,如图5所示,脊线提取的具体步骤为:从U*大小为零处开始,在U*场中提取等值线;沿各条等值线对线上各点做断面分析,本实施例中以U*梯度的垂直面作为断面,断面上曲线的极大值点即为脊点,提取各条等值线上的脊点;需要指出的是,在断面与脊线的方向垂直时,断面分析法提取脊线的误差较小,因为脊线是沿梯度的方向,因此本方法将与梯度的垂直面作为断面来分析极大值;顺次连接各脊点得到U*场的脊线,即基于刚度变化的主载荷路径,如图3中的黑色粗实线所示;
如图3所示,利用本发明提出的方法,最终可以得到基于刚度变化原理的载荷路径,从图中可以看出,从约束端到加载点,载荷路径逐渐汇集于主载荷路径;这也充分说明主载荷路径即U*分布云图的脊线在载荷传递过程中起着更大的作用;从图3可以看到本方法得到了能够明确反映载荷在结构中传递方向和分布形态的载荷路径。
Claims (2)
1.基于刚度变化原理的承载结构载荷传递路径可视化方法,其特征在于,包括以下步骤:
1)求解承载结构的U*场:通过公式(1)计算出承载结构中各点的U*值并通过插值方法得到U*场,承载结构中某点的U*值的求解方法如下:
式中,U为承载结构受原始约束和载荷作用下的应变能;U'为承载结构在原始约束和载荷的基础上,另行添加对目标点全约束后的结构应变能;
2)可视化载荷路径:由U*场确定承载结构的载荷路径和主载荷路径,具体为:
2.1)生成载荷路径:利用步骤1)中得到的U*场求出承载结构内各点处的U*梯度,利用四阶龙格库塔数值积分方法对承载结构内各点U*的梯度矢量进行数值积分,得到的曲线即为刚度线,刚度线为承载结构中载荷路径;
2.2)生成主载荷路径:从U*大小为零处开始,在U*场中提取等值线;沿各条等值线对线上各点做断面分析,以U*梯度的垂直面作为断面,断面上曲线的极大值点即为脊点,提取各条等值线上的脊点;顺次连接各脊点得到U*场的脊线,脊线即基于刚度变化的主载荷路径。
2.根据权利要求1所述的基于刚度变化原理的承载结构载荷传递路径可视化方法,其特征在于:所述的步骤1)中采用有限元方法求解承载结构的U*场,具体步骤如下:
1.1)有限元建模:利用有限元方法,对承载结构建模并划分网格,得到承载结构的有限元模型,进而得到承载结构的总刚度矩阵K;
1.2)求解承载结构受载时的应变能:根据实际情况对步骤1.1)得到的有限元模型施加约束和载荷,利用静力学分析得到承载结构的节点位移向量X,并求得节点位移向量的转置XT,按照公式(2)求解结构受载时的应变能:
1.3)依次全约束各节点并求解承载结构的应变能:在步骤1.2)基础上,依次全约束目标结构有限元模型的各节点,求解承载结构受载下的应变能,定义全约束第i个节点时的结构应变能为U'i;
1.4)求解承载结构的U*场:承载结构中第i个节点的U*i为:
依据公式(3),分别求解各节点处的U*值,并通过插值方法得到承载结构的U*场。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610756951.0A CN106354929B (zh) | 2016-08-29 | 2016-08-29 | 基于刚度变化原理的承载结构载荷传递路径可视化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610756951.0A CN106354929B (zh) | 2016-08-29 | 2016-08-29 | 基于刚度变化原理的承载结构载荷传递路径可视化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106354929A CN106354929A (zh) | 2017-01-25 |
CN106354929B true CN106354929B (zh) | 2019-04-09 |
Family
ID=57857055
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610756951.0A Active CN106354929B (zh) | 2016-08-29 | 2016-08-29 | 基于刚度变化原理的承载结构载荷传递路径可视化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106354929B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107563102A (zh) * | 2017-10-11 | 2018-01-09 | 燕山大学 | 一种承载结构的传力骨架可视化方法 |
CN108536946A (zh) * | 2018-04-02 | 2018-09-14 | 西安交通大学 | 基于Landweber迭代法的运行工况传递路径分析方法 |
CN110096756B (zh) * | 2019-04-08 | 2020-01-17 | 河海大学 | 一种考虑荷载不确定性的自由曲面结构形态创建方法 |
CN114595612B (zh) * | 2022-03-16 | 2024-05-10 | 上海核工程研究设计院股份有限公司 | 基于实体单元积分路径的受力构件的配筋计算方法及*** |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103065015A (zh) * | 2013-01-04 | 2013-04-24 | 西安交通大学 | 一种基于内力路径几何形态的承载结构低碳节材设计方法 |
CN103778278A (zh) * | 2013-12-31 | 2014-05-07 | 西安交通大学 | 一种机床模块式刚度表征与求解方法 |
-
2016
- 2016-08-29 CN CN201610756951.0A patent/CN106354929B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103065015A (zh) * | 2013-01-04 | 2013-04-24 | 西安交通大学 | 一种基于内力路径几何形态的承载结构低碳节材设计方法 |
CN103778278A (zh) * | 2013-12-31 | 2014-05-07 | 西安交通大学 | 一种机床模块式刚度表征与求解方法 |
Non-Patent Citations (3)
Title |
---|
Design of conducting paths based on topology optimization;Yongcun Zhang 等;《Heat Mass Transfer》;20080831;第44卷;第1217-1227页 |
Stiffness design of machine tool structures by a biologically inspired topology optimization method;BaotongLi 等;《International Journal of Machine Tools & Manufacture》;20140401;第84卷;第33-44页 |
大型薄壁件的多点支承/定位方法研究;王少锋 等;《西安交通大学学报》;20160630;第50卷(第6期);第122-129页 |
Also Published As
Publication number | Publication date |
---|---|
CN106354929A (zh) | 2017-01-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106354929B (zh) | 基于刚度变化原理的承载结构载荷传递路径可视化方法 | |
CN106227977B (zh) | 一种基于bim的钢结构网架液压提升时变有限元分析的方法 | |
CN109145446B (zh) | 一种基于模态应变能和卷积神经网络的结构损伤识别方法 | |
CN109992900A (zh) | 一种大体积混凝土多场实时在线协同智能仿真方法和*** | |
CN102254061A (zh) | 输电塔线塔耦合体系的有限元建模及受力方法 | |
CN106202799A (zh) | 一种多点吊装计算方法 | |
CN102222140A (zh) | 复杂桥梁结构分析建模的设计图纸信息提取方法 | |
CN109241562A (zh) | 基于多尺度有限元方法的微结构材料弹性性能测定方法 | |
Keck | Validation of the standalone implementation of the dynamic wake meandering model for power production | |
Peng et al. | Numerical investigation of the effects of structural geometric and material nonlinearities on limit-cycle oscillation of a cropped delta wing | |
Fang et al. | Wind-induced fatigue analysis of high-rise steel structures using equivalent structural stress method | |
CN103065015B (zh) | 一种基于内力路径几何形态的承载结构低碳节材设计方法 | |
Wang et al. | Research on the flow-induced stress characteristics of head-cover bolts of a pump-turbine during turbine start-up | |
CN107085646A (zh) | 一种输电线路下地面场强测量方法 | |
Ye et al. | Quantification of numerical uncertainties of CFD-predicted wind turbine performance with overset grids | |
Cai et al. | A robust algorithm for the generation of integration cells in Numerical Manifold Method | |
Wang et al. | Structural dynamic analysis of a tidal current turbine using geometrically exact beam theory | |
CN109635452A (zh) | 一种高效的多峰随机不确定性分析方法 | |
CN102567606A (zh) | 给水管网等水压面实时绘制方法 | |
Tusch et al. | Modeling of turbulent atmospheric flow around tubular and lattice meteorological masts | |
Banerjee | A software independent tool for mapping thermal results to structural model | |
CN106934160A (zh) | 海底管道整体屈曲数值模拟中管道长度的确定方法 | |
Kashefi et al. | An advanced technique for the analysis of box girder bridge superstructures | |
Sreenivas et al. | Higher-Order Accurate Simulations of Wind Turbine Flow Fields: A Poor Man’s Approach | |
Li et al. | The Fractional Diagnosis of Multi-Span Continuous Bridge’s Structural Damage Based on Neural Network and Genetic Algorithm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |