CN105677983A - 基于软硬件实时交互优化的计算方法 - Google Patents
基于软硬件实时交互优化的计算方法 Download PDFInfo
- Publication number
- CN105677983A CN105677983A CN201610011844.5A CN201610011844A CN105677983A CN 105677983 A CN105677983 A CN 105677983A CN 201610011844 A CN201610011844 A CN 201610011844A CN 105677983 A CN105677983 A CN 105677983A
- Authority
- CN
- China
- Prior art keywords
- granule
- time
- discrete element
- intelligence
- real
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于软硬件实时交互优化的计算方法,包括数值模拟颗粒的形成、智能颗粒传感设备制造、“离散介质实时交互优化法则”及其应用等步骤。本发明简便、可靠、可高效地提高离散元等数值模拟计算方法精度;可以被广泛地应用于土木、机械、化工、冶金、农业、制药、环境等各项领域,对复杂、昂贵、耗时、危险的很多实验、现象进行模拟,解决工程实际问题,预测结构、材料、系统性能。
Description
技术领域
本发明涉及一种基于软硬件实时交互优化的计算方法。
背景技术
有限元及离散元等数值模拟方法目前被广泛应用在土木、机械、化工、冶金、农业、制药、环境等各项领域,用于模拟规则或不规则颗粒或者材料(不连续介质)的运动受力行为。离散元法的基本原理是以离散单元之间的相对位移为基本变量,通过力与相对位移的关系得到单元间的法向和切向的作用力并得到外力合力和合力矩。再根据牛顿第二定律,获得单元的加速度,进而积分得到单元的速度、位移等物理量。循环反复,加以表征颗粒单元的实时运动轨迹和受力行为。
离散元的数值模拟是通过数以万计次数的迭代完成的,每一次迭代过程的准确度直接影响到整个离散元模拟结果的精度。同时,数值模拟的过程中不可避免地需要对各种各样的实际问题进行假设、简化、甚至是改变。当所需要解决的实际问题是一个比较复杂的过程的时候,离散元的数值模拟精度和累计误差往往成为研究人员难以克服的一个瓶颈。因此,在解决大型工程问题,预测材料和结构的性能时,传统的离散元等数值模拟方法还往往依赖于昂贵、复杂、甚至是危险的现场实验作为后期校核的手段。
发明内容
本发明的目的在于提供一种简便、可靠、高效地提高离散元等数值模拟计算方法精度的基于软硬件实时交互优化的计算方法。
本发明的技术解决方案是:
一种基于软硬件实时交互优化的计算方法,其特征是:包括下列步骤:
(1)数值模拟颗粒的形成
a)选择欲进行模拟的实际颗粒并通过成像技术获得其三维视图及侧面、顶面、正面图;
b)通过对三维视图的正交延展,形成一个由这三个面构成的正交多面体;
c)取正交多面体在三个方向上重叠的公共部分,获得和原颗粒形状完全相似的模拟颗粒;
(2)智能颗粒传感设备制造
(3)“离散介质实时交互优化法则”及其应用
a)实验室颗粒实验:在一个既定容器中放入大小和形状不同的颗粒,并随机掺入若干已知形状和大小的智能颗粒,通过加入纵向重复荷载力,颗粒的运动状态、受力状态和位置也将随之改变,加载过程中,智能颗粒实时输出运动状态、位置、和受力情况的数据;
b)建立离散元模型:在计算机离散元模拟中建立与实验室相同的虚拟实验模型,其中智能颗粒的形状大小位置完全采用智能颗粒的初始数据;同时,整个实验模型被网格化,以智能颗粒作为基准点;
c)离散元模拟虚拟实验过程:
将实验模型以智能颗粒为基准点进行网格化;一旦基准点的智能颗粒的运动轨迹参数可知,那么其他颗粒和整个结构的运动轨迹参数也可以通过形函数而被准确获得;
利用智能颗粒获得特定时间点的颗粒运动轨迹,并及时与离散元计算结果加以比较,通过“卡曼过滤器”对离散元的预测结果进行实时更新,优化离散元对结点运动轨迹的预测精度;
根据所建立的形函数,对离散元在该时间点所预测的所有实验颗粒的运动轨迹和受力状态等参数进行优化;下一次迭代运算在优化后结果的基础上进行。
本发明简便、可靠、可高效地提高离散元等数值模拟计算方法精度;可以被广泛地应用于土木、机械、化工、冶金、农业、制药、环境等各项领域,对复杂、昂贵、耗时、危险的很多实验、现象进行模拟,解决工程实际问题,预测结构、材料、系统性能。
附图说明
下面结合附图和实施例对本发明作进一步说明。
图1是数值模拟颗粒形成过程示意图。
图2是掺入了智能颗粒的实验室颗粒实验示意图。
图3是基于“离散介质实时交互优化法则”的离散计算示意图。
具体实施方式
一种基于软硬件实时交互优化的计算方法,包括下列步骤:
(1)数值模拟颗粒的形成
a)选择欲进行模拟的实际颗粒并通过成像技术获得其三维视图及侧面、顶面、正面图;
b)通过对三维视图的正交延展,形成一个由这三个面构成的正交多面体;
c)取正交多面体在三个方向上重叠的公共部分,获得和原颗粒形状完全相似的模拟颗粒;
(2)智能颗粒传感设备制造;具体方法可以按照中国专利ZL201410094018.1进行;
例如:根据获得的需模拟颗粒的形态数据,采用3D打印技术,得到智能颗粒的外壳;在智能颗粒的内部放置可以实时获得颗粒运动形态、位置和表面受力状态的传感设备;
(3)“离散介质实时交互优化法则”及其应用
a)实验室颗粒实验:在一个既定容器中放入大小和形状不同的颗粒,并随机掺入若干已知形状和大小的智能颗粒。通过加入纵向重复荷载力(力的大小速率固定,但次数不断增加),颗粒的运动状态、受力状态和位置也将随之改变。加载过程中,智能颗粒可以实时输出运动状态(转动平移等)、位置、和受力情况的数据;
b)建立离散元模型:在计算机离散元模拟中建立与实验室相同的虚拟实验模型,其中智能颗粒的形状大小位置完全采用智能颗粒的初始数据;同时,整个实验模型被网格化,以智能颗粒作为基准点;
c)离散元模拟虚拟实验过程:
将实验模型以智能颗粒为基准点进行网格化;一旦基准点的智能颗粒的运动轨迹参数可知,那么其他颗粒和整个结构的运动轨迹参数也可以通过形函数而被准确获得;此原理类似于有限元的网格划分和性能预测。
利用智能颗粒获得特定时间点(可以是一次或多次循环迭代等)的颗粒运动轨迹,并及时与离散元计算结果加以比较,通过“卡曼过滤器(KALMANFILTER)”对离散元的预测结果进行实时更新,优化离散元对结点运动轨迹的预测精度;
根据所建立的形函数,对离散元在该时间点所预测的所有实验颗粒的运动轨迹和受力状态等参数进行优化;下一次迭代运算在优化后结果的基础上进行,从而最大化的减少误差累积,提高预测精度。
整个模拟优化过程可以循环反复,从而最终达到对系统的运动受力方式的准确预测。
Claims (1)
1.一种基于软硬件实时交互优化的计算方法,其特征是:包括下列步骤:
(1)数值模拟颗粒的形成
a)选择欲进行模拟的实际颗粒并通过成像技术获得其三维视图及侧面、顶面、正面图;
b)通过对三维视图的正交延展,形成一个由这三个面构成的正交多面体;
c)取正交多面体在三个方向上重叠的公共部分,获得和原颗粒形状完全相似的模拟颗粒;
(2)智能颗粒传感设备制造
(3)“离散介质实时交互优化法则”及其应用
a)实验室颗粒实验:在一个既定容器中放入大小和形状不同的颗粒,并随机掺入若干已知形状和大小的智能颗粒,通过加入纵向重复荷载力,颗粒的运动状态、受力状态和位置也将随之改变,加载过程中,智能颗粒实时输出运动状态、位置、和受力情况的数据;
b)建立离散元模型:在计算机离散元模拟中建立与实验室相同的虚拟实验模型,其中智能颗粒的形状大小位置完全采用智能颗粒的初始数据;同时,整个实验模型被网格化,以智能颗粒作为基准点;
c)离散元模拟虚拟实验过程:
将实验模型以智能颗粒为基准点进行网格化;一旦基准点的智能颗粒的运动轨迹参数可知,那么其他颗粒和整个结构的运动轨迹参数也可以通过形函数而被准确获得;
利用智能颗粒获得特定时间点的颗粒运动轨迹,并及时与离散元计算结果加以比较,通过“卡曼过滤器”对离散元的预测结果进行实时更新,优化离散元对结点运动轨迹的预测精度;
根据所建立的形函数,对离散元在该时间点所预测的所有实验颗粒的运动轨迹和受力状态等参数进行优化;下一次迭代运算在优化后结果的基础上进行。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610011844.5A CN105677983B (zh) | 2016-01-11 | 2016-01-11 | 基于软硬件实时交互优化的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610011844.5A CN105677983B (zh) | 2016-01-11 | 2016-01-11 | 基于软硬件实时交互优化的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105677983A true CN105677983A (zh) | 2016-06-15 |
CN105677983B CN105677983B (zh) | 2018-09-25 |
Family
ID=56299549
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610011844.5A Active CN105677983B (zh) | 2016-01-11 | 2016-01-11 | 基于软硬件实时交互优化的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105677983B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113720382A (zh) * | 2021-08-20 | 2021-11-30 | 中铁十四局集团大盾构工程有限公司 | 一种基于动态反分析计算和融合算法与智能化监测系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011145943A (ja) * | 2010-01-15 | 2011-07-28 | Japan Agengy For Marine-Earth Science & Technology | 粒子シミュレーション装置及び粒子シミュレーション方法 |
CN103850153A (zh) * | 2014-03-14 | 2014-06-11 | 黄炳勋 | 智能道渣 |
CN104504258A (zh) * | 2014-12-15 | 2015-04-08 | 广西科技大学 | 基于颗粒流方法的数值试样模拟方法 |
CN104978452A (zh) * | 2015-06-11 | 2015-10-14 | 昆明理工大学 | 一种模拟预测颗粒物料悬浮速度的方法 |
CN105138773A (zh) * | 2015-08-25 | 2015-12-09 | 北京建筑大学 | 基于离散元仿真平台的行人仿真方法及装置 |
-
2016
- 2016-01-11 CN CN201610011844.5A patent/CN105677983B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011145943A (ja) * | 2010-01-15 | 2011-07-28 | Japan Agengy For Marine-Earth Science & Technology | 粒子シミュレーション装置及び粒子シミュレーション方法 |
CN103850153A (zh) * | 2014-03-14 | 2014-06-11 | 黄炳勋 | 智能道渣 |
CN104504258A (zh) * | 2014-12-15 | 2015-04-08 | 广西科技大学 | 基于颗粒流方法的数值试样模拟方法 |
CN104978452A (zh) * | 2015-06-11 | 2015-10-14 | 昆明理工大学 | 一种模拟预测颗粒物料悬浮速度的方法 |
CN105138773A (zh) * | 2015-08-25 | 2015-12-09 | 北京建筑大学 | 基于离散元仿真平台的行人仿真方法及装置 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113720382A (zh) * | 2021-08-20 | 2021-11-30 | 中铁十四局集团大盾构工程有限公司 | 一种基于动态反分析计算和融合算法与智能化监测系统 |
CN113720382B (zh) * | 2021-08-20 | 2024-05-03 | 中铁十四局集团大盾构工程有限公司 | 一种基于动态反分析计算和融合算法和智能化监测系统 |
Also Published As
Publication number | Publication date |
---|---|
CN105677983B (zh) | 2018-09-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102207988B (zh) | 一种多自由度机械臂高效动力学建模方法 | |
CN106021828B (zh) | 一种基于格子‑玻尔兹曼模型的油层流体模拟方法 | |
CN110750933B (zh) | 一种耦合Lagrange质点和Euler方法的精确界面追踪处理方法 | |
US9311745B2 (en) | Systems and methods of analysis of granular elements | |
KR102014120B1 (ko) | 유체 역학 시스템에서의 경계 | |
JP7324932B2 (ja) | 動的計画コントローラ | |
CN107169236B (zh) | 一种基于有限元与离散元耦合的虚拟三轴试验仿真方法 | |
WO2017150626A1 (ja) | 粒子シミュレーション装置、粒子シミュレーション方法及び粒子シミュレーションプログラム | |
CN102865866A (zh) | 基于双星敏感器的卫星姿态确定方法及定姿误差分析方法 | |
CN104802167A (zh) | 基于MAPLESIM的Delta机器人快速建模方法 | |
CN108984829A (zh) | 堆石混凝土堆石体堆积过程的计算方法和系统 | |
CN108828959A (zh) | 一种新型的桥式起重机防摆与定位控制方法与装置 | |
CN105677983A (zh) | 基于软硬件实时交互优化的计算方法 | |
Lee | Developments in large scale discrete element simulations with polyhedral particles | |
CN104715133B (zh) | 一种待辨识对象的运动学参数在轨辨识方法和装置 | |
CN103294522A (zh) | 舰船雷达天线设备数值风洞仿真流程的简化方案 | |
CN110992475B (zh) | 水利水电工程大体积混凝土快速计算任意部位工程量的方法、系统及介质 | |
CN104657595A (zh) | 一种单颗粒曳力模型系数标定方法 | |
Huang et al. | Parallel Performance and Optimization of the Lattice Boltzmann Method Software Palabos Using CUDA | |
Fu et al. | Numerical model of three-dimensional motion of plate-type wind-borne debris based on quaternions and its improvement in unsteady flow | |
CN115329647B (zh) | 一种基于sph的筒形件旋压成形过程分析方法 | |
Masajedi et al. | Verification of bee algorithm based path planning for a 6DOF manipulator using ADAMS | |
CN109033665A (zh) | 一种基于cae仿真计算球阀转动扭矩的方法 | |
Li et al. | Sensitivity analysis approach to multibody systems described by natural coordinates | |
WO2024103241A1 (zh) | 一种基于物质点法的软体机器人仿真方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20191216 Address after: 210019 Room 209-38, 2nd floor, No. 18 East Street, Jialing River, Jianye District, Nanjing City, Jiangsu Province Patentee after: Nanjing Stardell Intelligent Technology Co., Ltd. Address before: Qidong City, Jiangsu province 226264 Nantong Huifeng Miao Jiao Industrial Park Nantong wandilai electromechanical Co Ltd Co-patentee before: Huang Hai Patentee before: Shen Shihui |