[数学][编程] 蒙特卡洛模拟法求圆周率π (Monte Carlo method to find PI)

实干、实践、积累、思考,创新。 蒙特· 卡罗方法是一种通过概率来得到问题近似解的方法。其实之前也测试过,用在这个投针求PI([数学][概率] Buffon’s Needle problem [蒲丰投针问题]),最近又搞起算法研究,顺带再整体各种东西都测算一次吧。 蒙特· 卡罗方法求圆周率PI 网上很多例子。直接上图吧。   微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[下载][软件][地震工程] Spectr_Evolution: Seismic response spectrum evolution [地震波反应谱演化][反应谱随积分时间长短的变化]

实干、实践、积累、思考、创新。 写在前面: 受 徐自国 博士 的idea的启发,写了这个地震波反应谱演化工具,idea来源于徐总。主要用于分析地震波反应谱随选取的地震波计算时间长短的一个演化。可以有不同的思考。感兴趣的朋友可以下来看看。以下是软件简介。 程序图标( Program Icon ) 程序介绍 ( Program Introduction) 主要用于分析地震波反应谱随选取的地震波计算时间长短的一个演化。还有一些其他思考。后续再陆续放上来。 程序视频:     程序截图: 程序下载 ( Program Download ):  下载: Spectr_Evolution(地震波反应谱演化).rar ( PS. 程序免费下载 ,欢迎大家给我提意见 ,敬请尊重劳动成果 ) 相关话题 ( Related Topics) [01]. [Tool] …

[下载][软件][结构设计] BeamWebRebarNum: RC梁的构造腰筋计算工具 [Calculation tool of waist reinforcement for RC beam]

实干、实践、积累、思考,创新! 之前写的小工具,整理分享给大家。 程序图标 ( 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 ) (  如果需要这个软件,请在评论区留下您的评论,软件会发送到您的评论邮箱,敬请尊重劳动成果 !!)   (大家注意评论邮箱要填写正确,后台自动发送软件,如果邮箱填写错误,会发送不到。) 相关软件资料 ( …

[下载][软件] BeamRebarNum: RC梁纵向钢筋单排最大根数[Maximum Number of Beam Longitudinal Reinforcing Bars]

[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.

[软件][动力学][Dynamics] NSDOF算例6——设置黏弹性阻尼器单自由度体系动力时程分析 (Dynamic time history analysis of a single degree of freedom system with viscoelastic dampers)

实干、实践、积累、思考、创新! 粘弹性阻尼器是一种有效的减震控制装置,主要依靠粘弹性材料的滞回耗能,为结构提供附加刚度和阻尼,减小结构动力反应,从而实现减震目标。粘弹性阻尼器既能提供刚度,也能提供阻尼,其典型的滞 回曲线为椭圆形,具有良好的耗能性能。下图 所示是粘弹性阻尼器的滞回环形状示意。 粘弹性阻尼器,恢复力由两部分构成,即速度相关部分和位移相关部分,具体如下 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 …

[工具][试验][编程] DataSmoothing试验滞回曲线平滑修正——案例4 (隔震支座竖向荷载-位移滞回曲线)

实干、实践、积累、思考、创新。 利用小伙伴的隔震支座竖向滞回曲线数据,做一个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] …

[工具][软件]OSFSV v2021: OpenSees Fiber Section Viewer [OpenSees纤维截面可视化工具]

实干、实践、积累、思考、创新。 一直想更新这个软件,这款软件最早在2017年开发出来的时候,算是全网唯一的opensees的纤维截面工具,最近终于可以抽时间进一步对这个软件进行更新。 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) OpenSees纤维截面可视化工具 (OpenSees Fiber Section Viewer)。程序基本功能是,解析OpenSees的.tcl文件,提取其中的纤维截面,显示具体的纤维截面剖分情况。OpenSees Fiber Section Viewer is a visualization program for fiber secion in OpenSees. 编写这个程序的主要目的是,由于很多.tcl文件编写时,纤维截面的定义采用的是 OpenSees 提供的path和layer命令,这些命令容易编写却不方便查看,难以检查具体的纤维截面剖分是否正确。尤其是在阅读别人编写的 .tcl …

[下载][软件]GB-SPECT V2021: 中国规范反应谱生成程序 [Chinese Code’s Design Response Spectrum]

实干、实践、积累、思考、创新。 有网友问可否更新一下这个程序:[下载][软件][规范]GB-SPECT: Chinese Code’s Design Response Spectrum[中国规范反应谱生成程序],在原来程序基础上增加《建筑隔震设计标准》 GB/T 51408-2021 的反应谱。这是个很简单小程序,2013年写的,那时候刚学编程不久,主要是练习编程用。既然网友提到,趁着中秋假期,把程序翻出来更新一下,有需要的可以看看。 程序图标( Program Icon ) 程序介绍 ( Program Introduction): 该程序用于生成中国规范的设计反应谱。包括 GB 50011-2010 《建筑抗震设计规范》 和、GB 50011-2001《建筑抗震设计规范》 和 GB/T 51408-2021《建筑隔震设计标准》。程序可以输出反应谱曲线的文本文件和图片。 This program is used for …

[软件][动力学][Dynamics] NSDOF算例2——单自由度体系非线性动力时程分析 ( NSDOF Example 2 — Nonlinear dynamic time history analysis of single degree of freedom system)

实干、实践、积累、思考、创新。 最近小伙伴做非线性粘滞阻尼器的参数分析,于是我们在 NSDOF (http://www.jdcui.com/?p=13947)软件上加了非线性粘滞阻尼器的分析功能。顺便做些测算例子。 这个例子与前面[软件][动力学][Dynamics] NSDOF算例1——单自由度体系弹性动力时程分析 例子的模型基本一致,不同之处在于此例结构为非线性,取二折线非线性本构,进行非线性动力时程分析。 单自由度体系参数: 质量m: 1.0 N-s2/m (kg); 阻尼比: 0.05; 初始刚度k0: 10 N/m; 屈服强度Fy: 0.75N 相应的屈服位移为 0.075m 重力加速度g: 9.807 m/s2; 加速度时程: A9OL 对应的阻尼系数c: 0.31623 N-s/m, 单自由度体系的自振周期为  …

[试验][软件][YPD] 等效屈服点和延性计算程序 YPD 提示 “Floating point division by zero” 的处理办法 (案例2)

实干、实践、积累、思考、创新。 小伙伴采用YPD( [软件][科研][试验] 2021版 YPD: Yield Point and Ductility [2021版 YPD: 等效屈服点和延性分析程序] )分析数据,遇到 “Floating point division by zero” 错误,这里记录一下处理方法。 (1)打开软件,导入数据,图形显示数据能导入 (2)点击“运行分析”,弹出报错 “Floating point division by zero” 错误对话框。 (3)后续检查发现,主要原因是骨架曲线数据是倒序的,是从大到小排序的。YPD应该读入从小到大的骨架曲线。 (4)将骨架曲线数据改为从小到大。 (5)重新导入修改后的骨架曲线数据,并点击“运行分析”,可正确计算并获得计算结果了。 …

[笔记][非线性][OpenSees] 多个单轴二折线非线性本构并联后的滞回曲线是什么样子的?

坚持实干、坚持实践、坚持积累、坚持思考,坚持创新。 两个单轴二折线随动强化本构并联叠加后滞回曲线是什么样子的呢?在讨论别的问题的时候突然遇到。一开始以为组合后还是二折线,后面仔细想想发现不对。那还是动手测试吧,看看多个单轴二折线随动强化本构并联叠加后滞回曲线是什么样子? 首先用 OSMatTest ( [Tool][软件] OSMatTest: OpenSees Material Test [OpenSees单轴材料测试工具] )生成几个二折线滞回,如下图所示。 2个二折线本构并联: 3个二折线本构并联: 4个二折线本构并联: 由上面的结果可见:大概率上,多个二折线并联后的滞回曲线不再是二折线,除非并联的各个本构的屈服位移相同;其中2个二折线本构并联后成为3折线,3个二折线本构并联后成为4折线,4个二折线本构并联后成5折线。 材料测试相关博文( Related Topics about Material Test) [01] [科研][软件][OpenSees] OSBWTest: OpenSees BoucWen Material Test ( OpenSees …

[LaTeX][Web][笔记] 利用MathJax 在网站上显示LaTeX公式

实干、实践、积累、思考、创新。 MathJax是一个开源的JavaScript库,专门用于在浏览器中显示LaTeX、MathMl 和 AsciiMath,适用于所有浏览器。 这个是MathJax的官网:https://www.mathjax.org/ 之前在博客上的公式基本都是图片,看以后有没有可能尝试用LaTeX在博客上写公式,估计是有点难。。。。 以下是几个MathJax显示的LaTex公式,效果还是可以的。 $$E=mc^2$$ $$E=mc^2$$ $$E=mc^2$$ \(E=mc^2\) \(E=mc^2\) \(E=mc^2\) \(x = {-b \pm \sqrt{b^2-4ac} \over 2a}\) 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[选波][地震波][工程] 某大跨结构的地震波选取(GMS选波系统-选波应用案例23)

实干、实践、积累、思考、创新。 用GMS选波系统( http://www.jdcui.com/?page_id=6118 )做的一个大跨结构的选波案例。 结构设防烈度为 7度 0.10g,设计地震分组一组,场地类别为 II 类。 前三阶周期主要分布在3~4s。由于是大跨空间结构,需要进行三向地震波的选取。 采用GMS系统选波,如下: 对于这类结构而已,结构的影子周期不仅仅是前三周期,从下面三个方向的有效质量参与系数可以看到,结构在前300阶周期范围内,累计有效质量参与系数依然有较大幅度的增加。 此外,对于需要进行多点激励分析的结构,地震波需要进行基线修正。 所选地震波主方向反应谱与规范反应谱的对比情况如下图: 相关案例 ( Related Examples) [01]. [工程][选波][地震波] 某超高层选波案例(GMS选波系统-选波应用案例1) [02]. [工程][选波][地震波] 某框筒高层建筑结构选波案例(GMS选波系统-选波应用案例2) [03]. [工程][选波][地震波] 某多层框剪建筑结构(短周期)选波案例(GMS选波系统-选波应用案例3) [04]. [工程][选波][地震波] 某钢筋混凝土框架-核心筒高层建筑结构选波案例(GMS选波系统-选波应用案例4) [05]. [工程][选波][地震波] 某大底盘-多塔-高位连体高层建筑结构选波案例(GMS选波系统-选波应用案例5) [06]. [工程][选波][地震波] 某8度区大底盘-多塔高层建筑结构选波案例(GMS选波系统-选波应用案例6) [07]. [工程][选波][地震波] …

[下载][软件][规范] GD_SPECTRUM: 广东省标准《高层建筑混凝土结构设计规范》(GD DBJ 15-92-2020) 反应谱计算工具

实干、实践、积累、思考、创新。 20200608列的题目,之前测算广东规范的时候写的工具。现在广东高规出来了可以更新了。有需要的可以下载。 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 程序下载 ( Program Download ) Download Program: GD_Spectrum (2021广东高规反应谱) ( PS. 程序免费下载 ,欢迎大家给我提意见 ,敬请尊重劳动成果 ) 相关博文( Related Topics) [01]. [Tool] SPECTR – A program for …

[MATLAB][笔记] EEEP等效屈服点及延性计算

坚持实干、实践、积累、思考,创新。 小伙伴让写个EEEP(Equivalent Energy Elastic Plastic) 等效屈服点及延性分析程序,用MATLAB测试编写。改天找时间添加到 YPD([软件][科研][试验] 2020版 YPD: Yield Point and Ductility [2020版 YPD: 等效屈服点和延性分析程序])软件。 相关博文( Related Topics) [01] [软件][科研][试验] 2020版 YPD: Yield Point and Ductility [2020版 YPD: 等效屈服点和延性分析程序] …

[试验][软件][科研] LoopModifier 2021: Hysteretic Loop Modification Program[2021版: 滞回曲线修正工具]

新版软件已发布,移步这个页面:[试验][软件][科研] LoopModifier 2024: Hysteretic Loop Modification Program[滞回曲线修正工具 2024版] 实干、实践、积累、思考、创新。 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 更新这个Loop Modifier 科研小工具,之前读博士的时候做试验的时候写的小工具,一个辅助科研的小工具,该程序主要用于根据试验的规律局部处理试验数据中的错误离散点。 关于旧版软件可以在这个链接查看(v2019之前版本 http://www.jdcui.com/?p=9260 ,2019版 http://www.jdcui.com/?p=12201)。 如今根据之前网友的反馈,重写了软件,新版软件导入数据更加方便,修改数据也更加简洁。 另外网站还提供了其他几款用于修正试验数据的工具: NoiseRemoval:http://www.jdcui.com/?p=15046 该程序主要作用是,修正试验数据中的噪声,适合那些许多波动试验数据曲线。 OutlierRemoval:http://www.jdcui.com/?p=14365 该程序主要作用是,剔除数据中的异常点,毛刺点,跳跃点。数据中的这些异常点及毛刺点通常是因为采集仪器信号不稳定引起的。 SawtoothRemove: …

[抗震][减震][笔记] 黏弹性阻尼器的减震性能曲线 Response Reduction Curve of Viscoelastic (VE) System

实干、实践、积累、思考、创新。 学习减隔震知识,研究各类阻尼器的减震性能曲线。做个笔记。 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号

[编程][计算机] Windows Installer [Msiexec] 命令行参数

坚持实干、坚持实践、坚持积累、坚持思考,坚持创新。 笔记,备忘!!! Windows ® Installer. V 5.0.18362.1 msiexec /Option <Required Parameter> [Optional Parameter] 安装选项 </package | /i> <Product.msi> 安装或配置产品 /a <Product.msi> 管理安装 – 在网络上安装产品 /j<u|m> <Product.msi> [/t <Transform List>] [/g …

[工具][试验][更新] CFHLA(v1.1): Cyclic/Fatigue Loading Hysteresis Loop Analysis [单向往复与疲劳滞回曲线分析工具][单向往复加载]

软件已更新了,移步:http://www.jdcui.com/?p=17854 实干、实践、积累、思考、创新。 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 疲劳试验滞回曲线的分析,或者单拉或者单压往复滞回曲线的分析,提取各个滞回环结果,结果曲线输出Excel等。 具体结果包括:位滞回环的位移,力,面积,阻尼比,割线刚度。 程序案例 (Examples) 程序界面 ( Program Interface ) 疲劳试验滞回曲线 往复受压(受拉)试验滞回曲线 程序操作教程 (Tutorials) [01]  [工具][教程] CFHLA单向往复滞回曲线分析 使用操作 下载 (  Download ) (  如果您需要使用这个软件,请在这给网站捐助:http://www.jdcui.com/?page_id=4813,捐助不少于【200RMB】,捐助后在评论区留下您的评论,软件将通过网站管理员QQ:3014479529给您发送,敬请尊重劳动成果 !!!) 相关博文( Related Topics) …

[地震][结构] 双向地震作用效应,【先振型组合,再方向组合】及【先方向组合再方向组合】的差异?(实际案例测算)

实干、实践、积累、思考、创新。 近日许多公众号分享了王亚勇大师提出的双向地震作用算法,许多群都进行了讨论。抱着好奇心,趁热打铁,这里也对这块内容做个测算研究。振型分解反应谱法,在计算双向地震作用时,涉及振型组合及方向组合。对于双向地震作用效应,我们是先振型组合,再方向组合?还是先方向组合,再振型组合?不同的组合顺序对结果有什么影响?以下通过算例做些探索。 1.1 测算目的 (1)了解振型分解反应谱法双向地震作用计算时,【先振型组合、再方向组合】与【先方向组合,再振型组合】的差异。先振型组合、再方向组合的结果大,还是先方向组合,再振型组合的结果大。 (2)测算王大师提出的双向地震作用计算方法与目前规范的双向地震作用计算方法的差异。这里贴一下王大师提出的算法的公式。原文见: 1.2 测算说明 (1)振型组合方式主要有SRSS方法、CQC方法和ABS法等。本文测算时,不管是【先振型组合、再方向组合】还是【先方向组合,再振型组合】,振型组合统一按CQC组合,这也是规范建议的方法。 (2)对于方向组合,考虑两种情况进行测算,分别是SRSS组合及ABS组合。对于SRSS及ABS组合,次方向的效应折减均取0.85。则,王大师建议的方法即,先方向组合,后振型组合,且方向组合采用ABS组合的方法。规范方法为,先振型租后,后方向组合,且方向组合采用SRSS组合的方法。 (3)测算的效应。由于王大师的ppt对新旧方法构件层次的内力响效应给出了较多例子。这里主要测算结构的宏观效应:楼层剪力及扭矩(累积扭矩)。结构的楼层剪力及累积扭矩值也是一种效应。 (4)测算的算法 假定,Sj(x),Sj(y)分别为X向及Y向单向地震作用下振型分解反应谱法获得的结构的地震效效应。对于本文的测算,指的是结构的楼层剪力或累积扭矩。 A.对于先振型组合,后方向组合,方向组合采用SRSS组合的验算过程: Step1:进行CQC振型组合,获得的两个方向地震作用振型组合后的效应S(x)及S(y) 其中, Step2: 进行SRSS方向组合,获得的X向为主方向及Y向为主方向的地震作用效应S(EX0.85EY)及S(EY0.85EX)。其中S(EX0.85EY) = sqrt( S(x)*S(x)+0.85*0.85*S(y)*S(y)) B.对于先方向组合,后振型组合,方向组合采用SRSS组合的验算过程: Step1:先进行SRSS方向组合,X方向为主方向的地震效应为Sjmx =sqrt( Sj(x)*Sj(x)+0.85*0.85*Sj(y)*Sj(y));Y方向为主方向的地震效应为Sjmy =sqrt(0.85*0.85* Sj(x)*Sj(x)+Sj(y)*Sj(y)) Step2:分别对Sjmx 及Sjmy 进行CQC振型组合,获得的X向为主方向及Y向为主方向的地震作用效应,同样命名为S(EX0.85EY)及S(EY0.85EX),其中, …

[工具][试验][编程] DataSmoothing V2020: A Program for Data Smoothing [试验数据曲线平滑+降噪工具]

新版软件已更新,请移步:[工具][试验][软件] DataSmoothing v2024: A Program for Test Data Smoothing [试验数据曲线平滑+降噪工具 2024版] 实干、实践、积累、思考、创新。 小伙伴让做的一个小工具,用于修正试验数据中的噪声,适合那些许多波动试验数据曲线。 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 主要用于平滑试验数据,可处理滞回曲线的波动情况,同时可一定程度处理试验数据中的锯齿问题、噪声问题。 另外网站还提供了其他几款用于修正试验数据的工具: NoiseRemoval:http://www.jdcui.com/?p=15046 该程序主要作用是,修正试验数据中的噪声,适合那些许多波动试验数据曲线。 OutlierRemoval:http://www.jdcui.com/?p=14365 该程序主要作用是,剔除数据中的异常点,毛刺点,跳跃点。数据中的这些异常点及毛刺点通常是因为采集仪器信号不稳定引起的。 Loop Modifier: http://www.jdcui.com/?p=12201  该程序主要用于根据试验的规律局部处理试验数据中的错误离散点。 …

[编程][算法][优化] 编程训练:函数的CONLIN线性化

坚持实干、坚持一线、坚持积累、坚持思考,坚持创新。     2018 12 03 挖的坑,现在来更新…… 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号