2D Truss by Matlab [Matlab 2D桁架有限元]
MATLAB 2D Truss 用MATLAB实现2D桁架
实干、实践、积累、思考、创新! Tag: 结构工程博士 结构工程师 伪程序员 结构抗震 地震工程 超限设计 软件定制 环评减振 振动控制 减隔震 施工过程模拟 小品钢结构 有限元研发 参数化设计 大震弹塑性
MATLAB 2D Truss 用MATLAB实现2D桁架
实干、实践、积累、思考、创新。 分享一个小程序。 A small program to share with you. 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 平时做构件抗震性能分析时经常需要进行构件的低周往复试验模拟。对于一般的低周往复试验,文献仅给出加载过程的分级点,如力分级点和位移分级点。而做数值模拟我们经常根据这些分级点生成一定间距的加载点进行低往复模拟,该程序可以将分级点按一定的间距生成加载历程,用于数值模拟。 Usually, the low-cycle reciprocating test simulation of the component is often required when analyzing the …
这个软件更好的替代软件为HLA,HLA连接:http://www.jdcui.com/?p=22143 实干、实践、积累、思考、创新。 分享一个小程序。 A small program to share with you. 程序图标 ( Program Icon ) 程序作者 ( Author) JiDong Cui (崔济东) 1, XueLong Shen (沈雪龙)2 1. 华南理工大学 土木与交通学院; 2. 华南理工大学 建筑设计研究院。 程序介绍 ( Program Introduction) …
分享一个小程序。 A small program to share with you. 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 地震反应谱是地震工程中重要的概念,我们经常需要将地震加速度时程转换为反应谱。SPECTR 是一个简单易用的反应谱计算程序。 Seismic Response Spectra is an important concept in earthquake engineering, we often need to evaluate seismic response …
反应谱的基本概念(Introduction to Response Spectra) 地震动反应谱:单自由度弹性系统对于某个实际地震加速度的最大反应(可以是加速度、速度和位移)和体系的自振特征(自振周期或频率和阻尼比)之间的函数关系。工程中常用到反应谱的概念,如,我们目前采用的是基于反应谱的抗震设计方法、在静力弹塑性分析的时候我们还需要用到反应谱、在动力时程分析选择地震波的时候,我们还常需要以规范的设计反应谱为参考,要求选择的地震波的反应谱和设计反应谱相比不能偏差过大。 反应谱的理论公式(Basic Equations of Response Spectra) 单自由度地震作用下的动力方程: 从公式(2-5)可知,结构的地震动响应只和 3 个变量有关:地面加速度()、结构的阻尼比(ξ)和无阻尼结构体系的圆频率(ω)。因此对于一个给定的地震动,结构的位移响应ur 是圆频率(ω)和阻尼比(ξ)的函数。根据反应谱的概念,对于特定的地震波,只要选定阻尼比ξ,结构的最大响应(包括,速度、位移和加速度)是结构频率ω的函数,对一条特定的地震波,圆频率(ω)取不同的值进行动力时程分析,将获得响应最大值作为纵坐标,圆频率(ω)为横坐标绘图,就得到反应谱曲线。 常用的地震反应谱(Commonly used Response Spectra) (1)反应谱(Response Spectra) 当最大响应分别取相对位移、相对速度、绝对加速度时,我们可获得以下反应谱: 相对位移反应谱: 相对速度反应谱: 绝对加速度反应谱: (2)伪反应谱(Pseudo-Response Spectra) 另外,工程中还常用到伪反应谱。所谓伪反应谱是指,当我们获得相对位移位移反应谱后,通过以下公式获得伪相对速度反应谱、伪绝对加速度反应谱: 伪相对速度反应谱: 伪绝对加速度反应谱: …
做研究做分析我们经常需要跟各种数据打交道,为了清晰的表达数据之间的关系,我们往往需要用到各种图,如散点图,折线图等等。首先,图是表达各种变量之间关系的最直接的方法,一个好的图往往使人理解问题事半功倍。因此,我经常觉得,一篇好的文章,如果有图,那么它的图肯定是很清晰的。因此,你不妨可以通过文章中的图的质量来判断文章的质量,如果一篇文章的图做得很丑,表达又不清晰,通常这不是一篇好文章 🙂 🙂 ,当然事情也不绝对! 因为图的重要,所以我们尝尝需要各种绘图软件,如Excel,Origin,Sigmaplot,Matlab等 除此之外,图还有另外一个功能:寻找变量之间的关系。由于图可以表达变量之间的关系,因此,我们常常通过变量之间的图来寻找变量之间的关系,从图中获取变量信息是一个十分重要的学问。对于从图获得信息这个功能,上面的一般软件都可以做,但是往往不方便,尤其是当变量特别多、同时需要不断变换各种坐标轴的变量、需要利用各种信息筛选数据并进行绘图时,上述提到的软件就特别不方便。最近处理数据,为此,我写了个简单的小软件,并且命名为: Multiparameter Plot,软件基本功能是:导入数据,并快速选取坐标轴变量,绘制2D 和 3D 坐标图,同时能够根据不同的变量对数据进行分类,并且每一个变量可以设置多个区间,按不同区间交叉绘图等。 Multiparameter Plot 基本满足了我的要求,节省了大量的时间,软件适合去做前期参数分析、观察和找规律,然后可以导出数据,用前面提到的那些软件绘制漂亮的图形。 软件界面: 2D绘图分析: 3D绘图分析: Go on!
整理了Saatcioglu & Razvi混凝土应力应变模型的参数计算方法,并编制Excel表格,希望对大家有用。Saatcioglu & Razvi混凝土本构模型可以考虑圆形截面、矩形截面、长方形截面,可考虑多种配箍形式,并可通过等效侧向力叠加的方式考虑截面中存在多种配箍形式的情况。
钱稼茹混凝土应力-应变模型是由钱稼茹等人[1]提出的,通过配置普通箍、拉筋复合箍和井字复合箍RC柱的中心受压试验,给出了配箍特征值λv在我国抗震设计规范[2]规定范围内的约束混凝土的应力应变关系。整理了钱稼茹混凝土应力应变模型的参数计算方法,并编制了Excel表格,希望对大家有用。如有错误或遗漏,欢迎大家批评指正。
软件图标 (Program ICON) 引言(Introduction) 上一篇博文《Mander Confined Concrete Model 资料整理》总结了Mander模型参数的计算方法,其中该模型十分关键的一个参数为约束混凝土的峰值强度fcc’。关于fcc’的确定,Mander模型原文(Ref. 1)给出的方法为:首先确定两个方向的约束系数,然后通过查图得到约束强度系数(Confined Strength Ratio),进而得到fcc’(博文《Mander Confined Concrete Model – Confined Strength Ratio [Mander混凝土本构约束强化系数]》也提到),不方便使用。 原文给出的约束强度系数(Confined Strength Ratio – fcc’/fco’)图: 为此,编制小程序计算Mander模型的约束强度系数,程序如下图所示。 程序功能(Program Features): 程序通过输入两个方向的有效约束系数,计算Mander混凝土模型的约束强度系数,省去查图的过程,方便平时做数值模拟。 This …
箍筋的配置对核心混凝土具有侧向约束作用,可以显著提高混凝土的强度和延性。工程实践中,常采用约束混凝土本构模型来考虑箍筋对核心混凝土的约束作用。目前已提出多种模型,其中Mander混凝土本构模型是较为常用的约束混凝土模型之一,该模型由Mander、Priestley和Park于1988[1]年提出。将Mander约束混凝土模型(Confined Concrete Model)的参数取值方法整理成文档,给出实例,并编制Excel表格,方便使用。
元包自动机(Cellular Automata)这里就不介绍了,可以在我的其他博文中了解到,有兴趣的朋友也可以去百度一切。这里利用元包自动机进行森林火灾蔓延模拟,当然这肯定不是基于物理的模拟了。最简单的基于元包自动机的森林火灾模型(Cellular Automata Forest Fire Model)大家可以在这个网站 Cellular Automata Forest Fire Model 了解到,作者用Java写了小程序进行模拟。 森林火灾模型介绍(The rules of Cellular Automata Forest Fire Model): 燃烧的树(红色格子)会在下一步变成空地(白色格子)。A burning tree (red cell) becomes an empty site (white cell)。 如果一颗非燃烧的树(绿色格子)周围有树燃烧,那么它会在下一步变成燃烧的树(白色格子)。A tree (green cell) becomes …
Kent & Park和修正Kent & Park 混凝土模型 ( Kent & Park vs Modified Kent & Park Concrete Model ) Kent & Park混凝土应力应变模型由Kent和Park[1]于1971年提出,该模型既可以考虑约束混凝土又可以考虑非约束混凝土。但对于约束混凝土,该模型仅考虑箍筋对混凝土的延性提高作用,并不考虑箍筋对混凝土强度的提高作用。因此,Scott、Park和Priestley[2]于1982发表的文章中又对改模型进行了修正,通过引入强化系数K,同时考虑了箍筋对混凝土的延性和强度的提高作用,这一模型被称为修正的Kent & Park模型。由于Kent & Park类模型的简单而不失精度的特性,使得模型应用十分广泛。OpenSees中的Concrete 01和Concrete 02单轴混凝土材料受压部分的骨架曲线就是基于Kent & Park类模型。 Kent …
趁着闲暇的时间,将以前学过的一些知识整理成文档,一方面可以巩固自己的知识,另一方面也希望这些资料能够帮助别人。今天这篇文档主要介绍 PERFORM-3D中低周往复试验模拟的两种方法(A discussion of two methods of conducting low-frequency cyclic test modeling in PERFORM-3D)。文档由我和我的小伙伴——沈雪龙(XueLong Shen,华南理工大学建筑设计研究院)一起整理。 摘要(Abstract) PERFORM-3D中不存在具体针对构件低周往复位移加载的工况,然而我们在试验或者参数校准的过程中通常需要进行低周往复试验的数值模拟并与试验结果进行对比。本文介绍了在PEROFRM-3D中可以实现构件的低周往复试验模拟的两种方法:一种是基于Push-Over工况的方法,一种是基于Dynamic Force工况的方法。并以PEER柱子数据库(http://nisee.berkeley.edu/spd/)中的一根柱子(SPD-28)的低周期往复试验为参考,采用两种方法在PERFORM-3D中实现该柱子的低周往复试验模拟,并对两种分析方法的优缺点进行探讨。 截图(Screenshots) 希望后续可以有闲暇时间陆续整理了一些学习资料上来与大家分享。 如有错误,希望大家批评指正。 😎 😎 我的邮箱: jidong_cui@163.com
最近在微信上有一个很火的小游戏——Core Ball,又名见缝插针游戏。游戏有很多个不同的版本,不同人编写。游戏玩法很加单,就是用一个轮盘不断的转,轮盘上插满了球,然后点击鼠标,发射球(针头)到轮盘上,但是球不能和轮盘上已有的球发生碰撞,碰撞则重新玩该关卡。游戏难点是:寻找合适的时机发射球,考验你的耐性,考验你是否眼疾手快。趁着休息时间,我也来编写一个:CJD Core Ball。大家快来玩。
生命游戏其实是一个零玩家游戏(百度百科),这个小游戏是一个简单的二维元包自动机(Cellular Automaton),1970年由英国的数学家John Horton Conway提出,关于生命游戏的介绍参考Wikipedia的Conway’s Game of Life。
基于1D元包自动机实现Sierpinski三角形。Sierpinski Triangle by 1D Cellular Automata。
著名的几何概率问题 —— 蒲丰投针问题(Buffon’s Needle problem ),最初由数学家Georges-Louis Leclerc, Comte de Buffon于18世纪提出。问题可表述为:假定长度为L的的针,随机投到画满间距为T的平行线的纸上,求针和平行线相交的概率。同时有趣的是,该概率值和圆周率(PI)有关系,因此,我们可以利用投针试验来计算圆周率(PI)的值。实际上,这个试验有价值的地方在于:我们可以设计一个试验,它的概率与我们感兴趣的一个变量有关,然后通过大量试验来估算这个量,这其实就是目前使用得十分广泛的蒙塔卡罗法(Monte Carlo method)的主要思想。
终于可以稍微静下来几天了,前一排一直都在实验室做试验,短短半个来月连续做了9个柱子的低周往复试验,实在是累趴了。做试验真心是一个十分辛苦的事情,因为绝对不是简单地浇筑完试件就去推完那么简单,期间遇到的各种奇奇怪怪的问题都需要自己想办法解决。作为师兄的我,自然压力很大,然而庆幸的是“构件组”的师弟妹们十分团结和互助,大家分工合作,一步步将问题解决。从梁试验到现在柱子试验,我们遇到了大大小小的各种困难,同时我们也努力去克服这些困难,为这份难得的经历而感动 。 😎 😎 。其中,龚涣钧师弟认真负责的态度十分值得肯定 😛 😎 ,虽然这些试验是他论文的一部分。 接下来还有两个未做完的柱子和一批剪力墙试件需要进行试验,困难还会有很多,相信经过大家的团结努力,苦难一定会被克服。 其实,生活就是这样子,很多时候事情不管你计划得如何周密,它也不一定会按照你的计划进行。我们需要做的是积极去面对它、积极地去解决其中的问题、真切地去经历和感受这个过程,这样你得生活才有厚度。因为生活不是一个点,生活是一个过程。 由于师弟(LEE)毕业论文的需要,这几天趁着闲暇的时间,帮师弟写了个关于ABAQUS的数据处理软件。这一次,决定尝试点新意,做了个酷炫的界面,还是很激动的。(PS. 每次打开这个软件,我都感觉自己不是在做论文,而是在打游戏,在干大project!!) 生活很需要变化,一层不变总是最容易的,但是一层不变的东西也是最无趣的;新东西是最有趣的,然而很多时候我们却因为自己的懒惰将新东西杜绝了。我们要勇于而且敢于尝试新东西,杜绝懒惰,因为只有新东西和新思想才能使人类不断进步和不断发展。 PS. 能够有能力去帮助到别人是一件很开心的事情。继续加油! 😎
分享一个小程序。 A small program to share with you. 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 针对PERFORM-3D软件的工字形梁纤维截面剖分小工具。程序通过导入文本参数(.csv),直接生成纤维截面的参数,并导出PERFORM-3D需要的二进制文件(.PF3CMP)。然后通过PERFORM-3D导入.PF3CMP文件完成繁琐的纤维截面输入工作,节省你的时间。 This program is used for the data input of the “Inelastic Fiber Beam Section” in PERFORM-3D. Through the import …
分享一个小程序。 A small program to share with you. 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 针对PERFORM-3D软件的圆形柱纤维截面剖分小工具。程序通过导入文本参数(.csv),直接生成纤维截面的参数,并导出PERFORM-3D需要的二进制文件(.PF3CMP)。然后通过PERFORM-3D导入.PF3CMP文件完成繁琐的纤维截面输入工作,节省你的时间。 This program is used for the data input of the “Inelastic Fiber Column Section” in PERFORM-3D. Through the import of …
针对PERFORM-3D软件的箱形柱纤维截面剖分小工具。程序通过导入文本参数(.csv),直接生成纤维截面的参数,并导出PERFORM-3D需要的二进制文件(.PF3CMP)。Through the import of section properties in text format (.csv), the program can generate the fiber section properties automatically. And then, the program export the binary data input file (.PF3CMP) for “Inelastic Fiber Column Section” which PERFORM-3D can then read to complete the fiber section definition.
程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 针对PERFORM-3D软件的RC剪力墙纤维截面剖分小工具。程序通过导入文本参数(.csv),直接生成纤维截面的参数,并导出PERFORM-3D需要的二进制文件(.PF3CMP)。然后通过PERFORM-3D导入.PF3CMP文件完成繁琐的纤维截面输入工作,节省你的时间。 This program is used for the data input of the “Inelastic Fiber Shear Wall Section” in PERFORM-3D. Through the import of section …
程序图标 ( Program Icon) 程序介绍 ( Program Introduction) 针对PERFORM-3D软件的RC矩形柱纤维截面剖分小工具。程序通过导入文本参数(.csv),直接生成纤维截面的参数,并导出PERFORM-3D需要的二进制文件(.PF3CMP)。然后通过PERFORM-3D导入.PF3CMP文件完成繁琐的纤维截面输入工作,节省你的时间。 This program is used for the data input of the “Inelastic Fiber Column Section” in PERFORM-3D. Through the import of section properties in text format (.csv), …
程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 针对PERFORM-3D软件的RC矩形梁纤维截面剖分小工具。程序通过导入文本参数(.csv),直接生成纤维截面的参数,并导出PERFORM-3D需要的二进制文件(.PF3CMP)。然后通过PERFORM-3D导入.PF3CMP文件完成繁琐的纤维截面输入工作,节省你的时间。 This program is used for the data input of the “Inelastic Fiber Beam Section” in PERFORM-3D. Through the import of section properties in text …