(1) 如果你模拟结构体中裂缝扩展过程的模拟，在Ansys中可以用全解耦损伤分析方法来近似模拟裂缝扩展，我曾用Ansys软件中提供的可以定义10,000个材料参数和单元ekill/alive 功能完成了层状路面体中表面裂缝和反射裂缝在变温作用下的扩展过程的模拟。我模拟的过程相对来说比较简单，模拟过程中我们首先要知道裂缝的可能扩展方向，这样在裂缝可能扩展的带内进行网格加密处理，加密到什么程度依据计算的问题来确定。
(2) 如果采用断裂力学理论计算含裂缝结构体的应力强度因子，建模时只需在裂尖通过命令kscon生成奇异单元即可。Ansys模块中存在的断裂力学模块可以计算I、II、III型应力强度因子(线弹性断裂力学)和J积分(弹塑性断裂力学)，在Ansys中verification里面有一个计算I型应力强度因子的例子vm143，参见该例子就可以了。
(3) 如果通过断裂力学模拟裂缝的扩展过程，需要采用动态网格划分，这方面我没有做，通过Ansys的宏命令流应该可以实现。技术参考可参阅文献：杨庆生、杨卫.断裂过程的有限元模拟.计算力学学报，1997,14(4).
(4) 我现在做动荷载作用下路面结构体中应力强度因子的分布规律，我是通过位移插值得到不同时间点处的应力强度因子。如果想这样做，可参阅理论参考中关于应力强度因子计算说明。
1. 讨论两种Ansys求极限荷载的方法
（1）力加载
可以通过对应的方法（比如说特征值屈曲）估计结构的极限荷载的大致范围，然后给结构施加一个稍大的荷载，打开自动荷载步二分法进行非线性静力分析，最后计算会因不收敛终止，则倒数第二个子步对应的就是结构的极限荷载；另外，也可以选择弧长法，采用足够的子步（弧长法可以一直分析到极限承载力之后的过程）同样可以从绘制的荷载位移曲线或计算结果中找出结构的极限荷载。
（2）位移加载
给结构施加一个比较大的位移，打开自动荷载步二分法进行非线性分析，保证足够的子步数，这样也可以分析到极限荷载以后，通过绘制荷载位移曲线或查看相应结果文件也可知道结构的极限荷载。
希望众高手讨论一下
（1）弧长法求极限荷载的收敛性问题，如何画到荷载位移曲线的下降段？
（2）位移法求极限荷载的具体步骤？
2. 需要注意的问题
1. 由于SOLID 65单元本身是基于弥散裂缝模型和最大拉应力开裂判据，因此在很多情况下会因为应力集中而使混凝土提前破坏，从而和试验结果不相吻合，因此，在实际应用过程中应该对单元分划进行有效控制，根据作者经验，当最小单元尺寸大于5cm 时，就可以有效避免应力集中带来的问题；
2. 支座是另一个需要注意的问题。在有限元分析中，很多时候约束是直接加在混凝土节点上，这样很可能在支座位置产生很大的应力集中，从而使支座附近的混凝土突然破坏，造成求解失败。因此，在实际应用过程中，应该适当加大支座附近单元的尺寸或者在支座上加一些弹性垫块，避免支座的应力集中；
3. 六面体的SOLID 65 单元一般比四面体的单元计算要稳定且收敛性好，因此，只要条件允许，应该尽量使用六面体单元；
4. 正确选择收敛标准，一般位移控制加载最好用位移的无穷范数控制收敛，而用力控制加载时可以用残余力的二范数控制收敛。在裂缝刚刚出现和接近破坏的阶段，可以适当放松收敛标准，保证计算的连续性；
3. 关于下降段的问题
1）在实际混凝土中都有下降段，但是在计算的时候要特别小心下降段的问题。
2）下降段很容易导致计算不收敛，有时为了计算的收敛要避免设置下降段，采用rush 模型。
3）利用最大压应变准则来判断混凝土是否破坏。
4. Solid65单元中的破坏准则
1）采用Willam&Warnke五参数破坏准则
2）需要参数：
单轴抗拉强度，单轴，双轴抗压强度，围压压力，在围压作用下双轴，单轴抗压强度
5. 近来我对混凝土单元进行了一点思考，有一些想法，贴在下面，共同探讨：
1）分析混凝土结构，选择合理的材料特性是建立模型的关键，所以有必要弄清混凝土的材料特性。混凝土是脆性材料,并具有不同的拉伸和压缩特性。典型混凝土的抗拉强度只有抗压强度的8％-15%。
在ANSYS中，对于混凝土单元，材料特性ANSYS要求输入以下数据（为了清楚起见，我将几个系数均译为了中文）：弹性模量、泊松比、张开与闭合滑移面的剪切强度缩减系数、抗拉与抗压强度、极限双轴抗压强度、周围静水应力状态、静水应力状态下单轴与双轴压缩的极限抗压强度、断裂发生时刚度乘子。其中，1～6是必须输入的，7～11要么不输入，都采用默认值，如果输入其中一个，其他的都需要输入；另外，与在0～1之间取值，具体如何取值，是值得探讨的话题，但有一点是肯定的，不能将剪切缩减系数，取的太小，否则，就很难不收敛，据我体会，分析一个梁的极限荷载时，剪切缩减系数的取值影响也不是很大；
2）据我理解，如果定义：TB,concer,mat(mat是需要定义的材料号)，则混凝土定义了破坏准则，没有定义屈服准则，主要是W-W模型。由于混凝土材料的复杂性（太随机了），很难得到一个完全适合混凝土的屈服的材料模型；
3）如果考虑混凝土的压碎破坏，有限元模型会较早失效，得不到真实极限荷载，建议在研究钢筋混凝土结构极限荷载时，关闭混凝土压碎能力；材料模型的选取对荷载变形曲线路径影响不大，即模拟曲线与真实曲线相对应部分吻合较好；不考虑混凝土的压碎破坏，并不意味着不考虑混凝土的抗压能力，相反，为了得到较准确的极限荷载，采用受压混凝土模型是必需的，也只有采用受压混凝土模型才能得到整个荷载变形曲线；
6. do a mesh sensitivity analysis, which means you need analyze the model with three different mesh sizes in the span direction, and maybe you can find something.
7. 开裂荷载就是第一到裂纹出现的时候，所加荷载即为开裂荷载；至于如何检测，我也没有好的办法，就是在开裂荷载附近，将荷载不加的很小！观察第一条裂纹的出现.
8. 由于没有考虑混凝土的应力应变关系，所以荷载-位移曲线接近双折线；我的观点是：
１）：仅设置CONCR，不管是否设置压碎，极限荷载偏小，荷载位移曲线一般接近二折线；采用CONCR＋MISO则荷载位移曲线形状明显是曲线的，更接近试验结果，所以设置CONCR 后，还要需要具体定义混凝土的应力应变关系；
２）：不考虑压碎时，计算相对容易收敛；而考虑压碎则比较难收敛，即便是没有达到压碎应力时，且有限元模型会较早失效，得不到真实极限荷载。
3）关掉压裂，应该在考虑其他屈服准则，如von mises。
****米赛斯准则在特定的条件下是适合混凝土的．米赛斯准则不适合混凝土是因为它不能包括静水压力效应，对简单的简之梁摸拟，根本没有考虑静水压力效应，所以还是可以用米赛斯准则的；
****我很赞成你的观点：＂定义Concr时把压碎特性去掉（抗压强度为“-1”）肯定是不合理的＂；我也在相关的帖子，发表过向类似的观点；比如，不考虑混凝土压碎，如何观察混凝土结构的裂缝发展？我也一直在考虑此问题，只是高我们这一行的人太少，连找一个讨论的人都很难！！
我的观点是：如果必设压碎检查，则要通过大量的试算，并设置不同的网格密度、子步数，以达到目的；
***另外，在ANSYS软件的帮助里，明确提到混凝土可以将抗压强度设为－１，以不考虑混凝土的压碎功能！原叙述如下：
Absence of the data table removes the cracking and crushing capability. A value of -1 for constant 3 or 4 also removes the cracking or crushing capability，
注：constant 3:Uniaxial tensile cracking stress.
constant 4:Uniaxial crushing stress (positive).
即混凝土开裂、压碎均可去掉。
***有时我想是不是＂ANSYS在混凝土这一块功能还有待于改善＂？
9. 有限元方程求解
采用Newton-Raphson迭代法求解，求解时采用多荷载步，每荷载步又设若干子步。为便于加载，将整个求解过程分为四个阶段：（1）加载到开始出现裂缝；（2）很多裂缝出现；（3）钢筋屈服；（4）一直到破坏。第一阶段为线性阶段，可以加比较大的荷载步，当出现裂纹后，就逐步减小荷载步。进入非线性后，方程较难收敛甚至不收敛，在裂缝刚刚出现和接近破坏的阶段，可以适当放松收敛标准，保证计算的连续性
10. 我在用ansys计算沈聚敏教授做过的试验（《钢筋混凝土构件的变形和延性》清华大学，建筑结构学报1980.2）总共10根梁。
发现我不考虑混凝土压碎的时候，得到的结果基本上与试验相同。
但是考虑压碎的时候（把-1改为轴心抗压强度，命令流相同），计算得到的结果与试验值相去甚远。这是为何呢？
在考虑压碎的前提下，我把收敛准则里面的0.05改为0.2之后，得到的最大承载力得到了显著的提高，可以和试验值相当这样的话，我就感到很不解了，应该说收敛准则，只是决定是否收敛，而不能改变结果阿？？？（而且我感觉0.2是不是不行啊，最大应该就是0.05）
子步数和迭代次数应该是越大越好吧，（在不考虑计算时间的前提下），我好像感觉有人说，大了也不好，不好收敛。表示怀疑？？？
11. 如何显示SOLID 65单元内部配筋方向？
操作步骤：
1）使实体透明化：
Utility Menu>PlotCtrls>Style>Translucency
注意：1——完全透明；0——完全不透明
2）显示与实常数相关的单元形状：
Utility Menu>PlotCtrls>Style>Size and Shape
将ESHAPE开关打开。
3）消影
Utility Menu>PlotCtrls>Style>Edge Options
对GLINE进行选择。
4）将单元缩减以更好地观察
Utility Menu>PlotCtrls>Style>Size and Shape
对SHRINK开关进行选择缩减比例0％～50％。
12. 为什么ANSYS对WW破坏面进行改动？
1) 混凝土材料本身是一种非常复杂的材料，受压破坏和受拉破坏表现完全不同；
2) WW破坏曲面中没有区分这两种破坏形式，但实际程序计算中却需要加以区分；
3) 人为强制在压-压-拉分区和拉-拉-压分区中按开裂处理，体现了理论研究和实用计算程序编写之间的差异。
13. Solid 65的本构关系
弹塑性行为
1) 受拉弹性
2) 受压可以选择一些塑性模型：
a) Von Mises屈服面
b) Drucker-Prager屈服面
c) 理想弹塑性
可以使用以下本构关系：
1) 等强硬化模型（Multilinear Isotropic Hardening）
2) 随动硬化模型（Multilinear Kinematic Hardening）
3) Drucker-Prager模型(DP模型)，DP模型在混凝土中为理想弹塑性模型
14. Solid65单元中的破坏准则
采用Willam&Warnke五参数破坏准则
需要参数：
1) 单轴抗拉强度，
2) 单轴，双轴抗压强度，
3) 围压压力，
4) 在围压作用下的双轴，单轴抗压强度
15. 混凝土与钢筋组合模型的选取
1. 在条件允许的情况下，优先使用整体式模型；
2. 滑移影响可以通过折减钢筋弹模加以模拟；
3. 在类似节点受往复荷载作用的问题，由于滑移严重，必须使用界面单元；
4. 预应力考虑预应力损失，必须使用界面单元。
16. ANSYS裂缝模型的说明
1. ANSYS的裂缝模型相对比较粗糙，使用时应加以慎重考虑；
2. 受拉软化给的是基于固定极限应变的软化模型，比较适用于配筋合适的钢筋混凝土构件，对素混凝土构件有很强的网格依赖性，效果一般不好；
3. ANSYS的裂面受剪模型为恒定的裂面剪力传递系数模型，在模拟受剪破坏时，剪力传递系数选择不当往往会高估构件受剪承载力。
17. ANSYS建模中需要注意的问题
!!单元尺寸大小
基于最大开裂应力准则，单元越细，应力集中越严重，开裂出现越早
解决方法
1) 使用半脆性裂缝模型，减小单元尺寸影响
2) 控制网格大小，单元尺寸不宜小于5cm
3) 控制网格划分，在容易出现应力集中的部位要避免过小的单元出现
!!施加支座
支座是个非常严重的应力集中部位，尽量避免把约束直接施加在支座上
解决方法
1) 加弹性垫块，利用圣维南原理减小应力集中
2) 加大支座部位单元尺寸，减小应力集中
18. 具体的系数及公式
1) 定义tb,concr时候的两个系数如何确定？
一般的参考书中，其值建议先取为0.3~0.5(江见鲸)，原话是“在没有更仔细的数据时，不妨先取0.3~0.5进行计算”，足见此0.3~0.5值的可用程度。根据我的经验和理由，建议此值取大些，即开裂的剪力传递系数取0.5,（定要>0.2）闭合的剪力传递系数取1.0。支持此说法的还有现行铁路桥规的抗剪计算理论，以及原公路桥规的容许应力法的抗计剪计算。
2) 定义混凝土的应力应变曲线
单向应力应变曲线很多，常用的可参考国标混凝土结构规范，其中给出的应力应变曲线是二次曲线＋直线的下降段，其参数的设置按规范确定即可。当然如有实测的应力应变曲线更好了。
19. 关于收敛的问题
ANSYS混凝土计算收敛(数值)是比较困难的，主要影响因素是网格密度、子步数、收敛准则等，这里讨论如下。
1）网格密度：网格密度适当能够收敛。不是网格越密越好，当然太稀也不行，这仅仅是就收敛而言的，不考虑计算费用问题。但是究竟多少合适，没有找到规律，只能靠自己针对情况慢慢试算。
2 )子步数：NSUBST的设置很重要，设置太大或太小都不能达到正常收敛。这点可以从收敛过程图看出，如果F范数曲线在[F]曲线上面走形的很长，可考虑增大nsubst。或者根据经验慢慢调正试算。
3）收敛精度：实际上收敛精度的调正并不能彻底解决收敛的问题，但可以放宽收敛条件以加速吧。一般不超过5%(缺省是0.5%)，且使用力收敛条件即可。
4）混凝土压碎的设置：不考虑压碎时，计算相对容易收敛；而考虑压碎则比较难收敛，即便是没有达到压碎应力时。如果是正常使用情况下的计算，建议关掉压碎选项；如果是极限计算，建议使用concr＋MISO且关闭压碎检查；如果必设压碎检查，则要通过大量的试算（设置不同的网格密度、NSUBST）以达到目的，但也很困难。
5）其他选项：如线性搜索、预测等项也可以打开，以加速收敛，但不能根本解决问题。
6）计算结果：仅设置concr，不管是否设置压碎，其一般P-F曲线接近二折线；采用concr+miso 则P-F曲线与二折线有差别，其曲线形状明显是曲线的。
20. 如何用AUTOCAD如何出ANSYS的图?有什么好办法吗?
这两种方法：
1）用apdl生成数据文件，然后用excell打开并处理，然后到autocad中绘制；如post26中的变量曲线即用此法。
2）直接用apdl生成DXF格式的文件，如变形图等。
21. 用ansys计算钢筋混凝土，当第一条裂缝出现（此时荷载很小）时，往往很难收敛，你可以改变一下混凝土的参数，或调整一下收敛准则，或改变网格划分，可能有用。
22. 我想问什么时候是开裂荷载呢，当有第一道裂纹是吗？那么在ansys中怎么检测了，用etable吗？还是别的什么命令呢？
开裂荷载就是第一到裂纹出现的时候，所加荷载即为开裂荷载；至于如何检测，我也没有好的办法，就是在开裂荷载附近，将荷载不加的很小！观察第一条裂纹的出现。有文献在开裂后取到每个子步5N。真是有耐心。
23. 请教，在混凝土的压碎图中，绿色小圆圈跟红色小圆圈的表示有什么区别
红色是第一次开裂的，绿色是第二次开裂的。
24. 混凝土裂缝显示
PLCRACK displays circles at locations of cracking or crushing in concrete elements. Cracking is shown with a circle outline in the plane of the crack, and crushing is shown with an octahedron outline. If the crack has opened and then closed, the circle outline will have an X through it. Each integration point can crack in up to three different planes. The first crack at an integration point is shown with a red circle outline, the second crack with a green outline, and the third crack with a blue outline.
Symbols shown at the element centroid (LOC = 1) are based on the status of all of the element's integration points. If any integration point in the element has crushed, the crushed (octahedron) symbol is shown at the centroid. If any integration point has cracked or cracked and closed, the cracked symbol is shown at the element centroid. If at least five integration points have cracked and closed, the cracked and closed symbol is shown at the element centroid. Finally, if more than one integration point has cracked, the circle outline at the element centroid shows the average orientation of all cracked planes for that element.
25. 关于ANSYS分析效果
1) 受弯——强度一般都比较准，刚度要差一些
2) 受剪——精度依赖于裂面剪力传递系数，要调整。
3) 受轴压——高轴压比，高围压，很难算准
4) 往复荷载——效果一般都不是很理想，除非很精细地调整参数
26. LNSRCH, Key — Activates a line search to be used with Newton-Raphson
You cannot use line search [LNSRCH], automatic time stepping [AUTOTS], or the DOF solution predictor [PRED] with the arc-length method [ARCLEN, ARCTRM]. If you activate the arc-length method after you set LNSRCH, AUTOTS, or PRED, a warning message appears. If you choose to proceed with the arc-length method activation, ANSYS disables your line search, automatic time stepping, and DOF predictor settings.
27. 钢筋混凝土整体模型计算小结
1) 单元太小易开裂，难收敛，且计算结果不稳定；
2) 子步数太少，计算结果不精确但曲线较光滑。
28. ANSYS文件指南
最大程度地减小由于误操作引起的文件覆盖等，我们建议您培养以下习惯：
1) 针对每个分析项目，设置单独的子目录；
2) 每求解一个新问题使用不同的工作文件名，在AYSYS启动对话框中设置工作文件名；
3）ANSYS的Output文件在交互操作中并不自动被写出，在交互操作中，您必须用Utility Menu: File > Switch Output to > File把output写到一个文件中；
4）分析完成后，您必须保存如下文件: log 文件( .log), 数据库文件( .db), 结果文件( .rst, .rth等), 载荷步文件(.s01, .s02, ...), 输出文件( .out), 物理环境文件(.ph1, .ph2, ...)。
注意：log 文件只添加，不会覆盖.
29. 制订分析方案是很重要的。一般考虑下列问题：
a) 分析领域
b) 分析目标
c) 线性/非线性问题
d) 静力/动力问题
e) 分析细节的考虑
f) 几何模型对称性
g) 奇异：是有限元模型中由于几何构造或载荷引起弹性理论计算应力值无限大。
h) 单元类型
i) 网格密度
j) 单位制
k) 材料特性
l) 载荷
m) 求解器
30. 在应力奇异处:
1）单元网格越是细化，越引起计算应力无限增加，并且不再收敛。
2）网格疏密不均匀时网格离散误差也大小不一（自适应网格划分结果是失败的或者网格错误）。
31. “热点”
1）对于面或体，热点为图形中心；
2）对于线，有三个热点。
为什么这一点非常重要? 需要在图形窗口拾取取图元时，应该点取图形的热点，确保拾取所需要的图元。这对于有多个图形重叠的情况非常重要.
32. 应力上下限
应力上下限可以帮助确定由于网格离散误差对模型的应力最大值的影响，显示或列出的应力上下限包括：
估计的上限- SMXB
估计的下限- SMNB
应力上下限限并不是估计实际的最高或最小应力.对于有些情况，SMXB 过于保守. 而有些情况比实际的要小.
应力上下限定义了一个”确信范围.” 如果没有其他的确凿的验证，就不能认为实际的最大应力低于SMXB.
33. 接触单元
注意：点对点接触只能用于低次单元.
接触12单元和接触52单元既能用直接生成法创建, 也能在重合节点处创建单元.
前处理器-> 创建-> 单元-> 在重合节点
( Preprocessor -> Create -> Elements -> At Coincid Nd)
接触12单元应该在重合节点间创建. 然而接触52单元要求1E-6的距离来定向单元.
34. 接触刚度
点对点(接触12单元和接触52单元)和节点对表面(接触48单元和接触49单元)接触单元都要求给出罚刚度.
罚刚度越大, 接触表面的侵入量越小. 然而, 若此值太大,则会由于病态条件而引起收敛困难.
可以通过一些实验来确定一个合适的接触刚度, 使求解收敛,而且侵入量可以接受.
35. 选择接触刚度
接触刚度是接触面的相对刚度的函数.
对于块状实体, 通常赫兹接触刚度(Hertz contact stiffness)适用于罚刚度, 可以这样估算: k = fE
式中f 是介于0.1~10之间的系数, E是较软的接触材料的弹性模量. 设f=1 通常是一个较好的起始值.
应该设置罚刚度比例系数(FKN)和拉格朗日乘子侵入比例系数(FTOLN). FKN 通常介于0.01~10之间. 对于体积变形问题, 用值1.0(默认), 对于以弯曲变形为主的问题, 用值0.1.
FTOLN 默认为0.1. 可以改变此值, 但若容差太小, 会使迭代数过多或不收敛.
36. 面对面接触处理
对于面对面接触单元, 一个面指定为“ 目标”面, 另一个面为“ 接触”面.
对于刚体对柔体接触, 刚体表面总是指定为目标面.
对于柔体对柔体接触, 接触面与目标面都与变形体相关联.
接触单元被约束, 不能侵入目标面. 然而, 目标单元能够侵入接触面.
37. 目标/接触面的指导方针
如果凸面与平面或凹面接触, 那么平面或凹面应该是目标面.
如果一个表面网格粗糙, 而另一个表面网格较细, 那么网格粗糙的表面应该是目标面.
如果一个表面比另一个表面的刚度大, 那么刚度大的表面应该是目标面.
如果一个表面划分为高次单元, 而另一个表面划分为低次单元, 那么划分为低次单元的表面应该是目标面.
如果一个表面比另一个表面大, 那么更大的表面应该是目标面.
38. 接触算法
选择一个接触算法关键字选项(2)
增强的拉格朗日法(关键字选项(2)=0)是缺省选项, 推荐于一般应用. 它对罚刚度不太敏感, 但是也要求给出一个侵入容差.
能够用罚函数法(关键字选项(2)=1)这个选项. 它推荐应用于单元非常扭曲、大摩擦系数和/或用增大的拉格朗日法收敛行为不好的问题.
39. 确定罚刚度
对于面对面接触单元, ANSYS基于单元类型、材料性质和它下面的单元尺寸确定接触刚度. 可以用实常数FKN给出接触刚度的一个比例系数或绝对值.
惩罚刚度(FKN)应该足够大, 使接触侵入量小。同时也应足够小, 使问题没有病态矩阵.
FKN值通常在0.01~10之间. 对于体积变形问题, 用值1.0 (默认), 对于弯曲问题, 用值0.1.
40. 对称接触
对称接触不如不对称接触有效. 然而, 许多分析需要用它(典型用于减少侵入). 对称接触增加了接触检查点的数目.
对称接触的准则:
目标面和接触面没有明显的区别.
目标面和接触面的网格都粗糙.
注意：用对称接触时, 后处理更困难. 接触压力是两个接触单元对的平均值.
41. 自接触
对于自接触, 使用不对称接触更有效, 但是难于预测接触面和目标面. 对于自接触, 用对称接触时, 只要简单地将目标单元和接触单元放在相同的表面上即可.
42. 刚体位移
如果在求解中的任一时刻, 两个物体没有联系, 刚度矩阵就会奇异. ANSYS将会发出一个负主元警告信息. 由于物体初始时没有联系, 要克服刚体位移有几个选项:
在“ 恰好碰上” 的位置建立几何模型
?动力学
?位移控制
?软弹簧
?用不分离接触(关键选项(12), 在后面讨论)
?调整初始接触条件
43. 接触刚度FKN
但对于大多数情况而言，最好用一个合理而不过度精确的刚度值进行第一次求解，然后用10倍于该值的刚度进行第二次求解，如果两者结果相差很小，而迭代数增加很多，那么我们则正好取得了曲线上的突变点，从而获得相当好的结果。
当接触单元的刚度为10e6时，可获得合理精确的结果。任何大于该值的刚度对下梁的偏移量没有什么影响，而求解所需的迭代数却显著的增加。对于这个题目，10e6的刚度是很适合的。但是，如果改变边界条件、网格密度、两梁之间的相对位置、材料特性或梁的几何形状，能获得满意结果的接触刚度值将是不同的。比如，如果网格密度增加，则接触单元数将
增加，每一个单元上的载荷将降低。如果接触单元数增加两倍，一个合适的接触单元刚度值应为原来的一半。
为了确定一个较好的接触刚度值，可能需要一些经验。用户可以按照下面的步骤来进行尝试：
1、开始时取一个较低的值。低估值要比高估值好，因为由一个较低的接触刚度导致的穿透问题，比过高的接触刚度导致的收敛性困难，要容易解决。
2、对前几个子步进行计算分析，直到最终载荷的一个比例（刚好完全建立接触）。
3、检查每一子步中的穿透量和平衡迭代次数。如果总体收敛困难是由过大的穿透引起的，那麽可能低估了FKN的值，或者是将FLOLN的值取得大小。如果总体的收敛困难是由于不平衡力和位移增量达到收敛值时需要过多的迭代次数，而不是由于过大的穿透量引起的，那麽FKN的值可能被估高。
4、按需要调整FKN和FTOLN的值，重新进行完整的分析。
我理解的接触问题求解过程，是一个调整接触刚度的过程，不知理解得是否对。接触分析，是要通过大量的结果画出一条曲线，选取曲线的最低点，作为最终结果。
44. 建立刚－柔接触对
如果想建立一个柔－柔接触对的话，可以将接触的两个面都划分网格，然后再用接触向导建立接触对。如果想建立一个刚－柔接触对的话，可以先将想定义成柔体的部分划分网格，然后定义接触对就可以了，如果将两个面都划分网格再定义接触对的话就又成了定义柔－柔接触对了。
45. ANSYS坐标系总结
工作平面（Working Plane）
工作平面是创建几何模型的参考(X,Y)平面，在前处理器中用来建模(几何和网格)
总体坐标系
在每开始进行一个新的ANSYS分析时，已经有三个坐标系预先定义了。它们位于模型的总体原点。三种类型为:
CS,0: 总体笛卡尔坐标系
CS,1: 总体柱坐标系
CS,2: 总体球坐标系
数据库中节点坐标总是以总体笛卡尔坐标系，无论节点是在什么坐标系中创建的。
局部坐标系
局部坐标系是用户定义的坐标系。局部坐标系可以通过菜单路径Workplane>Local CS>Create LC来创建。
激活的坐标系是分析中特定时间的参考系。缺省为总体笛卡尔坐标系。当创建了一个新的坐标系时，新坐标系变为激活坐标系。这表明后面的激活坐标系的命令。菜单中激活坐标系的路径Workplane>Change active CS to>。
节点坐标系
每一个节点都有一个附着的坐标系。节点坐标系缺省总是笛卡尔坐标系并与总体笛卡尔坐标系平行。节点力和节点边界条件（约束）指的是节点坐标系的方向。时间历程后处理器/POST26 中的结果数据是在节点坐标系下表达的。而通用后处理器/POST1中的结果是按结果坐标系进行表达的。
例如: 模型中任意位置的一个圆，要施加径向约束。首先需要在圆的中心创建一个柱坐标系并分配一个坐标系号码(例如CS,11)。这个局部坐标系现在成为激活的坐标系。然后选择圆上的所有节点。通过使用"Prep7>Move/Modify>Rotate Nodal CS to active CS", 选择节点的节点坐标系的朝向将沿着激活坐标系的方向。未选择节点保持不变。节点坐标系的显示通过菜单路径Pltctrls>Symbols>Nodal CS。这些节点坐标系的X方向现在沿径向。约束这些选择节点的X方向，就是施加的径向约束。
注意：节点坐标系总是笛卡尔坐标系。可以将节点坐标系旋转到一个局部柱坐标下。这种情况下，节点坐标系的X方向指向径向，Y方向是周向（theta)。可是当施加theta方向非零位移时，ANSYS总是定义它为一个笛卡尔Y位移而不是一个转动（Y位移不是theta位移）。
单元坐标系
单元坐标系确定材料属性的方向（例如，复合材料的铺层方向）。对后处理也是很有用的，诸如提取梁和壳单元的膜力。单元坐标系的朝向在单元类型的描述中可以找到。
结果坐标系
/Post1通用后处理器中(位移, 应力，支座反力）在结果坐标系中报告，缺省平行于总体笛卡尔坐标系。这意味着缺省情况位移，应力和支座反力按照总体笛卡尔在坐标系表达。无论节点和单元坐标系如何设定。要恢复径向和环向应力，结果坐标系必须旋转到适当的坐标系下。这可以通过菜单路径Post1>Options for output实现。/POST26时间历程后处理器中的结果总是以节点坐标系表达。
显示坐标系
显示坐标系对列表圆柱和球节点坐标非常有用(例如, 径向，周向坐标)。建议不要激活这个坐标系进行显示。屏幕上的坐标系是笛卡尔坐标系。显示坐标系为柱坐标系，圆弧将显示为直线。这可能引起混乱。因此在以非笛卡尔坐标系列表节点坐标之后将显示坐标系恢复到总体笛卡尔坐标系。
(1) 如果你模拟结构体中裂缝扩展过程的模拟，在Ansys中可以用全解耦损伤分析方法来近似模拟裂缝扩展，我曾用Ansys软件中提供的可以定义10,000个材料参数和单元ekill/alive 功能完成了层状路面体中表面裂缝和反射裂缝在变温作用下的扩展过程的模拟。我模拟的过程相对来说比较简单，模拟过程中我们首先要知道裂缝的可能扩展方向，这样在裂缝可能扩展的带内进行网格加密处理，加密到什么程度依据计算的问题来确定。 (2) 如果采用断裂力学理论计算含裂缝结构体的应力强度因子，建模时只需在裂尖通过命令kscon生成奇异单元即可。Ansys模块中存在的断裂力学模块可以计算I、II、III型应力强度因子(线弹性断裂力学)和J积分(弹塑性断裂力学)，在Ansys中verification里面有一个计算I型应力强度因子的例子vm143，参见该例子就可以了。 (3) 如果通过断裂力学模拟裂缝的扩展过程，需要采用动态网格划分，这方面我没有做，通过Ansys的宏命令流应该可以实现。技术参考可参阅文献：杨庆生、杨卫.断裂过程的有限元模拟.计算力学学报，1997,14(4). (4) 我现在做动荷载作用下路面结构体中应力强度因子的分布规律，我是通过位移插值得到不同时间点处的应力强度因子。如果想这样做，可参阅理论参考中关于应力强度因子计算说明。 1. 讨论两种Ansys求极限荷载的方法 （1）力加载 可以通过对应的方法（比如说特征值屈曲）估计结构的极限荷载的大致范围，然后给结构施加一个稍大的荷载，打开自动荷载步二分法进行非线性静力分析，最后计算会因不收敛终止，则倒数第二个子步对应的就是结构的极限荷载；另外，也可以选择弧长法，采用足够的子步（弧长法可以一直分析到极限承载力之后的过程）同样可以从绘制的荷载位移曲线或计算结果中找出结构的极限荷载。 （2）位移加载 给结构施加一个比较大的位移，打开自动荷载步二分法进行非线性分析，保证足够的子步数，这样也可以分析到极限荷载以后，通过绘制荷载位移曲线或查看相应结果文件也可知道结构的极限荷载。 希望众高手讨论一下 （1）弧长法求极限荷载的收敛性问题，如何画到荷载位移曲线的下降段？ （2）位移法求极限荷载的具体步骤？ 2. 需要注意的问题 1. 由于SOLID 65单元本身是基于弥散裂缝模型和最大拉应力开裂判据，因此在很多情况下会因为应力集中而使混凝土提前破坏，从而和试验结果不相吻合，因此，在实际应用过程中应该对单元分划进行有效控制，根据作者经验，当最小单元尺寸大于5cm 时，就可以有效避免应力集中带来的问题； 2. 支座是另一个需要注意的问题。在有限元分析中，很多时候约束是直接加在混凝土节点上，这样很可能在支座位置产生很大的应力集中，从而使支座附近的混凝土突然破坏，造成求解失败。因此，在实际应用过程中，应该适当加大支座附近单元的尺寸或者在支座上加一些弹性垫块，避免支座的应力集中；
ANSYS热应力分析例题 实例1圆简内部热应力分折： 有一无限长圆筒，其核截面结构如图13—1所示，简内壁温度为200℃，外壁温度为20℃，圆筒材料参数如表13．1所示，求圆筒内的温度场、应力场分布。 该问题属于轴对称问题。由于圆筒无限长，忽略圆筒端部的热损失。沿圆筒纵截面取宽度为10M的如图13—2所示的矩形截面作为几何模型。在求解过程中采用间接求解法和直接求解法两种方法进行求解。间接法是先选择热分析单元，对圆筒进行热分析，然后将热分析单元转化为相应的结构单元，对圆筒进行结构分析；直接法是采用热应力藕合单元，对圆筒进行热力藕合分析。 /filname,exercise1-jianjie /title,thermal stresses in a long /prep7 $Et,1,plane55 Keyopt,1,3,1 $Mp,kxx,1,70 Rectng,0.1,0.15,0,0.01 $Lsel,s,,,1,3,2 Lesize, all,,,20 $Lsel,s,,,2,4,2 Lesize,all,,,5 $Amesh,1 $Finish /solu $Antype,static Lsel,s,,,4 $Nsll,s,1 $d,all,temp,200 lsel,s,,,2 $nsll,s,1 $d,all,temp,20 allsel $outpr,basic,all solve $finish /post1 $Set,last /plopts,info,on Plnsol,temp $Finish /prep7 $Etchg,tts Keyopt,1,3,1 $Keyopt,1,6,1 Mp,ex,1,220e9 $Mp,alpx,,1,3e-6 $Mp,prxy,1,0.28 Lsel,s,,,4 $Nsll,s,1 $Cp,8,ux,all Lsel,s,,,2 $Nsll,s,1 $Cp,9,ux,all Allsel $Finish /solu $Antype,static D,all,uy,0 $Ldread,temp,,,,,,rth Allsel $Solve $Finish /post1 /title,radial stress contours Plnsol,s,x /title,axial stress contours Plnsol,s,y /title,circular stress contours Plnsol,s,z /title,equvialent stress contours Plnsol,s,eqv $finish
一学习ANSYS需要认识到的几点 相对于其他应用型软件而言，ANSYS作为大型权威性的有限元分析软件，对提高解决问题的能力是一个全面的锻炼过程，是一门相当难学的软件，因而，要学好ANSYS，对学习者就提出了很高的要求，一方面，需要学习者有比较扎实的力学理论基础，对ANSYS分析结果能有个比较准确的预测和判断，可以说，理论水平的高低在很大程度上决定了ANSYS使用水平；另一方面，需要学习者不断摸索出软件的使用经验不断总结以提高解决问题的效率。在学习ANSYS的方法上，为了让初学者有一个比较好的把握，特提出以下五点建议：（1）将ANSYS的学习紧密与工程力学专业结合起来 毫无疑问，刚开始接触ANSYS时，如果对有限元，单元，节点，形函数等《有限元单元法及程序设计》中的基本概念没有清楚的了解话，那么学ANSYS很长一段时间都会感觉还没入门，只是在僵硬的模仿，即使已经了解了，在学ANSYS之前，也非常有必要先反复看几遍书，加深对有限元单元法及其基本概念的理解。 作为工程力学专业的学生，虽然力学理论知识学了很多，但对许多基本概念的理解许多人基本上是只停留于一个符号的认识上，理论认识不够，更没有太多的感性认识，比如一开始学ANSYS时可能很多人都不知道钢材应输入一个多大的弹性模量是合适的。而在进行有限元数值计算时，需要对相关参数的数值有很清楚的了解，比如材料常数，直接关系到结果的正确性，一定要准确。实际上在学ANSYS时，以前学的很多基本概念和力学理论知识都忘得差不多了，因而遇到有一
定理论难度的问题可能很难下手，特别是对结果的分析，需要用到《材料力学》，《弹性力学》和《塑性力学》里面的知识进行理论上的判断，所以在这种情况下，复习一下《材料力学》，《弹性力学》和《塑性力学》是非常有必要的，加深对基本概念的理解，实际上，适当的复习并不要花很多时间，效果却很明显，不仅能勾起遥远的回忆，加深理解，又能使遇到的问题得到顺利的解决。 在涉及到复杂的非线性问题时（比如接触问题），一方面，不同的问题对应着不同的数值计算方法，求解器的选择直接关系到程序的计算代价和问题是否能顺利解决；另一方面，需要对非线性的求解过程有比较清楚的了解，知道程序的求解是如何实现的。只有这样，才能在程序的求解过程中，对计算的情况做出正确的判断。因此，要能对具体的问题选择什么计算方法做出正确判断以及对计算过程进行适当控制，对《计算方法》里面的知识必须要相当熟悉，将其理解运用到ANSYS的计算过程中来，彼此相互加强理解。要知道ANSYS是基于有限元单元法与现代数值计算方法的发展而逐步发展起来的。因此，在解决非线性问题时，千万别忘了复习一下《计算方法》。此外，对《计算固体力学》也要有所了解（一门非常难学的课），ANSYS对非线性问题处理的理论基础就是基于《计算固体力学》里面所讲到的复杂理论。 作为学工程力学的学生，提高建模能力是非常急需加强的一个方面。在做偏向于理论的分析时，可能对建模能力要求不是很高，但对于实际的工程问题，有限元模型的建立可以说是一个最重要的问题，而后
ANSYS建模分析 报 告 书 课题名称ANSYS建模分析姓名 学号 院系 专业 指导老师
问题描述 在ANSYS中建立如图一所示的支承图，假定平面支架沿厚度方向受力均匀，支承架厚度为3mm。支承架由钢制成，钢的弹性模量为200Gpa，泊松比为0.3。支承架左侧边被固定，沿支承架顶面施加均匀载荷，载荷与支架共平面，载荷大小为2000N/m。要求：绘制变形图，节点位移，分析支架的主应力与等效应力。 图1 GUI操作步骤 1、定义工作文件名和工作标题 （1）定义工作文件名：执行Utility Menu> Jobname命令，在弹出【Change Jobname】对话框中输入“xuhao144139240174”。选择【New log and error files】复选框，单击OK按钮。 （2）定义工作标题：执行Utility Menu> Title命令，在弹出【Change Title】对话框中输入“This is analysis made by “xh144139240174”，单击OK按钮。 （3）重新显示：执行Utility Menu>Plot>Replot命令。 （4）关闭三角坐标符号：执行Utility Menu>PlotCtrls>Window Options命令，弹出【Window Options】对话框。在【Location of triad】下拉列表框中选择“Not Shown”选项，单击OK按钮。 2、定义单元类型和材料属性
（1）选择单元类型：执行Main Menu>Preprocessor>Element Type>Add/Edit/Delete命令，弹出【Element Type】对话框。单击Add...按钮，弹出【Library of Element Types】对话框。选择“Structural Solid”和“Quad 8node 82”选项，单击OK按钮，然后单击Close按钮。 （2）设置材料属性：执行Main Menu>Preprocessor>Material Props>Material Models命令，弹出【Define Material Models Behavior】窗口。双击【Material Model Available】列表框中的“Structural\Linear\Elastic\Isotropic”选项，弹出【Linear Isotropic Material Properties for Material Number1】对话框。在【EX】和【PRXY】文本框中分别输入“2e11”及“0.3”。单击OK按钮，然后执行Material>Exit命令，完成材料属性的设置。 3、创建几何模型 （1）生成两个矩形面：执行Main Menu>Preprocessor>Modeling>Create>Areas>Rectangle>By Dimensions命令，弹出【Create Rectangle by Dimensions】对话框。输入第一个矩形的坐标数值：“X1=0,X2=-0.25,Y1=0.025,Y2=-0.025”,单击Apply按钮。输入第二个矩形的坐标数值：“X1=-0.2,X2=-0.25,Y1=0.025,Y2=-0.225”，单击OK按钮关闭该对话框。 （2）生成两个半圆：执行Main Menu>Preprocessor>Modeling>Create>Areas>Circle>Solid Circle
Ansys有限元分析实例[教学] 有限元分析案例:打点喷枪模组(用于手机平板电脑等电子元件粘接)，该产品主要是使用压缩空气推动模组内的顶针作高频上下往复运动，从而将高粘度的胶水从喷嘴中打出(喷嘴尺寸,0.007”)。顶针是这个产品中的核心零件，设计使用材料是:AISI 4140 最高工作频率是160HZ(一个周期中3ms开3ms关)，压缩空气压力3-8bar, 直接作用在顶针活塞面上，用Ansys仿真模拟分析零件的强度是否符合要求。 1. 零件外形设计图:
2. 简化模型特征后在Ansys14.0 中完成有限元几何模型创建:
3. 选择有限元实体单元并设定，单元类型是SOILD185，由于几何建模时使用的长度单位是mm, Ansys采用单位是长度:mm 压强: 3Mpa 密度:Ton/M。根据题目中的材料特性设置该计算模型使用的材料属性:杨氏模量 2.1E5; 泊松比:0.29; 4. 几何模型进行切割分成可以进行六面体网格划分的规则几何形状后对各个实体进行六面体网格划分，网格结果: 5. 依据使用工况条件要求对有限元单元元素施加约束和作用载荷:
说明: 约束在顶针底端球面位移全约束; 分别模拟当滑块顶断面分别以8Bar，5Bar，4Bar和3Bar时分析顶针的内应力分布，根据计算结果确定该产品允许最大工作压力范围。 6. 分析结果及讨论: 当压缩空气压力是8Bar时: 当压缩空气压力是5Bar时:
当压缩空气压力是4Bar时: 结论: 通过比较在不同压力载荷下最大内应力的变化发现，顶针工作在8Bar时最大应力达到250Mpa，考虑到零件是在160HZ高频率在做往返运动，疲劳寿命要求50百万次以上，因此采用允许其最大工作压力在5Mpa，此时内应力为156Mpa，按线性累积损伤理论[3 ]进行疲劳寿命L-N疲劳计算，进一部验证产品的设计寿命和可靠性。
1. 讨论两种Ansys求极限荷载的方法 （1）力加载 可以通过对应的方法（比如说特征值屈曲）估计结构的极限荷载的大致范围，然后给结构施加一个稍大的荷载，打开自动荷载步二分法进行非线性静力分析，最后计算会因不收敛终止，则倒数第二个子步对应的就是结构的极限荷载；另外，也可以选择弧长法，采用足够的子步（弧长法可以一直分析到极限承载力之后的过程）同样可以从绘制的荷载位移曲线或计算结果中找出结构的极限荷载。 （2）位移加载 给结构施加一个比较大的位移，打开自动荷载步二分法进行非线性分析，保证足够的子步数，这样也可以分析到极限荷载以后，通过绘制荷载位移曲线或查看相应结果文件也可知道结构的极限荷载。 希望众高手讨论一下 （1）弧长法求极限荷载的收敛性问题，如何画到荷载位移曲线的下降段？ （2）位移法求极限荷载的具体步骤？ 2. 需要注意的问题 1. 由于SOLID 65单元本身是基于弥散裂缝模型和最大拉应力开裂判据，因此在很多情况下会因为应力集中而使混凝土提前破坏，从而和试验结果不相吻合，因此，在实际应用过程中应该对单元分划进行有效控制，根据作者经验，当最小单元尺寸大于5cm 时，就可以有效避免应力集中带来的问题； 2. 支座是另一个需要注意的问题。在有限元分析中，很多时候约束是直接加在混凝土节点上，这样很可能在支座位置产生很大的应力集中，从而使支座附近的混凝土突然破坏，造成求解失败。因此，在实际应用过程中，应该适当加大支座附近单元的尺寸或者在支座上加一些弹性垫块，避免支座的应力集中； 3. 六面体的SOLID 65 单元一般比四面体的单元计算要稳定且收敛性好，因此，只要条件允许，应该尽量使用六面体单元； 4. 正确选择收敛标准，一般位移控制加载最好用位移的无穷范数控制收敛，而用力控制加载时可以用残余力的二范数控制收敛。在裂缝刚刚出现和接近破坏的阶段，可以适当放松收敛标准，保证计算的连续性； 3. 关于下降段的问题 1）在实际混凝土中都有下降段，但是在计算的时候要特别小心下降段的问题。 2）下降段很容易导致计算不收敛，有时为了计算的收敛要避免设置下降段，采用rush模型。 3）利用最大压应变准则来判断混凝土是否破坏。 4. Solid65单元中的破坏准则 1）采用Willam&Warnke五参数破坏准则 2）需要参数： 单轴抗拉强度，单轴，双轴抗压强度，围压压力，在围压作用下双轴，单轴抗压强度 5. 近来我对混凝土单元进行了一点思考，有一些想法，贴在下面，共同探讨： 1）分析混凝土结构，选择合理的材料特性是建立模型的关键，所以有必要弄清混凝土的材料特性。混凝土是脆性材料,并具有不同的拉伸和压缩特性。典型混凝土的抗拉强度只有抗压强度的8％-15%。 在ANSYS中，对于混凝土单元，材料特性ANSYS要求输入以下数据（为了清楚起见，我将几个系数均译为了中文）：弹性模量、泊松比、张开与闭合滑移面的剪切强度缩减系数、抗拉与抗压强度、极限双轴抗压强度、周围静水应力状态、静水应力状态下单轴与双轴压缩的
《有限元基础教程》作业二：平面薄板的有限元分析 班级：机自101202班 姓名：韩晓峰 学号：201012030210 一．问题描述： P P h 1mm R1mm 10m m 10mm 条件：上图所示为一个承受拉伸的正方形板，长度和宽度均为10mm ，厚度为h 为1mm ，中心圆的半径R 为1mm 。已知材料属性为弹性模量E=1MPa ，泊松比为0.3，拉伸的均布载荷q ＝ 1N/mm 2。根据平板结构的对称性，只需分析其中的二分之一即可，简化模型如上右图所示。 二．求解过程： 1 进入ANSYS 程序 →ANSYS 10.0→ANSYS Product Launcher →File management →input job name: ZY2→Run 2设置计算类型 ANSYS Main Menu: Preferences →select Structural → OK 3选择单元类型 ANSYS Main Menu: Preprocessor →Element Type →Add/Edit/Delete →Add →select Solid Quad 4node 42 →OK → Options… →select K3: Plane Strs w/thk →OK →Close 4定义材料参数 ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX: 1e6, PRXY:0.3 → OK 5定义实常数以及确定平面问题的厚度 ANSYS Main Menu: Preprocessor →Real Constants …→Add/Edit/Delete →Add →Type 1→OK →Real Constant Set No.1,THK:1→OK →Close 6生成几何模型 a 生成平面方板 ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Rectangle →By 2 Corners →WP X:0,WP Y:0,Width:5,Height:5→OK b 生成圆孔平面 ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Circle →Solid Circle →WPX=0,WPY=0,RADIUS=1→OK b 生成带孔板 ANSYS Main Menu: Preprocessor →Modeling →Operate →Booleans → Subtract →Areas →点击area1→OK →点击area2→OK 7 网格划分 ANSYS Main Menu: Preprocessor →Meshing →Mesh Tool →(Size Controls) Global: Set →SIZE: 0.5 →OK →iMesh →Pick All → Close 8 模型施加约束
5、ANSYS输出mnf文件 模型单位要统一，最好都适用国际单位米制的，那么弹性模量、密度也要统一单位。然后进行单元添加：solid45、beam4、mass21给beam4设置实常数（real constant）：基本都是1e-12（米制单位，毫米要相应改变） 给mass21设置实常数（real constant）：基本都是1e-12（米制单位，毫米要相应改变） 添加材料设置：包括两种材料，一种是实体需要的材料，即为应该模型材料。 一种就是需要刚度大但是质量轻的材料，一般用的是密度为1e-12，弹性模量比模型实体的高出5个数量级（这个数值对能否导成功有直接影响，可以进行试算，用高5个数量级保证了稳定输出）。 在attachpoint铰链位置添加两个keypoint，然后用mass21去划分网格。可以得到node 1、node2，然后对模型整体用solid45划分。现在要把这两个孔刚化，就需要用到刚性梁单元。 用beam4单元连接孔上每一个节点与孔中心节点（需要成为attachpoint的点）。 6、ansys中的add、glue、overlap的区别及联系 1、相加（add）：相加是指对所有图元进行叠加，包含原是个图元的所有部分，生成一个新图元，各个原始图元的公共边界将被清除，形成一个单一的整体。在ansys的面相加中只能对共面的图元进行操作.
对两个已经存在的面进行相加操作 命令：aadd,na1,na2,na3,na4,na5,na6,na7,na8,na9 2)对两个已经存在的体进行相加操作命令： vadd,nv1,nv2,nv3,nv4,nv5,nv6,nv7,nv8,nv9 3）对两条已经存在的线进行操作 命令：lcomb,nl1,nl2,keep keep表示保留进行相加操作的图元，deleted表示进行相加操作后删除原始图元。 2、搭接（overlap）：搭接食指将分离的同阶图元转变为一个连续体，其中图元的所有重叠区域将独立成为一个图元。搭接与相加操作类似，但相加操作是由几个图元生成一个图元整体，而搭接则是由几个图元生成更多的图元，相交的部分则被分离出来。 1）、线和线之间进行搭接操作 命令：lovlap,nl1,nl2,nl3,nl4,nl5,nl6,nl7,nl8,nl9 2）、面和面之间进行搭接操作 命令：aovlap,na1,na2,na3,na4,na5,na6,na7,na8,na9 3）、体和体之间进行搭接操作 命令：vovlap,nv1,nv2,nv3,nv4,nv5,nv6,nv7,nv8,nv9 3、粘结（glue）粘结操作是将多个图元组合成一个连续体，图元之间仅在公共边界处相连，其公共边界的维数低于原始图元一维。粘结操作与加操作类似，但不同的是这些图元之间仍然相互独立，只是在边界上连接。粘结操作通常还与搭接操作配合使用。
有 限 元 分 析 作 业 作业名称 输气管道有限元建模分析 姓 名 陈腾飞 学 号 3070611062 班 级 07机制（2）班 宁波理工学院
题目描述： 输气管道的有限元建模与分析 计算分析模型如图1所示 承受内压：1.0e8 Pa R1=0.3 R2=0.5 管道材料参数：弹性模量E=200Gpa；泊松比v=0.26。 图1受均匀内压的输气管道计算分析模型（截面图） 题目分析： 由于管道沿长度方向的尺寸远远大于管道的直径，在计算过程中忽略管道的断面效应，认为在其方向上无应变产生。然后根据结构的对称性，只要分析其中1/4即可。此外，需注意分析过程中的单位统一。 操作步骤 1.定义工作文件名和工作标题 1.定义工作文件名。执行Utility Menu-File→Chang Jobname-3070611062，单击OK按钮。 2.定义工作标题。执行Utility Menu-File→Change Tile-chentengfei3070611062，单击OK按钮。 3.更改目录。执行Utility Menu-File→change the working directory –D/chen 2.定义单元类型和材料属性 1.设置计算类型 ANSYS Main Menu: Preferences →select Structural →OK
2.选择单元类型。执行ANSYS Main Menu→Preprocessor →Element Type→Add/Edit/Delete →Add →select Solid Quad 8node 82 →apply Add/Edit/Delete →Add →select Solid Brick 8node 185 →OK Options…→select K3: Plane strain →OK→Close如图2所示，选择OK接受单元类型并关闭对话框。 图2 3.设置材料属性。执行Main Menu→Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic，在EX框中输入2e11，在PRXY框中输入0.26，如图3所示，选择OK并关闭对话框。 图3 3.创建几何模型 1. 选择ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints →In Active CS →依次输入四个点的坐标：input:1(0.3,0),2(0.5,0),3(0,0.5),4(0,0.3) →OK
《大型结构分析软件的应用及开发》 学习报告 学院：建筑工程学院 专业班级：工程力学141 姓名：付贤凯 指导老师：姚激 学号：201411012111
1．模型介绍 如下图所示的一桁架结构，受一集中力大小为800N的作用，杆件的弹性模量为200GPa，泊松比为0.3。杆件的截面为正方形达长为1m，横截面面积为1m2。现求它的变形图与轴力图。 图1 桁架模型与受力简图（单位：mm） 2．建模与划分网格 利用大型有限元软件ANSYS，采用Link，2Dspar 1的单元进行模拟，通过网格的划分得到如图2所示的有限元模型。 图2 有限元模型
结合有限元模型中的约束条件为左侧在X与Y方向铰支固定，荷载条件为最右侧处施加向下的集中力P=800N。施加约束与荷载后的几何模型如图4所示。 图3 施加荷载与约束的几何模型 3．位移与轴力图 因在Y方向受力，所以主要做Y方向的位移图，又因为杆件在轴线方向有变形，故在X 方向仍有一定的位移。则图5为变形前后的板件形状。图6为模型沿Y方向的位移图，图7为模型沿X方向的位移图，图8为模型的总位移图。 图4 桁架变形前后形状图
图5 Y方向位移图 图6 X方向位移图
图7总位移图 分析所有的位移图可以看出从以看出左端变形最小，为零，右端变形最大。从总位移图可以看出最大的位移在左下点处，大小为0.164×10?5m。从X方向位移图可以看出，左下点处在X方向位移最大为0.36×10?6。从Y方向位移图可以看出最大位移在左下点处为0.164×10?5。都符合实际情况，图9为模型的轴力图。 图8 轴力图
学习ansys的一些心得 学习ansys的一些心得（送给初学者和没有盟币的兄弟） 1 做了布尔运算后要重画图形（删除实体）时：需拾取Utility Menu>Plot>Replot 2 标点的输入是在英文状态下，―,‖。 3 线段中点的建立：Modling>Creat>Keypoints>Fill between kps 4 还不会环形阵列。 5 所谓杆系结构指的是长度远远大于其他方向尺寸（10：1）的构件组成的结构，如连续梁，桁架，钢架等。 6 静力学分析的结果包括结构的位移，应变，应力和反作用力等，一般是使用POST1处理（普通后处理器）和查看这些结果。 7 干系结构的静力学分析—平面桁架的建模，用NODE（节点），ELEMENT(元素)创建。复杂体积的建模一般用KPS(关键点)，LINE（Straight line—直线），再生成面，再生成体。 8 如果输入的数据单位是国际单位制单位，则输出的数据单位也是国际制单位。 9 创建正六边形：Creat>Areas>Polygon>Hexagon.指定中心和半径。 10 由面沿线挤出体：Modling>Operate>Extrude>Areas>Along Lines. 11 Ansys中没有Undo命令.需及时保存数据库文件. Def Shape Only:只显示变形图.Def + Undeformed:显示未变形的图.Def + Udef egde:显示未变形的图形的边界. 13 用等高线显示：Plot Results>Contour Plot>Nodal Solu.
14 模态分析用于分析结构的振动特性，即确定结构的固有频率和振型，它也是谐响应分析，瞬态动力学分析以及谱分析等其他动力学分析的基础。 15 Ansys的模态分析是线型分析。任何非线型分析，例如，塑性，接触单元等，即使被定义了也将被忽略。 16 平面桁架：Beam(2D elastic 3) 厚壁圆筒：Solid(8 node 13)>Options(K3—Plane strain) 17 一般材料的弹性模量（EX）：2e11.泊松比(PRXY)：0.3.密度：7800 18 做完静力学分析后，再做模态分析时，要再次求解，同时预应力效果也应该打开（PSTRES,on）.可以在命令行中输入：pstres,on 也可以用菜单路径：Solution>Analysis Type>Analysis Options. 19 弹簧阻尼器单元：Combination-Spring damper 14. 20 接触问题属于状态非线性问题，是一种高度非线性行为，需要较多的计算资源。接触问题有两个基本类型：刚体-柔体的接触，柔体-柔体的接触（许多金属成型的接触问题）。在刚体-柔体的接触问题中，有的接触面与它接触的变形体相比，有较大的刚度而被当做刚体。而柔体-柔体的接触，是一种更普遍的类型，此时两个接触体具有近似的刚度，都为变形体。 21 1 点-点接触：过盈装配问题是用点点接触单元模拟面面接触的典型例子。 2 点-面接触：不必预先知道准确的接触位置，接触面之间也不需要保持一致的网格，并且允许有较大的变形和相对滑动。典型实例：模拟插头插入插座里。 3 面-面接触：刚性面作为目标面，柔性面作为接触面。 22 打开自动时间步长：Solution>Load Step Opts>Time Frequenc>Time And Substps.
ansys有限元分析工程实例大作业
————————————————————————————————作者：————————————————————————————————日期：
辽宁工程技术大学 有限元软件工程实例分析 题目基于ANSYS钢桁架桥的静力分析专业班级建工研16-1班（结构工程）学号 471620445 姓名 日期 2017年4月15日
基于ANSYS钢桁架桥的静力分析 摘要：本文采用ANSYS分析程序，对下承式钢桁架桥进行了有限元建模;对桁架桥进行了静力分析，作出了桁架桥在静载下的结构变形图、位移云图、以及各个节点处的结构内力图（轴力图、弯矩图、剪切力图），找出了结构的危险截面。 关键词：ANSYS；钢桁架桥；静力分析；结构分析。 引言：随着现代交通运输的快速发展，桥梁兴建的规模在不断的扩大，尤其是现代铁路行业的快速发展更加促进了铁路桥梁的建设，一些新建的高速铁路桥梁可以达到四线甚至是六线，由于桥面和桥身的材料不同导致其受力情况变得复杂，这就需要桥梁需要有足够的承载力，足够的竖向侧向和扭转刚度，同时还应具有良好的稳定性以及较高的减震降噪性，因此对其应用计算机和求解软件快速进行力学分析了解其受力特性具有重要的意义。 1、工程简介 某一下承式简支钢桁架桥由型钢组成，顶梁及侧梁，桥身弦杆，底梁分别采用3种不同型号的型钢，结构参数见表1，材料属性见表2。桥长32米，桥高5.5米，桥身由8段桁架组成，每个节段4米。该桥梁可以通行卡车，若只考虑卡车位于桥梁中间位置，假设卡车的质量为4000kg，若取一半的模型，可以将卡车对桥梁的作用力简化为P1，P2，和P3，其中P1=P3=5000N，P2=10000N,见图2,钢桥的形式见图1，其结构简图见图3。
有限元与CAE分析报告 专业： 班级： 学号： 姓名： 指导教师： 2016年 1 月 2 日
简支梁的静力分析 一、问题提出 长3m的工字型梁两端铰接中间1.5m位置处受到6KN的载荷作用，材料弹性模量E=200e9，泊松比0.28，密度7850kg/㎡ 二、建立模型 1.定义单元类型 依次单击Main Menu→Preprocessor→Elementtype→Add/Edit/Delete，出现对话框如图，单击“Add”，出现一个“Library of Element Type”对话框，在“Library of Element Type”左面的列表栏中选择“Structural Beam”，在右面的列表栏中选择3 node 189,单击“OK”。
2设置材料属性 依次单击Main Menu→Preprocessor→MaterialProps>Material Modes，出现“Define Material ModelBehavior”对话框，在“Material Model Available”下面的对话框中，双击打开“Structural→Linear→Elastic→Isotropic”,出现对话框，输入弹性模量EX=2E+011,PRXY=0.28，单击“OK”。 依次单击Main Menu→Preprocessor→MaterialProps>Material Modes，出现“Define Material ModelBehavior”对话框，在“Material Model Available”下面的对话框中，双击打开“Structural→Density”弹出对话框，输入DENS为7850 3.创建几何模型 1)设定梁的截面尺寸
ANSYS使用心得体会 本次结构力学课程设计是学习使用ANSYS软件对框架结构内力进行计算，在未学习该软件前，对于此类问题，通常会采用力矩分配法来进行计算，计算过程繁复，计算量大。导致过程缓慢。 通过对ANSYS软件的学习和了解，知道了它的一些明显的优点。 相对于其他应用型软件而言，ANSYS作为大型权威性的有限元分析软件，对提高解决问题的能力是一个全面的锻炼过程，是一门相当难学的软件，因而，要学好ANSYS，对我们提出了很高的要求，一方面，需要我们有比较扎实的力学理论基础，对ANSYS分析结果能有个比较准确的预测和判断，可以说，理论水平的高低在很大程度上决定了ANSYS使用水平；另一方面，需要我们不断摸索出软件的使用经验不断总结以提高解决问题的效率。 刚开始接触ANSYS时，没有限元，单元，节点，形函数等的基本概念没有清楚的了解话，会感觉还没入门，只是在僵硬的模仿，即使已经了解了，必要先反复看几遍书，加深对有限元单元法及其基本概念的理解。 ANSYS在对结构力学的静力学分析非常方便，用来求解外载荷引起的位移、应力和力。静力分析很适合求解惯性和对结构的影响并不显著的问题。ANSYS 程序中的静力分析不仅可以进行线性分析，而且也可以进行非线性分析，如塑性、膨胀、大变形、大应变及接触分析。 但是学习的过程是充满烦恼和惊喜的，因为总是会碰到许多的新问题，需要较好的耐心去解决这些问题，这是在学习过程中遇到的最大的难题。然而，在解决问题之后，就会有恍然大悟的喜悦，可以说是痛苦和快乐并存的。所以对于初学者，缺乏经验是非常难的。必须保持良好的心态，对于不断出现的ERROR提示要坚定自己的信心，坚信自己可以解决这些问题。所有困难都会迎刃而解。 本次的学习让我认识到了提高建模能力是非常急需加强的一个方面。在做偏向于理论的分析时，可能对建模能力要求不是很高，但对于实际的工程问题，有限元模型的建立可以说是一个最重要的问题，而后面的工作变得相对简单。建模能力的提高，需要掌握好的建模思想和技巧。 ANSYS软件是一款在建模等方面非常实用的软件，本次的学习我其实并没有完全熟练地掌握它的应用，以后还要加强对它的学习，相信在以后的学习和工作中会带来巨大的便利。
ansys有限元分析大作业
有限元大作业 设计题目: 单车的设计及ansys有限元分析 专业班级: 姓名: 学号: 指导老师: 完成日期: 2016.11.23
单车的设计及ansys模拟分析 一、单车实体设计与建模 1、总体设计 单车的总体设计三维图如下，采用pro-e进行实体建模。 在建模时修改proe默认单位为国际主单位(米千克秒 mks) Proe》文件》属性》修改
2、车架 车架是构成单车的基体，联接着单车的其余各个部件并承受骑者的体重及单车在行驶时经受各种震动和冲击力量，因此除了强度以外还应有足够的刚度，这是为了在各种行驶条件下，使固定在车架上的各机构的相对位置应保持不变，充分发挥各部位的功能。车架分为前部和后部，前部为转向部分，后部为驱动部分，由于受力较大，所有要对后半部分进行加固。
二、单车有限元模型 1、材料的选择 单车的车身选用铝合金（6061-T6）T6标志表示经过热处理、时效。 其属性如下： 弹性模量：） .6+ 90E （2 N/m 10 泊松比：0.33 质量密度：） 3 2.70E+ N/m （2 抗剪模量：） 60E .2+ N/m （2 10 屈服强度：） .2+ （2 75E 8 N/m 2、单车模型的简化 为了方便单车的模拟分析，提高电脑的运算
效率，可对单车进行初步的简化；单车受到的力的主要由车架承受，因此必须保证车架能够有足够的强度、刚度，抗振的能力，故分析的时候主要对车架进行分析。简化后的车架如下图所示。 3、单元体的选择 单车车架为实体故定义车架的单元类型为实体单元（solid）。查资料可以知道3D实体常用结构实体单元有下表。 单元名称说明 Solid45 三维结构实体单元，单元由8个节点定义，具有塑性、蠕变、应力刚化、 大变形、大应变功能，其高阶单元是 solid95
有限元分析 一个厚度为20mm的带孔矩形板受平面内张力，如下图所示。左边固定，右边受载荷p=20N/mm作用，求其变形情况 P 一个典型的ANSYS分析过程可分为以下6个步骤： ①定义参数 ②创建几何模型 ③划分网格 ④加载数据 ⑤求解 ⑥结果分析 1定义参数 1.1指定工程名和分析标题 (1)启动ANSYS软件，选择File→Change Jobname命令，弹出如图所示的[Change Jobname]对话框。 (2)在[Enter new jobname]文本框中输入“plane”，同时把[New log and error files]中的复选框选为Yes，单击确定 (3)选择File→Change Title菜单命令，弹出如图所示的[Change Title]对话框。 (4)在[Enter new title]文本框中输入“2D Plane Stress Bracket”，单击确定。 1.2定义单位
在ANSYS软件操作主界面的输入窗口中输入“/UNIT,SI” 1.3定义单元类型 (1)选择Main Menu→Preprocessor→Element Type→Add/Edit/Delete命令，弹出如图所示[Element Types]对话框。 (2)单击[Element Types]对话框中的[Add]按钮，在弹出的如下所示[Library of Element Types]对话框。 (3)选择左边文本框中的[Solid]选项，右边文本框中的[8node 82]选项，单击确定，。 (4)返回[Element Types]对话框，如下所示 (5)单击[Options]按钮，弹出如下所示[PLANE82 element type options]对话框。
关于网格的几个误区 尽管当前出现了不少使用无网格方法的FEA及CFD代码，但是网格划分依然是大多数CAE工作者们最重要的工作任务，对于高质量网格生成的重要性怎么强调都不过分。 但是如何生成高质量的或更精细的网格呢？查看网格生成软件所输出的网格质量报告是最基本的方式，使用者还需要对网格是否适用于自己的物理问题做出自己的判断。 不幸的是，使用者对于“好网格”存在很多的误区。如今已经很难在工程学科中找到关于网格划分方面的课程，数值算法在大多数工程学科中成了选修课程。因此，新生代CAE使用者对于网格在CAE系统中的工作机理方面的欠缺也不足为怪了。这里有5个最主要的误区： 误区1：好的网格必须与CAD模型吻合 越来越多的CAE使用者来自于原来的设计人员，他们在CAD方面受到了良好的培训，因此他们倾向于CAE模型体现所有的几何细节特征，他们认为更多的细节意味着计算结果能够更加贴近于真实情况。 然而这种观点是不正确的，好的网格是能够解决物理问题，而不是顺从CAD模型。 CAE仿真的目的是为了获取物理量：应力、应变、位移、速度、压力等。CAD 模型应当是从物理对象中提取的。大量与物理问题不相干的或对于仿真模型影响较小的细节特征在建立CAD模型之前就应当进行简化。因此，了解所仿真的系统中的物理细节是最基本的工作任务。好的网格应当简化CAD模型并且网格节点是基于物理模型进行布置。
这意味着：只有在充分了解所要仿真的物理系统前提下才可能划分出好的网格。误区2：好的网格一直都是好的 我们经常看到CAE使用者花费大量的心血在改变网格尺寸、拆解几何及简化几何上，以期能够获得高质量的网格。他们仔细的检查网格生成软件输出的网格质量报告，这是很有必要的。但是这事儿做得太过也不一定好，因为好的网格也不一定永远都好，网格的好与坏，还取决于要仿真的物理问题。 例如，你生成了一套非常好的网格，其能够很好的捕捉机翼的绕流，能够很精确的计算各种力。但是当你将流动攻角从0°调整到45°，试问这网格还是好的网格吗？很可能不是了。 好的网格总是与物理问题相关。当你改变边界条件、改变载荷、改变分析类型、改变流动条件，好的网格也可能变成坏网格。 误区3：六面体网格总比四面体网格好 很多老的书籍会说六面体（四边形）网格要比四面体（三角形）网格要好，同时告诉你说引入四面体（三角形）网格会造成很大的数值误差。一些情况下这种观点是正确的，特别是15~20年前。 历史上，人们热衷于六面体网格，主要有以下原因：1）在当时，CFD求解器仅能使用结构网格；2）计算条件不允许使用大量网格，为了节省内存和节省时间；3）非结构网格还不成熟。 在过去的几十年里，大部分商用FEA及CFD求解器技术获得的了极大的发展，对于绝大多数问题，利用六面体网格及四面体网格都能获得相同的计算结果。当然，四面体网格通常需要更多的计算资源，但是其能在网格生成阶段为使用者节省大量的时间。对于大多数工程问题，六面体网格在计算精度方面的优势已经不