声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 7376|回复: 5

[转帖]abaqus问答整理

[复制链接]
发表于 2006-4-28 16:46 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
本帖最后由 wdhd 于 2016-4-6 11:00 编辑

  Q:预拉钢筋怎样施加预应力,请各位指点~~~~

  Q:我在文档里看到要在inp文件定义一个rebar,但是rebar只能用于shell, membrane, and solid elements 。我现在想做的是一个预应力拉索,不是镶嵌在shell, membrane, and solid 这些单元里的,而是独立的一根拉锁。拉索单元打算用truss,但是怎样在truss上使用rebar啊?请高手指点

  还有个问题,我看到别人的inp文件,如下:

  *rebar,element=continuum,material=rebar2,name=ubar

  top1,1.005e-4,0.15,0.0,0.5,1

  第二行第一个是setname(top1),第二个是rebar的截面面积(1.005e-4),那第三、第四、第五是指什么?(0.15,0,0.5),最后一个应该是方向,是1方向。 哪位高人指点下第三、四、五项分别代表什么?

  A:施加预应力

  *initinial conditions,type=stress,rebar

  elset,rebar name,所施加预应力的值 ,另prestress hold 为保持所施加的预应力的值不变,我的理解是防止别的构件吃掉所施加的预应力,造成所施加预应力的损失。使用了这个命令之后就避免了这种损失,保证所施加的预应力都施加到了钢筋上。

  A:谢谢指点,你所说的应该是把预应力加在rebar上面,但我发觉truss单元不能定义成rebar,其实是我多想了,truss本来就可以当拉索,实际工程中加预应力只是为了使钢绞线拉紧,起到张拉作用,而在abaqus里,truss本身就是拉紧的,不用施加预应力

  A:我知道模拟加强筋的时候需要用rebar,但 钢筋确实可以直接用truss来模拟 ,而lz所说的预应力其实其实只是施工时的张力而并不是真正意义上的预应力,比如螺栓预应力之类的。如果是索的话可能是要施加预应力的,仅个人看法。

  Q:请教:做一个空间钢框架结构,梁柱用梁元,板采用壳元,打算采用tie命令(共用节点),但不知该如何实现?

  A:我想可以用*equation实现,共用节点的约束情况自己在这一命令下定义。

  A:我因为用命令比较多,但是用cae我想一样, 在CAE里进入命令编辑器,然后编辑就是了,写入*equation命令, 指定约束的自由度(这个看一下标准手册,写得很清楚)

  Q:“Response spectrum analysis(响应谱分析)与Modal dynamic analysis(模态动力分析)区别在什么地方? 如Response spectrum analysis可以进行结构设计? 但Modal dynamic analysis是用来干什么的阿?

  A:就我知道的,modal dynamic analysis应该是振型分解法做动力解析。分解为单自由度体系再取有限个进行组合求反应。

  Q:abaqus如何施加地震荷载?

  A:可以参考abaqus 6.3的例子,Seismic Analysis of a Concrete Gravity Dam

  可以使用:

  1。*amplitude, name=amp, input=seismicdata.dat输入地震波

  2。*boundary, type=acceleration, amplitude=amp 来施加荷载。

  在的2.1.15 Seismic analysis of a concrete gravity dam

  A:这是一个相对的问题,你可以推导一下那个动力方程,

  结果是:ANSYS是取基础固定,解出结构相对基础的相对时程,而ABAQUS是在边界上施加加速度时程,解出的是结构的绝对时程。因而两者相差一个所施加加速度时程。

  在边界上施加加速度时程用显式计算速度较快(主要取决于你的CPU)ABAQUS的STANDERD里面也是显式计算,另外它的EXPLICT模块也可以,具体例子可参见Seismic Analysisof a Concrete Gravity Dam 的分析; ANSYS用的是隐式求解(主要取决于你的硬盘速度和CPU),并且结果暴大(应为它默认输出一大堆数据),可用ANSYS/LSDYNA进行显式计算。

  偶曾经作过共同作用的时程分析,有一点点经验。

  Q:我老板要求我用的是concrete damage plasticity model

  1。在做非线性分析的时候,怎么知道结构破坏了,如何得知什么时候获取了结构的极限承载力

  2。如何确定concrete compression damage和concrete tension damage?二者似乎都无法通过试验得到。如果说可以忽略compression damage,那么由于tension stiffening的存在(下降段相对而言变长了),tension damage的定义就显得很重要了。在ABAQUS关于重力坝的例子中,tension stiffening和tension damage的数据,似乎有一些关系。

  3。定义concrete tension stiffening的时候,我看一些论文上面关于tesion stiffening说明的时候,说如果tensiong stiffening 定义得太小,收敛会很难。但是tensionstiffening的大小在用*concrete tension stiffening定义如何体先呢,在前两个数据是remaining direct stress after cracking和direct cracking strain.

  4。tension stiffenig 与配筋率是否有关系呢?如果有的话,关系是怎么样的呢? 谢谢大家了。

  A:我认为damage plasticity model主要用于混凝土受循环往复荷载的情况,concrete compression damage和concrete tension damage也只在这时需要定义,tension stiffening好像就是定义混凝土单轴受拉本构关系,与tension damage应该没有关系。与配筋率好像也没有关系。

  damage plasticity model既可以模拟循环拟静荷载的情况,也可计算结构在动力荷载(地震、风以及冲击荷载的作用),通过concrete compression damage和concrete tension damage考虑混凝土的拉压异性,tension stiffening则可考虑拉伸屈服后的软还阶段本构行为,同样是混凝土损伤全过程的一部分,它与配筋率及网格划分密度的相互关系很大程度上影响到计算的收敛性。

  A:损伤主要针对循环加载问题,如果你不是循环加载的话可以不考虑损伤 ,我做过比较

  A:楼上的大侠好象问题考虑的不是很全面吧.如果在单向加载时同时存在拉压,就应该考虑.而其中的恢复系数是循环加载下考虑损伤恢复

  A:我现在也在搞concrete damage plasticity model,希望能和各位大虾交流交流

  A:。在做非线性分析的时候,怎么知道结构破坏了,如何得知什么时候获取了结构的极限承载力 ?

  这个问题牵涉到结构在损伤到多少时即认为是损坏.也就是通常说的怎么把结构算倒.不知道我理解的是否正确.如果是这个问题就牵涉到损伤系数的取法.因为是整个结构,就须要把所有的损伤系数加权平均.加权的取法也很有讲究,有人认为达到0.7即认为倒塌,可以参考王光远的文章.现有有限元程序还不能完成,是个急待解决的难题.

  当然是有关系的了concrete compression damage和concrete tension damage分别是受压损伤系数和受拉损伤系数函数,在单项受力时是直接乘入混凝土刚度的.在往复荷载下要乘损伤恢复系数,在三轴应力下还要乘加权因子。

  考虑了该参数后,开裂应变就不等于塑性应变。开裂应变是考虑弹性模量折减的非弹性应变。因此此时塑性应变不等于非弹性应变。这是关键啊

  A:谢谢各位的回复, 帮助很大. 现在riks 去求解的时候,可以获取极限何在,结果与实验结果差别不大,但是在极限荷载下的位移与实验结果差别很大,无论用shell还是solid 来模拟钢筋混凝土的话, 位移误差都很大. 请问大家有没有碰到着个问题呢?

  Q:钢筋混凝土梁中的箍筋如何加入?用CAE操作可以直接加入吗?

  A:cae中可以用truss单元模拟箍筋和纵筋,再采用embed约束实际到混凝土实体中.钢筋单元的位移由他所在的混凝土单元的节点位移来决定.abaqus中的帮助文件里面有详细的说明,abaqus加钢筋的办法有几种:

  1。采用rebar layer 的办法,在part里面画一个面,然后在property里面定义一个surface为rebar layer,把这个surface的属性赋给前面的part里面的那个面。然后在interation中embed中把钢筋层embed到混凝土实体中去。

  2。采用桁架单元的办法,在part里面建好纵筋和箍筋的钢筋骨架,在property中分别赋予截面和属性,在interation中的embed把钢筋骨架embed到混凝土的实体中去。

  3。如果是作构件的话,第二种办法建的比较精确,而且后处理比较方便,查看钢筋单元的应力比较直观,如果是作结构的话,第一种钢筋层的办法比较好,但是个人觉得钢筋层的办法纵筋和箍筋的位置定义的不是很明确。

  楼上高手多谢指点。你说的建钢筋骨架是指在part中用wire画吗?一根一根的画?我做一根混凝土梁的话,我先创建了一个实体代表混凝土梁,然后再create---part画钢筋,画钢筋的坐标系与前面的坐标系应该是相同的?在part模块中创建两个part,对吗?这两个part的坐标系完全相同,要不然在interation中没法embedded.。不知我的想法对不?

  liuxinmei2005兄见笑,我学abaqus也不到一个月,绝对不是什么高手.这是我自己想的办法,我自己的办法是这样的,part里面画一个纵筋,画一个箍筋,在assembly中逐个加入纵筋和箍筋并且定位,形成钢筋骨架,然后再embed.我的办法很笨的,好处就是模型中钢筋的位置非常明确,纵筋箍筋的后处理也都很方便,这种办法对于建一个构件来说还可以,一般10分钟就建完了.对于建结构来说,钢筋层可能好些,很方便,但是后处理想直观的看钢筋的应力就不方便了.

  箍筋的建立的办法画图说明一下:

  1.在part中画一个wire,

  2.赋予wire截面和材料属性

  3.在asemnly中插入,平移或者旋转来布置箍筋的位置

  下图就是箍筋的布置,同样的办法也可以建立纵筋,这样就形成了钢筋的骨架,这么建模比较仔细,后处理比较方便。

  对于剪切破坏的钢筋混凝土构件,箍筋的作用比较关键,必须要建立箍筋单元,对于受弯曲破坏的钢筋混凝土构件,可以不建箍筋。

  对于剪切破坏的混凝土构件,箍筋的作用在任何一本混凝土教材里已经叙述的很清楚啦,由于箍筋的作用对该类型构件影响非常大,当然一定要考虑。

  对于轴向受压构件,由于箍筋对侧向膨胀的约束作用在后期非常明显,所以当然也要考虑。

  而对于受弯破坏的构件,起主要作用的是受压区混凝土以及纵向配置的钢筋,其性能受箍筋的影响很小,从计算成本及建模的方便与否等因素考虑,可以忽略。

  请问ilxy,assembly中的插入操作怎么实现?俺在菜单中只找着了平移,旋转

  所谓的"插入"只是形象的说法,并不是说constrain里有相应的操作,都是通过那几个基本命令实现的.我是用点重合的命令插入钢筋的.在混凝土的相应位置定义一个datum point,再让该点与钢筋相应的关键点重合就好啦.比较苯的方法,供参考.

  钢筋嵌入单元可以直接横穿嵌入钢管混凝土吗?也就是,即嵌入了混凝土又嵌入了钢管?

  请问另外一个问题,不同材料属性的两个instant可以用merge geometry 在一起, 而属性不变吗?我要把型钢梁和栓钉焊接在一起,所以想用这个命令,总说有属性错误,有办法吗?还是要用merge mesh吗?还不会用merge mesh,请高手指点。

  Q:怎么在属性里面定义应力与应变的关系?

  A:可以直接在材料定义的地方输入数据。自己把本构关系的曲线分段,

  然后输入每一段的应力和塑性应变数据对。

  看一下入门手册关于材料非线形那章说明就知道了。

  可以在inp文件中直接输入,*PLASTIC,HARDENING=(强化方法) , 就可以了。
回复
分享到:

使用道具 举报

 楼主| 发表于 2006-4-28 16:46 | 显示全部楼层

回复:(TNC)[转帖]abaqus问答整理

本帖最后由 wdhd 于 2016-4-6 11:02 编辑

  Q:我以前没接触过ABAQUS,有一点ANSYS的基础,但也只是刚入门,没做过具体的项目。

  目前在做竖井施工对软粘土的扰动效应研究的硕士论文,其中一部分工作是用有限元软件模拟竖井施工的挤土效应,考虑土体自重,得到竖井周围土体的塑性区范围,以便与理论解互相验证。应该是一个比较简单的圆孔扩张问题:轴对称建模,在土体中给定一个孔半径大小的位移,不需要考虑接触问题。

  理论解采用的是摩尔库仑屈服条件,而ANSYS里面没有这个模型,所以打算采用ABAQUS。由于只是论文的其中一部分,不可能花大量的时间来深入学习ABAQUS,这两天在版上看了许多贴子,发现在短时间内掌握ABAQUS并不是一件容易的事情,甚是着急,希望高手能给小弟指点一下迷津,如何针对我这个课题做一下分析?3x

  A:把abaqus documentation——abaqus/cae user's manual里面的第二个例子做一下,就了解abaqus cae的基本过程了。例子里面没有而你需要解决的问题是初始地应力的形成。

  过程如下:(假定基本操作你已经会了)

  1、建好模型,设置好材料参数,只施加重力荷载,计算一遍,点view results进入visualization模块,点report菜单,将mises应力的s11,s22,s33,s12,s13,s23输出至文件,计算方式可以选择average,将文件命名为xxx.inp。

  2、将文件中的非数据行删掉,再在每个数据之间加一个逗号(用excel很容易的),每行的最后一个数据不用加,保存。

  3、在cae中重新进入你原先建好的模型,在step模块里面再新建一个step,进入load模块,沿着圆孔径向加位移(和加重力没什么区别),这个步应该发生在新建的step里面,注意原来加的重力荷载不要删了!

  4、进入job模块,在job manage里面点write input,于是inp文件被刷新了,用文本编辑器打开inp文件,找到*step区域和*material区域,在这两个区域之间加入*initial conditions, type=stress, input=xxx.inp

  5、在开始——程序里面将abaqus command界面打开,假设你原先的工作文件名为yyy.inp,打入abaqus job=zzz.inp input=yyy.inp,所以路径都应该一致。

  6、打开log文件,不管有错误还是顺利都会有提示,计算完成后,在command界面打入abaqus viewer database=zzz.odb查看结果

  忘了说了,abaqus可以形成两种inp文件,应该用没有part和assembly的那种,这样单元和结点的编号不会重复,*initial conditions, type=stress, input=xxx.inp才能正确执行,参考"cae_no_parts_input_file"

  Q:就是一根三点弯曲的普通钢筋混凝土梁,我想看它的荷载-位移曲线,如何操作?

  A:利用riks分析,记录位移和反力,再在xy plot operation里面作图。

  *step,INC=100, NLGEOM

  *static,RIKS

  0.001,1.0,,0.04,,4,3,-0.080

  riks指的是弧长法

  也可以利用*NODE PRINT命令,把节点位移和支座反力输出到.dat文件,用excel画图,很方便,修改也容易。

  *node print只能在ABAQUS/standard中起作用,在ABAQUS/explicit中不能使用

  Q:有个问题想请教大侠:

  我在做一个双材料界面的斜裂纹问题,模型见图。裂纹两侧是不同材料。如果不定义局部坐标,输出的裂纹附近的正应力和剪应力应该都是整体坐标系下的,显然这时候的S22,S12并不满足连续条件。

  现在就是想定义局部坐标系,并且输出局部坐标系定义下的S22',S12',不知道abaqus具体怎么样实现,请指教。

  看到过帮助里提到过orient子程序,不知道怎么用

  A:在simwe上坛友的帮助下(当时只记了他的mail:asia_999@sina.com,忘了他的ID,很抱歉),解决了这个问题。

  1。在打开ODB文件的时候,前面的read_only前面的勾去掉。

  2。create coordinate system选择一种方法建立局部坐标。

  3。coordinate system manager选择建立的csys-->move to odb

  4。重新打开odb,main menu bar-->result-->options...-->transformation-->user-spcified-->choose your csys-->apply

  Q:adaptive mesh是不是只能用在explicit里头呀

  A:如果想要在standard里面用,在cae里该如何操作?

  在step模块,在other菜单里设定

  我昨天凌晨三点终于把弹簧加上了,幸亏我是一个人一个宿舍,否则别人会有意见的。在我学习ABAQUS的过程中,从中SIMWE论坛中获益不少,有了办法,不敢独享,特传上来。

  要点:1、弹簧单元只能在INP文件中手工加入。2、*ELEMENT,----和*SPRING,-----这两个语句要挨着。3、新加的节点号和单元号,不要和已有的冲突。4、有不明白的,我愿意解释,最好是晚上,白天我和老师在一个办公室,包教包会

  QQ号:328895637,请注明是ABAQUS 同行

  Q:我想将abaqus cae中的图形导成jpeg,或bmp 文件格式。找了好几天,还是不知道怎么弄?知道的,说一下具体操作。谢谢!

  A:File -> Print -> to file, select tif or png or ..., you will get it.

  Q:section那项里没有point section这一项啊,不知道如何定义,望大虾们指教,谢谢

  A:create section对话框,

  Category一栏选other,然后在type中有point。

  在定义属性的模块下,工具兰有一个定义特殊属性的菜单可以定义惯性属性

  Q:?有哪位大侠知道怎样在ABAQUS中输入初应力文件吗:)

  a:ABAQUS Example Problems Manual

  1.1.10 Stress-free element reactivation里面有部分相关内容,可以查阅之。

  ** --------------------------------------

  *STEP

  step 1: add initial stress state & remove liner

  *STATIC

  *DLOAD

  SOIL, BY, -20.0E3

  *MODEL CHANGE,REMOVE

  LINER,

  *CLOAD,AMPLITUDE=RELAX

  100,1, 5.4086E+04

  101,1, 4.3918E+04

  Q:如何让abaqus以double precision运行?

  A:job manager->precision

  Q:无意间发现abaqus里输出的都是积分点的应力,如果想得到节点处的应力,该怎么办?

  A:这是节点处的应力输出命令:

  *EL PRINT, POSITION=NODES

  *EL FILE, POSITION=NODES

  这是单元质心处的应力输出命令:

  *EL PRINT, POSITION=CENTROIDAL

  *EL FILE, POSITION=CENTROIDAL

  Q:我现在在分析ABAQUS EXAMPLE PROBLEM MANUAL的1.1.7例子,但是当把inp文件导入时出现这个问题:

  WARNING: The following keywords/parameters are not yet supported by the input file reader:

  ----------------------------------------------------------------------------

  -----

  *ELFILE

  *ELPRINT

  *INITIALCONDITIONS, TYPE=HARDENING

  *NODEFILE

  *NODEPRINT

  为什么硬化初始条件导不进去了?如果导不进去又该怎么处理了?请各位高手帮忙

  Q:请问如何实现,用户子程序里如何实现分层填土荷载

  A:

  SUBROUTINE DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,

  1 COORDS,JLTYP,SNAME)

  C

  INCLUDE 'ABA_PARAM.INC'

  C

  DIMENSION TIME(2), COORDS (3)

  CHARACTER*80 SNAME

  H=KSTEP*2

  WIDTH=12

  IF(COORDS(1).LE.WIDTH)THEN

  PLOAD=(WIDTH-COORDS(1))*2/3*20000

  ELSE

  PLOAD=0

  END IF

  IF(COORDS(1).LE.(WIDTH-1.5*H))THEN

  NLOAD=-(WIDTH-1.5*H-COORDS(1))*2/3*20000

  ELSE

  NLOAD=0

  END IF

  F=PLOAD+NLOAD

  RETURN

  END

  Q:在ABAQUS/EXPLICIT计算中,请问:

  *OUPUT ,FIELD,NUMBER INTERVAL=n,TIME MARKS=YES (OR NO)

  *RESTART,WRITE,NUMBER INTERVAL=n,TIME MARKS=YES (OR NO)

  中,这里的n是什麽?取多少合适?MARKS=YES或NO有何区别?谢谢!!

  上述的n在ABAQUS/CAE中怎样设置?

  A:n是你取的输出结果的间隔数,取n为4是指每计算4个迭代步骤输出一次结果。time mark 取yes就行了。

  Q:在abaqus中做钢筋混凝土的分析,钢筋怎么样加进去,能不能象ansys中直接画进去?

  还是要在inp文件中自己输进去!请大侠指点!

  A:使用rebar element

  Q:abaqus 的cae中有没有弹簧元件, 怎么调用?

  A:You should use editor to edit *.inp. Add spring element and the property of spring.

  Q:请问各位高手,如何在viewer里将contour最大值与最小值所在处以max及min文字符号表示?

  A:U will see the max. & min. node # in legend.

  Q:我在做一个例子时发现,在一个input文件中不能定义两种不同的材料属性。后来我在keywords manual中关于*material的用法中看到了以下的文字:

  Required parameter:

  NAME

  Set this parameter equal to the name to be used for the material when it is referenced in the element property options. Material names in the same inputfile must be unique. Material names adhere to the naming convention for labels (see ``Input syntax rules,'' Section 1.2.1 of the ABAQUS/Standard User's Manual and the ABAQUS/Explicit User's Manual), except that they cannot begin with a number.

  这是不是就是说ABAQUS在一个算例中没法定义一种以上的材料?这个局限性好像有点大了吧?很多问题都需要定义多种材料属性的!我印象里好像其CAE中还可以分别给不同的部分施以不同的材料属性,难道反而在程序流中不能适用吗?还是我的理解有偏差? 万望高手能指点迷津!

  A:Material name is unique. You could definie many different materials by different name.

  Q:abaqus里做静力分析时怎么加重力加速度?在inp文件里的card是什么?

  如:

  *DLOAD

  ELSET, GRAV, 9.81, 0., 0., 1.

  A:有疑问查一些*DLOAD卡就清楚了

  *DLOAD和*CLOAD有什么不同?

  *CLOAD : concentrated load, 集中負載

  *DLOAD : distributed load, 分佈負載

  Q:请问如何设置ABAQUS的默认路径?

  A:我是这样做的:

  比如,工作目录在d:\abaqusdata

  从开始运行abaqus command。

  c:>d:

  d:>cd abaqusdata

  d:\abaqusdata>abaqus job=**** interactive

  or

  d:\abaqsudata>abaqus viewer/cae

  所有工作文件都保存在d:\abaqusdata中。

  对于command方式,一般就是在文件的目录下运行abaqus就可以!生成的文件全部在运行的文件夹下面。

  而对于cae方式也可以通过另存为...而指定路径的
发表于 2008-12-7 09:53 | 显示全部楼层
真是好东西啊,3X:handshake
发表于 2010-8-29 14:41 | 显示全部楼层
怎么看不完整啊?
发表于 2010-8-29 14:42 | 显示全部楼层
如何查看完整版本,文章里面好多符号,顺序很乱,我登录了啊。
发表于 2010-8-30 17:37 | 显示全部楼层
让人看了眼更晕,哈哈,不过整理到是件好事。
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-25 03:53 , Processed in 0.063062 second(s), 19 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表