实干、实践、积累、思考、创新！ Tag: 结构工程博士 结构工程师 伪程序员 结构抗震 地震工程 超限设计 软件定制 环评减振 振动控制 减隔震 施工过程模拟 小品钢结构 有限元研发 参数化设计 大震弹塑性
实干、实践、积累、思考、创新！ 随后更新…… 相关博文 ( Related Topics) [00] [数学][地震动][软件] FOUR_TRAN: Fourier Analysis Tool [傅里叶分析工具] [01] [数学][软件] FOUR_TRAN Example 1: Filtering [FOUR_TRAN傅里叶分析工具使用案例1: 滤波] [02] [数学][软件] FOUR_TRAN Example 2: Square Wave …
实干、实践、积累、思考、创新。 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) GMP ( A tool for Calculating Ground Motion Parameters for Seismic Analysis of Structures) 是一款结构抗震分析地震动强度指标/地震动参数计算工具，程序将地震波积分+反应谱分析+地震动参数分析等功能集合于一身，使用方便。GMP的主要功能包括： (1) 丰富的地震波数据格式导入支持 (2) 地震动加速度、速度、位移时程 (3) 地震加速度谱计算功能，包括伪加速度谱、速度谱、位移谱 (4) …
实干、实践、积累、思考、创新。 最近采用YJK转模型到midas Gen，发现两个软件的地震力结果对不上，且差异较大，最后发现是阻尼模型参数设置有问题，而且可能存在一些bug，在此做个记录。 对比模型： 对比模型如下，下部4层材料为混凝土，上部6层为钢结构。其中YJK采用的是2.03版，midas Gen测试的是2020版。 两个软件模型的质量和周期结果一致，周期结果对比如下图所示。 具体模型： 在此基础上，进行了两大组模型对比测算，具体如下表所示： 编号 模型 备注 1 YJK 5%振型阻尼 阻尼比统一5% 2 Gen 5%振型阻尼(不修改阻尼比) 采用单条5%阻尼反应谱曲线，不勾选修改阻尼比 3 Gen 5%振型阻尼(修改阻尼比) 采用单条5%阻尼反应谱曲线，勾选修改阻尼比 1 YJK 材料阻尼(砼5%钢2%) 采用应变能阻尼，分别定义钢材和混凝土的材料阻尼 2 Gen …
实干、实践、积累、思考、创新。 以下内容摘选自《结构地震反应分析——编程与软件应用》一 书的 第14章，关于书本更多信息可查看链接：http://www.jdcui.com/?page_id=16529 在地震作用下，能量不断输入到结构体系中，其中一部分能量以动能和可恢复弹性应变能的形式存储，另一部分能量则被结构体系的阻尼和结构构件产生的非弹性变形耗散。当结构停止震动时，体系动能和可恢复弹性应变能归零，地震输入到结构体系的能量全部被结构体系的阻尼和结构构件产生的非弹性变形耗散。结构在地震作用下的反应过程，是地震输入能量在结构体系中以多种形式不断转化和耗散的过程。本章从能量平衡方程出发，给出各类耗能的定义，在此基础上给出逐步积分法时程分析时各类能量的求解方法，并给出具体的MATLAB编程代码。 14.1 能量平衡方程 在地震动持续过程中的任意时刻，结构体系储存和耗散的总能量等于地震动输入到结构体系中的能量[1]，即 \[{E_{In}} = {E_k} + {E_s} + {E_d} + {E_p} (14.1‑1)\] 其中EIn表示地震输入的总能量，Ek表示体系的动能，Es表示结构的可恢复弹性应变能，Ed表示结构阻尼耗能，Ep表示结构非弹性耗能。其中动能Ek和弹性应变能Es是瞬时变量，阻尼耗能Ed和非弹性耗能Ep是累积的。 以下讨论上述公式中各项能量的计算公式。 14.1.1 单自由度体系能量平衡方程 水平地震作用下单自由度体系的地震动力方程： 以相对位移表示的动力方程： \[m\ddot u\left( t \right) + …
新版已推出，请移步：[软件][编程][动力学] NSDOF v2023: A Tool for Nonlinear Dynamic Analysis of SDOF System (NSDOF单自由度系统动力非线性分析工具 v2023) 实干、实践、积累、思考、创新。 趁着国庆，来个大更新，再2020版NSDOF的基础上，增加了非线性黏滞阻尼器、摩擦阻尼器、黏弹性阻尼器、金属阻尼器，适合做减震研究。 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) NSDOF 是一个基于微软的windows窗口程序，用于单自由度结构的动力非线性分析。结构可是弹性也可以是弹塑性。动力荷载可以是施加在结构基座的地震加速度，也可以是施加在结构顶部的动力荷载。程序使用逐步积分法求解增量非线性运动方程。可以输出结构的各种响应结果，包括抗力，阻尼力，参考惯性力，位移，速度，加速度，耗能，滞回曲线等。 NSDOF is a …
实干、实践、积累、思考、创新！ NMDOF是一个多自由度剪切层模型动力非线性分析工具。NMDOF软件链接：[软件][动力学][编程] NMDOF v2022: A Tool for Nonlinear Dynamic Time History Analysis of Shear-Type MDOF System (多自由度剪切层模型系统动力非线性时程分析工具 v2022) 本篇博文用NMDOF做一个单自由度非线性动力时程分析案例，采用模型是Plastic Wen模型，也就是常说的Bouc-Wen模型，并同时利用SAP2000进行对比验证。 （1）模型参数 质量【MASS：100】 材料模型【MATERIAL：Plastic Wen】 Stiffness: 80000 Yield Strength: 640000 Post …
实干、实践、积累、思考、创新。 一个遗留了好几年的问题，起因是一位网友找了条十分特殊的自己合成波测试 SPECTR 2015 (当时搞选波工具，SPECTR是仿Seismosignal的界面做的)，结果发现Newmark法计算的反应谱结果是发散的，之所以说“特殊”，是因为其他波都没出现这个问题（其他波不同方法计算都是重合的），由于一般波测试没问题，而且默认用的逐步精确法，十分稳健，当时没太在意，就去研究其他东西去了，但是问题倒是记录到了word，一直放在硬盘的某个角落，电脑换了几轮，但记录的这个word每次都舍不得删掉，虽然极其碍眼，今天又整理硬盘，又看到了这个word，这次决定把这个word删掉，一次性解决问题，于是打开以前的代码，果然看回好几年前写的代码，写得的确丑，还有点乱，但仔细检查公式倒是没问题，那问题到底出在哪？最后灵机一动，发现问题在加速度的递推公式上，当时采用了增量法的NEWMARK递推公式，各step的加速度均通过上一步的加速度加上加速度增量逐步累加得到，而可能对于这条波，恰好出现了浮点数精度丢失。于是把加速度改为平衡方程反推，果然问题解决了！！！改天更新SPECTR 2022!!! PS：浮点数精度丢失，整型数溢出是编程中很难发现的bug，隐藏最深。Floating point precision lost and integer overflow is a hard-to-find bug in programming. 相关话题 ( Related Topics) [01]. [Tool] SPECTR – A program for Response …
实干、实践、积累、思考、创新。 介绍 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 …
实干、实践、积累、思考、创新。 在网友建议下，NSDOF ( [Tool][软件][Dynamics] NSDOF v2020: A Tool for Nonlinear Dynamic Analysis of SDOF System (NSDOF单自由度系统动力非线性分析工具 v2020) ) 增加了非线性粘滞阻尼器。可以在考虑或者不考虑结构阻尼的情况下，考虑非线性粘滞阻尼器进行动力时程分析。下面算例测试NSDOF同时设置非线性粘滞阻尼器+材料非线性的动力时程分析功能，同时用SAP2000进行同样的分析，并对比验证。其中非线性阻尼器阻尼指数取0.1，整个动力方程高度非线性。 算例参数 质点质量： 1kg 体系弹性刚度：100N/m 屈服强度： 0.75N 屈服后刚度强化系数：0 结构的粘滞阻尼系数c：0N-s/m 粘滞阻尼器的阻尼系数cvd: 1.0 N-s/m …
实干、实践、积累、思考，创新。 在网友建议下，NSDOF ( [Tool][软件][Dynamics] NSDOF v2020: A Tool for Nonlinear Dynamic Analysis of SDOF System (NSDOF单自由度系统动力非线性分析工具 v2020) ) 增加了非线性粘滞阻尼器。可以在考虑或者不考虑结构阻尼的情况下，考虑非线性粘滞阻尼器进行动力时程分析。下面做两个算例测试一下NSDOF的非线性粘滞阻尼器计算功能，同时用SAP2000进行同样的分析，并对比验证。 算例1：线性粘滞阻尼 具体参数： 质点质量： 1kg 体系弹性刚度：100N/m 结构的粘滞阻尼取 ：0 粘滞阻尼器的阻尼系数c: 1.0 N-s/m 粘滞阻尼器的阻尼指数alpha：1.0 …
实干、实践、积累、思考、创新。 最近小伙伴做非线性粘滞阻尼器的参数分析，于是我们在 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, 单自由度体系的自振周期为 …
实干、实践、积累、思考、创新。 最近小伙伴做非线性粘滞阻尼器的参数分析，于是我们在 NSDOF软件 （http://www.jdcui.com/?p=13947）上加了非线性粘滞阻尼器的分析功能。顺便用NSDOF做些测算例子。这是一个单自由度体系动力分析的例子。 单自由度体系参数： 质量m: 1.0 N-s2/m (kg); 阻尼比: 0.05； 弹性刚度k: 10 N/m; 重力加速度g: 9.807 m/s2; 加速度时程： A9OL 对应的阻尼系数c: 0.31623 N-s/m, 单自由度体系的指针周期为 1.98692s 采用NSDOF进行计算，设置参数，并分析 同时采用NONLIN进行计算，并对比验证。 时程结果对比 NONLIN的结果 NSDOF的结果 滞回曲线结果对比 …
实干、积累、思考、创新。 今天和小伙伴讨论问题。突然说到绝对加速度和相对加速度。小伙伴潜意思认为绝对加速度大于相对加速度。 因为，一致激励地震动力方程分析的时候，有动力方程可知，实际计算获得是结构的相对加速度，结构的绝对加速度等于相对加速度加上地面加速度。这么一听，似乎绝对加速度比相对加速度要大。 是不是这样呢？以下用NSDOF（ [Tool][软件][Dynamics] NSDOF v2020: A Tool for Nonlinear Dynamic Analysis of SDOF System (NSDOF单自由度系统动力非线性分析工具 v2020) ）做几个线性单自由度系统的时程分析案例。 例子1： 例子2: 从上面两个例子看，结构的相对加速度和地面加速度不总是同向的，绝对加速度可以比相对加速度大，也可以比相对加速度小，与结构的刚度、阻尼等参数有关。 其实我们可以这么想：当结构无限刚的时候，结构相对加速度为0，绝度加速度等于地面加速度，当结构无限柔的时候，结构的相对加速度与地面加速度反向，绝对加速度等于0。因此，当结构刚度K介于0~∞之间时，绝对加速度可能大于相对加速度也可能小于相对加速度。如下图所示： 相关博文( Related Posts ) [01]. [Tool] SPECTR – …
实干、实践、积累、思考、创新。 《建筑抗震设计规范》（GB50011-2010）中给出了采用振型分解反应谱法计算地震作用时的地震力计算公式：\({F_{ji}} = {\alpha _j}{\gamma _j}{X_{ji}}{G_i}\)，其中\({\gamma _j} = \frac{{\sum\limits_{i = 1}^n {{X_{ji}}{G_i}} }}{{\sum\limits_{i = 1}^n {X_{ji}^2{G_i}} }}\)，\({F_{ji}}\)为j振型i质点的水平地震作用标准值；\({\alpha _j}\)为相应于j振型自振周期的地震影响系数；\({X_{ji}}\)为j振型i质点的水平相对位移；\({\gamma _j}\)为振型的参与系数。以下根据结构动力学的相关理论，给出上述公式的一种推导。 1多自由度体系振型分解法 Mode Superposition Method 对于多质点体系，地震动力方程为： $${\left[ M \right]\left\{ {\ddot u} \right\} …
实干、实践、积累、思考、创新。 最近研究舒适度，做些算例测算。对两端铰接、两端固接、一端固接一端铰接、悬臂等截面梁进行振型分析，获得各类梁的前三阶振型，并对振型向量进行最大位移值归一化，并利用归一化后的振型向量求解前3阶振型的振型质量。测试算例梁截面统一为，梁截面为200X200，梁长度为1000mm，沿梁长划分80个单元，振型的质量通过公式 \({M_n} = \int_0^L {m(x)\phi _n^2(x)dx} \) 进行计算。 1 简支梁 1.1 振型形状 一阶振型 二阶振型 三阶振型 1.2 振型质量 振型 节点质量 总质量 振型质量 振型质量/总质量 1 0.0009815 0.314 0.1570 0.500 2 0.0009815 …
坚持实干、坚持积累、坚持思考，坚持创新。 题目如题，结论肯定是没有影响的，因为振型向量本来就是不定的，振型元素之间只有相对关系，要求解振型向量元素的具体值，必须对振型向量进行标准化。简单说即先假定某个元素的值，然后才能求解出其余元素的值。 最近在研究舒适度，顺便把相关东西整理一下，正好还有小伙伴问，同时正好测试一下在网站上用LATEX写公式，看看是不是会专业点。 基本公式 结构的运动方程： \[[M]\{ \ddot u\} + [C]\{ \dot u\} + [K]\{ u\} = \{ P\} （公式1） \] 将位移向量\(\{ u\} \) 用振型展开， \[\{ u\} = [\phi ]\{ q\} …
实干、实践、积累、思考、创新。 最近研究舒适度，涉及振型向量的标准化，顺便测试一下SAP2000默认的振型向量标准化方法。 大家都知道，振型向量是不定的，振型向量的参数之间只有比值关系。为了求解振型向量的元素绝对值，必须对振型向量进行标准化。 我们接下来测试SAP2000中振型向量的标准化方法，在SAP2000中建立一根简支梁模型，如下图： 梁的前三阶Z向振型形状如下： 振型形状是与理论分析结果一样的。 将软件输出的振型变形的平方乘以节点质量，可获得各振型的广义质量，结果均为1。即SAP2000默认输出的振型是满足关于质量矩阵内积为1的条件的。即采用的是关于质量矩阵的正交归一化方法。 相关博文( Related Posts ) [01] [Structural Dynamics][Mode superposition] 振型参与质量系数(Participating Mass Ratio) [02] [动力学][振型分解][Mode Superposition] 振型向量与振型参与系数的乘积公式推导 [03] [结构设计][地震作用][规范] 振型分解反应谱法的一些概念总结 (Basic Concepts of Response Spectra …
实干、实践、积累、思考、创新。 如题，一看是一个莫名其妙的想法。实际上也是一个错误的想法，不过，当时脑子一热，就测算一下。具体看看结果。 模型： 20层的剪切层模型。 阻尼矩阵： （1）模态阻尼，20阶振型计算 （2）模态阻尼，20阶振型计算，把非对角线元素取为0。 采用MATLAB编程，采用Newmark-β积分法进行弹性时程分析，两种阻尼模型的计算结果对比如下。 其中，参考阻尼为完整的模态阻尼，对比阻尼为去掉阻尼矩阵对角线元素后的矩阵。 由图可见，采用仅保留对角线元素的阻尼矩阵，结构的位移及剪力响应远小于完整的阻尼矩阵，而楼层加速度响应似乎相差不大！！十分诧异！！ 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号
实干、实践、积累、思考、创新。 来自小伙伴 刘骥 的分享。这个是继《 [抗震][动力学] 对于整体结构，X向地震作用下有Y向剪力吗？有！！ 》后对对称结构进行的简单测试。直接看测试结果吧。 算例1： 算例2： 可见，剪力作为矢量，是满足平行四边形法则的。同时，对于对称结构，往哪个方向输入地震加速度，总剪力均一致。 相关博文 ( Related Topics) [01] [地震][动力学] 对称结构的地震剪力规律 [02] [地震][结构] 双向地震作用效应,【先振型组合,再方向组合】及【先方向组合再方向组合】的差异？(实际案例测算) [03] [地震计算][反应谱][动力学][CQC] 振型叠加法随着组合振型数量的增加各种响应量是怎么变化的？ [04] [抗震][动力学] 对于整体结构，X向地震作用下有Y向剪力吗？有！！ 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号
实干、实践、积累、思考、创新。 对于整体结构，X向地震作用下，结构有Y向的剪力吗？以下通过两个简单的时程分析案例进行测算。 案例1： 振型结果如下： 振型1为Y向平动 振型2为X向平动 振型3为绕Z轴扭转 案例2：案例1模型逆时针旋转45度 周期值与模型1是一样的，只是因为结构转了一个角度，振型方向不同了。 振型1为135度方向平动 振型2为45度方向平动 振型3为绕Z轴的扭转 分别对两个模型沿X向施加地震动加速度时程，进行直接积分动力时程分析。所选的地震波如下图所示。 模型1沿X方向与Y方向的基底剪力结果如下图所示。由图可知，对于模型1，沿X向输入地震，Y向剪力几乎为0。 模型2沿X方向与Y方向的基底剪力结果如下图所示。由图可知，对于模型2，沿X方向输入地震，结果Y向会产生剪力，且剪力大小不可忽略。 粗看似乎有点难理解，外力和内力不是应该平衡的吗？为何施加X方向加速度，结构有Y向的剪力？对于静力情况下，结构受到到的外力与的内力平衡，比如，当沿结构X方向施加力F时，结构总的剪力必然是沿X方向，且大小为F，Y方向不存在剪力。为何到了动力情况，就不满足这个规律了？不妨看一下两种情况下结构的平衡方程。 结构静力平衡方程： 其中，为结构的外力，为结构的抗力，其中，当只有X向力作用时，即 ，，即结构的抗力也只有X反向的力，y及z向的力为0 结构动力平衡方程（地震）： 其中，当只有x向地震时，，即向量在非X向自由度上为0，此时结构的抗力，假设忽略阻尼，结构的抗力为，对比静力平衡下的公式（）可见，尽管向量在非X向自由度上为0，但是抗力的右边项不是，而是，即所谓的绝对加速度，其中相对加速度在非X向自由度上不一定为0，当在非X向自由度上存在非0值时，抗力就可以能存在非X向自由度上的力，即对于整体结构，X向地震作用下，结构整体在其他方向也可能存在抗力，包括Y向的剪力。 对于算例1，在X向地震加速度作用下，由于结构基本只有X向的位移，因此Y向的抗力很小，进而Y向剪力很小。对于案例2，由于结构扭转了45度，在X向地震加速度作用下，结构不仅有X向的位移，也有Y向的位移，有Y向的位移，就可能有Y向的抗力及Y向剪力。 由以上分析也可发现，引起动力与静力概念上不同的错觉的原因是，把动力情况下结构的抗力当成了，实际上动力情况下结构的抗力等于（忽略阻尼情况下）。 平衡还是满足的！！！。 相关博文 ( Related Topics) [01] …
实干、实践、积累、思考、创新。 7月拖到现在11月，终于可以更新了．．．．．． 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 基于FEMA 440等效线性化法 Pushover分析方法 的性能点求解程序。( A program for the Solution of Pushover Performance Point based on FEMA 440 Equivalent …
实干、实践、积累、思考、创新。 测算10组模型，模型1为2*3网格的5层框架结构，第一周期为X向平动，第二周期为Y向平动，第三周期为扭转。 模型2~模型10，由模型1沿X向不断拉长形成，层数不变。 看看往两端加宽结构的各周期如何变化。 采用YJK MultiModel Compare提取各模型的结果。 绘制各模型的周期变化曲线。 第一周期 X向平动 第二周期 Y向平动 第三周期 扭转 可见，往X向不断加宽，第一周期不断减小，第二周期不断增加，扭转周期不断增加。 为什么会这样？其实 周期的增加或者减小其实与质量、刚度增加的相对比例有关。 这里给出了基本理论解释：[动力学][Structure Dynamics] 线性增加刚度K与质量M下单自由度(SDOF)结构的周期变化 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号
坚持实干、坚持一线、坚持积累、坚持思考，坚持创新。 单自由度体系，线性增加刚度K与质量M，结构周期的变化规律。 现性递增k与m 单自由度结构的周期变化关系：周期可增可减，与k、m各自增幅有关、即 k/m有关。增、减构件类似在原有基础上递增k、m，如果原先K、M的基数已经很大，曲线已经在平滑段，简单线性增减对周期结果影响小。 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号
实干、实践、积累、思考、创新。 采用 SPECTR反应谱分析软件 （下载链接： http://www.jdcui.com/?p=1875 ） 对几个地震波进行位移谱的求解，结果如图1~图5所示，可以发现，随着阻尼比的增大，大部分地震波在各周期范围内反应谱值减小。但是部分地震波在某些周期范围内，阻尼比增大，但是位移谱值不一定减小。 在图形上表现为不同阻尼比的反应谱曲线在某些周期范围内出现相交的现象（见图1及图3）。同样，拟加速度反应谱也存在这个现象，因为拟加速度反应谱是由位移谱转换过来的。 初看这个现象感觉很奇怪，其实仔细想也十分正常。 以图1的反应谱为例，在周期2.5s左右出现明显的相交现象。把2.5s周期各阻尼比单自由度结构的位移时程绘制出来，见图6. 由图6可见，各阻尼比下，位移时程曲线的整体趋势比较一致，该凸的地方大家一起凸起来，该凹的地方一起凹下去，即趋势是一致，但是随着阻尼比的增大，出现绝对值最大值的时间不同了。 如图，阻尼比为20%的位移最大值出现在30s左右，其他阻尼比下位移绝对值最大值出现在40s左右。且阻尼比增大到20%后绝对值比其他阻尼比在40s左右出现的绝对值大。因此不同阻尼比的反应谱曲线就出现了交点。 因此，不同阻尼比的反应谱曲线可能出现交点。因为，反应谱纵坐标是绝对值。阻尼比发生改变，可能整个响应时程的整体趋势没改变，但是最大值出现的位置会不同，大小的增大或减小规律也不同，而反应谱记录的是绝对值。 所以，千万别闭着眼睛说，阻尼比越大，位移越小。瞎说！！！ 
图 1 图 2 图3 图 4 图 5 图 6 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号
已经更新至 V2022版，连接：http://www.jdcui.com/?p=16697 实干、实践、积累、思考、创新。 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) NSDOF 是一个基于微软的windows窗口程序，用于单自由度结构的动力非线性分析。结构可是弹性也可以是弹塑性。动力荷载可以是施加在结构基座的地震加速度，也可以是施加在结构顶部的动力荷载。程序使用逐步积分法求解增量非线性运动方程。可以输出结构的各种响应结果，包括抗力，阻尼力，参考惯性力，位移，速度，加速度，耗能，滞回曲线等。 NSDOF is a Microsoft Windows based application for the dynamic analysis of single degree of freedom structural systems. …
[01] Ground Motion Selection (选波) 服务
[02] 著:《PERFORM-3D原理与实例》
[03] 著:《有限单元法-编程与软件应用》
[04] 著:《结构地震反应分析-编程与软件应用》
[05] 著:《有限单元法 Python编程》
[06] 著:《结构地震动力响应Python编程》（New!!!）
[07] 土木工程试验数据处理软件汇总
[08] 自编程序 [Software Box]
