坚持实干、坚持实践、坚持积累、坚持思考、坚持创新！ Tag: 结构工程博士 结构工程师 伪程序员 建筑结构 抗震设计 性能设计 弹塑性 编程 软件开发 有限元 超高层 超限设计
坚持实干、坚持实践、坚持积累、坚持思考，坚持创新。 在 BBCA 2020（ [科研][工具][软件] BBCA v2020: 滞回环骨架曲线提取工具 [BBCA: Extract Backbone Curve from Hysteretic Loop] ）的基础上，增加这个 BBCA_BAT软件，用于进行多个构件的批处理分析。方便做构件参数研究的朋友。 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 基本目标是：从滞回曲线提取骨架曲线，方便做试验的小伙伴. Extract Backbone Curve from Hysteretic …
坚持实干、坚持一线、坚持积累、坚持思考，坚持创新。 在 YPD 2020（ [软件][科研][试验] 2020版 YPD: Yield Point and Ductility [2020版 YPD: 等效屈服点和延性分析程序] ）的基础上，增加这个 YPD_BAT软件，用于进行多个构件的批处理分析。方便做构件参数研究的朋友。 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 软件根据构件或结构的力-位移骨架曲线计算的等效屈服变形、极限变形和延性系数。软件提供了3种常用的等效屈服位移计算方法：几何作图法、等能量法、R.Park、最远点法。方便研究而编写的小工具，特别适合于做实验研究的同学。 This program is used for evaluation …
坚持实干、坚持实践、坚持积累、坚持思考，坚持创新。 在小伙伴们的建议下，在 HLA2020（[科研][试验][工具] HLA v2020: Hysteretic Loop Analysis Program[HLA滞回环分析工具 2020版]）的基础上，增加了一个HLA_BAT，用于进行多个构件的批处理分析。 主要是给做大量构件参数化分析研究的小伙伴使用。 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 主要是针对做构件试验的研究生及科研人员做的，软件可以对对称的滞回曲线，不对称的滞回曲线进行分析，获得各圈滞回曲线的割线刚度、等效黏滞阻尼系数、耗能系数，等常用的做构件滞回性能评估的参数。程序非常方便做构件试验的同学使用，只要将试验的滞回曲线导入程序进行分析，可输出常用的分析参数，可直接输出 EXCEL图表，用户导入数据，进行分析后可选择输出EXCEL图表，就是我们连EXCEL图都懒得做了，输出EXCEL后可以对格式稍微做些调整，即可直接用于论文；另外还可以可将各圈滞回环曲线输出到文本，更加方便使用者进行数据分析。（HLA_BAT 支持批量分析勾选的曲线，输出各个曲线的详细分析结果，输出各个曲线指定参数的汇总结果，方便进行构件进行参数分析研究。） 程序使用 ( Steps to use) 程序包含以下 HLA 的功能： （1）导入滞回曲线（Import hysteretic …
坚持实干、坚持实践、坚持积累、坚持思考，坚持创新。 根据小伙伴的建议，制作了这个程序。 程序图标 ( Program Icon )： 程序介绍 ( Program Introduction) 程序主要主要作用是，剔除数据中的异常点，毛刺点，跳跃点。 数据中的这些异常点及毛刺点通常是因为采集仪器信号不稳定引起的。 比如我们做试验的时候，仪器采集的位移计数据可能就不稳定，通常就会出现这些毛刺点，而且一般仪器采集的频率比较搞，采集的数据比较密集，后期人工处理比较困难，因此可采用程序去处理。 程序界面 ( Program Interface ) 程序案例 ( Program Example ) 如下图所示，为低周往复试验的数据。其中第一幅图为位移计的位移数据，第二图为作动筒力传感器的数据，第三图为位移-力滞回曲线。 由图可见，试验力传感器的数据是比较稳定的，而位移数据则存在较多的异常点，是由于采集仪器信号不稳定引起的。进而导致绘制出的滞回曲线数异常，根本无法使用。 此外还可以发现，采集的位移计数据仅仅是一些零散的位置出现了波动，波动外的数据是稳定的，因此，只要处理好这些波动的位移数据，即可将数据恢复。 将位移数据导入OutlierRemoval 进行处理，修正这些异常数据点，处理后的位移历程数据及滞回曲线数据如下图所示。 可见异常点处理后，滞回曲线的形状出来了。 进一步对滞回曲线进行调整，基本的数据就出来了，如下图所示。(PS. 当然你可以调好看点，这里仅简单演示) 如何获得软件 ( …
坚持实干、坚持一线、坚持积累、坚持思考，坚持创新。 采用GMS选波系统（ http://www.jdcui.com/?page_id=6118 ）做的一个选波案例，一个实际工程的选波案例。 某长周期结构，前三周期分别为 8.3s、6.7s、4.3s，跨度较大，设防烈度7度0.1g，进行大震弹塑性时程分析选波。 由于第一周期大于6s，因此涉及到6s后反应谱的取值问题，本项目6s后按斜直线的方式选。 因为第一周期大于6s太多了，如果按拉平处理，结果太过保守。 以下是GMS选波过程： 以下是所选地震波的 YJK 时程分析结构 各组天然波的反应谱与规范反应谱对比结果， 相关案例 ( Related Examples) [01]. [工程][选波][地震波] 某超高层选波案例（GMS选波系统-选波应用案例1） [02]. [工程][选波][地震波] 某框筒高层建筑结构选波案例（GMS选波系统-选波应用案例2） [03]. [工程][选波][地震波] 某多层框剪建筑结构(短周期)选波案例(GMS选波系统-选波应用案例3) [04]. [工程][选波][地震波] 某钢筋混凝土框架-核心筒高层建筑结构选波案例（GMS选波系统-选波应用案例4） [05]. [工程][选波][地震波] 某大底盘-多塔-高位连体高层建筑结构选波案例（GMS选波系统-选波应用案例5） [06]. [工程][选波][地震波] 某8度区大底盘-多塔高层建筑结构选波案例（GMS选波系统-选波应用案例6） [07]. [工程][选波][地震波] 某7度区框架核心结构选波案例（GMS选波系统-选波应用案例7） …
坚持实干、坚持一线、坚持积累、坚持思考，坚持创新。 采用GMS选波系统（ http://www.jdcui.com/?page_id=6118 ）做的一个选波案例，主要用于科研。 小伙伴做大跨刚构桥的地震易损性分析，需要选取一定数量的近场地震波，但最初选的波计算结果离散性比较大，后初步分析发现是因为选的波没有很好的与反应谱匹配，导致IDA分析结果较为离散。 为此，重新选波。选波基本参数，场地地震基本烈度为Ⅵ度，设计地震动峰值加速度为 0.05g，地震动反应谱特征周期为 0.40s，结构阻尼比为 0.05。断层距控制在0~30km。 结构主要周期分布在2s~4s。由于要分析地震波的脉冲效应对场地的影响，因此，选波分两大组，一组是近场有脉冲，一组是近场无脉冲波。 如下是选波过程： 近场 不控脉冲 近场 脉冲波 如下是选取的地震波的反应谱，可见选取的地震波反应谱与目标反应谱十分吻合。 近场 不控脉冲 近场 脉冲波 如下是选取的部分脉冲波的波形，可以看见，有十分强的脉冲效应。 相关案例 ( Related Examples) [01]. [工程][选波][地震波] 某超高层选波案例（GMS选波系统-选波应用案例1） [02]. [工程][选波][地震波] 某框筒高层建筑结构选波案例（GMS选波系统-选波应用案例2） [03]. [工程][选波][地震波] 某多层框剪建筑结构(短周期)选波案例(GMS选波系统-选波应用案例3) …
坚持实干、坚持实践、坚持积累、坚持思考，坚持创新。 随后更新…… 微信公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号
坚持实干、坚持一线、坚持积累、坚持思考，坚持创新。 采用GMS选波系统（ http://www.jdcui.com/?page_id=6118 ）做的一个有趣的选波案例，小伙伴找我们帮忙做的，这个主要用于科研。 选波需求参数： （1）近场（有脉冲）(Near-fault)，远场长周期（无脉冲）(Far-fault），普通地震波（用于对比分析用尽量突出前几个波的特点）各10条。断层距要求：近场 2~20km，远场 >60km.（GMS选波系统可以实现这样的特殊参数控制要求！！） （2）同一地震波不超过2条 （3）结构基本周期2.5s左右，结构阻尼0.05。三类场地，地震分组第一组，设防烈度8度。根据我国规范《抗规》和《高规》，反应谱的特征周期Tg=0.45s，中震得影响系数最大值为0.45，时程分析的地震加速度最大值（PGA）为200cm/s2。 （4）反应谱控制要求 强调按美国规范ASCE7-10的方法进行反应谱验算，即要求在0.2T~1.5T范围内，地震波的均值反应谱与目标反应谱比较吻合，且不应该低于规范的目标反应谱。本项目目标设计反应谱：中震规范设计谱。由于结构基本周期为2.5s左右，即控制结构在0.5s~4s（0.2T~1.5T）范围内，选取的时程波的反应谱与目标反应谱比较吻合，且不应低于目标反应谱。（GMS选波系统也可以实现这样的特殊要求！！） GMS选波结果： 【近断层】【Near-fault】 【远断层】【Far-fault】 【普通断层距】【moderate fault distance】 可见，各组波反应谱的平均值与目标反应谱较吻合，且在0.2~1.5T周期范围内反应谱均值不低于目标反应谱。 相关案例 ( Related Examples) [01]. [工程][选波][地震波] 某超高层选波案例（GMS选波系统-选波应用案例1） [02]. [工程][选波][地震波] 某框筒高层建筑结构选波案例（GMS选波系统-选波应用案例2） [03]. [工程][选波][地震波] 某多层框剪建筑结构(短周期)选波案例(GMS选波系统-选波应用案例3) [04]. [工程][选波][地震波] …
坚持实干、坚持一线、坚持积累、坚持思考，坚持创新。 用GMS选波系统（ http://www.jdcui.com/?page_id=6118 ）做的一个有趣的选波案例，小伙伴找到，帮忙做的，主要用于进行IDA分析研究。 基本信息如下： 三类场地，地震分组第二组，设防烈度8度，结构阻尼0.05；一共有多组结构，结构的前三阶周期分布在0.4s~1.1s；选25组地震波，多组模型公用。 GMS选波信息： 选波按中震选，按中震反应谱进行匹配，所给出的地震波峰值均按规范调整为200cm/s2，IDA分析进行其他强度地震分析时候，可以自行对中震地震波进行缩放来使用，超大震就按超大震得PGA来调整即可。规范要求时程分析时有效持时不小于15s及5倍T1，对于本结构，有效持续时间不小于15s，所选地震波控制时间均在15s以上。 选波兼顾0.5s~1.1s范围的反应谱，及1.1s后的反应谱。0.5s~1.1s权重大，务求平均谱与规范反应谱相差不大于20%，1.1s后，尽量满足，优先选最满足的。达到各组结构均适用。这样，做IDA分析也方便不同模型之间进行对比分析。如下是本次选波的反应谱结果。可见，在结构主要周期范围内（0.5s~1.5s）内，反应谱吻合很好。 GMS选波 反应谱对： 相关案例 ( Related Examples) [01]. [工程][选波][地震波] 某超高层选波案例（GMS选波系统-选波应用案例1） [02]. [工程][选波][地震波] 某框筒高层建筑结构选波案例（GMS选波系统-选波应用案例2） [03]. [工程][选波][地震波] 某多层框剪建筑结构(短周期)选波案例(GMS选波系统-选波应用案例3) [04]. [工程][选波][地震波] 某钢筋混凝土框架-核心筒高层建筑结构选波案例（GMS选波系统-选波应用案例4） [05]. [工程][选波][地震波] 某大底盘-多塔-高位连体高层建筑结构选波案例（GMS选波系统-选波应用案例5） [06]. [工程][选波][地震波] 某8度区大底盘-多塔高层建筑结构选波案例（GMS选波系统-选波应用案例6） [07]. [工程][选波][地震波] 某7度区框架核心结构选波案例（GMS选波系统-选波应用案例7） [08]. [工程][选波][地震波] 某8度区框架-剪力墙结构选波案例（GMS选波系统-选波应用案例8） …
坚持实干、坚持一线、坚持积累、坚持思考，坚持创新。 用GMS选波系统（ http://www.jdcui.com/?page_id=6118 ）做的一个选波案例，帮助小伙伴选波进行近断层（Near-fault earthquake motions）、远断层（Far-fault earthquake motions）相关的结构抗震研究。 基本的选波参数如下： 场地类型三类，设防烈度8度，阻尼比4%，结构主要周期范围在0.5~3s，近断层控制震中距不大于20km，尽量包含脉冲效应，远断层控制震中距大于60km。地震动记录震级应大于6.5级，且地震波PGA不应小于0.1g，同一组地震事件尽量不要超过2组。按中震选波，地震波反应谱与目标反应谱在解耦股主要周期范围内误差不小于20%…..等等，选波控制，还有一些其他参数，此处不一 一列举。由于小伙伴主要选波主要用于科研，因此参数较多。GMS选波软件可以控制这些参数选波。 采用GMS按上述要求进行地震波选取，如下图： 近断层结果： 远断层选波结果： 相关案例 ( Related Examples) [01]. [工程][选波][地震波] 某超高层选波案例（GMS选波系统-选波应用案例1） [02]. [工程][选波][地震波] 某框筒高层建筑结构选波案例（GMS选波系统-选波应用案例2） [03]. [工程][选波][地震波] 某多层框剪建筑结构(短周期)选波案例(GMS选波系统-选波应用案例3) [04]. [工程][选波][地震波] 某钢筋混凝土框架-核心筒高层建筑结构选波案例（GMS选波系统-选波应用案例4） [05]. [工程][选波][地震波] 某大底盘-多塔-高位连体高层建筑结构选波案例（GMS选波系统-选波应用案例5） [06]. [工程][选波][地震波] 某8度区大底盘-多塔高层建筑结构选波案例（GMS选波系统-选波应用案例6） [07]. [工程][选波][地震波] …
坚持实干、坚持实践、坚持积累、坚持思考，坚持创新。 根据网友反馈，再次抽时间更新了这个软件。 旧版 2019版软件可以在这里查看： [软件][科研][更新][试验] 2019版 HLA: Hysteretic Loop Analysis Program [2019版 HLA: 滞回环分析软件] 目前2020版再一次做了更新，在2019版的基础上，改进了软件界面风格，更加易用，更加酷炫，增加了骨架平均曲线的输出，同时根据其他网友反馈问题做了多出小修正和改进。 程序图标 ( Program Icon )： 程序介绍 ( Program Introduction) 基本目标是：从滞回曲线提取骨架曲线，方便做试验的小伙伴. Extract Backbone Curve from Hysteretic Loop.（以前做论文时候 为处理试验数据 …
坚持实干、坚持实践、坚持积累、坚持思考，坚持创新。 在多个网友的建议下，更新了这个 HLA: Hysteretic Loop Analysis Program [HLA: 滞回环分析软件]。 2019版之前的版本：[科研][Tool][软件][试验] HLA: Hysteretic Loop Analysis Program [HLA: 滞回环分析工具] 2019版： [软件][科研][更新][试验] 2019版 HLA: Hysteretic Loop Analysis Program [2019版 HLA: 滞回环分析软件] 目前2020版再一次做了更新，在2019版的基础上，改进了软件界面风格，更加易用，更加酷炫，同时增加了正负半圈各类响应的平均值曲线及其他网友反馈的问题。 程序图标 ( …
坚持实干、坚持一线、坚持积累、坚持思考，坚持创新。 忙完项目，根据网友要求，更新一下这个YPD程序，YPD: Yield Point and Ductility [YPD: 等效屈服点和延性分析程序]。 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 软件根据构件或结构的力-位移骨架曲线计算的等效屈服变形、极限变形和延性系数。软件提供了3种常用的等效屈服位移计算方法：几何作图法、等能量法、R.Park、最远点法。方便研究而编写的小工具，特别适合于做实验研究的同学。 This program is used for evaluation of yield deformation , ultimate deformation and ductility …
坚持实干、实践、积累、思考，创新。 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 随后更新…… 欢迎关注 “结构之旅” 微信公众号
坚持实干、实践、积累、思考，创新。 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 马上更新…… 欢迎关注 “结构之旅” 微信公众号
坚持实干、实践、积累、思考，创新。 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 一个数学题，红酒杯与白酒杯的一道数学题。从结构求真群听来的。主要是这么一个问题，有一杯白酒和一杯红酒，从白酒杯里蒯一勺白酒倒入红酒杯里，然后再从红酒杯里蒯一勺倒入白酒杯里，问：白酒杯里的红酒多，还是红酒杯里的白酒多? 于是写这个小软件测试一下，结果也十分有趣，看完结果也恍然大悟。 首先，相互蒯一勺后，白酒杯里的红酒和红酒杯里的白酒是一样的多的。 而且，只要不断相互倒酒，那么白酒杯里的红酒和红酒杯里的白酒就会一直是一样多。 其实与每次倒酒的量无关，不管倒几次，不管每次倒的量是否相等，只要最终状态还是两杯一样多，那么就是一样。最终状态控制，以白酒杯为例，出去多少白酒，就会进来多少红酒。 而且，不断相互倒酒，当倒酒次数无限多厚，白酒杯里面的红酒和白酒，红酒杯里面的红酒和白酒就一样多了。 其实也很好理解，相互不断倒酒，无限多次地倒，相当于酒杯倒均匀了。 十分有意思。感兴趣的朋友可以下载来玩玩： Download: 红酒与白酒 White_RedWine.rar 欢迎关注 “结构之旅” 微信公众号
坚持实干、实践、积累、思考，创新。 网友使用CAD2XTRACT出现问题，讨论后解决。以下是网友使用该插件进行的一些复杂核心筒的建模分析案例。 程序下载 ( Program Download ) Link：CAD to XTRACT 插件[AutoCAD转任意XTRACT截面] 相关博文（Related Posts） [01] CAD to XTRACT 插件[AutoCAD转任意XTRACT截面] [02] CAD to XTRACT插件应用的一些问题 [03] XTRACT PMM 3D屈服面的绘制(Origin版) [04] XTRACT PMM 3D屈服面的绘制(MATLAB版) …
坚持实干、实践、积累、思考，创新。 用GMS选波系统（ http://www.jdcui.com/?page_id=6118 ）帮小伙伴做的一个的选波，这里简单记录一下。 选波基本信息如下。项目为7度罕遇，二类场地，第二分组，按照 铁路工程抗震设计规范 进行选波，结构周期0.25s，在平台段。 反应谱形状如下： 采用GMS系统进行选波 地震波的反应谱与规范反应谱的对比如下，可见，在结构主要周期范围（平台段）吻合较好。 相关案例 ( Related Examples) [01]. [工程][选波][地震波] 某超高层选波案例（GMS选波系统-选波应用案例1） [02]. [工程][选波][地震波] 某框筒高层建筑结构选波案例（GMS选波系统-选波应用案例2） [03]. [工程][选波][地震波] 某多层框剪建筑结构(短周期)选波案例(GMS选波系统-选波应用案例3) [04]. [工程][选波][地震波] 某钢筋混凝土框架-核心筒高层建筑结构选波案例（GMS选波系统-选波应用案例4） [05]. [工程][选波][地震波] 某大底盘-多塔-高位连体高层建筑结构选波案例（GMS选波系统-选波应用案例5） [06]. [工程][选波][地震波] 某8度区大底盘-多塔高层建筑结构选波案例（GMS选波系统-选波应用案例6） [07]. [工程][选波][地震波] 某7度区框架核心结构选波案例（GMS选波系统-选波应用案例7） [08]. [工程][选波][地震波] 某8度区框架-剪力墙结构选波案例（GMS选波系统-选波应用案例8） [09]. [工程][选波][地震波] …
坚持实干、实践、积累、思考，创新。 20200205放假期间做的，现在整理出来。整个暑假，因为疫情都待在家中。 简单的傅里叶分析例子，方波信号分解。用前面自己写的 FOUR_TRAN 傅里叶分析工具进行分析。 准备如下的方波信号数据。 打开FOUR_TRAN，导入方波数据，进行傅里叶分析。 将前2阶频率的数据累加，获得下图粉色部分曲线。曲线类似一个sin函数，与方波差异较大。 将前5阶频率的数据累加，获得下图灰色部分曲线，曲线相对有之前多，与方波接近了。 将前20阶频率的数据累加，获得下图蓝色部分曲线，发现曲线以方波为中心进行震荡，并且与方波十分接近了。 接着将前50阶频率的数据累加。 以下是，将前300阶频率的数据累加的结果，可见曲线在离散点上与方波吻合了。方波可以由这么多三角函数累加而成。 绘制成EXCEL图形如下。 关于软件( About the Program) [01] [数学][地震动][软件] FOUR_TRAN: Fourier Analysis Tool [傅里叶分析工具] 案例 ( Examples ) …
坚持实干、实践、积累、思考，创新。 滤波，傅里叶分析的一个应用。用前面自己写的 FOUR_TRAN 傅里叶分析工具进行分析。 函数 y1： y1=10*SIN(2*PI()*0.5*x)，振幅为10，频率为0.5； 函数2： y2=0.5*SIN(2*PI()*10*x)，振幅为0.5，频率为10； 函数yw 为函数y1与函数y2的叠加： y=y1+y2 从图可见，由于y2振幅小，频率大，y1振幅大，频率小。叠加后y整体趋势与y1类似，相当于在y1的基础上加上了一些波动。 我么记下来通过傅里叶变化，想办法把从y中把这些波动的y2过滤掉。 打开FOUR_TRAN，导入离散后的函数y，并进行傅里叶分析。 从上图可见，傅里叶变换，可准确解析出两个主要频率成分，一个是0.5，一个时10，其中频率为0.5的傅里叶振幅谱大，而频率为10的傅里叶振幅谱小，能量成分低。与我们前面的函数一致。 接着我们把频率为0.5附近的成分累加起来（逆变换），如下图红色部分，红色曲线基本就是最初的y1，y2被成功过滤掉了。 如下，红儿曲线即是过滤掉频率10后的曲线，曲线与y1一致。 接着我们把频率为10附近的成分累加起来（逆变换），得到下图的红色部分，此时红色部分基本是函数y2，y1成功过滤掉了。 这就是简单的滤波的例子。 关于软件( About the Program) [01] [数学][地震动][软件] FOUR_TRAN: Fourier Analysis …
坚持实干、坚持一线、坚持积累、坚持思考，坚持创新。 有空学习一下减隔震相关的知识。 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 非线性粘滞阻尼器耗能阻尼指数 lambda1，按《JGJ297-2013建筑消能减震技术规程》规范编制。 程序界面 ( Program Interface ) 公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号
坚持实干、坚持一线、坚持积累、坚持思考，坚持创新。 网友反馈，HLA软件（ HLA: Hysteretic Loop Analysis Program [HLA: 滞回环分析软件] ）、 BBCA软件（[科研][工具][软件] BBCA: 滞回环骨架曲线提取工具 [Extract Backbone Curve from Hysteretic Loop]）软件界面太大，电脑桌面无法显示完整，如下图。缩放到最大，依然无法显示完整个软件界面。 因为不断更新，软件加入了很多变量输出，同时为了一目了然，都直接显示到软件上。但这个不是根本原因，电脑是不可能显示不完整的。 解决办法： （1）将电脑屏幕分辨率设置到最大。 （2）如果分辨率设置到最大，依然无法显示完整。那么则调整一下字体的缩放即可。 对于WIN7系统，在显示设置里面，设置为【较小】-100%，也就是不进行字体缩放。貌似默认一般尺寸的笔记本都是125%缩放，因为现在电脑屏幕分辨率都比较高，但是电脑屏幕的尺寸却没有增大。 对于WIN10系统，在显示模块，缩放与布局位置，将缩放改为100%。 关于程序 [01] HLA: [软件][科研][更新][试验] 2019版 HLA: …
随后更新……，敬请期待…… 程序图标 ( Program Icon ) 程序介绍 ( Program Introduction) 程序界面 ( Program Interface ) 公众号 ( Wechat Subscription) 欢迎关注 “结构之旅” 微信公众号
[1] Ground Motion Selection (选波) 服务
[2] 书:《PERFORM-3D原理与实例》
[3] 书:《有限单元法-编程与软件应用》
[4] 我的博士论文 My Ph.D. Paper
[5] Software Notes [结构软件笔记汇总]
