ANSYS中混凝土的计算问题(分离式)解析
- 格式:doc
- 大小:55.00 KB
- 文档页数:9
0 前言利用ANSYS分析钢筋混凝土结构时,其有限元模型主要有分离式和整体式两种模型。
这里结合钢筋混凝土材料的工作特性,从模型建立到非线性计算再到结果分析的全过程讲述了利用ANSYS进行钢筋混凝土结构分析的方法与技巧,并以钢筋混凝土简支梁为例,采用分离式有限元模型,说明其具体应用。
1 单元选取与材料性质1. 1 混凝土单元ANSYS中提供了上百种计算单元类型,其中Solid65单元是专门用于模拟混凝土材料的三维实体单元。
该单元是八节点六面体单元,每个节点具有三个方向的自由度( UX , UY , UZ) 。
在普通八节点线弹性单元Solid45 的基础上,该单元增加了针对于混凝土的材性参数和组合式钢筋模型,可以综合考虑包括塑性和徐变引起的材料非线性、大位移引起的几何非线性、混凝土开裂和压碎引起的非线性等多种混凝土的材料特性。
使用Solid65 单元时,一般需要为其提供如下数据:1)、实常数(Real Constants) :定义弥散在混凝土中的最多三种钢筋的材料属性,配筋率和配筋角度。
对于墙板等配筋较密集且均匀的构件,一般使用这种整体式钢筋混凝土模型。
如果采用分离式配筋,那么此处则不需要填写钢筋实常数。
2)、材料模型(Material Model) :在输入钢筋和混凝土的非线性材料属性之前,首先必须定义钢筋和混凝土材料在线弹性阶段分析所需的基本材料信息,如:弹性模量,泊松比和密度。
3)、数据表(Data Table) :利用数据表进一步定义钢筋和混凝土的本构关系。
对于钢筋材料,一般只需要给定一个应力应变关系的数据表就可以了,譬如双折线等强硬化(bilinear isotropic hardening)或随动硬化模型( kinematic hardening plasticity)等。
而对于混凝土模型,除需要定义混凝土的本构关系外,还需要定义混凝土材料的破坏准则。
在ANSYS中,常用于定义混凝土本构关系的模型有:1)多线性等效强化模型(Multilinear isotropic hardening plas2ticity ,MISO模型),MISO模型可包括20条不同温度曲线,每条曲线可以有最多100个不同的应力-应变点;2)多线性随动强化模型(Multilinear kinematic hardening plas2ticity ,MKIN 模型),MKIN 模型最多允许5个应力-应变数据点;3)Drucker2Prager plasticity(DP)模型。
钢筋混凝⼟梁ansys分析附命令流钢筋混凝⼟⾮线性分析2015⼤作业上海交通⼤学陈明1、参数选择梁的截⾯宽度为200mm,上部配置2Φ8受压筋,混凝⼟的净保护层厚度为25 mm(从纵向钢筋外边缘算起),箍筋两端区采⽤8@100的双肢箍,中间区取8@200 双肢箍1)梁的截⾯⾼度选300mm;2)两加载间的距离选1000mm;3)混凝⼟选C30;4)纵向受拉钢筋配筋选218;2、描述选⽤的有限元模型及单元的特点采⽤ansys软件进⾏模拟计算,钢筋混凝⼟模型采⽤分离式模型,不考虑钢筋与混凝⼟之间的相对滑移。
混凝⼟采⽤solid65单元模拟,solid65⽤于模拟三维有钢筋或⽆钢筋的混凝⼟模型。
该单元能够计算拉裂和压碎。
在混凝⼟应⽤中,该单元的实体功能可以⽤于建⽴混凝⼟模型,同时,还可⽤加筋功能建⽴钢筋混凝⼟模型。
另外,该单元还可以应⽤于加强复合物和地质材料。
该单元由⼋个节点定义,每个节点有三个⾃由度:节点坐标系的x,y,z⽅向的平动。
⾄多可以定义三种不同规格的钢筋。
钢筋单元采⽤link180单元模拟,link180是⼀个适⽤于各类⼯程应⽤的三维杆单元。
根据具体情况,该单元可以被看作桁架单元、索单元、链杆单元或弹簧单元等等。
本单元是⼀个轴向拉伸⼀压缩单元,每个节点有三个⾃由度:节点坐标系的x,y,z⽅向的平动。
本单元是⼀种顶端铰接结构,不考虑单元弯曲。
本单元具有塑性、蠕变、旋转、⼤变形和⼤应变功能。
缺省时,当考虑⼤变形时任何分析中LINK180单元都包括应⼒刚化选项。
3、描述选⽤的混凝⼟与钢筋粘结滑移本构关系的具体形式、参数等。
钢筋的应⼒应变关系曲线考虑到极限塑性应变最⼤值为0.01,钢筋本构模型采⽤多线性模型kinh,初始弹性模量为Es=200000Mpa,强化系数为0.001。
混凝⼟的应⼒应变关系曲线混凝⼟选⽤各向同性的miso模型,当计⼊下降端时,程序报错,所以只取了前⾯的上升段,⽤5段折线模拟混凝⼟应⼒应变曲线。
ANSYS混凝土计算问题引言ANSYS作为一款强大的有限元分析软件,在混凝土结构设计中扮演着重要的角色。
但在使用过程中,我们会遇到一些问题,本文主要分析在ANSYS中进行混凝土计算时可能出现的问题以及相应的解决方法。
问题一:材料属性的选择在进行混凝土计算时,材料的选择是非常重要的,而在ANSYS中,材料属性的选择却非常的繁琐。
首先需要在ANSYS中创建新的材料属性,并指定相应的弹性模量、泊松比以及混凝土的强度参数等。
在这个过程中,我们需要确保选择的材料属性符合我们所使用的混凝土标准,否则计算结果可能会存在误差。
解决方法:建议在材料属性的选择上,我们应该非常谨慎,并注意选择我们所使用的混凝土标准对应的材料属性。
同时,在进行计算时添加合适的材料力学模型和屈服准则,以获得更为准确的计算结果。
问题二:边界条件的设定在进行混凝土计算时,经常需要设置不同的边界条件以模拟实际的工程情况。
然而,在ANSYS中,边界条件的设定较为繁琐,需要用户自己手动输入边界条件参数。
这样很容易出现手误,导致计算结果的误差增大。
解决方法:可以采用ANSYS提供的图形化界面进行边界条件的设定,避免手动输入参数导致的误差。
同时,我们应该明确每个边界条件的物理意义,并根据实际情况进行合理的选择和设置。
问题三:网格剖分的影响在ANSYS中,网格剖分对于计算结果的精度有着直接的影响。
对于混凝土的计算而言,网格剖分的密度直接决定了计算结果的准确性和精度。
解决方法:建议在进行混凝土计算时,应根据所需精度和计算要求,对模型进行合理的网格剖分。
在进行初始计算前,可以采用自适应网格划分方法,确保计算结果的准确性和精度。
结论本文主要介绍了在ANSYS中进行混凝土计算时可能遇到的三个主要问题,包括材料属性的选择、边界条件的设定以及网格剖分的影响。
针对这些问题,我们提出了相应的解决方法,同时也提醒读者在使用ANSYS进行混凝土计算时,需要格外的谨慎,选择合适的材料属性并进行合理的模型设置和计算分析。
发信人: rubors (宝马), 信区: FEA标题: 混凝土单元的应用(solid65)[转载]发信站: 同舟共济站(2002年09月08日17:16:34 星期天), 站内信件ANSYS中混凝土的计算问题【精华】最近做了点计算分析,结合各论坛关于这方面的讨论,就一些问题探讨如下,不当之处,敬请指正。
一、关于模型钢筋混凝土有限元模型根据钢筋的处理方式主要分为三种,即分离式、分布式和组合式模型。
考虑钢筋和混凝土之间的粘结和滑移,则采用引入粘结单元的分离式模型;假定混凝土和钢筋粘结很好,不考虑二者之间的滑移,则三种模型都可以;分离式和分布式模型适用于二维和三维结构分析,后者对杆系结构分析比较适用。
裂缝的处理方式有离散裂缝模型、分布裂缝模型和断裂力学模型,后者目前尚处研究之中,主要应用的是前两种。
离散裂缝模型和分布裂缝模型各有特点,可根据不同的分析目的选择使用。
随着计算速度和网格自动划分的快速实现,离散裂缝模型又有被推广使用的趋势。
就ANSYS而言,她可以考虑分离式模型(solid65+link8,认为混凝土和钢筋粘结很好,如要考虑粘结和滑移,则可引入弹簧单元进行模拟,比较困难!),也可采用分布式模型(带筋的solid65)。
而其裂缝的处理方式则为分布裂缝模型。
二、关于本构关系混凝土的本构关系可以分为线弹性、非线性弹性、弹塑性及其它力学理论等四类,其中研究最多的是非线性弹性和弹塑性本构关系,其中不乏实用者。
混凝土破坏准则从单参数到五参数模型达数十个模型,或借用古典强度理论或基于试验结果等,各个破坏准则的表达方式和繁简程度各异,适用范围和计算精度差别也比较大,给使用带来了一定的困难。
就ANSYS而言,其问题比较复杂些。
1.ANSYS混凝土的破坏准则与屈服准则是如何定义的?采用tb,concr,matnum则定义了W-W破坏准则(failure criterion),而非屈服准则(yield criterion)。
ANSYS混凝土问题分析1.关于模型钢筋混凝土有限元模型根据钢筋的处理方式分为三种:分离式、整体式和组合式模型◆分离式模型:把混凝土和钢筋作为不同的单元来处理,即混凝土和钢筋各自被划分为足够小的单元,两者的刚度矩阵是是分开来求解的,考虑到钢筋是一种细长的材料,通常可以忽略起横向抗剪强度,因此可以将钢筋作为线单元处理。
钢筋和混凝土之间可以插入粘结单元来模拟钢筋与混凝土之间的粘结和滑移。
一般钢筋混凝土是存在裂缝的,而开裂必然导致钢筋和混凝土变形的不协调,也就是说要发生粘结的失效与滑移,所以此种模型的应用最为广泛。
◆整体式模型:将钢筋分布与整个单元中,假定混凝土和钢筋粘结很好,并把单元视为连续均匀材料,与分离式模型不同的是,它求出的是综合了混凝土与钢筋单元的整体刚度矩阵;与组合式不同之点在于它不是先分别求出混凝土与钢筋对单元刚度的贡献然后再组合,而是一次求得综合的刚度矩阵。
◆组合式模型组合式模型分为两种:一种是分层组合式,在横截面上分成许多混凝土层和若干钢筋层,并对截面的应变作出某些假设,这种组合方式在钢筋混凝土板、壳结构中应用较广;另一种组合方法是采用带钢筋膜的等参单元。
当不考虑混凝土和钢筋二者之间的滑移,三种模型都可以。
分离式和整体式模型使用于二维和三维结构分析。
就ANSYS而言,可以考虑分离式模型:混凝土(SOLID65)+钢筋(LINK单元或PIPE单元),认为混凝土和钢筋粘结很好。
如要考虑粘结和滑移,则可引入弹簧单元进行模拟,如果比较困难也可以采用整体式模型(带筋的SOLID65)。
2.本构关系及破坏准则◆本构关系混凝土本构关系的模型对钢筋混凝土结构的非线性分析有重大影响。
混凝土的本构就是表示在各种外荷载作用下的混凝土应力应变的响应关系。
在建立混凝土本构关系时一般都是基于现有的连续介质力学的本构理论,在结合混凝土的力学特性,确定甚至调整本构关系中各种所需的材料参数。
通常,混凝土的本构关系可以分为线性弹性、非线性弹性、弹塑性及其他力学理论等四类。
ANSYSsolid单元整体建模的实常数问题Ansys里的solid65单元可以用来模拟混凝土,像陆新征,王新敏等人已经做过一些算例。
最近本人也在使用这个单元,以前看过一些资料,但是自己没有亲自动手作分析。
哎,什么事情都是看着容易做起来难啊,千万不能眼高于顶,呵呵。
这次做了几个算例,用到的是solid65的整体式建模。
先说一下,solid65单元模拟混凝土有两种方式。
一种是分离式建模,这又有两种思路:其一是不考虑混凝土与钢筋的滑移,钢筋和混凝土可以耦合或者共用节点,钢筋一般采用link8或者pipe20来模拟——使用这两种单元可能会因具体问题而有取舍,有的问题两者的计算结果相差很大,有时间我会做一些具体实例跟大家共同讨论;其二是考虑钢筋和混凝土之间的滑移,钢筋单元和混凝土单元之间的滑移用界面单元来模拟,在ansys中其combin39单元就是一个不错的选择。
另一种就是整体式建模了。
只有在整体式建模中solid65的实常数才真正派得上用场。
你比如,在分离式建模中混凝土单元的实常数是这样:r,1而在整体式建模中,要用到这个实常数;先看实常数的格式:r,real number,mat1,vr1,theta1,phi1,mat2,vr2,theta2,phi2,mat3,vr3,theta3,phi3real number——实常数编号;mat1,mat2,mat3——三个方向的钢筋的材料编号;theta1,phi1——前者是rebar1在xoy面上的投影与x轴的夹角,后者是rebar1与xoy平面的夹角;theta2,phi2——前者是rebar2在xoy面上的投影与x轴的夹角,后者是rebar2与xoy平面的夹角;theta3,phi3——前者是rebar3在xoy面上的投影与x轴的夹角,后者是rebar3与xoy平面的夹角;解释一下:与x轴平行的钢筋,两个角度应该是:0,0与y轴平行的钢筋,两个角度应该是:90,0与z轴平行的钢筋,两个角度应该是:0,90give several examples as follows:r,1,2,0.001,0,0,2,0.01,90,02,0.1,0,90 !reinforcements in x,y,z direc.or,r,1,2,0.001,0,90,2,0.01,90,02,0.1,0,0 !reinforcements inz,y,x direc.or,r,1,2,0.001,90,0,2,0.01,0,02,0.1,0,90 !reinforcements in y,x,z direc.上面的三个例子说明x,y,z方面的钢筋方向的定义顺序可以随意,没有特别的限制;同样有下面:r,1,2,0.001,0,0,2,0.01,90,0,2,0.1,0,90 !reinforcements in x,y,z direc.or,r,1,2,0.1,0,90,2,0.01,90,0,2,0.001,0,0 !reinforcements inz,y,x direc.or,r,1,2,0.01,90,0,2,0.001,0,0,2,0.1,0,90 !reinforcements in y,x,z direc.上面三个例子中实常数的意义是一样的。
一、关于模型钢筋混凝土有限元模型根据钢筋的处理方式主要分为三种,即分离式、分布式和组合式模型。
考虑钢筋和混凝土之间的粘结和滑移,则采用引入粘结单元的分离式模型;假定混凝土和钢筋粘结很好,不考虑二者之间的滑移,则三种模型都可以;分离式和分布式模型适用于二维和三维结构分析,后者对杆系结构分析比较适用。
裂缝的处理方式有离散裂缝模型、分布裂缝模型和断裂力学模型,后者目前尚处研究之中,主要应用的是前两种。
离散裂缝模型和分布裂缝模型各有特点,可根据不同的分析目的选择使用。
随着计算速度和网格自动划分的快速实现,离散裂缝模型又有被推广使用的趋势。
就ANSYS而言,可以考虑分离式模型(solid65+link8,认为混凝土和钢筋粘结很好,如要考虑粘结和滑移,则可引入弹簧单元进行模拟,比较困难!),也可采用分布式模型(带筋的solid65)。
而其裂缝的处理方式则为分布裂缝模型。
二、关于本构关系混凝土的本构关系可以分为线弹性、非线性弹性、弹塑性及其它力学理论等四类,其中研究最多的是非线性弹性和弹塑性本构关系,其中不乏实用者。
混凝土破坏准则从单参数到五参数模型达数十个模型,或借用古典强度理论或基于试验结果等,各个破坏准则的表达方式和繁简程度各异,适用范围和计算精度差别也比较大,给使用带来了一定的困难。
就ANSYS而言,其问题比较复杂些。
1 ANSYS混凝土的破坏准则与屈服准则是如何定义的?采用tb,concr,matnum则定义了W-W破坏准则(failure criterion),而非屈服准则(yield criterion)。
W-W破坏准则是用于检查混凝土开裂和压碎用的,而混凝土的塑性可以另外考虑(当然是在开裂和压碎之前)。
理论上破坏准则(failure criterion)和屈服准则(yield criterion)是不同的,例如在高静水压力下会发生相当的塑性变形,表现为屈服,但没有破坏。
而工程上又常将二者等同,其原因是工程结构不容许有很大的塑性变形,且混凝土等材料的屈服点不够明确,但破坏点非常明确。
定义tb,concr matnum后仅仅是定义了混凝土的破坏准则和缺省的本构关系,即W—W破坏准则、混凝土开裂和压碎前均为线性的应力应变关系,而开裂和压碎后采用其给出的本构关系。
但屈服准则尚可另外定义(随材料的应力应变关系,如tb,MKIN,则定义的屈服准则是Von Mises,流动法则、硬化法则也就确定了)。
2 定义tb,concr后可否定义其它的应力应变关系当然是可以的,并且只有在定义tb,concr后,有些问题才好解决。
例如可以定义tb,miso,输入混凝土的应力应变关系曲线(多折线实现),这样也就将屈服准则、流动法则、硬化法则等确定了。
这里可能存在一点疑问,即ANSYS中的应力应变关系是拉压相等的,而混凝土材料显然不是这样的。
是的,因为混凝土受拉段非常短,认为拉压相同影响很小,且由于定义的tb,concr中确定了开裂强度,所以尽管定义的是一条大曲线,但应用于受拉部分的很小。
三、具体的系数及公式1 定义tb,concr时候的两个系数如何确定?一般的参考书中,其值建议先取为0.3~0.5(江见鲸),原话是“在没有更仔细的数据时,不妨先取0.3~0.5进行计算”,足见此0.3~0.5值的可用程度。
根据我的经验和理由,建议此值取大些,即开裂的剪力传递系数取0.5,(定要>0.2)闭合的剪力传递系数取1.0。
支持此说法的还有现行铁路桥规的抗剪计算理论,以及原公路桥规的容许应力法的抗计剪计算。
ANSYS 用命令tb,concr matnum及随后的TBDATA,,C1,C2,C3,C4,C5,C6,C7,C8,C9共同定义了混凝土的破坏准则。
其中,C1为开口裂缝的剪切系数,一般为0.3-0.5,对于剪力墙取0.125,深梁0.25,一般梁取0.5;C2为闭口裂缝的剪切系数,一般为0.9-1.0;C3为单轴抗拉强度ft;C4为单轴抗压强度fc;C5为双轴抗压强度;C6为围压大小;C7为围压下双轴抗压强度;C8为围压下单轴抗压强度;C9为拉应力折减系数。
其中,C3和C4分别定义了混凝土的开裂大小和压碎破坏大小,但若C3=-1,则混凝土不考虑开裂,若C4=-1,则混凝土不考虑压碎;C5-C9可取默认值,否则若C5-C9有一个输入,在所有C1-C9都需要输入。
详见ANSYS之SOLID65帮助文件Table 65.1 SOLID65 Concrete Material Data。
SOLID65本构关系的判断在破坏准则之前。
一般情况下,不使用混凝土的压碎判据(即令C4=-1,后面的5个参数C5-C9不要设定,取默认值),而使用塑性或塑性下降段(即命令TB,MISO,1,1,5及TBPT,, ...所定义的非线性材料数据表)。
2 定义混凝土的应力应变曲线单向应力应变曲线很多,常用的可参考国标混凝土结构规范,其中给出的应力应变曲线是二次曲线+直线的下降段,其参数的设置按规范确定即可。
当然如有实测的应力应变曲线更好了。
四、关于收敛的问题ANSYS混凝土计算收敛(数值)是比较困难的,主要影响因素是网格密度、子步数、收敛准则等,这里讨论如下。
1 网格密度:网格密度适当能够收敛。
不是网格越密越好,当然太稀也不行,这仅仅是就收敛而言的,不考虑计算费用问题。
但是究竟多少合适,没有找到规律,只能靠自己针对情况慢慢试算。
2 子步数:NSUBST的设置很重要,设置太大或太小都不能达到正常收敛。
这点可以从收敛过程图看出,如果F范数曲线在[F]曲线上面走形的很长,可考虑增大nsubst。
或者根据经验慢慢调正试算。
3 收敛精度:实际上收敛精度的调正并不能彻底解决收敛的问题,但可以放宽收敛条件以加速吧。
一般不超过5%(缺省是0.5%),且使用力收敛条件即可。
4 混凝土压碎的设置:不考虑压碎时,计算相对容易收敛;而考虑压碎则比较难收敛,即便是没有达到压碎应力时。
如果是正常使用情况下的计算,建议关掉压碎选项;如果是极限计算,建议使用concr+MISO且关闭压碎检查;如果必设压碎检查,则要通过大量的试算(设置不同的网格密度、NSUBST)以达到目的,但也很困难。
5 其他选项:如线性搜索、预测等项也可以打开,以加速收敛,但不能根本解决问题。
6 计算结果:仅设置concr,不管是否设置压碎,其一般P-F曲线接近二折线;采用concr+miso则P-F曲线与二折线有差别,其曲线形状明显是曲线的。
*******************************************************************************例题1!----------------------------------------------------!题目:钢筋混凝土简支梁模拟计算!方法:分离式;solid65和link8!材料:混凝土采用concr和钢筋为弹性材料,但不考虑压碎!---------------------------------------------------!为方便,假定钢筋置于梁底两侧.!===================================================/config,nres,2000/prep7!定义单元及其材料特性等rd0=20.0 !钢筋直径et,1,solid65et,2,link8mp,ex,1,33e3mp,prxy,1,0.20r,1hntra=28hntrl=2.6tb,concr,1tbdata,,0.7,1.0,hntrl,-1mp,ex,2,2.1e5mp,prxy,2,0.30r,2,acos(-1)*0.25*rd0*rd0!定义梁体即单元划分blc4, , ,100,200,3000/view,1,1,1,1/ang,1gplot!定义网分时边长控制lsel,s,loc,z,1,2999lsel,r,loc,y,0latt,2,2,2lesize,all,,,20 !钢筋网格数目lmesh,alllsel,s,loc,z,0lesize,all,,,4 !截面上的网格数目4x4vsel,allvatt,1,1,1mshape,0,3dmshkey,1vmesh,allallsel,allfinish/solu!施加约束lsel,s,loc,z,0lsel,r,loc,y,0dl,all,,uydl,all,,uzlsel,all lsel,s,loc,z,3000lsel,r,loc,y,0dl,all,,uylsel,all ksel,s,loc,x,0ksel,r,loc,y,0dk,all,uxallsel,all!施加荷载qmz=0.3 asel,s,loc,y,200 sfa,all,1,pres,qmz allsel,allnsubst,40outres,all,alltime,qmz*10neqit,40solvefinish/post1pldisp,1 etable,zxyl,ls,1 plls,zxyl,zxyl,1finish/post26nsol,2,33,u,yprod,3,1,,,,,,1/100prod,4,2,,,,,,-1xvar,4plvar,3*******************************************************************************例题2!----------------------------------------------------!题目:钢筋混凝土简支梁模拟计算!方法:分离式;solid65和link8!材料:混凝土采用concr+Miso和钢筋为弹性材料,但不考虑压碎!增加网格密度!---------------------------------------------------!为方便,假定钢筋置于梁底两侧.!===================================================/config,nres,2000/prep7!定义单元及其材料特性等rd0=20.0 !钢筋直径et,1,solid65et,2,link8mp,ex,1,33e3mp,prxy,1,0.20r,1hntra=28hntrl=2.6tb,concr,1tbdata,,0.7,1.0,hntrl,-1tb,miso,1,,14tbpt,,0.0002,hntra*0.19tbpt,,0.0004,hntra*0.36 tbpt,,0.0006,hntra*0.51 tbpt,,0.0008,hntra*0.64 tbpt,,0.0010,hntra*0.75 tbpt,,0.0012,hntra*0.84 tbpt,,0.0014,hntra*0.91 tbpt,,0.0016,hntra*0.96 tbpt,,0.0018,hntra*0.99 tbpt,,0.0020,hntra*1.00 tbpt,,0.0025,hntra*0.95 tbpt,,0.0030,hntra*0.90 tbpt,,0.0035,hntra*0.85 tbpt,,0.0040,hntra*0.80 mp,ex,2,2.1e5mp,prxy,2,0.30r,2,acos(-1)*0.25*rd0*rd0 !定义梁体即单元划分blc4, , ,100,200,3000/view,1,1,1,1/ang,1gplot!定义网分时边长控制lsel,s,loc,z,1,2999lsel,r,loc,y,0latt,2,2,2lesize,all,,,20lmesh,alllsel,s,loc,z,0lesize,all,,,4vsel,allvatt,1,1,1 mshape,0,3dmshkey,1vmesh,allallsel,allfinish/solu!施加约束lsel,s,loc,z,0lsel,r,loc,y,0dl,all,,uydl,all,,uzlsel,all lsel,s,loc,z,3000lsel,r,loc,y,0dl,all,,uylsel,all ksel,s,loc,x,0ksel,r,loc,y,0dk,all,uxallsel,all!施加荷载qmz=0.3 asel,s,loc,y,200 sfa,all,1,pres,qmz allsel,alloutres,all,alltime,qmz*10nsubst,40neqit,40solvefinish/post1pldisp,1 etable,zxyl,ls,1 plls,zxyl,zxyl,1finish/post26nsol,2,33,u,y prod,3,1,,,,,,1/100 prod,4,2,,,,,,-1xvar,4plvar,3。