CN108108244A - 一种边坡强度折减系数多线程并行计算方法 - Google Patents

一种边坡强度折减系数多线程并行计算方法 Download PDF

Info

Publication number
CN108108244A
CN108108244A CN201711353469.3A CN201711353469A CN108108244A CN 108108244 A CN108108244 A CN 108108244A CN 201711353469 A CN201711353469 A CN 201711353469A CN 108108244 A CN108108244 A CN 108108244A
Authority
CN
China
Prior art keywords
reduction coefficient
max
slope
limit value
reduction
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
CN201711353469.3A
Other languages
English (en)
Other versions
CN108108244B (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.)
Central South University
Original Assignee
Central South 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 Central South University filed Critical Central South University
Priority to CN201711353469.3A priority Critical patent/CN108108244B/zh
Publication of CN108108244A publication Critical patent/CN108108244A/zh
Application granted granted Critical
Publication of CN108108244B publication Critical patent/CN108108244B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F9/00—Arrangements for program control, e.g. control units
    • G06F9/06—Arrangements for program control, e.g. control units using stored programs, i.e. using an internal store of processing equipment to receive or retain programs
    • G06F9/46—Multiprogramming arrangements
    • G06F9/50—Allocation of resources, e.g. of the central processing unit [CPU]
    • G06F9/5005—Allocation of resources, e.g. of the central processing unit [CPU] to service a request
    • G06F9/5027—Allocation of resources, e.g. of the central processing unit [CPU] to service a request the resource being a machine, e.g. CPUs, Servers, Terminals
    • 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
    • G06F9/00—Arrangements for program control, e.g. control units
    • G06F9/06—Arrangements for program control, e.g. control units using stored programs, i.e. using an internal store of processing equipment to receive or retain programs
    • G06F9/46—Multiprogramming arrangements
    • G06F9/48—Program initiating; Program switching, e.g. by interrupt
    • G06F9/4806—Task transfer initiation or dispatching
    • G06F9/4843—Task transfer initiation or dispatching by program, e.g. task dispatcher, supervisor, operating system
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2209/00—Indexing scheme relating to G06F9/00
    • G06F2209/50—Indexing scheme relating to G06F9/50
    • G06F2209/5018—Thread allocation

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种边坡强度折减系数多线程并行计算方法，将多线程并行运算应用于边坡强度折减系数计算，同时将FLAC3D命令流内嵌于Python脚本中，通过Python脚本同时调用多个FLAC3D应用程序进行多线程数值模拟计算。通过具体算例分析可以得到，本发明提供的边坡强度折减系数多线程并行计算方法与传统的二分法强度折减法相比，能够有效减少迭代次数及计算时间。

Description

一种边坡强度折减系数多线程并行计算方法
技术领域
本发明属于边坡稳定性分析技术领域，具体涉及一种边坡强度折减系数多线程并行计算方法的设计。
背景技术
强度折减法将边坡的安全系数定义为边坡刚好达到临界破坏状态时的折减系数，该过程一般采用二分法对边坡内聚力和内摩擦角进行迭代计算。随着FLAC3D编制强度折减法计算程序在边坡稳定性分析中的广泛采用，研究边坡强度折减法计算时间具有重要经济效益和实际意义。
随着计算机的迅速发展，有限元强度折减法在边坡工程中受到关注。其中较多学者选择FLAC3D软件进行边坡强度折减法数值模拟试验研究。薛雷以FLAC3D应用软件为计算平台，基于计算收敛性准则利用内嵌FISH语言二次开发出能够自动搜索安全系数的整体强度折减代码和局部折减强度代码。周元辅依托珍珠坝边坡采用FLAC3D对位移突变判据及塑性区判据进行分析，并与ANSYS计算结果进行比较得出：在滞滑带明确的情况下采用塑性区贯通率增量突变判据，在滞滑带未知的情况下采用位移增量突变判据。李跃以攀钢兰尖铁矿山第7排土场边坡为实例，采用FLAC3D软件进行三维强度折减法，数值模拟试验结果与现场调查及平面计算结果较为一致。张昊通过理论分析和FLAC3D强度折减法计算探讨了强度折减法边坡稳定性过程中安全系数和滑动面之间的关系。敬静利用FLAC3D软件建立节理岩质边坡桩基加固的数值模型，通过强度折减法计算边坡安全系数。采用FLAC3D软件进行边坡强度折减法稳定性分析时，判断边坡是否稳定是基于模型计算最终不平衡力比率满足10-5的要求。强度折减法对折减系数的折减一般采用普通的二分法进行迭代计算，采用该方法进行数值模拟计算需要较多的迭代次数及计算时间。因此，研究如何减少FLAC3D边坡稳定强度折减法过程中迭代次数及计算时间具有一定的实际意义。
发明内容
本发明的目的是为了解决现有技术中传统二分法强度折减算法进行数值模拟计算需要较多的迭代次数及计算时间的问题，提出了一种边坡强度折减系数多线程并行计算方法。
本发明的技术方案为：一种边坡强度折减系数多线程并行计算方法，包括以下步骤：
S1、根据边坡岩土体性质确定边坡强度的折减系数上限值Kmax和折减系数下限值Kmin。
折减系数上限值Kmax初始设置为20，折减系数下限值Kmin初始设置为0。
S2、通过线性插值形成折减系数等差序列(Kmin,K1,K2,...,KT,Kmax)；T为并行计算线程数，T为正整数且T＝N-1，N为计算机处理器的线程数。
S3、将FLAC3D命令流内嵌于Python脚本中，通过Python脚本同时调用T个FLAC3D应用程序，分别对折减系数为K1,K2,...,KT的数值模拟计算。
S4、根据数值模拟计算结果确定折减系数收敛区域，折减系数收敛的准则为不平衡力比率小于或等于10-5。
S5、根据折减系数收敛区域分别对折减系数上限值Kmax和折减系数下限值Kmin进行更新，得到更新后的折减系数上限值K′max和折减系数下限值K′min。
若K1不收敛，则K′min＝Kmin,K′max＝K1；
若Kt-1收敛，Kt不收敛，则K′min＝Kt-1,K′max＝Kt,t＝2,3,...,T；
若KT收敛，则K′min＝KT,K′max＝Kmax。
S6、判断是否满足K′max-K′min<η，若是则进入步骤S7，否则返回步骤S2，将Kmax和Kmin更新为K′max和K′min后进行迭代；η为给定计算误差，η根据边坡安全系数的设计值设置为0.001。
S7、根据更新后的折减系数上限值K′max和折减系数下限值K′min计算得到边坡强度折减系数K，计算公式为：
本发明的有益效果是：本发明将多线程并行运算应用于边坡强度折减系数计算，能够有效减少迭代次数及计算时间。同时本发明创新性地将FLAC3D命令流内嵌于Python脚本中，然后通过Python脚本调用FLAC3D应用程序，实现Python脚本与FLAC3D软件的交互式运算，为FLAC3D数值模拟计算研究提供了一种新的研究途径。
附图说明
图1所示为本发明实施例提供的一种边坡强度折减系数多线程并行计算方法流程图。
图2所示为本发明实施例提供的不同线程所需迭代次数曲线图。
图3所示为本发明实施例提供的FLAC3D数值模拟计算模型示意图。
图4所示为本发明实施例提供的4线程并行运算与二分法Kmax-Kmin变化趋势比较示意图。
具体实施方式
现在将参考附图来详细描述本发明的示例性实施方式。应当理解，附图中示出和描述的实施方式仅仅是示例性的，意在阐释本发明的原理和精神，而并非限制本发明的范围。
本发明实施例提供了一种边坡强度折减系数多线程并行计算方法，如图1所示，包括以下步骤S1-S7：
S1、根据边坡岩土体性质确定边坡强度的折减系数上限值Kmax和折减系数下限值Kmin。
本发明实施例中，折减系数上限值Kmax初始设置为20，折减系数下限值Kmin初始设置为0。
S2、通过线性插值形成折减系数等差序列(Kmin,K1,K2,...,KT,Kmax)；T为并行计算线程数，T为正整数。
本发明实施例中，根据计算机的性能来确定T值。设置T＝N-1，N为计算机处理器的线程数，例如计算机处理器为16线程，可将T设置为(16-1)＝15。
S3、将FLAC3D命令流内嵌于Python脚本中，通过Python脚本同时调用T个FLAC3D应用程序，分别对折减系数为K1,K2,...,KT的数值模拟计算。
Python语言设计明确、简单，在科学计算领域应用广泛。本发明实施例将FLAC3D命令流内嵌于Python脚本中，并由Python脚本同时调用多个FLAC3D应用程序实现并行运算，分析计算模型收敛情况，决定下一步计算。
为实现Python脚本调用FLAC3D应用程序，先定义一个flac(reduction)函数，该函数主要实现折减系数为reduction的FLAC3D数值模拟计算，其中flac为函数名，括号内参数为强度折减系数。通过预先安装多线程包threading实现调用多个FLAC3D应用程序。为避免多线程并行运算过程中出现线程间计算资源相互竞争或死锁现象，采用如下配套命令流：
Thread1＝threading.Thread(target＝flac,args＝(K1,))
Thread2＝threading.Thread(target＝flac,args＝(K2,))
……
ThreadT＝threading.Thread(target＝flac,args＝(KT,))
Thread1.start()
Thread2.start()
……
ThreadT.start()
Thread1.join()
Thread2.join()
……
ThreadT.join()
上述命令流中Thread1,Thread2,...ThreadT为线程名；K1,K2,...,KT为线程名为Thread1,Thread2,...,ThreadT所对应的FLAC3D应用程序计算的折减系数。由于Python脚本中对字母大小写敏感，上述配套命令流中不得随意更改其大小写状态，避免程序运行过程中出现错误。
S4、根据数值模拟计算结果确定折减系数收敛区域。
本发明实施例中，折减系数收敛的准则为不平衡力比率(节点平均内力与最大不平衡力的比值)小于或等于10-5。
S5、根据折减系数收敛区域分别对折减系数上限值Kmax和折减系数下限值Kmin进行更新，得到更新后的折减系数上限值K′max和折减系数下限值K′min。
若K1不收敛，则K′min＝Kmin,K′max＝K1；
若Kt-1收敛，Kt不收敛，则K′min＝Kt-1,K′max＝Kt,t＝2,3,...,T；
若KT收敛，则K′min＝KT,K′max＝Kmax。
S6、判断是否满足K′max-K′min<η，若是则进入步骤S7，否则返回步骤S2，将Kmax和Kmin更新为K′max和K′min后进行迭代；η为给定计算误差，本发明实施例中η根据边坡安全系数的设计值设置为0.001。
S7、根据更新后的折减系数上限值K′max和折减系数下限值K′min计算得到边坡强度折减系数K，计算公式为：
采用传统二分法对边坡强度折减系数进行迭代是不断将折减系数上限值Kmax和折减系数下限值Kmin之间的范围缩小为上一次迭代时的1/2，当折减系数上限值Kmax与下限值Kmin小于给定计算误差η时则停止计算。根据上述折减系数二分法迭代算法，则折减系数迭代次数Imin可用下式表示：
由公式(2)可知，采用二分法进行折减系数迭代计算时，至少要进行Imin次迭代，折减系数才能达到所需精度η：
式中，Ceil(·)表示向上取整运算。
根据多线程并行运算原理，采用FLAC3D软件进行计算时，数值模拟计算共需要迭代次数IT满足下式：
采用T线程并行运算时，至少要进行IT-min次迭代才能达到所需的计算精度：
通过比较公式(3)和公式(5)可知，采用二分法所需要的迭代次数是采用T线程并行运算的一个特例，当公式(5)中的T＝1时，即采用单线程运算时所需要的迭代次数与采用二分法所需的迭代次数相等。
令公式(5)中的T分别等于1,3,5,7，得到不同线程所需迭代次数曲线图，如图2所示。
对图2分析可知，在较小的时候，采用二分法(T＝1)与多线运算(T＝3,5,7)所需迭代次数相差较少。随着的不断增大，采用多线程并行运算的优点逐渐显现，即采用多线程并行运算所需迭代次数相比于二分法迭代降低较为明显。当时，采用传统二分法(T＝1，即单线程运算)所需迭代次数14次才能达到预期的计算精度，而采用3、5、7线程并行运算时，迭代次数分别为7、6、5次，所需迭代次数仅为传统二分法的35％-50％。
多线程并行运算时，应先采用公式(5)进行试算，从而确定合适的线程数。如图2所示，当时，采用5线程与7线程并行运算所需迭代次数均为5次，实际应用过程中，过多线程进行应用程序并行运算可能造成计算机内部进程间计算资源相互竞争，导致计算时间延长，这种情况下建议采用5线程并行运算。因此，采用多线程强度折减法计算时，应根据迭代次数及计算机性能选择并行运算程序个数，不能盲目提高并行运行线程个数以求快速得到计算结果。
下面以一个具体实例对本发明实施例提供的一种边坡强度折减系数多线程并行计算方法针对于传统二分法强度折减算法的优势作详细描述。
本发明实施例并行运算采用的应用软件为FLAC3D V3.0，Python版本为Python2.7.5，所使用Python编译器为Pycharm Community 4.5.3(开源，无需付费)。并预先安装数值模拟试验中需要的Python软件包如表1所示：
表1
以上软件包均遵守GNU协议，即软件包开源无需付费且用户可随意更改这些软件包中的程序。计算机硬件对多线程并行运算执行效率影响较大，为此表2给出数值模拟实验的主要硬件设备及系统等参数指标：
表2
本发明实施例以均质边坡作为分析对象，根据平面应变建立计算模型(模型厚度方向取1m)。边坡的高度20m，坡角45°，其余形态参数如图3所示。模型边界条件设置为：边坡的左右和前后边界约束法向位移，底部边界为位移固定约束，其他面为自由边界。视边坡为理想线性弹性体，外荷载为岩土体自重，计算时先将岩土体视为弹性体，待自重应力施加平衡后，在将模型调整为Mohr-Coulomb模型进行弹塑性分析。
模型网格采用均匀划分，该边坡模型共816个单元，1850个节点。数值模拟采用Mohr-Coulomb模型(考虑拉伸截断tension cutoff)、非关联流动法则，边坡模型的主要力学参数如表3所示：
表3
其中密度ρ、体积模量K及剪切模量G通过公式(6)～(8)计算得到：
其中g为重力加速度。
数值模拟计算过程总运行30000步，其计算收敛准则为最终不平衡力比率小于或等于10-5。初始折减系数下限值Kmin和折减系数上限值Kmax分别为0、20，给定计算误差η为0.001。
综合考虑数值模拟试验计算平台及硬件设施的实际情况，本发明实施例计算机的内核是3.5GHz，向上取整采用4线程并行运算。采用4线程法与二分法数值模拟计算得到的边坡安全系数F均为1.034。因此，采用多线程数值模拟试验结果与传统的二分法数值模拟试验结果相同，不存在误差，4线程并行运算计算结果准确可靠。
由于本次数值模拟试验初始折减系数下限值Kmin和折减系数上限值Kmax分别为0、20，给定计算误差η为0.001，因此，根据公式(3)和公式(5)，采用传统二分法和4线程法所需迭代次数分别为15次和7次，这与数值模拟计算结果吻合。采用4线程强度折减法所需迭代次数是传统二分法强度折减法迭代次数的46％。因此，采用4线程并行运算相比于传统二分法能有效减少数值计算过程中的迭代次数。
图4给出数值模拟计算过程中分别采用二分法和4线程并行运算迭代的Kmax-Kmin变化趋势。由图4可知，采用4线程并行运算在计算过程中的收敛速度明显大于采用二分法的收敛速度，且所需迭代次数也较少。相比二分法将折减系数的计算范围不断划分为2个区域，采用4线程并行运算时，每进行一次迭代都将折减系数的范围划分为5个相等的区域。每次迭代计算结束后，折减系数的范围缩小到原来的1/5。采用4线程并行运算，折减系数缩小范围的速率要快于二分法。因此，采用4线程并行运算能加快折减系数的收敛。
同时，在数值模拟计算过程中，通过采用Python脚本中time模块对二分法和4线程法数值模拟运算时间分别进行统计，所需计算时间分别为3123.9秒和1735.6秒(不同硬件性能计算机相应计算时间可能不同)，采用4线程并行运算所需时间是传统二分法强度折减法所需时间的55％。采用4线程并行运算能有效减少强度折减法迭代计算过程中所需要的计算时间。
综上所述，相比于传统的二分法折减系数迭代，采用4线程并行运算能有效减少迭代次数及数值模拟分析计算所需要时间，其中迭代次数为传统二分法的46％，计算时间为二分法的55％。从迭代次数及总的计算时间角度考虑，4线程并行运算计算优势较为明显。
本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合，这些变形和组合仍然在本发明的保护范围内。

Claims (7)

1.一种边坡强度折减系数多线程并行计算方法，其特征在于，包括以下步骤：
S1、根据边坡岩土体性质确定边坡强度的折减系数上限值Kmax和折减系数下限值Kmin；
S2、通过线性插值形成折减系数等差序列(Kmin,K1,K2,...,KT,Kmax)；T为并行计算线程数，T为正整数；
S3、将FLAC3D命令流内嵌于Python脚本中，通过Python脚本同时调用T个FLAC3D应用程序，分别对折减系数为K1,K2,...,KT的数值模拟计算；
S4、根据数值模拟计算结果确定折减系数收敛区域；
S5、根据折减系数收敛区域分别对折减系数上限值Kmax和折减系数下限值Kmin进行更新，得到更新后的折减系数上限值K′max和折减系数下限值K′min；
S6、判断是否满足K′max-K′min<η，若是则进入步骤S7，否则返回步骤S2，将Kmax和Kmin更新为K′max和K′min后进行迭代；η为给定计算误差；
S7、根据更新后的折减系数上限值K′max和折减系数下限值K′min计算得到边坡强度折减系数K。
2.根据权利要求1所述的边坡强度折减系数多线程并行计算方法，其特征在于，所述步骤S1中折减系数上限值Kmax初始设置为20，折减系数下限值Kmin初始设置为0。
3.根据权利要求1所述的边坡强度折减系数多线程并行计算方法，其特征在于，所述步骤S2中并行计算线程数T＝N-1，N为计算机处理器的线程数。
4.根据权利要求1所述的边坡强度折减系数多线程并行计算方法，其特征在于，所述步骤S4中折减系数收敛的准则为不平衡力比率小于或等于10-5。
5.根据权利要求1所述的边坡强度折减系数多线程并行计算方法，其特征在于，所述步骤S5具体为：
若K1不收敛，则K′min＝Kmin,K′max＝K1；
若Kt-1收敛，Kt不收敛，则K′min＝Kt-1,K′max＝Kt,t＝2,3,...,T；
若KT收敛，则K′min＝KT,K′max＝Kmax。
6.根据权利要求1所述的边坡强度折减系数多线程并行计算方法，其特征在于，所述步骤S6中计算误差η根据边坡安全系数的设计值设置为0.001。
7.根据权利要求1所述的边坡强度折减系数多线程并行计算方法，其特征在于，所述步骤S7中边坡强度折减系数K的计算公式为：
<mrow> <mi>K</mi> <mo>=</mo> <mfrac> <mrow> <msubsup> <mi>K</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> <mo>&amp;prime;</mo> </msubsup> <mo>-</mo> <msubsup> <mi>K</mi> <mrow> <mi>m</mi> <mi>i</mi> <mi>n</mi> </mrow> <mo>&amp;prime;</mo> </msubsup> </mrow> <mn>2</mn> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>.</mo> </mrow>
CN201711353469.3A 2017-12-15 2017-12-15 一种边坡强度折减系数多线程并行计算方法 Active CN108108244B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711353469.3A CN108108244B (zh) 2017-12-15 2017-12-15 一种边坡强度折减系数多线程并行计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711353469.3A CN108108244B (zh) 2017-12-15 2017-12-15 一种边坡强度折减系数多线程并行计算方法

Publications (2)

Publication Number Publication Date
CN108108244A true CN108108244A (zh) 2018-06-01
CN108108244B CN108108244B (zh) 2021-09-28

Family

ID=62217344

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711353469.3A Active CN108108244B (zh) 2017-12-15 2017-12-15 一种边坡强度折减系数多线程并行计算方法

Country Status (1)

Country Link
CN (1) CN108108244B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109359368A (zh) * 2018-10-02 2019-02-19 东北电力大学 一种均质边坡稳定性评价的双折减系数强度折减法
CN111008465A (zh) * 2019-11-25 2020-04-14 中南大学 考虑粘聚力与内摩擦角贡献度的边坡安全系数计算方法
CN111931273A (zh) * 2020-07-16 2020-11-13 中南大学 边坡安全系数理论值的快速计算方法
CN113806099A (zh) * 2021-09-13 2021-12-17 北京计算机技术及应用研究所 一种基于二值计算的多核自旋锁设计方法
CN114398166A (zh) * 2022-01-27 2022-04-26 重庆大学 基于二分法的分布式计算任务调度方法及设备
CN116991679A (zh) * 2023-09-26 2023-11-03 四川华鲲振宇智能科技有限责任公司 一种内存压力测试时自动分配最佳线程数的方法

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090172773A1 (en) * 2005-02-01 2009-07-02 Newsilike Media Group, Inc. Syndicating Surgical Data In A Healthcare Environment
CN103543468A (zh) * 2013-10-28 2014-01-29 北京大学 一种地震正演模拟的方法和系统
CN103914265A (zh) * 2014-04-09 2014-07-09 江苏物联网研究发展中心 集群细粒度内存管理方法
CN104406605A (zh) * 2014-10-13 2015-03-11 中国电子科技集团公司第十研究所 机载多导航源综合导航仿真系统
CN104615826A (zh) * 2015-02-02 2015-05-13 同济大学 一种复杂岩土工程可靠性分析装置和方法
US9116524B1 (en) * 2014-04-21 2015-08-25 Trutrak Flight Systems, Inc. System and method for vertical navigation based on GPS waypoints and autopilot programming
CN105824987A (zh) * 2016-03-09 2016-08-03 浙江大学 一种基于遗传算法的风场特征统计分布模型建立方法
CN105893483A (zh) * 2016-03-29 2016-08-24 天津贝德曼科技有限公司 大数据挖掘过程模型总体框架的构造方法
CN106599444A (zh) * 2016-12-09 2017-04-26 山东大学 一种盐渍土路基上盐分隔断处置方法
CN107015861A (zh) * 2016-11-07 2017-08-04 珠江水利委员会珠江水利科学研究院 一种基于Fork/Join框架的梯级水库群优化调度多核并行计算设计方法
CN107015852A (zh) * 2016-06-15 2017-08-04 珠江水利委员会珠江水利科学研究院 一种大规模梯级水电站群多核并行优化调度方法
CN107220401A (zh) * 2017-04-12 2017-09-29 中国地质大学（武汉） 基于并行蒙特卡洛法的边坡可靠性参数获取方法及装置

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090172773A1 (en) * 2005-02-01 2009-07-02 Newsilike Media Group, Inc. Syndicating Surgical Data In A Healthcare Environment
CN103543468A (zh) * 2013-10-28 2014-01-29 北京大学 一种地震正演模拟的方法和系统
CN103914265A (zh) * 2014-04-09 2014-07-09 江苏物联网研究发展中心 集群细粒度内存管理方法
US9116524B1 (en) * 2014-04-21 2015-08-25 Trutrak Flight Systems, Inc. System and method for vertical navigation based on GPS waypoints and autopilot programming
CN104406605A (zh) * 2014-10-13 2015-03-11 中国电子科技集团公司第十研究所 机载多导航源综合导航仿真系统
CN104615826A (zh) * 2015-02-02 2015-05-13 同济大学 一种复杂岩土工程可靠性分析装置和方法
CN105824987A (zh) * 2016-03-09 2016-08-03 浙江大学 一种基于遗传算法的风场特征统计分布模型建立方法
CN105893483A (zh) * 2016-03-29 2016-08-24 天津贝德曼科技有限公司 大数据挖掘过程模型总体框架的构造方法
CN107015852A (zh) * 2016-06-15 2017-08-04 珠江水利委员会珠江水利科学研究院 一种大规模梯级水电站群多核并行优化调度方法
CN107015861A (zh) * 2016-11-07 2017-08-04 珠江水利委员会珠江水利科学研究院 一种基于Fork/Join框架的梯级水库群优化调度多核并行计算设计方法
CN106599444A (zh) * 2016-12-09 2017-04-26 山东大学 一种盐渍土路基上盐分隔断处置方法
CN107220401A (zh) * 2017-04-12 2017-09-29 中国地质大学（武汉） 基于并行蒙特卡洛法的边坡可靠性参数获取方法及装置

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
SHIDAZAI1951: ""边坡稳定性分析报告"", 《HTTPS://JZ.DOCIN.COM/P-1723683990.HTML》 *
刘晓: ""汶川地震区斜坡动力反应研究"", 《中国博士学位论文全文数据库 基础科学辑》 *
武瑞婵: ""基于Java多线程的预处理迭代并行求解器"", 《山西大同大学学报（自然科学版）》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109359368A (zh) * 2018-10-02 2019-02-19 东北电力大学 一种均质边坡稳定性评价的双折减系数强度折减法
CN111008465A (zh) * 2019-11-25 2020-04-14 中南大学 考虑粘聚力与内摩擦角贡献度的边坡安全系数计算方法
CN111008465B (zh) * 2019-11-25 2022-09-20 中南大学 考虑粘聚力与内摩擦角贡献度的边坡安全系数计算方法
CN111931273A (zh) * 2020-07-16 2020-11-13 中南大学 边坡安全系数理论值的快速计算方法
CN113806099A (zh) * 2021-09-13 2021-12-17 北京计算机技术及应用研究所 一种基于二值计算的多核自旋锁设计方法
CN114398166A (zh) * 2022-01-27 2022-04-26 重庆大学 基于二分法的分布式计算任务调度方法及设备
CN114398166B (zh) * 2022-01-27 2024-05-28 重庆大学 基于二分法的分布式计算任务调度方法及设备
CN116991679A (zh) * 2023-09-26 2023-11-03 四川华鲲振宇智能科技有限责任公司 一种内存压力测试时自动分配最佳线程数的方法
CN116991679B (zh) * 2023-09-26 2023-12-08 四川华鲲振宇智能科技有限责任公司 一种内存压力测试时自动分配最佳线程数的方法

Also Published As

Publication number Publication date
CN108108244B (zh) 2021-09-28

Similar Documents

Publication Publication Date Title
CN108108244B (zh) 一种边坡强度折减系数多线程并行计算方法
CN108920739B (zh) 一种考虑损伤累积效应的材料本构模型数值分析方法
CN103077006B (zh) 一种基于多线程的长事务并行执行方法
US20170168779A1 (en) Automated generation of memory consumption aware code
US20200117846A1 (en) Information processing device, information processing method, and program
Park et al. Parametric design process of a complex building in practice using programmed code as master model
US9176732B2 (en) Method and apparatus for minimum cost cycle removal from a directed graph
Chiang et al. Determinism and reproducibility in large-scale HPC systems
Bonenfant et al. Early WCET prediction using machine learning
Tillmann et al. Application-independent autotuning for GPUs
Dubey et al. A priority based job scheduling algorithm using IBA and EASY algorithm for cloud metaschedular
Wulf et al. Increasing the throughput of pipe-and-filter architectures by integrating the task farm parallelization pattern
He et al. New speedup algorithms for nonlinear dynamic time history analysis of supertall building structures under strong earthquakes
CN113961439A (zh) 一种基于smt方法的dag任务wcrt计算方法
Garcia et al. A dynamic load balancing approach with SMPSuperscalar and MPI
Cools et al. On rounding error resilience, maximal attainable accuracy and parallel performance of the pipelined Conjugate Gradients method for large-scale linear systems in PETSc
Racz et al. Parallelizing boundary surface computation of Chua's circuit
CN111274667A (zh) 一种跨尺度材料计算软件集成计算系统及方法
Abudalfa et al. Comparative study on behavior-based dynamic branch prediction using machine learning
Kalling et al. Accelerating the numerical simulation of magnetic field lines in tokamaks using the GPU
CN111124694A (zh) 一种基于petri网的可达图的死锁检测和解决方法
Chiu et al. An Efficient Task-Parallel Pipeline Programming Framework
He et al. Efficient response time bound for typed dag tasks
Chiu et al. Pipeflow: An efficient task-parallel pipeline programming framework using modern C++
Wu et al. Research on OpenMP model of the parallel programming technology for homogeneous multicore DSP

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