ANSYS中混凝土的计算问题(分离式)解析
ANSYS钢筋混凝土板受力分析
ANSYS分析钢筋混凝土结构技巧及实例详解
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命令流
lsel,r,loc,y,0
dl,all,,uy
dl,all,,uz
lsel,all
lsel,s,loc,z,3பைடு நூலகம்00
lsel,r,loc,y,0
dl,all,,uy
lsel,all
ksel,s,loc,x,0
ksel,r,loc,y,0
!方法：分离式；solid65和link8
!材料：混凝土采用concr和钢筋为弹性材料，但不考虑压碎
!---------------------------------------------------
!为方便，假定钢筋置于梁底两侧.
!===================================================
lesize,all,,,20
lmesh,all
lsel,s,loc,z,0
lesize,all,,,4
vsel,all
vatt,1,1,1
mshape,0,3d
mshkey,1
vmesh,all
allsel,all
finish
/solu
!施加约束
/config,nres,2000
/prep7
!定义单元及其材料特性等
rd0=20.0 !钢筋直径
et,1,solid65
et,2,link8
mp,ex,1,33e3
mp,prxy,1,0.20
r,1
hntra=28
hntrl=2.6
tb,concr,1
/post1
pldisp,1
钢筋混凝土梁ansys分析附命令流
钢筋混凝⼟梁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中进行混凝土计算时可能遇到的三个主要问题，包括材料属性的选择、边界条件的设定以及网格剖分的影响。
针对这些问题，我们提出了相应的解决方法，同时也提醒读者在使用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,ADINA在钢筋混凝土结构非线性分析中的应用与算例分析
表1跨中截面的应变和应力计算 结果对比
结论与展望
结论与展望
通过以上分析和算例对比，我们可以得出以下结论：ANSYS和ADINA在钢筋混 凝土结构非线性分析中均具有较高的准确性和可靠性，二者的计算结果非常接近。 这些软件提供了丰富的建模功能、强大的求解器和灵活的材料本构模型选择，使 得它们能够广泛应用于各种复杂的钢筋混凝土结构非线性分析问题中。
ANSYS和ADINA在钢筋混凝土结构非线性分析中的应用步骤和常见问题
钢筋混凝土结构非线性分析基本原理和方法
ANSYS和ADINA都是广泛用于工程模拟的有限元软件，它们在钢筋混凝土结构 非线性分析中具有相似的应用步骤。具体而言，可以分为以下几个步骤：
钢筋混凝土结构非线性分析基本原理和方法
1、建立模型：根据实际结构，建立相应的计算模型，包括几何形状、材料属 性、边界条件等。
谢谢观看
ANSYS,ADINA在钢筋混凝土结构 非线性分析中的应用与算例分
析
目录
01 引言
03 表1跨中截面的应变 和应力计算结果对比
02 算例分析 04 结论与展望 Nhomakorabea 内容摘要
ANSYS和ADINA在钢筋混凝土结构非线性分析中的应用与算例分析
引言
引言
钢筋混凝土结构非线性分析在工程实践中具有重要意义，它能够帮助工程师 们更准确地预测结构的性能和行为。在众多非线性分析软件中，ANSYS和ADINA是 最常用的两个工具。本次演示将介绍这两个软件在钢筋混凝土结构非线性分析中 的应用原理、方法和算例，并对其优势和不足进行分析，以期为工程师们提供有 益的参考。
ANSYS Solid 65分析混凝土讲义
坐标系举例
整体笛卡 儿坐标系
局部柱（球） 坐标系
清华大学研究生精品课程——《钢筋混凝土有限元》
Solid65单元的本构关系
可以使用以下本构关系： 等强硬化模型（Multilinear Isotropic Hardening） 随动硬化模型（Multilinear Kinematic Hardening） Drucker-Prager模型(DP模型)
0.0000 0.0001 5
清华大学研究生精品课程——《钢筋混凝土有限元》
Solid65中的裂缝模型
弥散固定裂缝模型(Smeared Fixed Crack Model) 采用Rankine最大拉应力准则，达到最大 拉力就开裂，裂缝方向不随主应变方向 变化 正交裂缝模型，一个高斯点上可以有3条 正交裂缝
Crush & Tension Cut Off
Crush
Ideal Failure Surface
Tension Cut Off
清华大学研究生精品课程——《钢筋混凝土有限元》
混凝土立方体模拟试验
5.00E+06 0.00E+00 -0.000 -0.000 -0.000 -0.000 -0.000 -5E-05 0 -5.00E+06 3 3 2 2 1 -1.00E+07 -1.50E+07 -2.00E+07 -2.50E+07 -3.00E+07
清华大学研究生精品课程——《钢筋混凝土有限元》
Solid 65的单元性态
每个单元有2×2×2＝8个高斯积分点 每个高斯积分点上最多可以有3条相互垂 直的裂缝 使用分布固定裂缝模型 可以使用整体式钢筋模型 在破坏准则之外，还可以考虑塑性
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
一、关于模型钢筋混凝土有限元模型根据钢筋的处理方式主要分为三种，即分离式、分布式和组合式模型。
考虑钢筋和混凝土之间的粘结和滑移，则采用引入粘结单元的分离式模型；假定混凝土和钢筋粘结很好，不考虑二者之间的滑移，则三种模型都可以；分离式和分布式模型适用于二维和三维结构分析，后者对杆系结构分析比较适用。
裂缝的处理方式有离散裂缝模型、分布裂缝模型和断裂力学模型，后者目前尚处研究之中，主要应用的是前两种。
离散裂缝模型和分布裂缝模型各有特点，可根据不同的分析目的选择使用。
随着计算速度和网格自动划分的快速实现，离散裂缝模型又有被推广使用的趋势。
就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。