[数学][编程][混沌] 虫口模型的数值模拟 [Numerical simulation of insect population model]
实干、实践、积累、思考,创新。 1976 年 , 美国数学家梅 (May .R ) 在美国 《 自然 》 杂志上发表的题为 “ 具有复杂的动力学的简单模型 ” 文章中指出 , 在生态学中一些非常简单的确定性的数学模型却能产生看似随机的行为。如虫口模型:xn+ 1= μxn ( 1 -x n ),其中 xn 是第n 年的虫口数,xn+ 1 …
实干、实践、积累、思考、创新! Tag: 结构工程博士 结构工程师 伪程序员 结构抗震 地震工程 超限设计 软件定制 环评减振 振动控制 减隔震 施工过程模拟 小品钢结构 有限元研发 参数化设计 大震弹塑性
实干、实践、积累、思考,创新。 1976 年 , 美国数学家梅 (May .R ) 在美国 《 自然 》 杂志上发表的题为 “ 具有复杂的动力学的简单模型 ” 文章中指出 , 在生态学中一些非常简单的确定性的数学模型却能产生看似随机的行为。如虫口模型:xn+ 1= μxn ( 1 -x n ),其中 xn 是第n 年的虫口数,xn+ 1 …
实干、实践、积累、思考,创新。 小伙伴试验数据出问题,找我们处理。 用小伙伴的试验数据做的 SawtoothRemove([软件][试验][研究] SawtoothRemove: Remove the Sawtooth Patterns in Your Test Data [剔除试验数据中的锯齿]) 案例。直接上图,看修正过程和结果。 (1)将滞回曲线数据导入SawtoothRemove,从下图可以看见,滞回曲线有非常多的局部跳动,这种跳动也可以理解为一种锯齿。 (2)点击Analyze Loop分析滞回环,分析完毕我们可以在LOOP No.下拉菜单中查看各个滞回环,各个滞回环会在右边的绘图面板中显示。 (3)可以单独显示某个滞回环曲线,可以发现单个滞回环曲线有非常多的锯齿。 (4)直接点击 “修正数据的锯齿” 按钮,可一键修正锯齿,如下图所示,修正的曲线如红色所示,变得非常平滑,原来的锯齿修正了。 (5)可单独显示修正后的曲线,如下图。 (6) 最后点击 “输出EXCEL曲线图” ,可将曲线图直接输出到EXEL表,方便使用者进行后续数据的处理。 另外网站还提供了其他几款用于修正试验数据的工具,感兴趣的小伙伴可以看看: NoiseRemoval:http://www.jdcui.com/?p=15046 …
实干、实践、积累、思考,创新。 如题 如何将n^m个可能枚举出来?通常枚举实现,直觉想到的就是用循环。 对于n和m很大的时候,写循环显然不可能。另外且当n和m都是运行时确定时,不管n和m是大还是小,都无法用循环了,因为用循环就意味着代码写死了。 最后想了一下,对于这种情况,可以用进制的方式进行实现。 n^m,相当于一个m维的数组,每个数组的取值可能有n个,那么m维数组,总的可能数就有n^m个。 举几个例子 2^m,可以表示为一个2进制位数 4^m,可以表示为一个4进制位数 具体的 2^2,一共有4种可能,具体情况如下: 1: 0 0 2: 1 0 3: 0 1 4: 1 1 3^2,一共有9种可能,具体情况如下: 1: 0 0 2: 1 0 …
实干、实践、积累、思考、创新。 博客渐渐成了自己的一个对外的笔记本。以后增加一种【图】的记录方式。 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号
实干、实践、积累、思考,创新。 最近研究了一下粒子群算法,上一篇博文《 [编程训练][软件] 粒子群优化算法求解无约束最优化问题 [Particle Swarm Optimization Algorithm for Solving Unconstrained Optimization Problems] 》进行了编程实现。再测试过程中用了几个Benchmark函数,都是一些复杂的高纬函数,通常有非常多的局部波峰,用来测试算法的全局最优搜索能力。这些函数都比较漂亮,截几个图欣赏一下。 0: Sphere 1: Eggcrate 2: Rosenbrock 3: Ackley 4: Schaffer 5: Rastrigrin 6: Cigar 7: Griewangk 8: …
实干、实践、积累、思考,创新。 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 假期花了几天研究粒子群优化算法,顺便写了这个小工具,简单测试无约束最优化问题,也顺便训练一下编程。 软件自带几个经典的PSO粒子群优化算法测试函数,具体包括: 0: Sphere 1: Eggcrate 2: Rosenbrock 3: Ackley 4: Schaffer 5: Rastrigrin 6: Cigar 7: Griewangk 8: Schwefel 几个函数的图像可以在这个博文查看:[笔记][智能算法]几个经典的PSO粒子群优化算法测试函数 (Benchmark …
实干、实践、积累、思考,创新。 其他程序用到,就研究一下。 导入BMP图像,设置70%缩放: 导入BMP图像,设置50%缩放: 导入BMP图像,设置20%缩放: 相关博文 ( Related Topics) [01] [编程][算法][几何] 编程训练: 凸包问题 [Convex-Hull Problem] [02] [编程][算法][几何] 编程训练: 离散点的闭合路径 [Closed Path for points] [03] [编程][算法][几何] 编程训练: 3次B样条曲线 [Cubic BSpline] [04] …
实干、实践、积累、思考,创新。 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 曲线图动画及Gif动图生成工具,可以用于制作一些数据演示。 程序界面 ( Program Interface ) 程序演示视频: 以下是测试的几个例子 生成的GIF动图。 (1)滞回曲线 (2)地震波时程 (3)SIN函数 (4)螺旋线 下载 ( Download ) ( 如果您需要使用这个软件,请在这给网站捐助:http://www.jdcui.com/?page_id=4813,捐助不少于【200RMB】,捐助后在评论区留下您的评论,软件将通过网站管理员QQ:3014479529给您发送,敬请尊重劳动成果 !!!) 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号
实干、实践、积累、思考、创新。 来自小伙伴 吴金诚 的分享 0 钢筋对筏板抗冲切承载力的贡献 前阵子做的一个项目,遇到一个比较有意思的问题,想跟大家分享一下。项目采用的是桩筏基础,在基础验算的时候,发现塔楼范围内有个别柱下的冲切算不过。我们常用的解决方法主要有两种:一是加厚塔楼范围内的筏板厚度;二是在柱下局部加厚,做成柱墩。但是业主觉得上述第一种方法增加成本比较多,第二种施工又不是很方便。因此,我们只能“另谋出路”。 其实,目前常用结构常用分析软件在进行冲切计算时,均只考虑混凝土对抗冲切承载力的贡献,不考虑钢筋的贡献。而《混凝土结构设计规范》(以下统称“《混规》”)6.5.3给我们提出考虑钢筋时的抗冲切承载力计算方法,而我们的项目最终也是通过配置抗冲切钢筋的方法解决局部抗冲切不满足的问题。下面以一个配置抗冲切钢筋的算例,介绍一下抗冲切钢筋的计算过程以及节点大样。 1.筏板抗冲切箍筋布置方式 我们提出了以下两种抗冲切箍筋的布置方式。 (1)根据《混规》9.1.11,我们设计了以下抗冲切钢筋的布置方式,抗冲切钢筋布置大样如图1所示,相关平面图图例说明如图2所示。 其中①号筋为抗冲切箍筋,计算时仅考虑冲切面以内的抗冲切钢筋对抗冲切承载力的贡献; ②号钢筋为架立筋。下述抗冲切钢筋的验算依据此类布置方式。 图1 抗冲切钢筋布置大样一 图2 抗冲切箍筋平面图图例说明 (2)第一种抗冲切钢筋的布置方法虽然有规范的依据,但在施工上可能不太方便。因此,我们提出了第二种布置方法,详图3。此种布置方法是在筏板底部与顶部通常筋之间布置抗冲切拉筋,布置范围为冲切面以内。抗冲切拉筋的计算原理与抗冲切箍筋类似,此处不再详述。 图3 抗冲切钢筋布置大样二 2.计算条件 (a)Z1范围底板平面图 (b)Z1冲切验算结果 图4 Z1抗冲切计算结果与抗冲切箍筋布置 由图4可知,未配置抗冲切箍筋时,Z1冲切系数为0.88,小于1,不能满足抗冲切承载力要求,因此配置抗冲切箍筋C18@150。下面对进行验算配置抗冲切箍筋后的柱下筏板冲切承载力进行验算,以下为计算参数(计算参数均可从盈建科柱冲切计算书中读取): ft = 1.57N/mm2;h0 = …
实干、实践、积累、思考、创新。 新写的一本书快要出版了,书的题目是《结构地震反应分析 编程与软件应用》,这是继《PERFORM-3D原理与实例》、《有限单元法 编程与软件应用》后,www.jdcui.com的第三本书,主题是结构地震动力计算,依然是一本实用的基础的书籍,将理论、编程及软件应用结合,通过典型的小例子手把手教大家编程及使用软件。终于在年前拿到了出版社的样稿,印刷效果还是不错。 查了一下已经可以在京东购买了。以下是几个购买连接: 中国建筑工业出版社官方旗舰店: https://item.jd.com/10043957730709.html 中国建筑书店有限责任公司图书专营店:https://item.jd.com/10043957749812.html 顺带也贴一下目录吧: 年后再详细介绍这本书!!! 😛 🙂 公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号
实干、实践、积累、思考,创新。 用C#编写程序,查询Access数据库时,出现异常“IErrorInfo.GetDescription 因 E_FAIL(0x80004005) 而失败。” 后检查发现,SELECT 语句查询语句的 表的名称Section有冲突,将 Section加上”[ ]”修饰即可。 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号
实干、实践、积累、思考、创新! 蒙特· 卡罗方法是一种通过概率来得到问题近似解的方法。其实之前也测试过,用在这个投针求PI([数学][概率] Buffon’s Needle problem [蒲丰投针问题]),最近又搞起算法研究,顺带再整体各种东西都测算一次吧。 蒙特· 卡罗方法求圆周率PI 网上很多例子。直接上图吧。 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号
实干、实践、积累、思考、创新。 程序图标 ( Program Icon ): 程序介绍 ( Program Introduction) 隔墙线荷载计算、隔墙面荷载计算。 Partition wall line load calculation, partition wall area load calculation 程序界面 ( Program Interface ) 下载 ( Download ) ( 软件是免费的,如果需要这个软件,请在评论区留下您的评论,软件会发送到您的评论邮箱,敬请尊重劳动成果 !!) …
实干、实践、积累、思考、创新。 写在前面: 受 徐自国 博士 的idea的启发,写了这个地震波反应谱演化工具,idea来源于徐总。主要用于分析地震波反应谱随选取的地震波计算时间长短的一个演化。可以有不同的思考。感兴趣的朋友可以下来看看。以下是软件简介。 程序图标( Program Icon ) 程序介绍 ( Program Introduction) 主要用于分析地震波反应谱随选取的地震波计算时间长短的一个演化。还有一些其他思考。后续再陆续放上来。 程序视频: 程序截图: 程序下载 ( Program Download ): 下载: Spectr_Evolution(地震波反应谱演化).rar ( PS. 程序免费下载 ,欢迎大家给我提意见 ,敬请尊重劳动成果 ) 相关话题 ( Related Topics) [01]. [Tool] …
实干、实践、积累、思考,创新! 直接看图。 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号
实干、实践、积累、思考,创新! 之前写的小工具,整理分享给大家。 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 根据 GB50010-2010 混凝土结构设计规范 计算钢筋混凝土梁的构造腰筋。 规范条文: 6.3.1 矩形、T形和I形截面受弯构件,截面腹板的高度hw:矩形截面,取有效高度;T形截面,取有效高度减去翼缘高度;I形截面,取腹板净高。 9.2.13 梁的腹板高度hw不小于450mm时,在梁的两个侧面应沿高度配置纵向构造钢筋。每侧纵向构造钢筋(不包括梁上、下部受力钢筋及架立钢筋)的间距不宜大于200mm,截面面积不应小于腹板截面面积(bhw)的0.1%,但当梁宽较大时可以适当放松。 此处,腹板高度hw按本规范第6.3.1条的规定取用。 程序界面 ( Program Interface ) 下载 ( Download ) ( 如果需要这个软件,请在评论区留下您的评论,软件会发送到您的评论邮箱,敬请尊重劳动成果 !!) (大家注意评论邮箱要填写正确,后台自动发送软件,如果邮箱填写错误,会发送不到。) 相关软件资料 ( …
实干、实践、积累、思考,创新! 最近换电脑,重新装个新点的abaqus。 fortran 及 vs都装好后,运行Abaqus Verification,发现与子程序的东西验证均报错:‘ifort’不是外部或内部命令,也不是可运行的程序。 百度查了很多解决方法似乎解决不了。 问题肯定是处在编译器路径问题上。 测试了一下,最后在 launcher.bat 文件开头补充一下 fortran的路径就行,改为如下: @echo off call ipsxe-comp-vars.bat intel64 vs2019 @echo off call “E:\SIMULIA\Commands\abq2020.bat” %* 相关话题 ( Related Topics) [01]. Plate With Hole Stress …
[Software][软件] 梁纵向钢筋单排最大根数。根据混凝土规范GB 50010-2010计算梁纵向钢筋单排最大根数。A Program for calculating the maximum number of beam longitudinal reinforcing bars based on Chinese concrete design code GB 50010-2010.
实干、实践、积累、思考、创新! 粘弹性阻尼器是一种有效的减震控制装置,主要依靠粘弹性材料的滞回耗能,为结构提供附加刚度和阻尼,减小结构动力反应,从而实现减震目标。粘弹性阻尼器既能提供刚度,也能提供阻尼,其典型的滞 回曲线为椭圆形,具有良好的耗能性能。下图 所示是粘弹性阻尼器的滞回环形状示意。 粘弹性阻尼器,恢复力由两部分构成,即速度相关部分和位移相关部分,具体如下 NSDOF提供了单自由度主体结构及黏弹性阻尼器并联的动力时程分析模型,以下做两个测算。 第一个算例,结果部附加任何阻尼器,结构刚度取10,阻尼比取0.05,对应的阻尼系数为0.15915,点击运行计算可获得对应的力、位移响应结果及能量图。由能量图可见,此时主要包含3种能量:动能+应变能+阻尼耗能 第二个算例,主体结构的阻尼取0,刚度取一个很小的值0.0001,几乎可忽略,同时设置黏弹性阻尼器的阻尼系数为0.15915,与算例1中的主体结构的阻尼系数一致,黏弹性阻尼器的刚度取10,鱼算例1中的主体结构的刚度一致。因此可知,算例2的计算结果应该与算例1是一致的,相当于用一个黏弹性阻尼器去等效一个单自由度系统。计算结果如下图所示。 由结果可见,算例2的位移、速度、加速度响应均是与算例1一致的。不同的是,算例2的能量图,此时主要包含2种能量:动能+黏弹性阻尼器的耗能,因为用黏弹性阻尼器等效算例1的主体结构,同时也可以发现,黏弹性阻尼器的耗能等于算例1主体结构的应变能及阻尼耗能的叠加。 另外从滞回曲线可以发现,黏弹性阻尼器的滞回曲线是椭圆。此时黏弹性阻尼器的滞回曲线也等于算例1主体结构他弹性恢复力滞回曲线及阻尼力滞回曲线的叠加。 相关话题 ( Related Topics) [01]. [Tool] SPECTR – A program for Response Spectra Analysis [反应谱计算程序] [02]. [程序][Tool] Ground Motion Selection [强震记录选取] [03]. [程序][软件]Ground Motion …
实干、实践、积累、思考,创新。 典型核心筒斜墙转换区的受力如上图所示。 在斜墙斜率取值合理的情况下,斜墙的竖向力的水平力分量通常主要通过起止层楼板、剪力墙及连梁自平衡抵抗。核心筒外梁板分担的力实质上是非常小的。 斜墙起始层,核心筒内受拉,核心筒外受压。 斜墙终止层,核心筒内受压,核心筒外受拉。 斜墙中间楼层,通常不分担水平力。 通常情况,斜墙起始层楼板受力较大,终止层受力较小。 由于起止层楼板的重要,因此,通常起止层楼板需要加厚,配筋需要加强。 沿斜墙方向,由于拉力需要靠内墙及连梁形成的自平衡体系抵抗,因此需要复核连梁及内墙的水平拉力,如果拉力较大,需要加强,比如在关键传力连梁上设置贯通剪力墙的钢骨。 斜墙主要是竖向荷载传递上特殊,因此需特别关注竖向荷载下斜墙区域的受力特性(楼板应力、斜墙位移、连梁拉力、内墙抗剪、斜墙面外受力)。 相关博文( Related Topics) [01] [工具][软件][规范] 广东省标准《高层建筑混凝土结构设计规范》反应谱计算工具 [02] [结构力学][结构设计] 两端固支梁弯矩为0点距端部的距离 [03] [抗震][结构设计][规范] 非抗震设计情况下混凝土柱的“轴压比”可达多大? [04] [结构设计][动力学] YJK中CQC振型组合地震力的复核 [05] [结构设计][楼梯] 混凝土楼梯施工图笔记 [06] [结构][设计][规范] …
坚持实干、坚持实践、坚持积累、坚持思考,坚持创新。 随后更新。。。。。。 相关博文( Related Topics) [01] [工具][软件][规范] 广东省标准《高层建筑混凝土结构设计规范》反应谱计算工具 [02] [结构力学][结构设计] 两端固支梁弯矩为0点距端部的距离 [03] [抗震][结构设计][规范] 非抗震设计情况下混凝土柱的“轴压比”可达多大? [04] [结构设计][动力学] YJK中CQC振型组合地震力的复核 [05] [结构设计][楼梯] 混凝土楼梯施工图笔记 [06] [结构][设计][规范] 关于结构倾覆力矩计算公式的另一种理解 …
实干、实践、积累、思考,创新。 小伙伴拿了两份OpenSees资料问,OpenSees中瑞丽阻尼( Rayleigh Damping)的刚度系数应该如何取,搞不清楚几个刚度系数的关系。 OpenSees中Rayleigh_Damping的命令在这:https://opensees.berkeley.edu/wiki/index.php/Rayleigh_Damping_Command 具体贴图如下: 由图可见,OpenSees的瑞丽阻尼可以考虑三种刚度矩阵,并且可以考虑不同的刚度系数。 三种刚度矩阵为: current stiffness maxtrix, initial stiffness maxtrix, committed stiffness maxtrix, 前两个比较好理解,如字面意思,分别为当前刚度矩阵和初始刚度矩阵,最后一个committed字面意思为 提交完成。其实说的是上一步收敛的刚度矩阵。在OpenSees中,非线性迭代收敛后,都是通过commit函数进行收敛后相关变量的保存或者传递工作,所以带committed的变量为上一步收敛状态的变量。这个可以在源代码中看到。 用不同的刚度矩阵建立阻尼模型,分析结果肯定是有所差异的,通常差异可能不大(具体情况得具体分析,可以看这条命令底下的参考文件),但可以确定的是,使用当前刚度矩阵或者收敛刚度矩阵,计算量肯定会比使用初始刚度矩阵大一些,一般常规的设计软件,主要是用初始刚度矩阵及质量矩阵建立瑞丽阻尼,此时betakinit及alphaM可取1.0,其他刚度系数可取0。 相关话题 ( Related Topics) [01]. [Tool][软件] OSMatTest: OpenSees Material Test [OpenSees单轴材料测试工具] [02]. OpenSees …
实干、实践、积累、思考、创新。 介绍 NSDOF ( [软件][更新][Dynamics] NSDOF v2021: A Tool for Nonlinear Dynamic Analysis of SDOF System (NSDOF单自由度系统动力非线性分析工具 v2021) ) 设置摩擦阻尼器的单自由度体系动力时程分析。简单来个step by step图片演示吧。 STEP 1: 导入一个震荡动力荷载 STEP 2: 假定主体结构为弹性,设置摩擦阻尼器的摩擦力及刚度。 STEP 3: 点击Run …
实干、实践、积累、思考、创新。 CFHLA软件连接:[工具][试验][更新] CFHLA(v1.1): Cyclic/Fatigue Loading Hysteresis Loop Analysis [单向往复与疲劳滞回曲线分析工具][单向往复加载] 借用小伙伴的数据,做个CFHLA简单操作。 STEP 1: 打开软件 STEP 2: 导入滞回曲线 (数据为两列的文本格式,一列为位移,一列为力) 数据较大,为了方便测试,可以指定导入最大数据行号 STEP 3: 运行分析 STEP 4: 滚动LOOP NO. 下拉菜单,可高亮显示各个滞回环,查看软件滞回环是否划分正确。 STEP 5:点击相关按钮,输出图形结果或者文本结果。 相关博文( Related …
实干、实践、积累、思考、创新。 利用小伙伴的隔震支座竖向滞回曲线数据,做一个DataSmoothing曲线修正的案例。 STEP 1:导入滞回曲线 (可以发现,滞回曲线非常抖动,主要是出现在Y轴,PS这里力和位移有点对调了,不过不影响分析,是位移出现了抖动) STEP 2:按方法1的初始参数对曲线Y数据进行修正,由图可见,初始参数不合适这个实验数据,修正效果很差,并没有剔除掉无用的凸起。 局部放大数据可发现红色修正数据依然非常多锯齿,需要调整参数进一步修正。 STEP 3:按方法1重新调整参数,进行修正。此时红色的修正数据效果较为理想。 局部放大数据可发现红色修正数据已经得到了平滑,而且较为精确的捕捉到了原始数据的趋势。 STEP 4:查看修正效果OK后,可以将数据输出为EXCEL绘图。 DataSmoothing 的软件案例 ( Application Examples) [01] [工具][试验][编程] DataSmoothing试验滞回曲线平滑修正——案例1 [02] [工具][试验][编程] DataSmoothing试验滞回曲线平滑修正——案例2 [03] [工具][试验][编程] DataSmoothing + OutlierRemoval 试验滞回曲线修正与平滑案例 [04] …