您所在位置： &
&nbsp&&nbsp&nbsp&&nbsp
动态系统的鲁棒障检测方法研究.pdf56页
本文档一共被下载：
次 ,您可全文免费在线阅读后下载本文档。
文档加载中...广告还剩秒
需要金币：180 &&
动态系统的鲁棒障检测方法研究
你可能关注的文档：
··········
··········
南京理工大学硕士学位论文 动态系统的鲁棒故障检测方法研究 摘要 随着生产水平和科学技术的不断发展，现代控制系统的规模日趋大型化、复杂化，
现代工业过程对系统安全性、可靠性的需求日益提高。最近三十年来，故障检测与诊
断 FDD 算法及其在工业过程中的应用受到了人们的日益重视。本文主要研究了含
有不确定性的线性系统和非线性系统的鲁棒故障检测与诊断方法。主要内容如下： 首先，综述了动态系统故障检测与诊断技术的研究背景与发展现状。研究了基于
Beard故障检测滤波器 BFDF 的故障诊断方法，该方法针对理想情况下的线性系
统可准确地实现系统故障的检测和分离；其次研究了含有不确定性的线性系统的鲁棒
故障检测方法，该方法改进了Beard故障检测滤波器方法，提高了故障检测的鲁棒性；
最后进一步研究了故障特性的辨识与跟踪问题，探讨了一种利用神经网络对非线性故
障模式进行在线逼近的鲁棒故障诊断方法。通过在线估计器和非线性滤波器实现了对
故障的实时跟踪和估计，并对算法的鲁棒性、灵敏性和稳定性问题进行了分析。并且
通过数字仿真验证了本文所研究方法的有效性。
关键词：故障检测与诊断，故障检测滤波器，非线性系统，鲁棒性，RBF神经网络
Alratract 硕士论文 Abstract WitIl ofthe andthe of the technologyimprovementproductivity,the development
modem are andmore control becominglarger complicated．11ledemand胁hi曲 systems and inmodem
increase detectionand
reliabilitysafety industry
rapidly,Fault their beenattractedmuchmore and have atcenfion
diagnosis FDD flgofithmsapplications
overthe this forlinear thesis，fault diagnosisstrategies thirty detection锄5l past
years．In with are
andnonlinear uncertainties syste
正在加载中，请稍后...ABAQUS常见问题汇总&-&2.0版（25-29）
云木斋开光佛珠手串奉请官网：
请教：关于abaus热分析
各位高手，+HH-wj1m(}y/Q
利用abaqus模拟一个颗粒增强复合材料的降温后产生的残余应力，因为是放在空气中自然降温，能通过定义边界热对流系数使其降温吗？不定义空气这个part可以吗？如何定义？用什么单元进行计算？用3dstress单元行吗？
不用定义空气
_2[?仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM'I%d(M#@f!nN
可以定义表面的热膜换热,这个在contact property里film condition
直接在CAE底下操作，注意选择 coupled
temperature-displacement类型的step,仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM,gB[1]_*T|
定义好film参数就可以分析了。
&R]#_!vD‑~||仿真|设计|有限元|虚拟仪器还有，注意时间要设置的足够长，默认的1秒时间肯定是不行的。5|
Z5w V0hv.ea
b2~c*?4P7Z&Z&M;t
我在做一个机械零件淬火（放在水或油中间冷却）模拟，用该方法可以运行的。
2w$mNJ0Tc,HXSimWe仿真论坛
& M)Z7~(W^8K我的模拟难在要自定义材料属性。。。 &
谁说film只能用于2D模型？3D底下你定义好 *surface不就可以了？
我用coupled
temperature-displacement的step试了一下，是可以计算，只是再请教一下。h
我在计算时film coefficient这个参数设置为1000，sink
temperature设置为293，模型内部温度为811，step时间是10秒，结果计算了1秒，整个模型温度就都降到了293。是不是与定义的film
coefficient的大小有关，这个参数主要影响什么？
这个参数当然很重要了，它是与外界热，能量交换的一个速率，manual当中有例子，可以选他们的试一下，不过注意单位，例子很有可能用的是英制单位。仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM@:ke$MO
ISimWe仿真论坛还有，好像这里温度不用absolute
temperature吧，用celsius就可以了。具体查一下convection。
还有要强调的一个问题是，你可以将模型剖开，很可能表面温度已经全部降下来了，但是里面还是高温。SimWe仿真论坛l&
R!e2w& & &
仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,FluentJ1C|'j!KXu‑K
这里你要正确设置， conduction coefficient。
再请教一个问题，我试了一下，热分析中film
coefficient定义成500和1000结果没什么差别，这是什么原因？这里面有什么临界值存在吗？这个参数是根据什么选取的？一般取值是多少？
稳态传输的话,这个系数自然没什么关系SimWe仿真论坛z,~&Ci;~T
SimWe仿真论坛D(S C
时间无限长,啥都和外部相同了,有限元计算,重要的是物理模型的正确建立
我在作热传导时碰到的一个奇怪问题！
我做一个简单的非稳定场的简单计算，外荷载只有一个面热流。结果发现作用在面上的热流只能影响最外面一层单元的温度分布。我在论坛中也找到了相关的例题，但是我具体比对之后没有发现我跟他的设置有什么不同，但我的模型就是不出正常结果。
修改step的参数，把 response 由默认的 transient 改为 steady-state 就可以了。见ABAQUS
Analysis User's
Manual&&&&
6.5.2 Uncoupled heat transfer analysis
但是，你那样处理后算得是稳定温度场，从结果中可以看到：只要热流一加上去瞬间就影响了整个模型，这是不符合实际情况的。下面是我以前下载的用于瞬态温度场的计算输入文件，他用的就是transient。明显可以看到热流对模型温度影响的区域逐步增大。
仔细比较，总会看出区别。
你的单元尺寸比例题的大1000倍，可能是你们的长度单位不同，相应地， 密度、热传导率、 比热容、heat
flux等单位也应该不同，见
我把你的模型尺寸缩小1000倍，所有密度、热传导率、 比热容、heat
flux等都改为和例题一样，结果温度变化的区域是逐步增大的，如果把step的time period
由0.1改为100，温度变化的区域就会继续增大。
问：我一直用的都是国际单位，但是由于我的模型尺寸很大所以以前算的结果总是不明显。今天听到您介绍了time period
中时间的实际意义才恍然大悟，对我那么大的模型几秒钟怎么能够有什么结果呢 ，后来我将时间大幅度提高就得到了的较好的结果。
答：你怎么改的？我试了只把time period 改为100000,
别的都不改，也还是只有三层单元温度变了，这也太慢了吧。见附件。
问：我跟你差不多，但是由于我以前的模型是两米一个单元并且混凝土的传导率只有1.5~2.0，并且我们用的还是以前例题对于钢材的热荷载。在这种情况，你在一天左右的时间里影响了6米多深的一层混凝土，我觉得比较符合现实。仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM5R"w,fKRRV
我算的很你差不多，但是我是按例题的方法计算完后再加个空荷载步。好像在这个空荷载步里温度影响范围扩散的比较快。不知道例题是为力模拟现实情况还是给这个空时间步让温度的影响范围得以一定扩散？
在做热冲击的应力分析中也做过类似的计算.和ABA不同的是,我设定的是当温度到达稳态的条件下,自动中止计算.keyword和参数如下:
& j)},h+ra!C**Heat Transfer, end=SS,
deltmx=10.
A!c3a‑h
k仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent0.001,
100000., 1e-05, 1., 0.5
在静力学分析的时候一般使用,N,mm,Mpa单位制,设计能量,质量的时候换算确实很麻烦.现在研究所用一套材料属性的数据库,里面有单位换算的部分,方便了很多.
问：请问要怎么把已算出的温度场结果代入下个计算，得到热应力啊?
我是做热应力分析的，现在温度场算出来了但却不知道如何算出热应力，麻烦高手们帮忙解答||仿真|设计|有限元|虚拟仪器"e%V&{,{‑Fe
Load module
3F8h6C!|/o[1]create field
,r/W7Q q([& & &
D0az{x仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent在edit
field对话框的distribution中选择from results or output database file
[1]h[DEm‑N)j9x;C+|仿真|设计|有限元|虚拟仪器然后选择前面完成的温度场分析.
#Qwdl2t)TSimWe仿真论坛SimWe仿真论坛
P/W.Ga&F6_
需要注意的是,在热应力分析的step定义中,time period尽量和前面的传热计算的time period一致.
问：請問一下～～～有關定溫加載的問題
一開始試體於常溫下加熱，之後升溫到一定的溫度再開始加載，是不是在step那邊利用"coupled
temp-displacement"可以達到此效果。
I,m`7d^ +jH
還是要先分析熱作用後，將熱分析結果RESTART輸入到力學加載
coupled temp-displacement是说温度场和位移场相互影响（变形产生热），不适合你的问题。
Restart是对的。
问：怎么读入热传导的计算结果！
如果要将第一步计算所得的温度场作为下一次计算的初始条件，怎么操作？~(T_*mk'gl
&&&我的目的：第一次算出温度场&
&第二次 做应力计算
,\8Br9A.h$tSimWe仿真论坛&
&&&我是在第一计算完了温度场以后，在field
output里面讲节点温度写出来（*.rpt)!
Y#T.H7~ |仿真|设计|有限元|虚拟仪器&
&&&第二次计算应力（我做的是一般的应力计算）时，将第一步的写出的温度结果文件作为温度场定义（选择form
result file ).3z[1]g5C&
&&&但是提交job检测时提醒我说
我指定的文件（*rpt) 可能不在（could not be localed in
)我指定的文件中。然后就自动中止了程序。仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent
JN%w\6xr6Tx0D
&&&我不知道是不是我写出结果的形式或者是读入结果的设置错了，请各位大侠给点意见，告诉我该具体操作，怎么先谢谢了！
读入温度场时，只能是使用传热分析的结果文件 .odb或 .fil，不能用rpt.
生成.fil文件：&&
问：我的轧制模型需要用到前一步计算的温度场结果，可是怎样在程序中生成.fil文件呢？
gy_tian答：
在*end step前加上关键字：
*node file
答：ABAQUS/CAE 不支持*node file，见 ABAQUS/CAE User's Manual A.1节
"ABAQUS keyword browser table"：
*NODE FILE
U ABAQUS/CAE reads output from the output database file
我在作热传导时碰到的一个奇怪问题！
我做一个简单的非稳定场的简单计算，外荷载只有一个面热流。结果发现作用在面上的热流只能影响最外面一层单元的温度分布。我在论坛中也找到了相关的例题，但是我具体比对之后没有发现我跟他的设置有什么不同，但我的模型就是不出正常结果。敬请各位高手指点！仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM)\2c[1]E^0|0s2N
& 我的材料参数
只填了密度&&热传导率&&比热容，具体数值见模型。4{Z
L,u&P& & &
下面是这个算例的inp文件。仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,FluentI
P:T9^R[1]iN._#bC
修改step的参数，把 response 由默认的 transient 改为 steady-state 就可以了。见ABAQUS
Analysis User's
Manual&&&&
6.5.2 Uncoupled heat transfer analysis
问：step 1 钢水在铸模中冷却，step 2 冷却后的钢不再与铸模作用，而是受水冷作用。 我在step2
中对铸模与钢水之间的接触设成 inactive，其辐射也是inactive，但是从计算结果看，铸模还在升温，钢板在step2
中的冷却有点问题，怎样能在把铸模和板的接触取消呢?
看了你的模型。我的理解是： ABAQUS不是立刻去掉接触，而是在step 2中让heat
flux逐渐降至0，你用的是transient 分析，且step 2时间很长，因此铸模一直在升温。解决方法是：在step 2
之前添加一个时间很短的step（例如0.01秒），专门用来去掉接触。
我试了一下这种方法，slab温度总体上是一直在下降的，局部有些上升，不知是否正常现象。
参见ABAQUS Analysis User's Manual 的 21.2.4 Removing/reactivating
ABAQUS/Standard contact pairs：
问：请问用Abaqus/cae建模型的时候，怎么样可以使一个立体或平面在温度升高时伸长相应的值啊？
以下是我的做法：
1。在material 里选expansion,使Expansion coefficient Alpha
为0.03，然后reference temperature 为0（事实上我不知道这个reference
temperature有什么用的。）
2。在load里面建立一个field，范围是整个普通的长方形平面， 定义为Temperature. 在initial
时magnitude为0，在 step-1 时magnitude为1,性质是uniform。
因为Alpha=3beta, beta=1%, 所以我预期长方形的每条边都伸长1%，
但事实上是伸长了10%。我就多次改了temperature field， 令step-1 = 0.1, 0.001,
0.00001...但每次都是伸长10%.即使我是用不同的Expansion coefficient Alpha = 0.1,
0.00001... 结果还是伸长10%
请问我要它每条边都伸长1%,应该如何做啊？
注意到ABAQUS把增量步降到很小才收敛，而这样的简单问题应该一步就收敛的，所以怀疑模型有问题。
你的弹性模量 E 是 1.6291e+11，太大，导致了ABAQUS结果错误（可能是数值误差方面的原因）。将 E
改为1.6291e+5 （即单位为N/mm2）,
问题就解决了，两个方向各伸长1%，与你的热膨胀系数Alpha=0.01相吻合。你说的Alpha=3beta是什么？泊松比吗？热膨胀的两个方向都遵循ALPHA。
所以E 的单位应该用 N/mm2，相应地模型尺寸也应该是mm。你模型中的其他方面都正确。
ABAQUS中的量没有单位，用户自己保持统一即可。建议尺寸用mm，弹性模量用N/mm2。
Pa = N/m2, MPa = N/mm2, 162GPa = 162000 N/mm2.
ABAQUS中的量没有单位，只是纯粹的数字。比如模型长度是1m，你希望用mm做单位，就输入1000,
ABAQUS就知道长度是1000，它不关心这是m还是mm.
问：我从其它途径获得了模型所有节点的节点号和对应节点温度结果，并做成了Excel格式，现在想把这些温度结果导入到abaqus中，以云图的形式显示出来，不知道能不能办到？我试着以边界条件的形式添加，但没有成功。
现在手头上已经编写了inp文件的前一部分，包括节点和单元信息，后面的温度结果导入就不知道怎么写了，请帮忙。
*Initial Conditions, type=TEMPERATURE
节点1， 温度
节点2, 温度
*Step, name=......
问：用一个刚体做恒定的热源，用temperature field定义温度，可是出现错误：*temperature may not
be used with elements that possess temperature degrees of
不要用temperature field. 用边界条件，类型为temperature.
问：板料初始温度高于室温，模具温度始终保持室温，板料热冲压过程中会与模具发生热传导，那么该怎样定义边界条件呢？要定义边界的冷却曲线吗？
模具的温度在那里定义呢(因为将来还要计算模具的温度分布）?
板料的初始温度呢？
模具也要作为part出现在模型里。
温度场：load模块，菜单field / create, 选other, temperature
．动态分析
静态分析还是动态分析
想用abaqus解决一个 等速旋转圆盘 的问题,计算结果老是不对,困惑了好几天了,请高手帮帮忙.SimWe仿真论坛%s
内径a=1.0m，外径b=2.0m，厚度0.1m，材料为steel,
E=200e9,v=0.3,ρ=7800仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMfc%~#QQ#|A.p
角速度ω=100rad/s，求等速旋转时的应力分布
: [切向方向应力理论公式]
( 10:46, 43.02 K)
你的公式是静态问题的解，离心力作为静态载荷施加在物体上。静态或动态并不取决于物体是否运动，而是说是否考虑冲击和局部效应。Getting
Started with ABAQUS 图13.1、图13.4、图13.5 比较了静态和动态问题的区别。
你用explicit分析，相当于离心力是一个冲击载荷,
得到的应力比静态分析大很多，而且应力场不是轴对称地均匀分布，随时间推移而变化。
应该用static分析步，load type为 rotational body force.
只取1/4来建模，在三个面上分别施加关于x,y,z的对称边界条件（不模拟转动过程，只是施加离心力），或者直接建2D轴对称模型也可以。
要提高应力结果精度，单元类型不要用C3D8R, 而是用C3D8I，并且细化网格。附件是我的模型，和解析解很吻合。
请教问题：关于齿轮啮合过程分析的
aba_aba大侠及各位朋友！这几天试验了几次，同样不成功，主要是变形实体单元的角速度没有找到合理的施加方法！我都是在boundary
condition里面加的，今天试验了一下将主动齿轮部分设置为刚体在参考点上加角速度，结构刚体与变形体处发生严重扭曲，结构根本不能用！现在新的cae模型在邮箱里面的主题为“new”的下面，该模型我是将主动轮转化为柱坐标，在T方向加的速度209.44，还用了一个幅值曲线，我自己都感觉不行，单位都不是radian/s,但是实体单元不能施加VR3呀！各位使用该模型的时候要施加正确的角速度，请指教指教吧！我学abaqus不久很多问题都还不懂！还有就是显式计算之后的判断是主要依靠总能量（ETOTAL)看其历史曲线是否不变吧？
可以用速度场：load模块，菜单field / Create （定义在主动轮的所有节点上，而不是参考点上）。见
/forum/viewthread.php?tid=710486
这是很典型的explicit和standard的区别，动态分析在突然加载时会出现部变形，见Getting Started
with ABAQUS 图13.1、图13.4、图13.5 的比较。参见
/forum/viewthread.php?tid=733906
如果不像考虑冲击载荷带来的局部变形，应该进行准静态分析（见Getting Started with ABAQUS
第13章），其关键问题就是如何选择合适的加载速度、分析步时间、质量放大等参数，使结果尽量接近静态分析的结果。或者干脆用standard做静态分析。
显式计算之后要看总能量（ETOTAL)是否不变。
可以先做个二维的试试，节省运算时间。
ABAQUS Benchmarks
Manual&&&&
1.3.11 Frictional braking of a rotating rigid body 的例子可以参考。
急切请教问题：齿轮啮合过程的（麻烦大家了）
最近一直在做这个问题，前面已经发过贴子。曾得到过aba_aba大侠的回复，非常感谢！但始终都没有解决好问题。我试验了好多次，居然出现了齿轮接触时一个齿穿过了另一个齿。我接触定义的是显式通用接触（all
self),阻碍力矩也不大，真是百思不得其解。现将avi放在了邮箱里面：账号：&&密码：cae000
主题为“动画（avi)”的那个文件。由于我是加了一个初始应力场，那个文件又太大，上传的那个cae模型可能不能算，可以看一下那些加载。或者是把初始场去了看大家有没有好的方法也能做出正确的结果，不甚感激！急着用呀！谢谢各位。
我这里计算还没结束，先说说我看到的问题：
- field和速度边界条件中的角速度单位都是 弧度 / 秒，你输入的好像是 度 /
秒。转得太快了，冲击效应大，可能导致异常。
- 边界条件和载荷都应该使用smooth类型的amplitude, 以避免冲击效应。见Getting Started with
&&&13.5 Example:
forming a channel in ABAQUS/Explicit：
“For the first attempt of this metal forming analysis, you will use
tabular amplitude curves with the default smoothing parameter for
both the application of the holder force and the punch stroke.
- 内壁和参考点的连接不要用rigid body约束，应该用coupling. 把coupling
type设为distributing, 内壁就可以发生变形，可能就不会只是内壁旋转而造成扭曲。
- 应力沿轴向是均匀的，在轴向划分5个单元估计就够了，节省计算时间。
我只算了你的caemodel.cae.&&穿透应该是因为应力应变太大：齿根应力高达10e5
MPa, 这是不可能的（钢淬火后的断裂强度不过是1500MPA）.
接触的齿也出现了明显的大变形（只显示一个齿轮就可以看得更清楚）。
应力应变太大的原因是受力太大。阻碍力矩 10000 Nmm 施加在半径40mm的齿轮上，相当于受力250N,
应该问题不大。再看转速，小齿轮的初始速度场是 -100 弧度/秒，但速度边界条件是
-209.44弧度/秒，相当于让小齿轮在一瞬间速度增加了一倍多，其冲击载荷估计是非常大的，所以应力应变太大。
只是内壁旋转而造成扭曲，我觉得原因还是速度变化太剧烈，并且amplitude没有用smooth。如同 Getting
Started with ABAQUS 图13.1所描述的，加载速度太大就会造成局部变形。
你现在在小齿轮上同时定义了initial field和速度边界条件，只要二者的速度值是一致的，就应该是可以的。
;P 附上我的二维模型，转得挺好看的。二维的好啊，8分钟就算完了。要点：
- 模型为平面应力，而不是shell (发生面内弯曲的才是shell). 单元类型为CPS4R.
- step属性中的nlgeom要设为on （大变形，几何非线性）。
- 接触类型为surface to surface （二维好像不能用general contact）
- contact property中没有设置任何参数，即默认的摩擦为0。没有用你模型中的geometric
properties.
- 内壁和参考点还是用kinematic coupling的好，如果用distribute
coupling，内壁变形太大。
- 载荷和速度边界条件都使用smooth类型的amplitude：
时间 0& ， && 幅值
时间 0.00001 ， 幅值 1
- Field output 的interval 设为400，以看清运动过程细节。interval
越大则变化过程显示得越清楚，但odb文件也会更大，计算时间更长。
比较二维模型和你的三维模型运行后生成的sta文件，看到里面的模型总质量基本相同。二维模型应该与三维模型的分析结果相差不大。
请教aba_aba,弹塑性体的旋转问题？
我是左旋压的，工件需要旋转，可是我把工件和模具tie后，工件的tie表面扭曲变形了。仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMV9R
考虑用field来计算，可是那是初始速度，随着旋压力的加大，速度会越来越小，甚至不转了。仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM#l4x4PL]‑f
有没有好的方法来实现弹塑性体的旋转？
如果工件和模具之间没有相对滑动，可以用tie，否则要用带摩擦的接触。
关于局部扭曲变形，见：
/forum/viewt ... 0823&pid=821705
第12，16楼
/forum/thread--1.html&&第7，9，11楼
如果工件在整个模拟过程中的转速都是已知的（而不是根据外载荷由ABAQUS的），我认为可以使用上面帖子里sqf0834使用的方法：&&
用field定义初始角速度，并且用速度边界条件定义整个模拟过程中的角速度（要使用smooth类型的amplitude，不要让速度在很短时间内发生太大的变化，否则会出现局部扭曲变形）。
請問這個錯誤...
請問前輩，這個錯誤應該如何改正...&
&仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,FluentUF@|J
&& & THE OPTION
*BOUNDARY,TYPE=DISPLACEMENT HAS BEEN USED; CHECK STATUS
仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM1Yf
O1B\)_1G3j
BETWEEN STEPS FOR WARNINGS ON ANY JUMPS PRESCRIBED ACROSS THE
$k-OO?8m仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM&
&& & STEPS IN
DISPLACEMENT VALUES OF TRANSLATIONAL DOF. FOR ROTATIONAL
Pt@,LG-`(T4B2mSimWe仿真论坛&
&& & DOF MAKE
SURE THAT THERE ARE NO SUCH JUMPS. ALL JUMPS IN
%n.P5aK$z3EC2I&
DISPLACEMENTS ACROSS STEPS ARE IGNOREDSimWe仿真论坛-w:xle)z'@
小弟把Model從STANDARD改成EXPLICIT之後就跑出來了...仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM_*L5^9V&_
可否請大大解答一下...
Tf$X仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM謝謝！～
shock111答：
你用指定displacement,
然后用的是STD默认的方式,这个和EXP不符合SimWe仿真论坛n2xq\,t&
修改amp就可以了
即在动态分析中，施加载荷的时候不要使用默认的Amplitude，而要自定义一个逐渐增大的Amplitude，并使用smoothing
parameter，见 Getting Started with
ABAQUS&& 7.5 Example: cargo
crane—dynamic loading，Table 7&1。
问：我建了JOB:K30-1运行如下:仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM
L/u(Y2`(^VUXk(h
The job "k30-1" has been created.
B_8A`9OThe job input file "k30-1.inp" has been submitted for
analysis.SimWe仿真论坛 so%M/N,R)D
Job k30-1: Analysis Input File Processor completed
successfully.SimWe仿真论坛V)TV[1]i9a0O
然后到这里就一直没有下一步了，是怎么回事,是由于电脑太慢了吗?
运行得挺好的，explicit分析的过程信息显示在sta文件里，msg文件什么都看不到。
我的机器很快，10分钟就正常运行结束了。
你为什么不用standard分析步？就不会那么慢了。接触部位的网格应该细化。
问：我做了一个椭圆物（多层构成）压平面的动态模拟，在椭圆中心加载荷，出现奇怪的变形结果。INP文件如下，诚恳各位高手的指点。
我的理解是，你的模型应该按照准静态问题来分析，见Getting Started with ABAQUS
13 Quasi-Static Analysis with
ABAQUS/Explicit，其中的13.1有一个电梯的例子，解释了为什么只在局部发生变形。
问：请大家讨论轧制模型！
轧制模型，不能正常模拟，请大家讨论是哪里有问题？
你是问为什么没有咬入吗？你的单元长度是1，初始速度是0.06, 几乎为0。把初始速度改为6，就咬入了。
另外你的密度数量级太大，是你故意做的质量放大吗？一般要保持单位统一，如果长度单位使用mm，那么质量的单位必须是“吨”，相应的密度单位就是吨/mm3，应力的单位是N/mm2（即MPa）。以钢为例：
210000., 0.3
如果用国际单位N, m, 以钢为例, 应该是：
2.1e11., 0.3
&amplitude不一致
问：我用＊amplitude定义了一条锯齿波曲线，然后在动态分析步中将此作为边界条件进行梁的动态分析，计算没有问题。计算结束以后，我通过xy-data提取施加锯齿波边界条件的节点的位移时间历程，却发现输出的波形被截波。如果我定义的曲线幅值为10米，则输出为8米，如果我定义为15米，输出为12米。如果将锯齿波的周期加大，则没有出现这样的问题。是不是我定义的边界条件周期太小，还是什么原因？可是我在定义曲线的时候已经设置smooth＝0.0，请大家帮我看看怎么回事。
时间增量步太大，改为0.001就行了：
0.001, 5.0
我的理解是： 如果对U1上影响较大的最高阶固有频率是f, 其相应的周期就是1/f, 你的时间增量步就要小于1/f。
Getting Started with ABAQUS&&
第7.5.2节介绍了如何在*Modal
dynamic分析步中确定这个最高阶固有频率：在*Frequency分析步中要保证提取了足够数量的模态，即在主要运动方向上的总有效质量要超过模型中可运动质量的90%
更正一下，不用改时间增量步，在step中加上下面一句就可以了：
*Output, field, variable=PRESELECT, frequency=1
因为*DYNAMIC默认的输出设置是frequency=10，正好跳过了位移的峰值。这也就可以解释为什么你把位移边界条件的周期加大就正常了。
问：我想請問一下我現在做的是一個一維自由度的彈簧和阻尼連接在一個物體的系統
我們知道當物體連接只有彈簧時連接所得到是自然頻率
當跟阻尼時得到的憑率會是另一種因為阻尼而改變的頻率
我結果檔都出現自然頻率
請問一下要加什麼指令才會出現另一種頻率
Frequency分析步中不考虑阻尼，要在它后面加上一个complex frequency
分析步，才能考虑阻尼的影响。ABAQUS Analysis User's Manual 6.3.6 Complex
eigenvalue extraction。另外，阻尼要很大，才看得出对频率的影响。
http://en.wikipedia.org/wiki/Damping
问：很奇怪是我有加damping
跟沒加damping的 u 結果是一樣的
frequency也是一樣的
應該要不一樣才對阿!!
*Frequency分析步中，有没有阻尼都是一样的。*Frequency分析步得到的位移只是一个相对值，最大位移总是1。
如果分析步类型为*Modal dynamic, 你在material里定义的阻尼不起作用，要在step里设置阻尼，例如：
*Step, name=Transient, perturbation
*Modal dynamic, continue=NO
0.005, 1.5
*Modal Damping, rayleigh
1, 30, 3., 0.
见Getting Started with ABAQUS 7.5 Example: cargo crane—dynamic
如果分析步类型为*Dynamic,Explicit, 像你那样在material里定义阻尼就可以了。见Getting Started
with ABAQUS 7.8 Comparison with direct time integration
wswu:& 一个突然反转的大变形问题——
这是一个突然反转的大变形的动力分析问题的练习（真实问题更复杂），在ABA-ABA的提示和指导下，已经达到了预想的结果（再次感谢ABA-ABA
！），只因忙于另外的一项工作而没有及时发此贴。
R~-d)q问题描叙如下：6J5GDQ,V:O$cF8o
某个汽车的钢板弹簧系统在某一载荷的作用下再施加一个加速度，模拟汽车在路面（带沟缝的水泥路）上运动时钢板弹簧的动强度问题。因为钢板弹簧在预载作用下，突然受到一个加速度的作用，因而它是一个突然反转的大变形问题。严重的非线形接触，使得问题的求解难以收敛。
R&}]0O仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent这是力学模型
问：半刚性节点分析时的警告信息处理问题（接触方面）
大家好，我在做一个钢结构半刚性节点的分析，在螺栓预紧力施加完成后的step里increment特别小。总共也就14000多个单元，算了14个小时了才算了那么一点点。我传上了inp和sta文件（cae文件太大）在楼上，请帮我诊断一下模型。谢谢！
答：动态分析就是计算时间很长。你一定要用动态分析吗？只有加载时间很短（例如0.2秒），或分析中不能忽略结构的惯性，才需要采用动态分析。
你的载荷虽然是变化的，但持续时间很长，我觉得可以用普通的静态分析步，模型中的Amplitude和边界条件同样可以用于静态分析步。这样计算时间就很短，各个step的initial
increment设为1即可。
位移好像很小，不用把nlgeom设为on.
螺杆和孔壁的接触有几个问题：
- 你用的是small sliding,
而螺杆和孔壁的距离很大，开始时不接触，运动后也不会接触。如果希望考虑运动后发生的接触，就要用finite sliding.
- 如果用finite sliding不收敛，可能是因为螺杆和孔壁的网格太粗，每90度的圆弧上应该划分10个单元。
．屈曲分析
热情奉献, 屈曲分析总结,适合结构人员
我给屈曲分析总结一下,希望支持!!
: 一起学习梁的屈曲分析.rar (
19:49, 10.93 K)
该附件被下载次数 21
: 一起学习板屈曲分析.rar (
19:49, 11.07 K)
该附件被下载次数 21
问：ANSYS里算非线性特征值时要用到线性算出的结果，也就是把线性算出的结果作为初始扰动加上，然后算非线性的加载、位移曲线，我想知道ABAQUS中怎样加这个初始扰动
chenwei0707答：6.2.4 Unstable collapse and postbuckling
Rik法用于跳越失稳问题的研究，也可以用于分支屈曲的后屈曲研究。分支屈曲的后屈曲分析不能直接在分支屈曲后面研究，而是要给一个初始缺陷，使力学响应呈连续状态（非线性）
7.6.1 Introducing a geometric imperfection into a model
定义初始缺陷
Abaqus用三种手段定义初始缺陷，根据分支屈曲模型取一个线性组合，根据静力分析结果，直接指定。
除非初始缺陷已经知道，一般采用第一种方法。
第一步，特征屈曲分析，Write the eigenmodes in the default global system to the
results file as nodal data (“Output to the data and results files,”
Section 4.1.2).
第二步，将这些特征屈曲模态添加到perfect 几何体中，作为初始缺陷，
where is the mode shape and is the associated scale factor.
一般来说系数w，在第一阶模态最大，而且w一般取结构几何参数的倍数，如壳的厚度的0.1倍，等。
第三步，用rik法进行分析。
ABAQUS通过节点标签来输入初始缺陷（imperfection）但是软件不会去确认两个模型的兼容性，所以要特别注意节点标号的一致性。
Imperfection 定义输入
Defining an imperfection based on eigenmode data
Input File Usage:
&&*IMPERFECTION,
FILE=results_file, STEP=step, NSET=name
Defining an imperfection based on static analysis data
可以基于前一步的静力分析输入初始缺陷，可以定义某个载荷步的结果，也可以不指定，默认输出文件里的最后一个载荷步
Input File Usage:
&&*IMPERFECTION,
FILE=results_file, STEP=step, INC=inc, NSET=name
Defining an imperfection directly
可以通过局部或整体坐标定义，也可以读入相应的文件
Input File Usage:
&&*IMPERFECTION, SYSTEM=name,
INPUT=input file
&&If no input file is specified,
ABAQUS assumes that the data follow the option.
缺陷敏感度的问题上，可以通过改变缺陷的系数（W），用分析结果变化大小来评估。
缺陷较小的结构初始位移变形较小，在极值点突变，而初始缺陷较大的结构，载荷位移曲线较平滑，可以说较容易分析。
ABAQUS计算结构屈曲问题
一个工程结构，受多项载荷作用，需要判断此结构是否结构失稳，看ABAQUS的文档，更多的是通过判断临界载荷来考察结构，或者施加一个小的imperfection，然后利用riks方法判断结构失稳载荷大小。仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM-]7KC:a)E!NO5hE
||仿真|设计|有限元|虚拟仪器9s7^& &
也就是说，ABAQUS更多是确定最终结构失稳载荷，假如对于结构只有一种载荷情况的话，那还好处理，但是现在有多种载荷共存，这种情况下该如何判断结构失稳情况？请予以指点，谢谢。
不必拘泥于以上方法，直接用Static加Stabilization算屈曲也不错，多少载荷都可以
我做了一个薄壳的屈曲和后屈曲问题，前边一步没问题，后屈曲总运行不了，请各位指点
*IMPERFECTION, FILE=fil文件名, STEP=2
其中的 “fil文件名” 不能以数字开头，比如你用的 5.fil 就不行。
问：用distributing coupling在参考点加载后发现计算量剧增
我计算曲壳在横向载荷作用下的非线形屈曲反映,由于不能直接加始终竖直的横向载荷,因此加载参考点,用用答：
看来是 distributing coupling 导致收敛得慢，可以用另一种方法，不需要参考点和distributing
coupling ，直接在shell上施加load，类型为surface traction,
把traction设为general，输入方向vector. 注意这种载荷类似于pressure，是单位面积上的力。
问：求助：收敛及restart问题
附件中的job-1收敛，job-2 和job-3都不收敛，这是哪的错呀？
(AM1[FU!pYSimWe仿真论坛另外，文献中说要对job-2重启动分析12次才能收敛，怎么分析12次？
你的文献P38中说了:
During the progress of top edge displacement it was necessary to
restart the analysis
when not achieving convergence.
The analysis had to be restarted twelve times to reach δ = 3 mm.
Concerning
restarting, there are two different ways to reach convergence:
decreasing the time
increment or increasing the damping factor.
就是说，不收敛时要修改模型参数（decreasing the time increment or increasing the
damping factor），然后从已经完成的increment处进行restart。
restart的方法见 Getting Started with ABAQUS& 11.4
Restart analysis
．材料失效破坏 / 裂纹
(chunyu) 材料失效及其模拟
在ABAQUS里，模拟材料失效主要有两种方法：一是使用已有的失效材料模型，在定义材料时指明最大剪切应变或最大拉应力失效准则，并给出相关的参数，详细内容可参见ABAQUS
Analysis User's Manual /11.2.8 Failure models
pq~WN仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM；另一种方法是使用VUMAT，在该子程序中定义一个状态变量（state
variable）来表征材料是否失效。为0，失效，为1，正常。然后根据自己定义的失效准则，来给这个变量赋值。子程序方法比较灵活，可以定义自己的失效准则和材料本构方程。详细内容参见ABAQUS
Analysis User's Manual&&/25.3.4
VUMAT。下面几幅截图是我自己的一个VUMAT产生的结果，用来模拟聚合物材料如PMMA的剪切及断裂破坏，本构方程及断裂准则均为自己定义，考虑了弹性，粘弹性，塑性，粘塑性。网格比较粗，示意而已。试件单向拉伸。||仿真|设计|有限元|虚拟仪器(e(o
谢谢大家的关注和鼓励。inp文件就不多说了，一个正方形的板，中间钻个洞。
*]‑D
\仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM至于材料破坏的算法，这里采用了B.P.
Gearing, L. Anand / International Journal of Solids and Structures
41 (&3150中推荐的方法。原理很简单:7uQH+Vmj7f
（1）根据应力状态，计算主应力。如果hydrostatic
stress为正（就是受拉了），最大主应力为正且大于一个值（材料参数），破坏开始形成（即initiation），并在最大主应力方向上产生一个塑性变形，文章里给出了一个公式计算塑性变形的大小（记为P1），当然自己也可以提出合适的公式。否则剪切变形占主导，采用一般的本构方程计算应力应变(等效塑性应变记为P2)。仿真分析,有限元,拟s#G&
& & & Rp)d
（2）当p1&p1_cr或p2&p2_cr时（p1_cr,p2_cr均为材料参数），积分点失效，在VUMAT中通知ABAQUS删除就行了（这是原方法），当一个单元内所有的积分点失效后，该单元自动失效。仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,FluentT‑Hj9jN!ij,KF
（3）为了计算稳定，并且更符合实际破坏的过程，我做了一些修改，不是在p1&p1_cr或p2&p2_cr时立即将单元删除，而是定义了一个我自己称之为质量因子（q_factor）的自变量，当p1&p1_cr或p2&p2_cr后，q_factor很快衰减，材料的模量等抵抗变形的参数也随之快速衰减，直至没有能力承受载荷，然后删除积分点，5X~,y~-_
（4）几点说明||仿真|设计|有限元|虚拟仪器7V4A t,Y6}[1]b"|
(a)破坏的initiation,growth直至材料breakdown，基本上没有统一的标准，所以可以根据实际材料的性能自己提出标准；仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM&d1egO)FeZ8h2z‑v_m
(b)这种failure model只能在abaqus/explicit中使用，所以计算效率可能不令人满意；
请问如何在abaqus中实现穿甲的模拟啊？
就是希望看到碎片（也就是失效的单元）飞溅的效果
f我在模型中定义了材料的失效（failure），可是不知道如何释放失效的单元（碎片），因此在穿甲的过程中看不到碎片飞溅的效果，失效的单元虽然变形很大，但弹丸穿过平面后这些单元仍然像橡胶膜一样包裹着弹丸，就是不碎裂开
*r(|V9{!Q+U仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM请众大侠指教，感激不尽！
自己定义的材料失效
|1yb%u!]XK#G在property模块中定义材料的elastic时，添加fail stress
P6o+dy$O&](U仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM不知这样的定义是否足够
我觉得你的材料模型好像有点简单，穿甲过程模拟应该会用到材料的塑性，听做过类似分析的人介绍过他用的是Johnson-Cook材料模型，而且穿甲过程涉及到大变形，所用的分析方法也得考虑好，比如ALE算法。
你要看到飞出去的效果,建议
^]n$rSimWe仿真论坛1:采用3dmax做|仿真|设计|有限元|虚拟仪器(y0lK,ec1A$j
2:采用cohesive单元
可是我看到ABAQUS Example Problems Manual中2.1.3 Rigid projectile
impacting eroding plate和2.1.4 Eroding projectile impacting eroding
plate 这两个例子的inp文件中都没有cohesive单元啊，那么它们的效果是通过什么单元模拟的呢？
E!u8]Y2B Z谢谢指教！！
这些例子是*shear failure的
SHEAR FAILURE在CAE中是不支持的,必须手动修改INP.
e-*SHEAR FAILURE
一般下一行跟的是塑性失效应变,其他参数可参看DOCUMMENT,例子里用到的也就只是*SHEAR FAILURE 1.0,
|+p r:@kY4b就是说塑性应该达到1.0,单元就失效,默认的是ELEMENT
DELETION=YES,不过实际情况下,失效单元仍然显示,只是失效了,不再参与承载等.
(林其米) ：ABAQUS裂纹的例,,互相交流啊
关于crack扩展的问题
在材料里定义了cohesive单元材料的破坏,但是分析的时候好象看不到裂纹的扩展,dava说在display
group里面可以定义，看了下没有发现，高手指点下啊,谢谢
oX8g!k dfSMSimWe仿真论坛display group--&
v‑cE$v!YT
Element value--&
o*Zp:[m[1]B(|仿真|设计|有限元|虚拟仪器SDEG--&SimWe仿真论坛'c)Q,W!nu
1.计算之前在field output里面选SDEG，这个参数在strains里面3L*~
2.Element value 那个选项两空都填 1
在进行裂纹扩展的过程中，我们需要使用*INITIAL
CONDITIONS将裂纹面bond在一起，使裂纹面在不满足断裂准则时，可以连接在一起不致断裂，但是我在模拟的过程中，在没有使用debond命令来启动裂纹面的扩展，按照手册上的说明，裂纹面应该不会分离，但在模拟的过程中，却出现了分离的现象。我已经研究了一个多月了，因为没有人指导和讨论，没有解决这个问题，还希望各位高手给点指导意见，谢谢！下面是我的inp文件。期待回复！！
根据ABAQUS Analysis User's Manual 7.10.3 Crack propagation
analysis：
“If a node set is specified, the initial conditions apply only to
the slave nodes in the node set.”
不属于node set 的部分就是unbonded portion of the slave surface ，“The
unbonded portion of the slave surface will behave as a regular
contact surface”，所以分开了。
同样，debond也只对node set 起作用。
*INITIAL CONDITIONS， TYPE=CONTACT
Slave surface ，Master surface ， node set
node set 是要产生裂纹的区域。slave surface上node
set之外的区域是普通的接触关系，可以随时分开。
怎么还是不明白呢，希望出现裂纹的地方应该定义为node set。
你现在把分界面的右侧定义为node set，这里就是将要出现裂纹的地方，用debond之前这里不会分开。 分界面的左侧不属于node
set，是普通的接触关系，可以随时分开。
所以应该把分界面的左侧定义为node set，右侧的正确处理方法见ABAQUS/CAE User's Manual 20.1
Modeling fracture mechanics 。
问：在帮助文件里发现一个例子，三点弯曲的断裂扩展，跟我所需建立模型一样，我就直接用它的model.inp，然后自己写了后面的，但倒入CAE里时，总说错误，大侠帮我看看哪！十分感谢！
错误信息：AbaqusException: in keyword *DEBOND, file "try1.inp", line
3825: The keyword is misplaced. It can be suboption for the
following keyword(s): step
意思是*DEBOND要放在step里面。你写了：
**STEP,NLGEOM,INC=500
**0.87, 2.87e-08, 0.00287
**BOUNDARY
两个星号后面的内容都不起作用，是注释行，应该为
*STEP,NLGEOM,INC=500
0.87, 2.87e-08, 0.00287
问：我刚刚接触ABAQUS，研究三维疲劳裂纹扩展问题，在建模时出现点问题。我计算时需要很小的网格（1/10微米），在网格细化时遇到了问题，请问三维网格如何做到局部细化？仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM'FFu
我曾经用tie约束来细化，但效果不理想。仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR
mesh模块，菜单 tools / partition
ABAQUS/CAE User's Manual&
&&&20.1 Modeling
fracture mechanics
Cohesive单元
:& 关于cohesive
model一些review
几个月前我做的总结。review of cohesive models in some journals
5wW0X{Ab6L*M9{不过是英文，说不定对学习cohesive相关的朋友们在理论上有些用。
定义cohesive element的厚度为0后出现的问题
大家好仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent4Vq'K!Db)EC
V6.5的文档中有复合材料加筋板拉伸破坏的例子，其中的破坏界面用cohesive element模拟，厚度为0。我也做答：
Re:定义cohesive element的厚度为0后出现的问题
粘结单元的厚度可以用本身的几何厚度、也可指定厚度。
S O J‑Z
|仿真|设计|有限元|虚拟仪器当几何厚度为0时，你可以直接指定厚度值作为计算中需要的厚度值
edit mesh， creat solid layers。
l.帮助文件中有。
关于cohesive单元受力分析中出现的问题，希望高手进来回答下，(有inp文件)
分析发现对于存在cohesive单元平板受力,(1)加载位移后只有加载位移的那一边受力,固定的另一边几乎不受力,或者力很小，而cohesive单元也就是起个胶层的作用，这跟实际不符合!
;Io"oPX(\1h0O
L-e仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM(2)而且其值很小，这个是因为cohesive单元部分的弹性摸量比较小，所以造成同样应变下应力比较小，可是增加弹性摸量后仍是如此，不知道什么原因?高手请回答下，谢谢了
C6[Z&}`8k4ZL3h.W)e&|$P8g
SimWe仿真论坛4Ks6^(g*p&D0j
&&&另外在分析的时候试着用了能量分散率,可是变化太大，掌握不了尺度不说，还是解决不了1问题的,跟实际情况不一样.是不是其中还需要一些特殊的设置,有知道的说下，困扰了半个月了，一直都在做这样的尝试,谢谢了.
_/lr-` |仿真|设计|有限元|虚拟仪器
q,v仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent&
& 我的qq是,希望能够指点下,
^'Q j3Y'D& &&
^T!b[1]R%i'C-w]&l仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent&
这个inp里面是按默认设定进行的，明显不符合实际的受力情况.&p&e*EBu!P
问题的原型是什么？不要太随意的给材料参数和模型。
4J:Y9aO9e/o0I1`&G&z(wSimWe仿真论坛
|仿真|设计|有限元|虚拟仪器inp文件给的材料参数很不合理，你两边材料用的是laminar，参数如下：
o,o:^do仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent*Elastic,
type=LAMINA4u!['}7] @*w5vr
1.,&&0.3,3894.,&
Oh(复合材料单层板纤维方向模量12.6GPa（模量较小的T300都至少要100多GPa），最后两个剪切模量怎么可能是1？至少和12方向剪切模量一个量级，而且应该G12＝G13。单层板的材料方向也不对，纤维方向和ce厚度方向一致，现实中不可能。两侧单元方向也不一致。‑]&Hvo
再说CE的材料，ce刚度还算合理，但还是有些小，这个要看具体情况。损伤起始的强度值给的是6MPa，太低，复合材料层间拉伸强度要高一个量级。损伤扩展剪切能量释放率给的太高，虽然你没用到，顶尖的树脂基体复合材料也少有达到这个数值。另外致命的是粘结单元的方向给错，也就是说你沿1方向的拉伸对于厚度方向为2方向的CE没有任何作用。仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,FluentF7k8~eL‑Y
建议仔细仔细看ABAQUS中的两个例子以及相关内容，如果研究的透彻，一篇硕士论文没问题，再往深里做，一篇博士论文也是可能的。-hy6Vz
ABAQUS两个材料例题中的材料参数堪称完美，直接拿过来放自己模型里就行。
Verifying the stack direction
b(H@&` G$A o]/_The stack direction can be verified visually in
ABAQUS/CAE by using the stack direction query tool (see
“Understanding the role of the Query toolset,” Section 45.1 of the
ABAQUS/CAE User's Manual). For three-dimensional elements
ABAQUS/CAE colors the top face purple and the bottom face brown.
For two-dimensional and axisymmetric elements, arrows indicate
the orientation of the element. In addition, ABAQUS/CAE
highlights any element faces and edges that have inconsistent
orientations.
?9LxtT仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM
Ue/cT仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAMAlternatively,
the material axes can be plotted in the Visualization module of
ABAQUS/CAE to verify that the 3-axis points in the desired normal
direction. If the element is oriented improperly, one of the
in-plane axes (either the 1- or 2-axis) will point in the normal
direction.SimWe仿真论坛kN-fyF0?1yAitB
{仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM以上摘自“ABAQUS
Analysis User's Manual&&18.5.4
Defining the cohesive element's initial
geometry”仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent7O%\w^&j
意思如下（只对二维），CE的厚度方向可以使用CAE中的Stack direction query tool来查询，具体参见CAE
45.1。对于二维的CE，箭头指示的方向就是CE的厚度方向。6f�‑{/BX9m}
或者，可以在CAE
Visualizaiton中通过查看材料方向的方式来验证CE方向，对于二维CE，2方向为其厚度方向。参见18.5.4一节Figure
18.5.4&7 Local directions for two-dimensional and axisymmetric
elements.仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM+P6Oj(E+T
h&Ud\呵呵，还是没认真看帮助。
．连接单元（CONNECTOR）
问：请问高手，如何输出CONNECTOR单元（连接两个部件）的节点所传递的作用力？
设置History Output, 将Domain设为Connector，选中CRF。
已投稿到：
以上网友发言只代表其个人观点，不代表新浪网的观点或立场。}