中国石油大学(北京)《高等工程力学》课件-高等工程力学(5)---弹性力学(3-1)_图文_百度文库
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
中国石油大学(北京)《高等工程力学》课件-高等工程力学(5)---弹性力学(3-1)
上传于||文档简介
&&中​国​石​油​大​学​(​北​京​)​《​高​等​工​程​力​学​》​课​件
阅读已结束，如果下载本文需要使用
想免费下载本文？
下载文档到电脑，查找使用更方便
还剩33页未读，继续阅读
你可能喜欢&& 查看话题
金属板材的弯曲试验中：三点弯、四点弯和纯弯曲有什么区别么？
是不是说三点弯、或者四点弯都是属于纯弯曲实验的一种呢？
只是说四点弯比三点弯更接近材料的真实性能？
三点弯曲不是纯弯曲，而四点弯曲在两个施力点之间的部分属純弯曲。纯弯曲由于弯矩恒定，且不受剪应力的影响，用来测试材料的抗弯性能应该更精确。 : Originally posted by peterflyer at
三点弯曲不是纯弯曲，而四点弯曲在两个施力点之间的部分属純弯曲。纯弯曲由于弯矩恒定，且不受剪应力的影响，用来测试材料的抗弯性能应该更精确。 大侠，能再讲讲三点弯和四点弯的区别么？三点弯是用来测什么的呢？ : Originally posted by spike_2014 at
大侠，能再讲讲三点弯和四点弯的区别么？三点弯是用来测什么的呢？
... 三点弯曲与四点弯曲都用来测材料尤其是脆性材料的抗弯强度和断裂韧性。三点弯曲没四点弯曲精确，但因为比四点弯曲简单，故经常得到应用。[转载]ABAQUS常见问题汇总&-&2.0版（18-24）
．边界条件
请教：第二步施加位移载荷
请教高手：SimWe仿真论坛(Y0P( d6|.O@(K9e
利用abaqus计算时，定义了两个分析步，第一步是一个热残余应力分析，第二步是施加的是位移载荷为10，结果由于第一步产生了残余应变，导致第二步加的位移载荷总是10减去第一步产生的残余应变。即第二步的位移载荷总是加不到10。这是什么原因？如何改进？
方法1： 用线性摄动分析步，见 Getting Started with ABAQUS 11.2.2 Specifying
loads in linear perturbation steps
方法2：用一般的分析步，先只作热残余应力分析，得到位移U0, 再用重启动分析加位移载荷，其大小为 U0 +
希望的位移增量。
边界约束施加问题
如图1所示，有一构件表面呈现周期性的起伏形貌，左右两侧尺寸可认为无限长。现取构件一部分进行有限元建模分析，如图2所示，那么图2左右两侧竖直边的边界约束该如何给呢？请各位高手指点一二，非常感谢。
感觉应该再多取这个的另一半吧,才是周期性的.完了左右两边定义位移协调应该可以.EQUATION可以定义位移协调.
同意 步行去远方 的说法。 另外，只有当载荷之类的其他参数也是周期性重复出现时
（最终得到的位移结果也是周期性的），才能这样只取一部分来模拟。
问：请教有关刚体的问题
对一薄板进行拉伸,给俩边TIE刚体时遇到问题,显示算法,位移控制
}5M@仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM不加刚体时,我就是给板的俩边各一个位移
^3[加刚体后,因为板和刚体tie连接,我把位移加在参考点上,但是对于initial步的定义我蒙了,
f*|{J}z;IW& & &
1:initial步不设边界条件的话,结果是侧面看板翘了,最大应力也没到屈服
:[;N*O+~;j)z
|仿真|设计|有限元|虚拟仪器2:把板设为set后,在initial步对其全约束然后在step中解除这个约束,给参考点加质量,其它条件同上,结果同上,如果不给参考点加质量的话,云图大概趋势一致,侧面看板是拱形4nt@-T[1]d)y)j7E
3:如果不管板,在 initial 中对俩个参考点全约束的话,在step中给参考点位移,结果是云图发生很大的变化
m不明白为什么会有这么大的差别,板和刚体tie&&的话,对板或对参考点施加约束不应该一样么,还有就是刚体到底需不需加质量,请高手解惑,谢谢
k/uhe+}$N(J
dy仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM期待各位大侠的指导,谢谢
"给参考点位移,侧面看板翘了": 参考点不在板的平面内，相当于板受到一个弯矩，所以会翘起来。
边界条件定义在 initial 或 step 1 中，作用是一样的。
板和刚体tie连接，板在tie的位置就不能再变形，所以和没有tie的时候结果是不一样的。
如果刚体很重，它运动的惯性对变形影响很大，就要定义刚体质量。
问题的关键不是在initial还是step 1中定义边界条件，而是在step 1时的边界条件是什么状态。例如：
- 如果在initial中约束了U2, 而在step 1中没有对U2的边界条件作修改， 则在step 1中仍然 U2=0；
- 如果在initial中约束了U2, 而在step 1中令U2=100，则在step 1中 U2=100。
你的模型1： 参考点上只定义了U1=4.1,
在其他5个自由度上可以自由运动（没有定义边界条件），所以板会翘起来（参考点在Y方向运动了）。
你的模型2： 在initial约束了板，在step 1又把这个约束deactivate了,
等于这个约束没起作用，所以结果和模型1一样。
你的模型3： 参考点上定义了U1=4.1, 其他5个自由度都为0，所以板不会翘（参考点在Y方向不能运动）。
前两个模型和第3个模型中参考点的运动不同，板受到的弯矩也就不同，所以云纹图不同。
如果实际问题中受力点位于板所在的平面内，板就根本不会受到弯矩，见附件。
问：让一个部件在10分钟内用定义边界条件的displacement向下移动一定距离，这能实现么？如果能，amplitude应如何定义？
静态还是动态分析？
创建边界条件时，把step由默认的initial改为某个step，然后输入位移值就行了，不用定义amplitude,
默认的amplitude就是从0开始线性增大。
．定义Set / 选取模型的局部
问：在part中划分单元，建立set，和在assembel中主要有什么区别
ABAQUS/CAE User's Manual
13.3.2 What is the difference between a dependent and an
independent part instance
问：关于创建几何/元素等集合
在建立好模型之后 需要建立元素和节点的 几何：按照tool－－set－－creat步骤点击之后
结果只能有”geomerty“& 一种选择， 请问
如何才能进行选择“node”& “element” ？&
选”geomerty“&
即可，ABAQUS/CAE会自动为所选的区域分别创建一个node集合和一个element集合。
问：请问，我在定义一个弹簧时，首先建立了两点，后来选择连接2点创建弹簧，但是cae中这两点显示是重合在一起的，请问该如何选择啊
在窗口底部提示区中出现“Select the first
point...”时，点击其后的图标，在弹出的Options对话框中，点击图标 （Select the Entity Closest
to the Screen），使其变灰。
这时再点击模型时，窗口底部会出现按钮"next","previous", 就可以选择，然后点ok确认。
问：请问spring2的弹簧刚度（受压而不受拉）的设定在CAE中怎么实现
CAE中恐怕不行，可以先建个线性的，再手工修改INP文件，例如:
*SPRING, ELSET=aaa, NONLINEAR
-90000 , -0.396
-18000 , -0.1315
-400 , -0.0048
(*SPRING 后面必须是一行空格)
我用的是SPRINGA:
*ELEMENT, TYPE=SPRINGA, ELSET=aaa
100, 23, 26
所以*SPRING 后面用一行空格，不用定义两点的自由度。
Data lines to define nonlinear spring behavior for SPRINGA or ITS
First line:
Enter a blank line.
Second line:
2。Relative displacement.
你用SPRING2的话，就要指定两点的自由度。
问：请问如何让弹簧不受压
我定义弹簧的语句后面加了
1~‑mn-U4a8d
Y&Q*Q*elastica[1]S`5k+q&
*no compression
wk:c!T.G+|仿真|设计|有限元|虚拟仪器结果出错：***ERROR: in keyword
*ELASTIC, file "Job-1.inp", line 190788: The keyword is
‑pxZ~J‑O3R仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent&
&&&misplaced. It
can be suboption for the following keyword(s): materialSimWe仿真论坛
j.T‑VL‑A:{
***ERROR: in keyword *NOCOMPRESSION, file "Job-1.inp", line 190789:
The SimWe仿真论坛2[
L,pE_1s*B#
&&&keyword is
misplaced. It can be suboption for the following
|)|仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM&
&&&keyword(s):
elasticSimWe仿真论坛iC*F:VTp'l
i仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM是不是我语句加的地方不对啊？
答：“ It can be suboption for the following keyword(s):
material”：
意思是 *elastic 要用在 *material 之后，不能用于 *spring.
．约束（coupling, tie）
请教：加转矩不成功问题
问：实体上一点与面coupling，限制了所以自由度，加转矩，提交后出错***ERROR: 1 nodes are missing
degree of freedoms. The
MPC/Equation/kinematic&&coupling
constraints can not be formed. The nodes have been identified in
ErrNodeMissingDofConstrDef.仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM-[2q
ba*u[1]r"[2p1o
此节点即为加转矩的控制点。不清楚为什么说它没有自由度。JD{#S#i
附INP文档，请各位指教，多谢
你把part的整个表面和它上面的一个节点coupling在一起，这个节点既是control point,
又是被couple的对象（整个表面的一部分），是自相矛盾的。
应该或者去掉coupling, 力直接作用在此节点上，或者定义一个reference point,
作为coupling中的control point。
另外，为什么要把整个表面coupling起来，应该只couple实际问题中的受力区域。
再有，用了kinematic coupling的区域就不能再变形，如果希望变形，可以用distribute
coupling。
问：请问梁单元和实体单元的自由度怎么耦合呢
RT,是不是要采用写约束方程的方法,如果是这样,按约束方程怎么写?
shock111答：
*tie或者 shell to solid coupling在interaction模块的constraint
问：注浆锚杆，由注浆体、锚固体、锚杆（钢筋）组成，请问各位前辈，如果钢筋只有在锚固体中的部分需要模拟、前部分不过是传力杆、且不和注浆体发生作用，应该如何做模型？钢筋的直径小、长度长，如果建立3d实体模型，那这个单元的数量……而如果用wire建模，又要如何处理传力杆部分和注浆体之间的关系？
不重要的那部分钢筋可以用analytical rigid，端部和正常的钢筋tie在一起。
问：钢筋的长度和截面相比实在是相差太大，我要做粘结界面，需要面积，难道只有用三维实体网格？
如果按照abaqus自己的mesh verify，区区一根钢筋就用掉数千甚至上万单元，尺寸太小计算时也容易出错~
可以试试沿钢筋的外围建一圈shell, 材料很软，作为粘结面。钢筋用beam，钢筋和shell之间tie在一起（slave
surface的类型为node region），设置tie的position
tolerance为大于钢筋的直径，以保证钢筋在各个方向上都和shell绑定。
问题同题。就是模型的一个面上的所有节点在计算过程中Z轴自由度位移是相同的，但具体是多少由计算才能知道，在abaqus中有什么命令可以实现？？？
Interaction模块，菜单constraint / create, 选 Equation
问：请教aba_aba版主及其他高手，薄膜和绳单元接触和建模问题(附CAE)？
因为CAE文件太大，应aba_aba版主要求已经删掉网格，但是没有定义接触。其它均定义完毕！&&
‑pE7N)G#^Jk} 我做的是薄膜和绳子大变形分析,薄膜和绳子是两个Parts。我用tie定义薄膜和绳子wire粘接部分接触,最后job出现问题。我认为不能用tie,
因为绳子远离薄膜边界不能满足tolerance。我用其它接触方法也不行!我应该怎么定义接触??SimWe仿真论坛's"E BY
用tie就可以（见附件），注意要把薄膜设为master surface, 伞绳选择node region, 让
tolerance 大于薄膜和伞绳的距离。
你模型中的问题：
1。 *No Compression 的用法是：
2.83e+10, 0.4
*No Compression
*No Compression
2.83e+10, 0.4
2。薄膜单元类型应该是membrane, 不是shell
3. 伞绳单元类型应该是truss, 不是beam
我附件中的模型仍不能收敛，因为truss单元没有旋转自由度，你的固支边界条件无法约束整个模型的刚体旋转。应该只建1/4模型，施加x和z方向的对称边界条件。可能还需要定义y方向的软弹簧，见
/forum/viewthread.php?tid=419588
aba_aba版主 你说的这一条:3. 伞绳单元类型应该是truss, 不是beam 小弟有点不同看法:
Q$W6Q‑i4X仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM我们都认为索不抗弯
即抗弯刚度为0&&但是你使用tie的时候
好像是整根索和薄膜tie&&个人觉得就已经考虑了索的弯曲(不知道是否正确理解版主的意图?)&&个人看法:是不是可以考虑用beam单元做下试试(因为索的截面很小,抗弯惯矩也相应很小.似乎青马大桥和汀九桥的索就是用beam单元做的)
如果用truss单元 在非线性情况下还应该对弹模进行折减.个人愚见 请aba_aba版主指教
讨论出真知！索我没做过，也许你是对的。庄教授的书《ABAQUS非线性有限元分析与实例》p409
讨论了索的建模，他用的是杆单元。
如果用beam, 则 property也要是beam类型, 不能是原来模型中的truss类型.
整根索和薄膜tie , 我并不认为考虑了索的弯曲。
truss单元没有弯曲刚度，如果薄膜弯曲，truss之间的铰接处会随之转动，truss单元本身并不弯曲。
当薄膜纯弯曲时（面内长度不变），意味着被tie的 truss单元要变短，而模型中定义了 truss单元不受压，所以
truss单元不会影响薄膜的纯弯曲。保险起见，可以把被tie的truss单元网格画得非常细
我前处理用的是直角坐标系，但是我想在后处理中输出关于柱坐标的位移分量
是不是要设计局部的坐标系 ？怎样设计？
如果建模时没用过局部坐标系，可以在后处理时点菜单tools / coordinates system / create,
创建柱坐标系（例如使用默认的名称csys-1）. 然后菜单result / options, 点tranformation,
点user-specified, 选中csys-1, 点OK. 窗口左上角显示的变量如果原来是U,
U1，现在就变为U,U1(CSYS-1).
如果建模时的边界条件或载荷中用过局部坐标系，可以在后处理时点菜单result / options,
点tranformation, 点nodal.
注意有时局部坐标系上的 histroy output 会不正常，应检查其正确性。
绘制曲线（X&Y
如果在分析结束后发现丢掉了一小项需要的结果，如何快速解决？
举个例子来讲，模型比较复杂，好不容易花了一个小时算完了，这时突然发现有个点的位移需要在history
output中反映出来，可是当时没有点选它。仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent)W-s5l
这是如何解决呢？我是指在不重新计算的情况下，如何快速解决。
菜单 tools / XY Data / Create, 选ODB Field Output
例如我想用odb文件建立这样一个曲线：
x y(自行指定)
currentmax("my-xy01") 1.0
currentmax("my-xy02") 3.3
………………
搜索了半天也找不到，在此向用过的前辈请教，或者有第三方软件也请指点。（虽然我手工excel也能做，但是那个时间就……-_-b）
后处理菜单 tools / XY Data / Create, 先创建你的"my-xy01"，再选择Operate on XY
Data, 来建立公式。
详见 Getting Started with ABAQUS 附录D.11 Operating on X&Y data
问：我在学习forming a channel例子的时候，发现帮助是将冲头力和位移画在一张图里了（图13-17，见getting
started with abaqus），并且将explicit 和standard的结果也都画在一张图中了，可是我只能将field
output的结果，也就是随时间变化的结果画出来，怎么修改x轴的输出量呢？
后处理菜单tools / xy-data / manager，点create, 选odb field
output，为每条曲线定义这样一个xy-data (两个不同的odb文件也没关系)，选中多个xy-data,
点plot，就可以把它们画在一起。
问：非线性分析，包括n-r，rik法，各增量步的载荷显示。我想画载荷位移曲线
还有，时间步长，弧长和载荷的关系是什么呢。
如果你没有为载荷设专门的Amplitude, 载荷大小就是你说得那样。
Field output中的CF就是集中载荷。如果你的载荷是pressure，恐怕不能直接输出，可以用distributing
coupling 把受力面和一个参考点耦合起来，然后在参考点上施加等效的集中载荷。
载荷位移曲线： Getting Started with ABAQUS 附录D.11 Operating on X&Y data
有一个应力应变曲线的例子。
field output
&/ history output
问：field variable
是什么概念，什么时候应用？请赐教
field output
的输出结果来自于整个模型或模型的大部分区域，被写入ODB文件的频率相对较低，用来在Visualization功能模块中生成云纹图、变形位移图、矢量图和XY图。例如，使用field
output 来在一个分析步结束时输出整个模型的位移场。
history output
的输出结果来自于模型的一小部分区域，被写入ODB文件的频率相对较高，用来在Visualization功能模块中生成XY图.
例如，使用history output 来输出某个节点在所有时间增量步上的位移。
简言之，field output是某个量随空间位置的变化，history output 是某个量随时间的变化。
问：桩土1/4模型，施加的symm边界约束，桩顶表面施加pressure,桩体内 embed 了
rebar，结果运行之后桩顶的四个顶点被疯狂地往下压，其他地方都没反应！
这是什么错误啊！
是显示的问题。你的field output只输出桩顶部的位移，不输出其他节点的位移，造成后处理显示异常。改为默认的field
output就正常了。
如果桩底面和侧面都用接触，等于在底面的圆弧上定义了两次接触，同一个节点有两个法线方向，会出现收敛问题。应该在此处倒一个小的圆角（土也倒圆角），圆角处不定义接触，就没有这种问题。
不倒角的话，土在拐角处的应力在理论上是无穷大的，可能也会造成收敛问题。
问：相对位移问题
我用ABAQUS做了一个楼的地震时程分析，但最后做动画的时候它总是向一个方向跑，我想请教一下，怎么样才能够得到楼体的相对位移值，并对做出相应的动画呢？
v6.6后处理有一个功能
eA3Qt$zda j}SimWe仿真论坛SimWe仿真论坛[1]C'BL:v'gk
set view camera to move with a node
X'h#h9mPv6.5需要你做个固定的坐标SimWe仿真论坛)y&N_zm
然后选择view-&view option里面SimWe仿真论坛't,u3itY9_/B5G
move camera with CSYS
如何设置云图在积分点、centroid、或者节点上显示？
比如，我定义了状态变量FV，如果是减积分单元，只有一个积分点，那么显示的结果就是FV最大为1，最小为0。但是如果我使用的单元内有多个积分点，比如C3D8，显示FV的结果就有可能最大值大于1，最小值小于0，看起来很不舒服，这应该是积分点插值显示云图造成的。abaqus如何设定云图显示在积分点、centroid、或者积分点？或者规范FV显示范围在0到1之间？
(六小) 答：
查帮助,*element output ,postion=?
如何输出拉力（内详）？
我现在做了个粘结剥离的模拟，如图所示，在模拟中剥离过程用A，B两端的位移控制。&Eq[1]]gm8W8M
现在的问题是如何输出整个过程中作用在A,B两端的拉力？谢谢！
把端面和一个参考点coupling起来，位移施加在参考点上，最后看参考点的反力。
查某几个节点的应力值，怎么其中有的节点输出好几个值啊？奇怪仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC
看了半天手册知道了，是averaging threshold的值的问题。但是不是越平均越好呢，比如设置为100％
不是月平均越好吧，要是节点之间差的太多，就把这种变化抹平了，相当于加了人工粘性什么的了，再者一些较细致的变化也就看不到了，类似于气动力学中用一阶迎风差分计算激波，就把激波给抹平了，后来发展了tvd等高阶模式才解决的，个人觉得用abaqus中的默认值0.75就行吧，一般不会出什么问题的
请教：删除失效单元的方法 （aba_aba大侠请进）
在ODB Display中如何删除失效的单元？具体怎么做啊？相信大家都很有兴趣的
提交分析前，在field
output中最下面的state/field/user/time中选中STATUS，然后在后处理中，用display group /
Element / Result value, 把field output 设为 STATUS.
例子见 ABAQUS Verification
Manual&& 2.2.20 Progressive
damage and failure。
不同的set能否用不同的颜色表示出来?同样被赋予不同材料属性的单元能否也用不同的颜色或其他方式相互区别开来?
后处理: 菜单tools / color code
答：云图叠加详细步骤
云图叠加详细步骤：（分别用于应力、变形位移、等）
t[[ u7mg?前提：同一几何模型建在同一cae文件， 在step1， 定义load1； step2,定义load2
(同时inactive load1), run
得到同一模型在两重载荷状态下的云图（分别存在step1和step2）.仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fmj$viG+e
叠加步骤：（应力叠加）
P|CJ9e%i仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM“visualization”&
“tools”& “create field outputs”& “from field” & in the opened
window: give name, for example “AA” &choose “step1” and last frame&
click “s1f0-s”&click “+”&choose “step1” and last frame& click
“s1f0-s”&click “apply”L9D&
To see the combined results, move on to “Result” &”field output”&
“AA”||有限元|模拟|计算|力学|航空|航天|ANSYS|MSC|ABAQUS|ALGOR|COMSOL|Matlab|Fluentj[^9f;|#W
U7wOIA`6Dt8u||有限元|模拟|计算|力学|航空|航天|ANSYS|MSC|ABAQUS|ALGOR|COMSOL|Matlab|Fluent应力叠加不叠加变形，若看变形叠加结果，重复以上叠加步骤，选择位移即可。
{||有限元|模拟|计算|力学|航空|航天|ANSYS|MSC|ABAQUS|ALGOR|COMSOL|Matlab|Fluent欢迎批评、补充。
问：abaqus得到的是积分点处的应力，但是由于不知道从属于那个面，这时候的应力方向如何确定呢？现在我非常着急
菜单plot / symbols, 可以画矢量图。
问：谢谢各位。现在明白了几点：R;i"[#rO2
仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM[1]c8g
1. 比较rpt文件输出的Element Nodal 和 Unique Nodal结果，Unique
Nodal的结果是同一节点相邻单元在此节点处应力结果的平均值，即后者通过前者在同一节点处平均所得，包括mises,s11,s22,s33,t12,t13,t23都是平均所得的。是否可以说Element
Nodal 的结果就是通过高斯点外推出的单元节点处的应力，Unique
Nodal就是那个加权平均值？8Z
X){J:KY:M&[
r!q V4Zx[1]A
2. 比较rpt文件输出的Unique Nodal结果和plot-&contours
view出的mises应力结果,在我的模型中，两者最大mises应力值差别很大（20%），最大mises应力值出现的位置（节点号）也不一样。前者是加权平均的话，后者是如何得出的呢？
y;JRx;^,|仿真|设计|有限元|虚拟仪器仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent:Q'V0J#v0t2y7O
3. 我猜想plot-&contours
view出的mises应力结果好像与积分点的应力值关系很大，但不知道如何把显示的mises应力结果输出到结果文件中保存，在rpt文件输出选项中没有这一项，rpt输出的只是积分点上应力。由于在强度校核时，经常要关心某些区域甚至某些点处的最大应力值，所以想把contours
view出的mises应力结果输出到文件中保存下来，不知能否办到？我觉得这个应力值和rpt输出的Unique
Nodal应力值有差别，但不知道差别在哪里？
看帮助吧。ABAQUS/CAE User's Manual&
24.4.1 Understanding how results are computed
24.4.2 Understanding result value averaging
Element-based field output variables are written to the output
database at the integration points, the element centroid, or the
element nodes, depending on the variable. For the display of
element-based field output variables, ABAQUS/CAE reads values from
the output database for all elements connected to all nodes
included in the plot. Computations are then applied to these values
to produce contours, nodal probe results, display groups or color
coding based on results, or X&Y data along a path.
问：要在后处理里面查看两个点的相对位移，应该如何操作？谢谢！
在后处理中可以直接看，菜单Tools / Query , 选 Distance。 其含义见：
问：question1:
用abaqus/cae, visualization, Query看形变后的长度，scaled 和 unscaled有什么不同
实际形变后的长度是Deformed distance (unscaled): 2.0.
变形前的长度 Base distance: 200
图形显示时的缩放系数 Scale: 8.8011
变形后的长度（不考虑缩放系数）Deformed distance (unscaled): 202.557
（所以变形量为 202.557 - 200 = 2.557 ）
变形后的长度 * 缩放系数 Deformed distance (scaled): 222.5
(即 200 + 2.557 * 8.8011 = 222.5)
注意：3个方向都有变形时，变形量是Deformed distance (unscaled)的magnitude 减去 Base
distance 的magnitude ，而不是Relative displacement (unscaled)的magnitude
问：在选项deformation scale factor里有三个选项
auto-compute
nonuniform
他们的用法是是什么
是不是这只是一个显示的选项 我选不同的选项 变形的结果不一样
是不是uniform=1才是真正的变形情况
auto-compute: ABAQUS自动确定变形的缩放系数。
uniform： 用户自己输入变形的缩放系数，x,y,z方向上的缩放系数相同。
nonuniform： 用户自己输入变形的缩放系数，x,y,z方向上的缩放系数可以不相同。
问：我把一個複合材料薄版
但是結果沒出現應力圖示
可否幫我看看哪裡出了錯
这种模型好像是不能在后处理中显示应力，例题也是这样，见ABAQUS Benchmarks
Manual&&&&
1.1.3 Composite shells in cylindrical bending中的
compositeshells_c3d20r.inp。
例题是把应力输出到dat文件：
&1,2,3,4,5,6,7,8,9,10,11,12,13,14,15
&S,MISES,PRESS
C3D8单元精度很差，例题中都是用C3D20R或C3D20, 有接触或塑性的话就要用C3D8I，
你的单元太薄，计算结果可能会有问题。
问：变形体用云图，而没有结果的刚体用透明，如何设置呢？
菜单options / contour, 把visible edges设为feature edges,
点other标签页，选translucency.
怎么把一连串ODB文件连接起来做成动画?
步行去远方& 答：
我只知道一个一个的输出，楼主分别输出动画以后然后去找个专门的视频处理软件让它们连续起来行不行啊？
MWINDOWNS自带有个可以把许多照片放在一起然后处理成视频格式文件的软件，处理连续动画应该也是可以的。名字好象叫MOVIE
MAKER。系统自带的，楼主不防找找看。
^仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,FluentABAQUS
V6.6 Analysis User's Manual&
[X'L9Jf7z4d仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM3.2.16
Execution procedure for joining output database (.odb) files from
restarted analyses
shell和beam的应力应变结果
关于shell单元的疑问？(附位移和应力分布图)
我手头有一个关于shell单元的模型，两端简支的圆管，用s4单元，给予圆管的右端一个延轴向的拉伸，为什么s33为0，很奇怪阿？附位移和应力分布图如下，谢谢！
"s1~s6分别为X，Y，Z和绕X，Y，Z的位移" 适用于3D
solid单元或平面应力/平面应变单元，shell和beam则不是：
- Shell单元：
S33是曲面的法向应力，S11，S22是曲面的两个切向应力。见
ABAQUS Analysis User's Manual& 15.6.7
Three-dimensional conventional shell element library
Element output：
If a local coordinate system is not assigned to the element, the
stress/strain components, as well as the section forces/strains,
are in the default directions on the surface defined by the
convention given in “Conventions,” Section 1.2.2.
以及& ABAQUS Analysis User's
Manual& 1.2.2
Conventions&&&
Figure 1.2.2&1 Default local surface directions.
- beam单元：
SF1是梁的轴向力，SF2和SM3分别是梁截面局部2方向（n2）和1方向（n1）上的剪力，SM1和SM2分别是梁截面局部1方向（n1）和2方向（n2）上的弯矩，SM3是梁的轴向扭矩。
问：我用EXPLICIT的SHELL作模拟
怎样输出每个积分点上的应力和应变呀？
學生& 答：
*node print
1,2,3,4,5,6
1~~6表示積分點由最下面一層算起
一層有三個(看你用什麼積分法就有幾個積分點我選用三個的)
所以在層與層的交界處會有兩點在同一平面重複
像是我舉兩層例子
1,2,3 4,5,6
3,4要分開但他們在同一平面上
henryzhang答：
在Field output 中可以选择所要输出的积分点
abaqus用梁单元计算时候遇到的两个问题
1.采用b33梁单元，想提取s11的力，结果发现总是出现四个截面点的力，我想提取的是梁单元积分点的力，这在应力云图上可以看出来，但总不能一个一个的用手记下来吧，所以求大家给个主意，怎么提取啊。||仿真|设计|有限元|虚拟仪器:E1b,v
2.在整个模型中，定义一个step，然后先加重力，后加集中力
[1]C1M;t#*DLOAD
XTOT,GRAV,0,0,-9.X4o&
{HT:d*h-`& & &
& q:a6ASimWe仿真论坛5,&
的计算结果怎么和只加一个集中力的计算结果相同呢，感觉好像重力并没有加上去，难道重力和集中力不能在同一个荷载步中出现吗？
重力的语句缺少一个参数，应该是
TOT, GRAV,&
1,&& 0, 0, -9.8
你密度的数量级正确吗？见&
&/forum/thread--2.html
ABAQUS常见问题汇总（1.0版）&&&
9.1& 各个量的单位要相互一致
你用probe values看到每个单元的S11有三个值，那是默认截面点上的三个值，如果点菜单Result / Section
Points / beam pipe, 选择不同的截面点，例如 angle = 180，
S11的值也就会变化。xy表或dat文件中显示了所有截面点上的应力（S11是应力，而不是力）。
如果希望得到准确的应力，可以在分析前在field output中添加SF,
以得到梁弯矩SM，然后根据材料力学的公式自己计算应力。
．多步骤分析
（Restart）
问：我学习了一下帮助文件anylysis user's manual中的7.7节，和Example
manual中的15.1节，大致知道关键词*import可以使一个计算过程的结果传递到下一个过程中，并作为起始条件进行后续的计算，但不知道初始计算的INP文件和后续计算的INP文件，怎么联系起来？也就是说怎么进行两个过程的计算。希望做过这个方面问题的大侠指点一下啊！
如果单元、节点等模型参数没变，只是想用新的step来施加载荷或边界条件，就不用*import，可以用RESTART,
见Getting Started with ABAQUS 11.4 Restart analysis。
（Submodel）
（个人经验）
***ERROR: AN ERROR OCCURED WHILE OPENING/READING THE .FIL FILE.
MAKE SURE YOU
&&&&&&&&&&
SPECIFIED THE CORRECT FILENAME WITH THE EXTENSION. FILENAMES
&&&&&&&&&&
CASE SENSITIVE ON CASE SENSITIVE OPERATING SYSTEMS.
需要读入的全局模型结果文件包括 .odb文件（或 .fil文件）和 .prt文件，
-&&&&&&&&&
这些文件都应该在ABAQUS的默认工作路径下
-&&&&&&&&&
这些文件的文件名不能太长
SUBMODEL的边界条件
SUBMODEL的边界条件b;`KV^0^)W6]
~.s"zl(YK,h/B*`
我的SUBMODEL是从母模型内部取出的一小部分，这个||仿真|设计|有限元|虚拟仪器0g5Wk@
*BOUNDY&&SUBMODEL&&STEP=1
DS~)p"zSimWe仿真论坛DRVEN，？？？
y I$|‑|$H*p1]边界条件该怎么加？如果不加就报错，（是SOLID TO
SOLID形式的）。
h#o AR)看例子是1-6自由度都被限制，这样的话，子模型的位移怎么办？
x仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMSimWe仿真论坛9}
还一个问题啊，就是在SHELL TO SOLID 里的
%p$Y3u& & &
]||仿真|设计|有限元|虚拟仪器*BOUNDY&&SUBMODEL&&STEP=|仿真|设计|有限元|虚拟仪器!?|m_sh&ko~k
DRVEN，仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM4J9l$?n&
& & & UP'l
这个限制的是什么自由度？我*PREINT之后看到DAT文件写的是限制SHELLTOS是什么意思？仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM+PR)BM4h
& &多谢大虾指点
*BOUNDY, SUBMODEL, STEP=1
节点集合, 1, 6
意思是: 这些节点各个自由度上的位移都将等于母模型在step 1中相应位置上的位移，而不是你理解的位移为0。
如果你没有shell单元，就不用SHELL TO SOLID 。
把计算结果传递到下一个过程中
问：请问各位大侠，显式分析后想继续用隐式分析，除了load外所有的其他定义都不改变，试过restart好像不行，请问应如何进行？CAE里可以实现么？谢谢！
ABAQUS/CAE User's
Manual&&& 16.6
Transferring results between ABAQUS analyses
ABAQUS Analysis User's
Manual&&&&
7.7.2 Transferring results between ABAQUS/Explicit and
ABAQUS/Standard
ABAQUS Example Problems
Manual&&&&
1.5 Import analyses
．岩土分析
nevercry_zju : ABAQUS岩土问题讲座（完全）
：[ABAQUS下载专区]&
abaqus讲座资料，岩土专题，没有积分限制&
asamurai : abaqus在基坑工程中使用小结
本人在使用abaqus解决基坑支护问题的过程中，得到了版上aba_aba，以及ustc的life大侠等人的帮助，写些自己使用abaqus的一些感受吧，（其实大多数东西都是从版上得到的相关参考）希望对其他新手能有所帮助。
**首先是建模方式。
对于初学者比较好掌握的是CAE界面操作方式。但是由于ABAQUS的CAE模块和输入文件两种建模方式是由两家不同的公司研制的，CAE模块推出比较晚，所以功能还不是很完全，对于CAE和INP交互操作需要注意几点：
1、一些功能只能通过编辑INP输入文件来实现，比如模拟开挖建造时需要用到杀死或者激活单元的语句model
change，此时需要直接在CAE下编辑关键字，或者生成INP文件后通过记事本等工具打开INP文件后加入关键字。
2、CAE对于INP输入文件中的一些关键字不支持，如果在CAE下编辑，则系统会给出不支持关键字的提示。对于这样的情况则只能通过文本编辑工具打开INP来操作。
3、对INP输入文件进行编辑的方式能支持所有关键字功能，因此应用更加广泛。在许多情况下都需要对INP文件进行操作。
**然后是初始应力场的施加。
土工有限元分析中初始应力场的施加是个关键。在ABAQUS中，施加初始应力场很方便，可以使用以下语句来完成。
initial conditions,type=stress,geostatic
可以有两种实现的办法：
1、如果地表水平，且初始应力就是自重应力时，可以在数据行中按照如下顺序填写数据
施加初始应力的集合，应力值1，应力值1对应的坐标值1，应力值2，应力值1对应的坐标值2，x方向侧压力系数，y方向侧压力系数（在二维问题中可以省略）
这样的顺序定义初始应力实际上默认了初始应力就是自重应力，只需要给定集合中2点的应力值及其对应的坐标值，其他部位的应力值就可以通过1，2两点的连线采用差值来得到，而水平向的地应力则乘一个侧压力系数得到。（这些内容在abaqus手册中都有）&&&
2、采用input文件的方式，即先通过计算一个重力荷载作用下没有初始应力的模型而得到各点的应力值，然后将这些值导入到新的模型中，在*initial
conditions,type=stress,后面加入input=XXX.XXX，让模型接受存有初始应力值的input文件而施加。这样的方法通用性非常强，但是比1稍微复杂些，具体可以见相关手册。
另外，对于施加初始应力，对于二维的情况，系统默认地应力方向为Y方向，对于三维情况，默认为Z方向，这一点在初学者中容易弄错。
**模拟开挖建造
ABAQUS中模拟开挖建造通过关键字model
change,type=element来实现，如果是杀死单元则在后面加上remove,如果激活单元则加上add。
事先需要将需要杀死或激活的土体单元定义为一个集合set。整个应用格式如下：
*Model change,type=element,remove
Set name1，set name2
其意义如下：首先Model
change说明是模型发生变化，后面是种类，因为除了杀死激活单元外，还常用到杀死激活接触对(type=contact)，此处需要指明是杀死激活单元导致的模型变化，后面remove表示杀死单元。重启一行写出需要杀死激活的单元集合列表。
在模拟建造过程中，通常一个单元的材料属性会发生变化，比如土钉支护中面层处单元，在没有开挖时是土体材料，当开挖建造之后就变成了面层的属性了。在ABAQUS中模拟材料属性随开挖步进行可以使用场量*Field关键字来实现。具体过程如下：首先是材料定义上要指定材料性质随着场变量而变化，如下面语句中，定义弹性时，利用dependencies=1表示材料根据1个场变量变化而变化，在接下来的属性数据行中每行后面的数字1，2表示了该行属性对应的这个场变量的值，如第一行的弹性模量和泊松比值是场变量为1的值，而第二行的弹性模量和泊松比值是场变量为2的值。
*Elastic, dependencies=1
&& 4e+07, 0.35,
&1.5e+10,& 0.2,
接下来就是在每个分析步中利用*field,variable=1这样语句指定分析步中使用的场变量数量，本例中值为1个场变量。接下来的数据行则是表示使用场变量的集合名称及其在该分析步中的场变量值。如果这个值是1，根据在定义材料时的场变量定义，值为1时对应于材料属性第一行值，这样就实现了材料属性随分析步变化的功能。
**模拟地连墙与土的接触ABAQUS中对于接触的模拟可以通过定义接触对来实现。系统可以允许面对面、点对面的接触等多种接触形式。在模拟开挖过程中，开挖土体与墙体的接触属性会丧失，根据经验，最好能将接触按照开挖步骤划分为许多对接触对（所谓的细化），根据开挖的进行，将接触对先杀死，这样系统比较容易找到准确的接触状态，使接触模拟顺利进行。模拟接触时，定义接触对要分清主控面和受控面。主控面和受控面的选择一般遵循以下原则：一、材质较软的为受控面，材质较硬的为主控面。二、网格较密的为受控面，网格稀疏的为主控面。应用到网格划分的时候就要先根据材质的相对软硬程度确定主控面，以便在网格划分中按照网格疏密关系正确划分网格。主控面和受控面定义一定要正确，这是影响接触模拟的关键。如果一个节点上同时定义了两个接触面，或者同时定义了一个接触面和边界条件，系统则会产生过约束(over
constraint)的提示，很容易影响到接触的模拟。
写的不清楚，请大家多包涵。
地应力平衡
问：如何施加土体初始应力！地应力
solong答：
我是新手、自学者，总结了一点经验，主要是给和我一样的新手
*initial condition,type=stress,geostatic 就不说了，啃啃手册就明白了，或者在论坛搜索下
我主要说一下*initial conditions,type=stress,input=XXX.XXX的方法：
1、先建立好你的模型，按照你的课题把边界条件、接触什么的都建好，然后使用命令mdb.models['你的模型名字'].setValues(noPartsInputFile=ON)，建立没有part的inp文件（abaqus不认带part的istress数据）
2、建立*geostatic的step，只施加重力荷载，然后求解
3、完成之后用viewer打开odb，选report-&field report,下拉菜单里面选择
centriod，然后依次把s11、s22、s33、s12、s13、s23点选上，去setup页面选择报告文件的名字（例如istres.rpt)把下面的和值之类的额外输出点掉，要不还要手动删除
4、编辑你的报告文件，把数据行之外的东西都删除。（正确的格式：单元号，s11，s22，s33，s12，s13，s23），使用excel或者ultraedit在每个数据间加上一个逗号（
, ），数据之间的空格不要紧，但是要注意不要有空行或者逗号之间是空白，否则会出现element 0错误
5、重新建立你的step，写入新的inp文件之后，编辑它，在材料属性之后、step之前加入*initial
conditions,type=stress,input=XXX.XXX
搞定，位移一般都被控制在10-6数量级，满足要求。
这是python命令，在ABAQUS/CAE底部标有&&&的窗口中输入。'你的模型名字'
的默认值是'Model-1'.
& 我的初始地应力平衡方法!
看了版上很多帖子,受益菲浅!
K;PoB?-X仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM关于地应力的平衡方法,综合了版上的一些意见,结合了自己的想法,对于初始地应力的施加,得到了e-6的效果,方法比较简单,与大家分享!仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,FluentR+sb6{.DP,VEB[1]M
1.先施加重力荷载的作用,可以在cae中实现;SimWe仿真论坛'g|H!f+
2.在inp文件中的output request中写上
"O‑@HPg5f3E)L%_*ws仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent*el
1Nh0tK& & &
Is,仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent3T)l"n$s_
这样就会将施加重力荷载后的应力输出到*.dat文件中了;
%D#^Vnp P1O
O-o3^t2PSimWe仿真论坛3.在*.dat文件中,将单元应力的序号及单元的应力拷出,例如仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM7h)c2b}M5E
ELEMENT&&PT
&&&仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM/X3DK:rA!L&cS(Q
SimWe仿真论坛 U(Uw ^/l1z&A
tS d1w p1H(T
&&&-1.2598E+05&&-1434.&
&&&-3.1852E+04&
u仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent&
&&&-1.2249E+05&&-6287.&
&&&-3.2194E+04&
9C8Pg R'j:n!puvSimWe仿真论坛&
&&&-1.1795E+05&&-497.7&
&&&-2.9611E+04&
&1664.仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM4B0_.H;Q7n
&&&-1.1210E+05&&-7240.&
&&&-2.9834E+04&
&1992.SimWe仿真论坛7z@9V5QD0R3C/}
&&&-1.0485E+05&
&&&-2.6068E+04&
&2600.仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent:qs@-cmm"a‑R.P
&&&-9.5803E+04&&-8272.&
&&&-2.6019E+04&
w7e8x%s,i仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM&
&&&-8.4709E+04&
&&&-2.0698E+04&
&4083.仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM_A&
&&&-7.0634E+04&&-9746.&
&&&-2.0095E+04&
GD||有限元|模拟|计算|力学|航空|航天|ANSYS|MSC|ABAQUS|ALGOR|COMSOL|Matlab|Fluent&
&&&-5.1088E+04&
&&&-1.1422E+04&
Pk仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM&
&&&-2.4353E+04
-1.1150E+04&&-8876.&
&1.2126E+04SimWe仿真论坛:VX+A/d B
&&&-1.2847E+05&
&&&-3.2050E+04&
T4SJ`||有限元|模拟|计算|力学|航空|航天|ANSYS|MSC|ABAQUS|ALGOR|COMSOL|Matlab|Fluent&
&&&-1.2786E+05&&-9868.&
&&&-3.4433E+04&
l.tj?5?‑j3r& &&
&&&-1.2938E+05&&-4224.&
&&&-3.3402E+04&
9er.{LLP(|9v仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent&
&&&-1.3039E+05&&-3458.&
&&&-3.3461E+04&
4Q5J'S,l$N‑~&YU仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent单独存为一个*.dat文件,
4.用excel打开该文件,将其中的1所在的列去掉,在每个单元号前面加上其instance.part.,即单元编号变为:instance名称.part名称.序号;注意不同的instance和part要都按照其所在的单元从小到大编号,而不是按照他们在整体单元编号来编号!
Mtf&I-[;Sd'V2xI||有限元|模拟|计算|力学|航空|航天|ANSYS|MSC|ABAQUS|ALGOR|COMSOL|Matlab|Fluent5.接下来就在excel把该文件另存为*.csv格式的文件(即带有逗号分隔符的格式),
"fDI%lP}$X
r6.最后在inp文件的step之前写上
U:HpnlTS;*initial
conditions,type=stress,input=文件名.csv即可!仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent&
这种方法不需要用python,比较简单,希望能对大家有用!
: solong和lfycsu提供的初始地应力解决问题.希望大侠跟进,把这个弄清楚.
我最近在搞初始地应力问题,只有一种方法*initial
condition,type=stress,geostatic应用起来十分顺手,简单好用,但是对于复杂地形,地面不水平的情况就有点麻烦,论坛里有另外两种方法:
0cP'y.O.W&F
B[1]第一种,solong在" [送给和我一样的新手]
如何施加土体初始应力!"文章中提到的方法,但是这个方法我在使用mdb.models['xxx].setValues(noPartsInputFile=ON)，建立没有part的inp文件(abaqus不认带part的istress数据)这个过程中,建立不了所需要的inp文件,我输入的命令没有任何错误提示,那就应该是没有语法错误的，但是在我的TEMP文件夹中找不到xxx.inp文件,我不知道这个过程我的操作有什么问题.请问哪位高手在使用这个方法运行成功了????希望提供点经验.
+Z6w2`g"|,]仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent第二种,lfycsu在"我的初始地应力平衡方法!"中提到的方法,但是在按我的理解提交后就是要出错.上面说到的这一步:用excel打开该文件,将其中的1所在的列去掉,在每个单元号前面加上其instance.part.,即单元编号变为:instance名称.part名称.序号;注意不同的instance和part要都按照其所在的单元从小到大编号,而不是按照他们在整体单元编号来编号!请问这个instance名称.part名称是什么,比如说我在ABAQUS6.5的model
tree的Assembly下的instance下面的名字diyingli2-1(我个人的例子)是不是就是instance名称,然后在Module的Part步后面的Part窗口里面的名字就是part名称??仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent
lGq5[4E:H1|
希望论坛里面的高手,大侠多多帮助,替大家解决一下这个一直很难搞明白的问题.我不是说上面两种方法是错的,而是我(很多人应该都碰到)在使用的过程中出现了好多问题,希望使用过这两个方法的高手们为大家提供点经验.在此谢谢大家了.以下有附图.仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent$D]c#x%k
A(UmIn‑Yck:C-R仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent[
本帖最后由 loveestboy1981 于
问：(请教)有谁做过对斜坡加初始应力
我在模拟斜坡上的路基,地基土的自重位移应该去掉,参考了*intial
condition的几个参数,没有一个适合的呀.仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMM%]M8ac+A7z
有过相同算例的xdjm,给点提示,谢谢
在*initial
condition选择*.inp作为初应力输入就可，要用inp编写Im‑h??:T0X2[
初应力可根据自重计算后导出，再用类如origin等的软件打开，修改成abaqus的应力格式输出，再在abaqus的inp里调用就可
刚才在傲雪论坛中，搜索了一下关于加载的问题，以下是引用：
{i仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent1、通过前一次的计算,通过report输出stress文件,并将与各分力的量不相关的信息删除,然后在excel可批量生成ISTRESS.INP所需要的格式,&#Q&
S[1]U[1]LL
点号,S11,S22,S33,S12,S13,S23 "["o~LN+U0L({A
1,-8.08E+05,-2.59E+06,-1.12E+06,-7.46E+04,819.572,-2.12E+04仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM9N
?+J"b& & &
2,-8.06E+05,-2.57E+06,-1.12E+06,-7.25E+04,-4.49E+03,-7.89E+03SimWe仿真论坛6[7`1[j
3,-8.25E+05,-2.62E+06,-1.11E+06,-7.48E+04,-8.39E+03,-2.82E+04
|‑NzSimWe仿真论坛4,-8.06E+05,-2.65E+06,-1.09E+06,-6.32E+04,1.82E+03,-2.57E+04
$W[1]ql't6l
|仿真|设计|有限元|虚拟仪器5,-8.07E+05,-2.72E+06,-1.08E+06,-6.51E+04,-1.58E+04,-3.53E+04
1d#h#~7l_&x‑I+H仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM6,-8.67E+05,-2.74E+06,-1.11E+06,-5.43E+04,-1.01E+04,-5.45E+03
I#o7V#?(z7b7,-9.93E+05,-3.28E+06,-9.98E+05,-1.03E+05,-1.22E+03,-8.85E+03
/_#R1J"{&c&Z#M'C.Bc仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM8,-9.90E+05,-3.27E+06,-9.99E+05,-1.03E+05,-6.20E+03,-1.81E+04
|仿真|设计|有限元|虚拟仪器9,-9.80E+05,-3.24E+06,-1.00E+06,-9.86E+04,-1.24E+04,-3.60E+04仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,A
10,-9.69E+05,-3.20E+06,-1.02E+06,-8.84E+04,-1.66E+04,-3.64E+04
b/q?){/`仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM11,-9.65E+05,-3.15E+06,-1.02E+06,-8.31E+04,-1.21E+04,-4.18E+04仿真分析,有y
12,-9.49E+05,-3.13E+06,-1.03E+06,-1.09E+05,-1.96E+04,-3.42E+04
-yzi[1]l xh4Y}6|g
tK仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent13,-9.53E+05,-3.06E+06,-1.09E+06,-8.38E+04,-1.65E+04,-8.31E+04www.simwe.coX
14,-4.22E+05,-1.69E+06,-7.02E+05,3.57E+04,6.76E+03,-2.57E+04
x"P[1]q"}/|仿真|设计|有限元|虚拟仪器15,-1.61E+05,-1.22E+06,-4.50E+05,-727.629,4.21E+03,-3.04E+04
Q4N%a.q&B仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent16,-1.03E+06,-2.96E+06,-1.30E+06,-1.95E+04,-4.05E+03,-2.04E+1MlZ+^bjP2W
17,-1.03E+06,-3.14E+06,-1.36E+06,-1.54E+04,2.20E+03,-1.04E+04w!
但最后的出错提示为：,X!er[1]G?yIR*R
***ERROR: in keyword *INITIALCONDITIONS, file "structure.inp", line
-F#Y8Kfg+_&
&&&Unknown
assembly id 1.SimWe仿真论坛 |!?{9k&N
要建立没有part的inp文件（有part的话，stress数据中的节点编号就不能写成1，2，3...,
而是要在前面加上instance名称，比较麻烦）, 方法见
/forum/thread--1.html
问：初始地应力平衡,
我的inp文件修改后二次读进后,ABAQUS里却并没有模型出现,是怎么回事?
ABAQUS/CAE User's
Manual&&&&
9.9.3 Writing input files without parts and assemblies
Some ABAQUS/Standard and ABAQUS/Explicit functionality (such as
*IMPORT) is not supported in a model that contains parts and
assemblies.
abaqus可以形成两种inp文件，应该用没有part和assembly的那种，这样单元和结点的编号不会重复，*initial
conditions, type=stress,
input=xxx.inp才能正确执行，参考"cae_no_parts_input_file"
对于有多个part的网格，在网格划分时候，inp文件中对应的网格和节点都是各自编号的，如何实现整体网格节点编号和单元编号？
不尽然, 其实CAE生成的INP文件完全可以不包括parts和assemblies信息,
&有两种办法可以实现此目标:
&(1) 运行CAE前修改ABAQUS环境文件(abaqus_v6.env)
中的cae_no_parts_input_file:
[color=red]cae_no_parts_input_file=ON[/color]
这种方法有个问题就是若要INP文件包括parts和assemblies信息就必须在&&&
CAE运行前改成OFF, 不是很方便.
&(2) 如果需要在CAE运行的时候实现上述目标, 可以在ABAQUS/CAE的
command line interface运行下述ABAQUS Python Script命令：
[color=red]mdb.models[[i]modelName[/i]].setValues(noPartsInputFile=ON)[/color]
此处[i]modelName[/i]是你的模型名称.
阅读了ABAQUS/CAE User's Manual&
&13.8.3 Writing input files without parts and
assemblies就很清楚了。
摘抄武大一位博士的论文中对初始地应力的阐述与大家分享
希望能对像我一样的新手有所启发。不过其中很多细节看了之后仍不得要领，希望大家各抒己见，共同提高~
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
在岩土工程数值计算中，初始地应力场是必须予以重视的问题。这主要是基于以下三个方面的原因:
(1)岩土工程的特点决定了分析手段多为增量分析，在增量分析中，分析域内的应力总是由应力增量加上初始应力而得，即初始地应力从一开始就影响了分析过程；
(2)岩土材料的刚度和应力状态有关；
(3)在开挖分析中，开挖导致的荷载是通过开挖之前的应力计算的，为了计算开挖荷载，必须首先知道初始应力状态。
虽然地应力场的测量技术取得了较大的发展，但在实际工程中往往由于测试条件、工程投资等方面的限制，只能对为数不多的点进行量测，难以获得对地应力场整体分布的完整认识。人们只能分析这些少数测点的资料，并采取一定的方法来推测整个地应力场的情况。目前主要是通过两种方法来获得初始地应力场，第一种是由一些位移实测值通过反分析得到，这种方法大多限于地下工程，且理论和方法都远未成熟；第二种是根据现场己测得的某些点的地应力值，通过简单回归而得，通常表现为在不同的区域，给出一组以深度为自变量的应力分布函数，这种方法虽然很粗糙，但仍然是目前使用最普遍的方法。
在获得了初始应力场的分布规律之后，接下来的问题就是如何在有限元模型中施加。不管采取何种方式施加初始应力场，始终要满足下面两个条件：
(1)平衡条件。指的是由应力场形成的等效节点荷载要和外荷载相平衡，如果平衡条件得不到满足，将不能得到一个位移为零的初始状态，此时所对应的应力场也不再是所施加的初始应力场;
(2)屈服条件。若通过直接定义高斯点上的应力状态的方式来施加初始应力场，常常会出现某些高斯点的应力位于屈服面之外的情况。超出屈服面的应力虽然会在以后的计算步中通过应力转移而调整过来，但这毕竟是不合理的。当大面积的高斯点上的应力超出屈服面之后，应力转移要通过大量的迭代才能完成，而且有可能出现解不收敛的情况。
基于以上两个条件，我们建议采用下面的方法来施加初始应力场。如果自重应力场就是初始应力场，直接将重力荷载施加于有限元模型，并施加相应的边界约束，计算得到在重力荷载下的应力场。再将得到的应力场和重力荷载一起施加于原始有限元模型，便可得到一个既满足平衡条件又不违背屈服准则的没有位移的初始应力场。
重力荷载的施加问题
假设分析过程中施加重力荷载，有以下两个不太确定的简单问题：
{A'H8Efo&ePqdSimWe仿真论坛1.&&created重力荷载的那个分析步的类型是不是应该选择geostatic类型？WM
^6[4Uq-]2c;x
2.&&load模块中创建重力荷载时，荷载类型是不是应该选择gravity类型？
y&] J请各位大虾对这两个问题做个判断，小弟感激不尽。
1.重力荷载不一定要施加在GEOSTATIC这一步，主要是看你需要步需要平衡地应力了，如果你要平衡地应力就要建立GEOSTATIC这个分析步，布需要平衡的话，随便施加在哪个类型的分析步都可以，只要你需要重力作用。Ibj.?0g&Y
2.对，重力荷载是要选择gravity。
求助：为什么算的结果这么小
埋深70m的隧道，为建模方便，我把上面55m等效为匀部荷载，整个岩层均一考虑，岩层密度2600kg/m3。最后的最大沉降只有零点几个毫米，麻烦帮忙看看，谢谢了！
你的模型错误很多。仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,FluentO-aBh!y
1、地应力不平衡；SimWe仿真论坛SL8K
2、上面55m的等效均布荷载施加的方向不正确；||仿真|设计|有限元|虚拟仪器#c LD
3、材料的单位不统一
D)] ln Z6p
apPz7q/s‑}H仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM*Material,Name=Steel仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMf"vIT6j`1G^?0t8yL
**U3tp+I J t
2.3,仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent"G
&4、努力吧！
andy7：自己做的一个基坑开挖的例子
这是弹性条件下的inp文件
( 23:08, 49.49 K)
基坑开挖隆起问题
觉得目前的模型都很难正确模拟好基坑问题,尤其位移情况可以发现与实际情况相差很远,个人觉得这主要是由于土的模型不完善导致的.4`B-n)wkzK%e
1.开挖过程中,我们输入的弹性摸量到底是什么?开挖实际上是土体泄荷的过程,但坑后的土体在竖向上又有压缩的趋势,而压缩摸量跟回弹摸量相差甚远,如何考虑?相对而言,修正剑桥模型更好的考虑了这一点.
E%|0s[1]jcg仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent2.土体基本不抗拉,出现了拉应力土体不能继续承担荷载,必然应力会相应迁移,但现有模型都不能考虑这一点,使得回弹力变得很大,导致回弹比实际大很多,而且这个回弹力使得挡土墙向后仰,而不是向前倾.
钢筋混凝土（Rebar
Layer和Embedded Element）
&[请教]钢筋混凝土本构模型与粘结滑移问题
钢筋混凝土本构模型有两种：concrete smeared model 和concrete damaged
model仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluentw&v'}`
两种模型中都有tension stiffening
部分，[手册]该部分可以模拟钢筋和混凝土之间的相互作用(bond-slip)。我看了些帖子，上面说考虑相互作用，要加弹簧单元。因此我有一个问题，就是加弹簧单元和两种本构中tension
stiffening 有什莫区别？我想考虑之间的相互作用关系，采用那种方式？谢谢
(紫菜)答：
tension stiffening
认为混凝土开裂后应力并未完全释放，仍有一部分抗拉能力，由于钢筋的存在，荷载在混凝土单元裂缝间传递
l&[)ZH~2[G& &
L仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent如果加弹簧的话，钢筋就不会被embed在砼中，如何定义弹簧刚度是一个小问题
L|xWW3B比较突出的问题是，如果你加弹簧，弹簧会在整个受力的分析过程中都起不变的作用，但实际上，粘结滑移只是在一定时间范围内的作用，在弹性段的时候，一定不会有粘结滑移作用，而后面的过程粘结滑移也不会是一成不变的，当然tension
stiffening也只是近似的模拟这种性质，也不好说就更好，但个人偏向于tension
stiffeningSimWe仿真论坛]A9S}
{:c5J"]‑g8V@
仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMQ/pwr‑G
另外好像cohesive单元能模拟粘结滑移，不过这样做的代价是钢筋都要用实体元，做局部分析还成，要是比较大一点的模型就。。。
YSimWe仿真论坛&q)h‑iqY!t;A:R
以上都是小弟拙见，难免错误，望大侠指正
lordonline： 一个简单的钢筋混凝土梁，供初学者参考 ，使用了*Rebar Layer和*Embedded
reinforcedconcrete beam.inp (16.22k)
混凝土采用实体单元，想在其中加钢筋，如何实现呢？请大家指教，谢谢了
zjhsimwe答：
*rebar和*embeded是两种完全不同的处理方法，如果要考察筋材的性能，用*rebar肯定是不行的，因为该关键词的含义是局部增强来模拟加筋，不能对筋材进行显示，只能考察加筋后的主体构建的性能。而*embeded可以把筋材和主体构建分别显示，方便对筋材性能的考察。
如果是单根钢筋，可以用*Rebar
如果是钢筋层，可以先建立一个Shell面Part，然后Embedded到实体Part中。
abaqus在混凝土中加钢筋的两个办法
lordonline答：
采用rebar layer 的办法，在part里面画一个面，然后在property里面定义一个surface为rebar
layer，把这个surface的属性赋给前面的part里面的那个面。然后在interation中embed中把钢筋层embed到混凝土实体中去。
2。采用桁架单元的办法，在part里面建好纵筋和箍筋的钢筋骨架，在property中分别赋予截面和属性，在interation中的embed把钢筋骨架embed到混凝土的实体中去。
3。如果是作构件的话，第二种办法建的比较精确，而且后处理比较方便，查看钢筋单元的应力比较直观，如果是作结构的话，第一种钢筋层的办法比较好，但是个人觉得钢筋层的办法纵筋和箍筋的位置定义的不是很明确。
&对于剪切破坏的钢筋混凝土构件,箍筋的对于抗剪作用比较关键,必须要建立箍筋单元,对于受弯曲破坏的钢筋
混凝土构件,可以不建箍筋。 &而在ABAQUS中，纵筋可以通过rebar
&layer的方式施加，箍筋采用rebar &layer
的话，位置定义的不适很明确，下边我自己想出了一个采用桁架单元定义纵筋和箍筋的办法。
CAE箍筋的建立的办法画图说明一下: &
1.在part中画一个wire, &
2.赋予wire截面和材料属性 &
3.在asemnly中插入,平移或者旋转来布置箍筋的位置 &
问：求助：大侠们帮忙看一下，附inp
我在实体单元里面定义了嵌入单元，然后加载
bm0b+F0]仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent如果嵌入和加载都没有问题的话，理想的膨胀状况应该是G5X2K|8U‑wT1gv
均匀的，但结构总是在某个直径方向膨胀的很厉害，不明白原因
我认为是rebar方向的问题，把你的
*REBAR LAYER
M1,2.1.00116,0.,BELT,&&
*REBAR LAYER
M2,2.1.00116,0.,BELT,&
*REBAR LAYER
M1,2.1.00116,0.,BELT,&
*REBAR LAYER
M2,2.1.00116,0.,BELT,&
变形就均匀了很多（增大变形显示系数还是会看到不均匀）。问题可能是，目前的rebar方向是基于默认的坐标系，可能应该定义一个局部坐标系（柱坐标系？），见ABAQUS
Analysis User's Manual& 2.2.3 Defining
reinforcement，Figure 2.2.3&8 。
另外可以试试细化网格是否能变形均匀一些。
你的material和set名称都用了BELT, 容易混淆，建议区分开。
inp文件中不要有空行，容易出错。需要空行的地方应写两个星号。
求助：急请高手帮忙看一下钢筋混凝土梁INP文件
一个钢筋混凝土简支梁，长3m，两个支座距两端各0.2m，两个集中力作用在距两端各1m处，梁截面为150&300，配筋为梁下部配5根直径为16mm的二级钢筋，上部配2根直径为6mm的一级钢筋，箍筋直径为6mm的已经钢筋，间距为200mm，混凝土弹性模量为42840MPa,强度为66.67MPa,我钢筋是用TRUSS来模拟的，将一个钢筋骨架那个PART，embedded混凝土中去。
K,X{2a L4MbSimWe仿真论坛&
这个文件是可以运行，但是出来跨中处的挠度跟实验结果差很多，差了一个数量级啊！！！郁闷！！！！
E.E)u0q$N& &&
&请各位高手指点一下，哪里需要修改啊？？？
~仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent&
&本人非常着急，求各位高手帮帮忙。
钢筋应该用beam单元吧，truss不能弯曲。
你载荷的单位是是牛顿么？不要用公斤。
网格要加密几倍才行。用C3D8I单元，不要用C3D8R或C3D8. 塑性问题最好不要用二阶单元。
INP文件中以下节点的z坐标3.0后面都多了一个小数点：
nodes: 31 131 231 331 431 531 631 731 831
nodes: 16 16 7016
实体单元的截面力/弯矩/转角
问：求助：请问abaqus里面怎样看一个构件截面（如：钢骨混凝土压弯柱）的内力啊
请问：SRC柱模拟后，如何提取截面内力:如某一截面处的轴力、弯矩、剪力等内容，谢谢。
你的 axialandsideforce1.inp中的
*surface,type=cutting surface,name=halfheight
*surface,type=cutting surface,name=halfheight
60.,50.,500.,0,0,1
意思是截面通过点60.,50.,500., 法向为0,0,1。分析后得到的
axialandsideforce1.dat文件中有
& THE FOLLOWING TABLES ARE PRINTED
SOFM&&&&&&&
SOF1&&&&&&&
SOF2&&&&&&&
SOF3&&&&&&&
SOMM&&&&&&&
SOM1&&&&&&&
SOM2&&&&&&&
1.0002E+04& 1.0000E+04&
-200.0&&&&&
1.6233E-09& 1.0052E+06&
2.0000E+04& 1..9879E+04
在dat文件中看到
&***ERROR: SURFACE ASSEMBLY_HALFHEIGHT COULD
NOT BE GENERATED BECAUSE THE GIVEN DATUM PLANE DOES NOT INTERSECT
THE ELEMENT SETS SPECIFIED UNDER THE *SURFACE, TYPE=CUTTING SURFACE
应该在 60.,50.,500.,0,0,1下面加一个空行，即
*surface,type=cutting surface,name=halfheight
60.,50.,500.,0,0,1
*End Assembly
见 ABAQUS Keywords Reference
Manual&&&&
*SURFACE：
Data lines to define a surface using a plane cutting through the
given element sets when the TYPE=CUTTING SURFACE parameter is
Second line:
List of elements or element set labels to be cut by the cutting
plane to generate an element-based surface that is an approximation
to the cutting plane. A blank data line can be specified to
generate a surface by cutting the whole model.
有两种方法：
方法1： 在ABAQUS/cae中，菜单model / edit keywords，添加上述*surface和*section
print语句，然后在ABAQUS/CAE中submit job.
方法2： 手工修改inp文件，然后在WINDOWS中点击 [开始] → [程序] → [ABAQUS 6.5-1] →
[ABAQUS Command]，在ABAQUS Command窗口中输入以下命令：
abaqus& job=inp文件名称
问：请教section print的问题
print可以求得实体单元上的弯矩。比如我求桩上的弯矩，inp文件的命令流如下：SimWe仿真论坛Kx&
*Section Print, name=xixi, surface=ding, axes=local
XS)^,j&e!H$D:A g8j
|仿真|设计|有限元|虚拟仪器仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent4k.H7F
||仿真|设计|有限元|虚拟仪器T J@qg4_*b
%MY"{QA3G'I?3OSimWe仿真论坛SOM SOF分别代表the tatal force和the tatal
moment。从帮助中可知：output variable SOM is computed with respect to the
origin of the coordinate system considered.
q"{:dw仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent我可以求出实体单元桩上的弯矩，但是不知道abaqus
是如何求出的。是如何从应力求出截面的弯矩来的？
应该是按照材料力学的基本原理吧：
弯矩 = 总和 （各个单元的截面力 * 受力点到局部坐标轴的距离）
问：进行了桩土分析计算，结果是应力与变形，不知如何反算桩身轴力和侧摩阻力？敬请各位版友赐教，谢谢！
在本版搜索"轴力"。
侧摩阻力： history output中有CFS（摩擦应力的合力），不知是不是你想要的。
solong答：我做桩身轴力曲线是先用
*surface，type=cuttingsurface在桩身上定义截面，然后用*section
print输出这些截面上的sof（具体请参考手册）
桩侧摩擦力我也是输出的整个摩擦界面的sof
试一下看看。
问：实体单元模拟桩，如何才能输出弯矩？还有土抗力？
先用 *surface，type=cuttingsurface在桩身上定义截面，然后用*section print
输出这些截面上的SOF（截面力）和SOM（截面力矩）.
问：我关注的模型是个环型的砼基础,里面配有环向的钢筋,混凝土用的是实体单元,钢筋用的是杆单元,现想知道杆单元沿环向的力以及实体径、环方向的力,请问大家怎么才能得到,谢谢了.
Beam的轴力: 在step的field output中加上SF。
实体的截面力: /forum/viewthread.php?tid=702302
问：想输出桩身的转角位移，特别是泥面处的转角？
梁单元： UR
实体单元： 可以试试在实体上tie一个很软的梁，然后看UR.
问：最近分析一个软土固结算例时遇到一个问题，提交任务后计算只完成每个step的第一个increment就结束了，结果明显有问题，我想问问怎样设置才可以完成全部step的所有increment？以前从来没遇到这样的问题，不知哪位大侠出手相助，不胜感激！
1。msg文件里有警告信息：
&***WARNING: SOLVER PROBLEM.&
ZERO PIVOT WHEN PROCESSING NODE 324 INSTANCE&
PART-1-1 D.O.F. 2
& Typically a zero pivot is associated with an
overconstraint in the model.
&However, ABAQUS was not able to indentify any
chains of constraints that might
&lead to an overconstraint at this node. This
suggests that the model might be
&in fact insufficiently constrained in this
particular degree of freedom.
&Please check the constraints used in this
意思是节点约束不足。后处理中要把scale factor 设为1e-35才能显示，可以看到那些约束不足的节点出现了very
very大的位移。
2。 只算一个increment可能是因为 *Soils, consolidation, end=PERIOD,
utol=8000 中的 utol, 见 ABAQUS Keywords Reference
Manual& *SOILS：
Set this parameter equal to the maximum pore pressure change
permitted in any increment (in pressure units) in a transient
consolidation analysis. ABAQUS/Standard will restrict the time step
to ensure that this value will not be exceeded at any node (except
nodes with boundary conditions) during any increment of the
3。 塑性问题不能用二阶单元，见 Getting Started with ABAQUS&
10.3 Selecting elements for elastic-plastic
problems。你的固结可能也不应该用二阶单元。
原创, ABAQUS/CAE划出无限元
先做一个几何体,然后采用swept方式划分网格,扫略方向自然是朝无限远处延伸仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM[1]@7iN*ba2I[
Z‑F@X
`zq仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM画好单元后生成mesh
V[1]g#X仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent把最后一层准备做无限元的做成个set,单元类型改成C3D8仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluentd0Jq5g
_k)y%K'[*E
其余的采用C3D8R
BJH k#J6z$Enj?然后写出inp,把单元类型更改,再读回去就可以了
楼 中说的节点编号顺序不同的问题吗？
用swept自然就这样给你编好号码了
jacobi& 答：
无限单元可以建立起来了，按大家的方法，在CAE中用C3D8R单元进行mesh，生成.inp文件后，
将需要设为CIN3D8单元的elements选出来，改其属性为CIN3D8
但是CIN3D8单元与C3D8R单元的node排列顺序不同，转变时必须 6] z)_&
手动修改，
原创：无限元建立方法，希望得到加分
这是我探索出来的关于无限元单元建立的过程，其实无限元建立最难的不是INP文件怎么通过*element去改变单元类型，而是如何把单元节点转换为与无限元单元节点一致的情况，也就是一个单元的最后几个节点要在同一侧（外侧）。
r^4RCQo~&zlSimWe仿真论坛&
以下是我建立的方法，如果能得到各位的同意的话，希望能得到加分肯定。
通过TOOLS-----Query看到你所选单元的NODE排列号码（首先在VIEW--PART DISPLAY
OPTIONS--MESH下选择show node labels）可以看到，所选单元node号为, 5434,
, , 5147.node 5148,
33这几个单元靠后的节点是在靠外的一侧。但是你需要在INP文件里把NODE靠后的单元号放在后面。
当不是这种情况，你想把单元node号反过来，即把nodalconectivit改为 , 5148,
, , 5433.则可以通过修改inp文件。首先建立JOB点write
input,然后用ULtraEdit打开建立的XXX.INP。如图。可以看到我想建立无限元单元集合为单元.则利用ULtraEdit把单元所有的节点编号复制下来新建一个文件，然后把NODE号改变，再粘贴覆盖原来这一部分，则已经把NODE号改变完成。可以看图7和图8的区别
最后别忘记了，对我想要建立的无限元部分进行单元类型修改。如下图。然后再通过以上查看节点的方法可以看到同样一个单元的NODE号把靠后大的号调了一个面。到此改变已经成功，单元类型也已经改成CIN3D8（三维无限单元）。谢谢，有不对的地方希望大家指正。5nqb-m
岩土本构模型
【讨论】在ABAQUS软件中关于岩土本构模型的选用
在ABAQUS中进行岩土工程问题的弹塑性计算时，鉴于常规地质报告所提供的参数，我们常常选用D-P或者M-C本构关系，而岩土材料中有两种典型的材料，一类是没有粘聚力或很小的砂土类，一类是有粘聚力的粘土类，对于后者我们采用任一本构关系均可（相对计算的收敛性而言），但对于前者，如果选用M-C，很多情况下难以收敛，这也是正常的，这是因为这两种本构关系在理论上都是针对有粘聚力的岩土材料，我认为这种情况下选用D-P模型可以解决这个问题（限于ABAQUS软件），因为ABAQUS软件中D-P模型参数的选用完全可以避开C值影响，在硬化参数选用时可以在单轴抗压强度、单轴抗拉强度、纯剪时粘聚力d的三者中任选一种办法，经过算例验证，我认为选用单轴抗压强度更好些，对于单轴抗压强度可以通过C、φ换算而来也可以以地质报告中提供的地基基本承载力作用参考，也就是说粘土类材料可以在M-C和D-P中任选，砂土类材料可以选用D-P，不知各位同行有何感想。
以上网友发言只代表其个人观点，不代表新浪网的观点或立场。
