清溪长河

对Midas/Civil做非线性分析的几点看法

长河 — Sun, 29 Apr 2012 20:13:15 +0800

前几日一直为midas做某桥的非线性分析而纠结，尝试了许多种方法，但始终得不到满意的阶段/步骤时程图形。阅读了一些文献、教程之后，虽然对midas做非线性分析的过程有所了解，但问题依然存在，只是略有改善而已。把一些想法记录在此，日后当个小问题研究一下。

下面这篇文档讲midas的非线性分析问题讲述得比较清楚了：

文档地址

对于非线性分析，我在一篇文献中看到的方法不错，先通过特征值分析得到结构的一阶模态特征值，然后再讲这些特征值作为放大倍数，将该工况的相关荷载放大、设置非线性分析选项，进行非线性分析。

但midas在非线性分析的时候，有几个问题：

对于梁单元，不能同时考虑几何非线性和材料非线性；
非线性分析中不能考虑温度荷载；
对于同一个结构，非线性分析的结果存在一定的离散性。结果一旦离散，就不好确定计算的准确性了。

等过些天有空的时候，做一个ansys和midas的对比，看看到底差在哪儿了。

DWG TrueView 2013 免费下载

长河 — Tue, 24 Apr 2012 19:29:27 +0800

刚刚从国外的网站上得到消息，AutoCAD DWG TrueView 2013已经可以免费下载了。

这个软件是AutoDesk公司出品，与CAD使用相同的内核，是官方出品的DWG文件浏览器。当然，如果仅仅是个文件浏览器的话，也就不怎么值得推荐了。它还有其他几个比较吸引人的功能：

1.文件版本切换，当把CAD2008存的文件转换成CAD2004能读取的格式，TrueView可以作为选择之一；

2.格式转换，将DWG文件转换成图片或者pdf文件，这或许能比“打印成pdf”效果更好一些。

更多的功能，参考这里。

软件下载地址，在这里http://www.autodesk.com/dwgtrueview

顺便问一下，哪位朋友有VPN或者在境外可以下载到那个免费软件的，可否发送到我邮箱dachanghe@163.com呢？先感谢了。

多参数敏感性分析的意义

长河 — Sat, 21 Apr 2012 06:29:07 +0800

结构参数化敏感性分析，往往是我们在尝试对结构计算模型调整之前需要进行的一个重要准备——毕竟如果先知道可以调整的几个参数中，哪个参数对目标结果的影响比较大，就会将其作为主要调整的对象，起到提高效率的作用。

但是，常规的敏感性分析方法，都是设定一个参数按照10%左右的幅度变化，其他各参数保持固定不变的前提。这样做是为了能够查看出单一变量对结构相应的影响，但却不符合基本事实——因为实际上所有的影响因素都是同时在变动的。当常规用有限元软件来操作的计算方法，又不能让各变量一起动，那结果就失去了意义。所以就只能舍弃一点精度，保证结果的可用性。

其实，完全可以让各个变量同时动起来的——那就需要借助人工神经网络技术来帮忙了。这个“人工神经网络”技术，并不是医学名词，而是一种数据处理方法，它本身其实没有任何物理意义。简单来说，工作流程就是我们先给它一组样本值和目标值，在一定的条件下(这些条件是可控的，比如神经元与隐含层数量，计算精度等等)，让程序自动建立二者之间的关系，当关系建立起来之后，也就是所谓的“神经网络训练完毕”，再输入新的样本值，就会输出我们需要的目标值了。

只需要几十组有限元分析的样本数据，就可以由神经网络技术来建立结构参数和结构响应之间的关系，由此，便可由神经网络的程序来替代有限元来进行大规模重复计算——我曾经用一套代码，在4个小时内进行了150万次计算，这要是用ansys来计算，单独一次就需要将近1个小时，等都算完，估计我这人都不在了。

神经网络除了可以提高有限元计算效率之外，还有个很强大的功能值得挖掘，就是可以输出网络中每个层的神经元之间的权值和阀值，这对我们研究“多参数敏感性分析”至关重要，这些权值是在网络训练过程中，由程序在样本和目标之间建立起来的数据转换关系，这从侧面可以得出输入变量对结果的影响权值——这就是我们需要的多参数敏感性分析的解答了。

至于权值之间的算法，有好多种，我个人比较倾向于Garson算法：

关于这个算法的详情，大家按照这个关键词查询文件就可以了，这里不多说。

个人以为，利用神经网络来做的多参数敏感性分析，至少有两个优点：一个是可以实现多参数联动，另一个是可以量化出每个参数的影响因子到底是多少。很多时候，当许多因素“看起来差不多”的时候，通过这种方法，得到因素A的敏感性是21%，因素B的敏感性是20%，这样再排序，就方便多了。

荷载分段增量法非线性分析的实现思路

长河 — Wed, 11 Apr 2012 12:35:47 +0800

有一种非线性的表现形式，是荷载-挠度关系的非线性，即结构的非线性主要体现在几何变性后其自身的位置变化。此时，利用荷载分段增量法的思路来求解，就会更加直观、适用。

这个方法的核心，是将前一次分析结果中的结构变形提取出来，并在下一次分析中作为初始状态。在每一次分析中，都是线性关系，即“将非线性的问题转化为许多个线性问题的组合”。大体的思路，可以这么来考虑：

1.非受荷状态下，建立结构模型，采用参数化建模方式，结构关键点坐标赋值变量；

2.在一次线性分析之后，用*get命令得到结构关键点的位移，比如*get,weiyi,node,1,u,x得到节点1的x方向位移，赋值给weiyi

3.在下一次分析中，更新关键点的坐标值，这样就相当于修改了模型的几何信息，在新的荷载(荷载增量)条件下重新分析。

4.每一步的分析，都将外力和关键点的坐标存入到数组或者输出成数据文件，以备后期使用。

5.当然，从1~4步，这需要一个循环语句来实现。

具体实例就不上了，其实很简单的，只是循环和增量怎么处理的问题而已。

ansys实体模型中查看任意点应力的方法

长河 — Fri, 06 Apr 2012 18:57:54 +0800

两年前在做论文的时候，涉及到查看实体模型中任意点应力的问题，当时我用“路径”这么一个比较讨巧的做法。但后来，marain老师分享给我一段很经典的apdl宏文件，直接按照提示输入坐标值就可以在cmd窗口中看到改点的应力值了。

后来，在一个国外论坛中，我又见到了这段代码，看来是公认的终极解决方案啊。下面就是这段代码的全文了，你可以将其复制粘贴到记事本里，另存为mac文件，再运行即可。

另外强调一下，这真不是我写的，以往有类似分享的时候闹出了不少误会，令我情何以堪啊。

/nopr
!
! sxyz.mac a macro which calculates stresses at a
!              given point in space
!
! usage: sxyz,X,Y,Z,avg
!
!    where: X= x coordinate
!            Y= y coordinate
!            Z= z coordinate
!            avg=0 for averaged result, 1 for unaveraged result
!
! written by: John Crawford date: 3-5-98
!
! make sure we are in /post1
*get,ar20,active,,rout
*if,ar20,eq,31,then
!
! turn off annoying warning messages
/uis,msgpop,3
/nerr,0
!
! set up the vector _str to hold the data after the
! macro is done running
_str=
*dim,_str,array,11
!
! delete any existing path with the same name
padele,path1
! set up the path
path,path1,2,,
!
! define two path points
ppath,1,,arg1-.00001,arg2-.00001,arg3
ppath,2,,arg1,arg2,arg3
!
! get sx
!
! define the result on the path for sx
*if,arg4,eq,0,then
pdef,sx,s,x,avg
ar32=_status
*else
pdef,sx,s,x,noav
ar32=_status
*endif
! check to see if the coordinate is in the model space
! ar32 will be zero if it is, and 2 if it isn’t
*if,ar32,eq,0,then
! get the result from the last point on the path
*get,ar21,path,,max,sx
! multiply by the scale factor and then put the
! value into the result vector _str
_str(1)=ar21
!
! get sy
!
! define the result on the path for sy
*if,arg4,eq,0,then
pdef,sy,s,y,avg
*else
pdef,sy,s,y,noav
*endif
! get the result from the last point on the path
*get,ar22,path,,max,sy
! multiply by the scale factor and then put the
! value into the result vector _str
_str(2)=ar22
!
! get sz
!
! define the result on the path for sz
*if,arg4,eq,0,then
pdef,sz,s,z,avg
*else
pdef,sz,s,z,noav
*endif
! get the result from the last point on the path
*get,ar23,path,,max,sz
! multiply by the scale factor and then put the
! value into the result vector _str
_str(3)=ar23
!
! get sxy
!
! define the result on the path for sxy
*if,arg4,eq,0,then
pdef,sxy,s,xy,avg
*else
pdef,sxy,s,xy,noav
*endif
! get the result from the last point on the path
*get,ar24,path,,max,sxy
! multiply by the scale factor and then put the
! value into the result vector _str
_str(4)=ar24
!
! get syz
!
! define the result on the path for syz
*if,arg4,eq,0,then
pdef,syz,s,yz,avg
*else
pdef,syz,s,yz,noav
*endif
! get the result from the last point on the path
*get,ar25,path,,max,syz
! multiply by the scale factor and then put the
! value into the result vector _str
_str(5)=ar25
!
! get sxz
!
! define the result on the path for sz
*if,arg4,eq,0,then
pdef,sxz,s,xz,avg
*else
pdef,sxz,s,xz,noav
*endif
! get the result from the last point on the path
*get,ar26,path,,max,sxz
! multiply by the scale factor and then put the
! value into the result vector _str
_str(6)=ar26
!
! get s1
!
! define the result on the path for s1
*if,arg4,eq,0,then
pdef,s1,s,1,avg
*else
pdef,s1,s,1,noav
*endif
! get the result from the last point on the path
*get,ar27,path,,max,s1
! multiply by the scale factor and then put the
! value into the result vector _str
_str(7)=ar27
!
! get s2
!
! define the result on the path for s2
*if,arg4,eq,0,then
pdef,s2,s,2,avg
*else
pdef,s2,s,2,noav
*endif
! get the result from the last point on the path
*get,ar28,path,,max,s2
! multiply by the scale factor and then put the
! value into the result vector _str
_str(8)=ar28
!
! get s3
!
! define the result on the path for s3
*if,arg4,eq,0,then
pdef,s3,s,3,avg
*else
pdef,s3,s,3,noav
*endif
! get the result from the last point on the path
*get,ar29,path,,max,s3
! multiply by the scale factor and then put the
! value into the result vector _str
_str(9)=ar29
!
! get sint
!
! define the result on the path for sint
*if,arg4,eq,0,then
pdef,sint,s,int,avg
*else
pdef,sint,s,int,noav
*endif
! get the result from the last point on the path
*get,ar30,path,,max,sint
! multiply by the scale factor and then put the
! value into the result vector _str
_str(10)=ar30
!
! get seqv
!
! define the result on the path for seqv
*if,arg4,eq,0,then
pdef,seqv,s,eqv,avg
*else
pdef,seqv,s,eqv,noav
*endif
! get the result from the last point on the path
*get,ar31,path,,max,seqv
! multiply by the scale factor and then put the
! value into the result vector _str
_str(11)=ar31
!
! print the result
*if,arg4,eq,0,then
*msg,info,arg1,arg2,arg3
%/ Averaged Stress Results at X= %g, Y= %g, Z= %g
*else
*msg,info,arg1,arg2,arg3
%/ Unaveraged Stress Results at X= %g, Y= %g, Z= %g
*endif
!
*msg,info,ar21
SX = %g
*msg,info,ar22
SY = %g
*msg,info,ar23
SZ = %g
*msg,info,ar24
SXY = %g
*msg,info,ar25
SYZ = %g
*msg,info,ar26
SXZ = %g
*msg,info,ar27
S1 = %g
*msg,info,ar28
S2 = %g
*msg,info,ar29
S3 = %g
*msg,info,ar30
SINT = %g
*msg,info,ar31
SEQV = %g
!
! if the coordinate isn’t in the model space then print a message
! and quit
*else
*msg,ui,arg1,arg2,arg3
Coordinate x= %g y= %g z= %g is not within the model. %/ &
Check your coordinates and try again.
*endif
!
! turn warning messages back on
/uis,msgpop,2
/nerr,defa
!
*else
!
! print warning message if user is not in /post1
*msg,ui
*** You need to be in /POST1 to run SXYZ ***
*endif
/gopr

采访本上的城市

长河 — Sat, 24 Mar 2012 20:40:18 +0800

第一次知道王军，是因为他的《城记》，那本书让我认识了这位正在思考着“正在逝去的北京城”现象的记者。这本《采访本上的城市》，我是用等班车和做班车的时间，在手机上读完的。虽然txt版本的电子书没有图片，但毕竟长安街、故宫、CCTV大楼、国家大剧院、香港汇丰银行……书中的建筑大多是我原本就熟悉，或者实地见过甚至在其施工现场实习过，所以即便没有影像也能用记忆来弥补。

王军并不是一位建筑师，但他的这部作品却值得城市规划专业的人读一读。这里面涉及到许多概念，比如“尺度”，对于建筑对于城市而言，怎么样的一个尺度是合理的、美的尺度，也许在城市中生活多年的我们并不了解。于是我们中的多数都喜欢大马路，喜欢高层建筑，以为这是现代化的象征——的确是现代化的象征，但不一定就是品质生活的象征。

有人说本书的内容看起来更像是新闻采访，但我觉得如贝聿铭等一流建筑师的访谈记录，里面传达出了建筑师对建筑和城市的思考，这远比新闻稿更有生命力。这并不是一本教程，但却很明确地在传递一种思考的习惯和谨慎的思维方式，这对我们路桥专业也很有用，起码许多失败的立交桥项目可以提前避免吧。

关于本书的精彩内容，我就不多透露了，看过之后，我忽然会觉得，建筑精品在我国能得以保存，是一件多么不容易的事情。

域名被“误伤”导致最近网站访问不正常

长河 — Fri, 16 Mar 2012 08:34:32 +0800

已经好几天，陆续有朋友向我反应，我的博客无法访问了。一开始，我还以为是他们的DNS被劫持或者污染，因为毕竟我和非联通的一些朋友访问还是正常的。但近三天，的确是遇到了个麻烦事，网站的确访问不了了——

问题现象

域名ping和解析都没有问题，在dnspod上做“域名健康检查”，结果现实网站可以访问。但事实上，的确不能。。。同时，虚拟主机空间的Cpanel后台和FTP均正常工作，分别联系过空间提供商和dnspod客服，他们都觉得问题很蹊跷——因为一切的数据显示都正常。

可能的答案

于是，我猜测这是godaddy域名的问题。前些天imharbin.com在哈尔滨网通和联通用户中访问不稳定，经衡天小张的提醒，知道godaddy最近在抽风，于是便使用了dnspod，问题得到解决。但同样使用dnspod解析的10kn.com却没那么幸运。

考虑到imharbin.com是已经有ICP备案号的域名，而10kn.com是曾经被删除备案信息的域名。我想问题大概就出现在这里吧，我们国家的防火墙可能会把它觉得不太靠谱的域名“定点清除”。10kn.com就是讨论技术性的问题，无毒无公害，此次大概是被误伤所致。

为什么现在可以访问了呢？

因为，两会结束了，谢谢

以后再不能访问怎么办？

请关注新浪微博 @长河

或者，通过备用域名dachanghe.com和邮箱dachanghe@163.com都可以找到我

博客可能会更新不稳定，但肯定不会关停的，这个大家放心好了

无巧不成书

长河 — Fri, 09 Mar 2012 09:36:59 +0800

这个世界真是小。

昨天下午，在爱建路跟朋友喝茶，聊一聊我们共同的理想。这原本就是一次短聚，大家都很忙，交流完毕意见，一个小时后就散了。在寻找公交站的时候，不忍还是想去看一看哈尔滨车辆厂旧址。这其实是哈尔滨的伤心地。哈尔滨车辆厂，建于1898年，1903年迁到现厂址并修建了东北地区最早也是规模最大的现代化重工业工厂——东省铁路哈尔滨总工厂。

这座将近1平方公里的百年老厂，有着非同一般的历史价值：

1905年11月，俄国布尔什维克在这个工厂成立，标志着马列主义最早传入哈尔滨
1905年成立中俄两国人参加的中国最早的工会，同年12月参与中东铁路全路大罢工
1907年，这个工厂的工人在太阳岛举行了中国首次庆祝“五一”的活动
1917年11月7日俄国“十月革命”后，工人与俄国士兵11月10日集会，用多次罢工支援布尔什维克
1946年4月28日，哈尔滨成为中国第一个解放的大城市以后，这个工厂成为解放区最大的工厂，生产军工产品，是对解放战争贡献最大的工厂
百年老厂经历了沙俄独占、中苏和办、日伪统制、中苏共管和新中国等不同的时期
直到2002年9月，厂区内仍保留着5座建于1903-1907年间的厂房和1座办公楼；机械加工分厂车间内，大部分20世纪初期产于俄、法、美、德等国的老机床仍在使用。

然而，2006年3月10日，上海爱建集团违法拆除铸铁车间，他们要在这里建住宅小区。为了利益，一家来自上海的公司，毁灭了哈尔滨唯一申请世界遗产的希望。

如今，保留下来的一座厂房，两个烟囱，在现代小区中显得尤为突兀，无知群众还以为这是保护利用的典型，成为若干所谓文艺青年们镜头下的模特。

昨天很偶然，看到厂房的门开着，我就进去，想看看厂房内部的现状如何了。碰巧的是，与里面做测绘的师傅聊天之后发现，我们竟共同做过黑河的同一座桥梁项目。只是他负责做下部结构定位放线，我负责做上部结构的施工监控。我们都很感慨，当年在做项目的时候无缘相见，却在多年后的哈尔滨见面了。世界就是这么小吧。他是个实在人，我们交流了不少技术问题，为我提供了另一种看待专业的思维方式，在专业领域，还能这样无意间地交到朋友，真是一件神奇的事。

愉快的聊天之后，我走出百年厂房，回头望去，又是伤心——2006年，正在上大学的我竟然对这场“车辆厂保卫战”毫无所知，愧为哈尔滨人啊！

王澍：赢得中国首个Pritzker建筑大奖

长河 — Tue, 28 Feb 2012 03:34:45 +0800

王澍

在看到芝加哥论坛报的这篇报道以前，很惭愧，我对王澍一无所知。他是一位杰出的艺术家和建筑师，这不是因为他获得了普利兹克建筑大奖，为中国建筑师赢得了自信——当然，获奖这点事儿，他只能代表自己，但从某种意义上来说，中国的建筑师集体不自信由来已久。这体现在建筑师们一味地跟随西方建筑学派再走，放弃了太多原本属于中国传统文化中的哲学思想和建筑元素。

建筑，我觉得是当今最接近博物学的一门学问。一座好的建筑，需要设计师从历史、人文、地理以及心理学、人体工程学等诸多方面综合考虑。并不是放上一个穹顶就是拜占庭式建筑，有个尖拱就可以说是哥特式建筑。

可悲的是我们如今的生活就是这样的，我们太过于讲究“元素”。比如妆点了盘长或者祥云，就说是有中国传统文化元素的设计，比如在建筑外立面粘贴上几处雕刻花纹，就说是巴洛克艺术。一座建筑需要有的是统一的思想，而无论窗、门、还是楼梯，这些都是表达建筑整体思想的“单词”，这些单词汇总到一起，成为语言。就好像果戈里所说的“当歌曲和传说已经缄默的时候，建筑还在说话”。

我觉得，王澍是一个符号。

中国势必有许许多多的建筑师在考虑，如何给与中国的建筑生命，让它们在迅速成长的城市中拥有属于自己的言语权。王澍做到了，这就可以给其他正在路上，或者正在犹豫中的建筑师一点信心。面对自然，解决人与建筑、人与社会的问题，这在过去几千年来一直是中国建筑师们的强项，如今是该复兴的时候了。

王澍的代表作，宁波历史博物馆，足以表达他想要表达的声音。[另外公认的代表作是中国美院象山中心校区]

宁波历史博物馆

这栋建筑的原址，曾经是一个美丽的村庄，城市化的疯狂蔓延，摧毁了这里曾有过的风貌。千篇一律的钢筋水泥，取代了婀娜多姿的亭台楼榭，马路越来越宽，城市中道路的功能从原来的以散步为主变为今天的以快速通行为主。借用王军在《采访本上的城市》一书中的观点，我们正在建设的，只是“郊区文化”罢了。

王澍在设计宁波历史博物馆的时候，采用了与众不同的构想。他在宁波搜集旧城改造中积累下来的旧砖瓦、陶片，并用这些“废弃”的建材，建造了24米高的 “瓦爿墙”。同时还运用具有江南特色的毛竹制成特殊模板清水混凝土墙，毛竹随意开裂后形成的肌理效果清晰地显现。

王澍这样谈到自己的设计初衷：“使用‘瓦爿墙’，大量使用回收材料，节约了资源，体现了循环建造这一中国传统美德，一方面除了能体现宁波地域的传统建造体系、其质感和色彩完全融入自然外，另一意义在于对时间的保存，回收的旧砖瓦，承载着几百年的历史，它见证了消逝了的历史，这与博物馆本身是‘收集历史’这一理念是吻合的。而‘竹条模板混凝土’则是一种全新创造，竹本身是江南很有特色的植物，它使原本僵硬的混凝土发生了艺术质变。” [来源]

这样的设计，体现了建筑师对城市的思考，对时代变迁的思考，以及对生命的思考。

他的工作室有个有趣的名字，叫“业余建筑工作室”，但他们的作品却得到了普利兹克大奖评委会的认可。

关于普利兹克大奖：

这是由美国芝加哥亿万富翁Pritzker家族设置的奖项，专门表彰对人类有重大贡献的建筑师(significant contributions to humanity)，奖品是十万美元、一枚青铜奖章和全球性的荣誉。这个奖项自1979年授予美国建筑师飞利浦·约翰逊(Philip Johnson)以后，这是第一次授予中国建筑师。

王澍说，他能得到这个奖项，“非常惊讶”。但被世界认可，对于他来说，完全在情理之中，因为他的作品之中蕴含着对文化和传统的思考。文化是我们的基因，希望有更多的人能认识到这一点，不要走其他三个文明古国的老路。

小链接

王澍，男，1963年11月出生，2000年获同济大学建筑学博士学位。现任中国美术学院建筑艺术学院院长、建筑系主任、教授、硕士生导师、建筑学学科带头人、浙江省高校中青年学科带头人。

P方法小结

长河 — Sun, 12 Feb 2012 21:01:26 +0800

在ansys分析过程中，若要获得更精细的解答，一般需要通过细分网格的方法实现。但也许有时候网格划分不是那么方便，此时就需要在网格密度不变的前提下，获得更精确的解答——这就需要P方法了。

P方法的基本原理，是提高单元位移形函数的多项式次数，以达到收敛需要的精度。

适用性

P方法用于线弹性结构分析—实体和壳体。P单元由以下5种单元：

2-D Quadrilateral (Plane145)
2-D Triangle (Plane146)
3-D Brick (Solid 147)
3-D Tetrehedron (Solid 148)
3-D Shell (Solid 150)

实现步骤

1.单元阶次的控制

PPRANGE,2,8
ET,1,PLANE145,3,8

et用来对某种类型的单元进行阶次设置，PPRANGE进行全局设置，两者兼有的时候，以et优先。

2.网格划分，不建议采用人工建立节点的方法，可自动划分网格。

3.为提高求解效率，求解器推荐条件共轭梯度(PCG)算法。

4.精度控制相关命令：PCONV；区域排除命令：PEXCLUDE；选择p 方法的计算模型中将包含/排除的单元：ESEL,,PINC/<ESEL,,PEXC；

下面小例来自清华大学的课件，非我原创，仅供参考。

/PREP7 !进入前处理
/TITLE, EX 4.3, p-Method Analysis of Planar Plate with hole (Zeng P)
ANTYPE,STATIC !设定为静力结构分析
ET,1,PLANE145,,,3 !设置 1 号单元(p 型单元)
MP,EX,1,2e5 $MP,PRXY,1,0.25 !设置 1 号材料弹性常数
L=150 $H=80 $R1=5 $R2=10 $t=0.5 !设置几何参数
Pressure=100 !设置载荷参数
RECTNG,0,L,0,H !建立一个矩形面
WPOFFS,L/3,, !平移工作平面(x1=L/3)
PCIRC,R2,0,0,360 !建立一个圆面(R2)
WPOFFS,L/3,, !再平移工作平面(x2=x1+L/3)
PCIRC,R1,0,0,360 !建立一个圆面
R,1,t !设置厚度
ASBA,1,2 !将面 1 减面 2，则生成面 4
ASBA,4,3 !将面 4 减面 3
SMRTSIZ,4 !设置自动网格划分指数
AMESH,1 !对面进行单元划分
NSEL,S,LOC,X,0 !选取 x=0 节点
D,ALL,UX,0 !施加对称约束
!DSYM,SYMM,X !施加对称约束(与上一条命令的效果相同)
NSEL,S,LOC,Y,0 !选取 y=0 节点
DSYM,SYMM,Y !施加对称约束
NSEL,S,LOC,X,L !选取右边界节点
SF,ALL,PRES,-Pressure !对所选择节点施加载荷
ALLSEL !选取所有对象
FINISH !退出前处理
/SOLU !进入求解模块
P1=NODE(2*L/3,R1,0) !获取位置(2*L/3,R1,0)处的节点号，赋给参数 P1
P2=NODE(L/3,R2,0) !获取位置(L/3,R2,0)处的节点号，赋给参数 P2
PMOPTS !设置 p 方法中固定单元 p 阶次的误差，默认 5%
PCONV,0.5,S,X,P1 !针对 P1 节点，采用 p 方法，设置应力 SX 的误差在 0.1%以内
PCONV,0.5,S,X,P2 !针对 P2 节点，采用 p 方法，设置应力 SX 的误差在 0.1%以内
SOLVE !进行求解
/POST1 !进入后处理
SET,1 !调出结果的 1 号数据集
PLNSOL,S,X !图形显示 x 方向的应力计算结果
PLDISP,1 !图形显示受力后的模型结构
*GET,P1_SX,NODE,P1,S,X !获取节点 P1 处的应力 SX，赋给参数 P1_SX
*GET,P2_SX,NODE,P2,S,X !获取节点 P2 处的应力 SX，赋给参数 P2_SX
PLCONV !图形显示 p 方法的收敛曲线
PPLOT !图形显示实际所用单元 p 阶次
*STATUS !列显参数的内容
FINISH !结束