ABAQUS中Truss单元预应力的两种施加方法
- 格式:docx
- 大小:20.54 KB
- 文档页数:6
ABAQUS模拟预应力筋的方法1.降温法这是目前很多人采用的方法。
即在预应力筋施加温度荷载(降温),使预应力筋收缩,从而使混凝土获得预应力。
2.ABAQUS自带的初始应力法直接用*Initial conditions, type=stress可以直接模拟先张法,能获得预应力筋和混凝土的后期应力增量,但无法获得预应力筋的真实应力。
3.Rebar element single 法利用ABAQUS提供的rebar功能,模拟预应力束,给出rebar与相关实体单元的信息,通过在rebar上施加初始应力即可模拟先张法和后张法。
4. MPC法分别定义预应力筋(比如truss单元)和混凝土,采用MPC将预应力筋与混凝土联系起来,对预应力筋施加初始应力,即可模拟预应力效应。
5.Rebar Layer法利用ABAQUS提供的rebar layer功能,将rebar layer定义到surface,membrane或shell基上,通过对rebar施加初始应力,即可模拟先张法和后张法。
经过一段时间的使用和尝试,发现实体内施加预应力还存在不少缺陷:1.无法模拟早期的预应力损失,如摩擦损失,锚具回弹损失等;2.无法准确模拟后张法中在张拉阶段净截面参与计算的问题,这在截面高度较小,预应力筋较多时,对计算结果影响会比较大;3.无法模拟换算截面的问题,尽管帮助文件中多次提到rebarlayer的刚度被添加到surface section等中,由于surface section没有内在刚度,多次测试发现rebar layer的刚度无法添加到结构中。
后尝试用shell section的方式来实现。
帮助文件中没有直接提到用shell section带rebar layer埋于solid 单元的方式可以模拟预应力。
经多次测试发现是可以考虑shell 和rebar layer的附加刚度,但结算结果不稳定。
几个要点:1>.shell section能自动采用换算截面,其但换算系数为N而不是N-1。
abaqus truss结构Abaqus Truss 结构引言:Abaqus 是一种广泛使用的有限元分析软件,它可以用于模拟和分析各种工程结构的行为。
在这篇文章中,我们将重点介绍Abaqus 中的 Truss 结构。
Truss 结构概述:Truss 结构是一种由直线杆件组成的简单结构,它们通过节点连接在一起。
这些杆件可以是理想化的,也可以是实际的。
Truss 结构常见于桥梁、建筑和机械结构中,其简单性和可靠性使其成为工程实践中常用的结构形式。
Abaqus 中的 Truss 结构建模:在 Abaqus 中,建立 Truss 结构模型需要进行以下步骤:1. 几何建模:首先,根据实际结构进行几何建模。
可以通过直线段来表示每个杆件,并通过节点将其连接在一起。
在建模过程中,需要考虑节点的位置、杆件的长度和方向等几何参数。
2. 材料属性:为了对Truss 结构进行分析,需要为杆件定义材料属性。
这包括杆件的弹性模量、截面积等材料参数。
这些参数将在后续分析中用于计算杆件的应力和应变。
3. 荷载和边界条件:在Truss 结构中,荷载可以通过节点施加。
可以在节点上定义不同的荷载类型,如集中力、分布力或位移约束等。
此外,还需考虑支座条件和约束条件,以确保结构的稳定性和平衡性。
4. 求解模拟:完成Truss 结构的建模和加载后,可以通过求解器执行模拟分析。
Abaqus 将根据所定义的材料属性、荷载和边界条件对结构进行求解,并给出相应的应力、应变和位移结果。
5. 结果分析:完成求解后,可以对Truss 结构的结果进行分析。
可以通过查看应力云图、应变云图和位移云图来了解结构的行为。
此外,还可以提取节点或杆件的特定结果数据,如最大应力、最大位移等。
Truss 结构的优势:Truss 结构具有以下几个优点:1. 高强度:由于Truss 结构杆件的特殊形状,它们可以在承载大荷载的情况下保持高强度和刚度。
2. 轻量化:相比其他结构形式,Truss 结构通常具有较低的重量。
[警告错误信息]【错误和警告信息汇总】(此贴为复件,请勿回复)[复制链接]zsq-w管理员CIO仿真币33975阅读权限255 电梯直达1#发表于 2009-5-7 17:08:16 |只看该作者|倒序浏览本帖最后由 zsq-w 于 2009-6-2 17:11 编辑*************************错误与警告信息汇总*************************--------------简称《错误汇总》***ERROR***WARNING***二次开发%%%%%%%%%%%%%%% @@@ 布局@@@ &&&&&&&&&&&&&&&&&&&&&&常见错误信息常见警告信息网格扭曲cdst udio斑竹总结的fortran二次开发的错误表%%%%%%%%%%%%%%%%% @@@@@@ &&&&&&&&&&&&&&&&&&&&&&&&&模型不能算或不收敛,都需要去monitor,msg文件查看原因,如何分析这些信息呢?这个需要具体问题具体分析,但是也存在一些共性。
这里只是尝试做一个一般性的大概的总结。
如果你看见此贴就认为你的warning以为迎刃而解了,那恐怕令你失望了。
不收敛的问题千奇万状,往往需要头疼医脚。
接触、单元类型、边界条件、网格质量以及它们的组合能产生许多千奇百怪的警告信息。
企图凭一个警告信息就知道问题所在,那就只有神仙有这个本事了。
一个warning出现十次能有一回参考这个汇总而得到解决了,我们就颇为欣慰了。
v1.0可编写可改正ABAQUS中预应力 Truss 单元的两种实现方法例题:100 m 长钢缆水平搁置从x = 0 到 x = 100 。
两头固定。
无初始拉力,计算下垂量。
截面; A = m 2,Density: r =7800 kg/m3, g = m/s 2,E=+11 N/m2Analytical solution of maximum displacement (u 2 ) at x = 50m :U2_max = -((3*r*g*L^4)/(64*E))^(1/3) = m方法一 .沿truss element加沿长度方向初始拉应力(see此文件中使用了initial condition, type = stress方法加初始拉应力。
因工程上无此初应力 , 更好的方法是使用降温法。
算完后再升温。
用降温法。
算完后再升温。
NOTE: 降温法施加预应力(激活钢绞线)。
温度=- 力 / (膨胀系数 * 弹模 * 钢绞线面积)1、第一步,在 truss单元中施加一个初始应力,让计算处于初始均衡状态;初始应力设置过小,可能不收敛,应多次试算,找到一个合理的应力值。
一般状况下,这个初始值对最后值的影响不大,能够忽视。
2、第二步,施加truss单元的自重荷载,翻开非线性开关(nlgeom=YES )考虑几何非线性问题;3、本例中初始值采纳。
自重作用下缆索的拉应力约为80Mpa。
最大位移为m,与理论计算值符合得很好。
*HeadingCable appling gravity load with initail stressThe maximum Analytical displacement without initail stress (at x = 50 m) U2 = meter****Method 1. Using * initial condition,type = stress method***Preprint, echo=NO, model=NO, history=NO,contact=NO *Node1,0.,0.101, 100., 0.*NGEN, NSET = NALL1, 101, 1*Element, type=T2D21, 1, 2*ELGEN, elset = ELALL1, 100, 1*ELSET,ELSET=EL_OUT1, 51, 100*Solid Section, elset=ELALL, material=steel,***Nset, nset=Left1,*Nset, nset=right101,*Nset, nset=mid51,****MATERIALS***Material, name=steel*Density7800.,*Elastic+11,*initial condition, type = stress**Note: the solution will not converge as the initial stress < 100,000N/m^2ELALL, 100000*Boundary Left,1, 2 Right, 1,2*Step, name=Step-0, inc=1000Initial stress equilibrium*Static1, 1., 1e-05, 1.*Output, field, variable=PRESELECT*Output, history, variable=PRESELECT*Node print, nset = mid, freq = 1000U,*EL PRINT, ELSET=EL_OUT, freq = 1000S*END STEP**---------------------------------------------------------------- ****STEP: Step-1***Step, name=Step-1, nlgeom=YES, inc=1000Apply gravity load*Static, 1., 1e-05, 1.** Name: GRAVITY-1Type: Gravity*DloadELALL, GRAV, , 0., -1.****OUTPUT REQUESTS***Restart, write, number interval=1, timemarks=NO *Output, field, variable=PRESELECT*Output, history, variable=PRESELECT*Node print, nset = mid, freq = 1000 U,*EL PRINT, ELSET=EL_OUT, freq = 1000S*End Step方法二 . 使用 STABILIZE parameter on the *STATIC.(see“stabilization ”在构造上附带 artificial viscous damping(粘滞阻尼),使得计算结果 to go beyond the instability point 。
问题:做了一个小模型,用降温模拟预应力,但是桁架在埋入砼后,我们开始计算分析时候报错,说没有给桁架赋予属性。
郁闷了。
不会了。
请教。
附上CAE
*************
今天回过头来看这个帖子,真的觉得好汗颜自己。
1 “说没有给桁架赋予属性”
lz在property步给part2赋予了truss的材料属性,然后在mesh步给划成beam单元,所以出此错误信息。
要么都用beam,要么都用truss。
改为truss单元即可收敛。
2 lz是采用降温法模拟钢筋预应力做得一个小例子。
经常有人求这个。
现贴出CAE文件和inp文件供大家讨论。
另外,我设置温度的时候是有一个公式可以换算的。
但是换算时候的压力值取什么值我还不是很明白,所以要请业内做桥梁的高人指点下。
其次,我发现在检测的时候有这个提示,说是没有赋予初始温度,我在预定义场里只定义了一个温度。
这个温度是末温度。
温度是从0度到这个预定义的值的。
也就是一个降温过程。
附图给大家看看。
这个不是错误。
若不定义初始温度,则缺省值取为0。
不过会有个warning提示你的。
当然你也可以在step之前设定初始温度。
ABAQUS建模如何施加预应力1.第一种方法是在模型建立阶段施加初始的预应力。
假设我们需要在一个弹性杆上施加预应力,可以选择一个适合的初始几何形状,然后在ABAQUS中建立一个非线性静力分析模型。
在模型中创建一个新的材料属性,并将该属性的引伸模量设置为预应力下的值。
然后,在加载步骤中施加适当的载荷以实现所需的预应力。
这种方法需要提前了解预应力的大小和方向。
2.第二种方法是使用一维元素(例如梁或弹簧元素)来模拟预应力的效果。
这种方法在模拟螺旋弹簧或拉索等情况下特别有用。
首先,在ABAQUS中创建一个线性静力分析模型,并将相应的材料属性分配给模型。
然后,在模型中创建一维梁或弹簧元素,并将其与表面节点相连。
在加载步骤中设置合适的载荷以利用这些元素施加预应力。
这种方法可以更加自由地控制预应力的大小和方向。
3.第三种方法是使用INITIALCONDITIONS卡来施加初始的预应力。
首先,在ABAQUS的输入文件中的合适位置添加一个INITIALCONDITIONS卡。
然后,在该卡中通过修改合适的变量和属性定义施加所需的预应力。
INITIALCONDITIONS卡的使用可以灵活地定义各种初始条件,但需要对ABAQUS的输入文件格式有一定的了解。
无论使用哪种方法,施加预应力时需要考虑一些因素。
首先,需要确定预应力的大小和方向,以便正确定义材料的属性和加载条件。
其次,需要注意预应力的影响范围,以便选择合适的单元类型和网格划分方式。
最后,需要进行适当的验证和调试,以确保模型的准确性和稳定性。
在使用ABAQUS进行建模时,建议先进行一些小规模的验证和参数敏感性分析,以确保所施加的预应力是可靠和合理的。
此外,在模拟中施加预应力时,需要根据具体情况和问题的要求选择合适的方法。
abaqus 粱单元预紧力
答:Abaqus是一款广泛使用的工程仿真软件,可以模拟各种材料和结构的力学行为。
在Abaqus中,预紧力通常是通过在载荷步中施加力或位移来实现的。
对于梁单元,可以通过以下步骤施加预紧力:
1. 打开Abaqus软件,创建或导入模型。
2. 在模型中创建梁单元,并指定其材料属性、截面尺寸等。
3. 在载荷步中施加预紧力。
可以通过以下几种方式实现:
a. 在载荷步选项卡中选择“力”或“位移”选项,并在右侧的输入框中输入预紧力值。
b. 如果需要分段施加预紧力,可以在载荷步选项卡中选择“多个”选项,并添加多个载荷子步。
在每个子步中输入相应的预紧力值。
c. 如果需要施加的预紧力与时间有关(例如蠕变预紧力),可以在载荷步选项卡中选择“时间”选项,并指定预紧力随时间变化的函数。
4. 确认载荷步设置后,执行模拟计算。
需要注意的是,施加的预紧力应该根据实际情况进行合理设置。
如果施加的预紧力过大,可能会导致梁单元在模拟过程中发生断裂或其他非预期结果。
因此,建议在进行模拟前进行必要的验证和校核。
ABAQUS建模如何施加预应力残余应力在ABAQUS中,可以通过几种方法施加预应力残余应力。
下面将详细介绍两种常用的方法:二层法和热加载法。
1. 二层法(Two Layer Method):二层法是一种模拟材料加工过程中产生的预应力方法。
基本思路是在模拟之前的一部分载荷历史之后,在初始状态下施加一些预应力。
其步骤如下:(1)准备一个加载步(Apply Load Step),在此步骤中定义预应力载荷的初始状态。
(2)定义载荷历史。
(3)在载荷历史的一部分之后,将模型还原为初始状态,并在此状态下施加预应力。
(4)在预应力载荷下继续加载模型。
(5)根据需要将数据保存。
例如,在ABAQUS/Standard中,可以在步骤中使用命令`STATIC`定义预应力载荷以施加预应力。
以下是一个使用二层法施加预应力的示例代码:```python*Step, name=initial step, nlgeom=yes*Static*End Step*Step, name=loading step, nlgeom=yes*Static*End Step```上述代码中,第一个步骤定义了预应力载荷的初始状态,并保持模型为非线性几何模型。
第二个步骤中的载荷历史定义了加载模型时施加的载荷。
预应力在两个步骤之间施加。
2. 热加载法(Thermal Loading Method):热加载法是一种在模拟焊接过程等应用中施加预应力的方法。
基本思路是通过施加热载荷引起温度梯度,从而产生预应力。
其步骤如下:(1)定义一个温度场,可以使用定义节点温度或通过导入温度场施加。
(2)应用热加载,在模型中引入相应的热载荷。
(3)在热载荷下,施加机械载荷以保持平衡。
(4)根据需要将数据保存。
例如,在ABAQUS/Explicit中,可以使用`*Temperature`和`*Amplitude`命令定义温度场。
下面是一个使用热加载法施加预应力的示例代码:```python*Amplitude, name=temperature, definition=SMOOTH STEP, smooth=ON*Initial Conditions, type=TEMPERATUREAll NSET*Temperatureall, type=AMP, amplitude=temperature, fixed=OFF*Step, name=loading step, nlgeom=yes*Static*End Step```上述代码中,首先定义了一个温度场,并将其应用于所有节点。
2.7索单元施加预应力的几种方法侯晓武索单元为只受拉的三维线性单元,只能传递轴向拉力,不能受压也不能受弯。
可按如下方式定义索材料。
图2.7.1 自定义索材料如图2.7.1所示,设计类型选择“用户定义”,而后分别输入各种参数。
泊松比、线膨胀系数、容重与钢材相同,弹性模量需单独定义。
线性分析时,索单元将被等效为桁架单元,其刚度由弹性刚度和下垂刚度两部分组成。
sagE K K K 111+=弹性刚度:LEA E K =下垂刚度:32312L w T K sag=因而,进行线性分析时,用索单元建模和用桁架单元建模其结果会有差异。
由于采用索单元建模时,考虑了下垂刚度后,会对结构的刚度进行折减。
因而采用索单元建模时,其位移结果要大于采用桁架单元建模的结果(如图2.7.2所示)。
(a)索单元(b)桁架单元图2.7.2 张弦梁恒载作用下线性分析结果非线性分析时,悬索结构的刚度并不是一次性计算得到。
而是通过多次重复计算,使结构达到平衡状态得到的。
在midas Gen中,有四种施加预应力的方法。
(1)建立单元同时定义预应力A.无应力索长LuLu<L时:引入初拉力F=EA∙∆L/L(∆L=L-Lu)L为单元长度,EA分别为弹性模量和横截面积;Lu>L时:索下垂,折减刚度;B.初拉力:输入沿单元坐标系X轴方向的预应力荷载;C.水平力:输入沿水平方向的预应力荷载。
说明:▪仅用于几何非线性分析,线性分析时将被忽略;如图2.7.3所示,在下方张弦梁中的324~332号索单元施加初拉力770KN,上方张弦梁中的索单元不施加初拉力。
对两个模型进行线性分析,得到恒载作用下的位移结果相同,如图2.7.4所示。
图2.7.3 建立单元时施加初拉力图2.7.4 恒载作用下位移结果▪对所有荷载工况结果都有影响;▪迭代计算时,第一步即产生初始刚度,该拉力对其它构件也有影响,如图2.7.6和图2.7.7所示。
菜单:荷载->静力荷载:如图2.7.5所示,定义一个空工况,在该工况下不定义任何荷载。
小弟最近做了一个加固的模拟,对于预应力钢筋的建模稍有体会,因此分享一贴分为两个大类进行阐述:1.CAE流首先建立一个分析步,对truss单元施加bolt那啥的力(忘记具体名字了,在寝室的电脑,没有装abaqus,抱歉),然后进行后续的分析布分析。
需要注意的是truss单元必须在中间分割,以便有一个中间点来施加这个力。
由于我不是采用这种方法,仅仅做了一个小东西验证下面的方法的正确性,因此也不好意思多说,具体参见分析手册29-5-1" p2 I2 L- F/ W- V; g 2.INP流(直接修改INP文件); O5 _' }# v/ y: ]: _7 |2 ~& ?首先强调:模型中应该有需要施加预应力的单元信息,无论是B31或者是T3D2,后面自行处理9 u# X( U! u- w/ k在*element章节与*node章节之间,加入*node1000XX,0,0,0(XX是编号,因为我一共有16个点需要施加预应力,因此XX就是1~16,而前面的1000是为了避免与已有node号码重复)+ S7 Q9 b/ u9 R8 c*PRE-TENSION SECTION, ELEMENT=XX, NODE=1000XX1 Z: Z2 p: x0 O8 {& C: c9 {; U' b4 @& e 1 z 然后在 \% X" g" p6 C- I" D# _*element章节把需要施加预应力的杆系单元改为T3D2单元,这个应该都会改吧$ f1 [$ S0 a: x/ G& i' w; R另外,如果模型中没有truss单元的截面信息,需要自行加入,如下6 j$ Z1 R5 A R7 y% z, _*SOLID SECTION, ELSET=EL-pretension, MATERIAL=rebar6 `9 O* L6 a/ z5 }5 o! I1 D0.002965然后加入预应力的分析工况,用*CLOAD写入即可。
ABAQUS施加预应力参考百度文库:/link?url=dt_VLOGCUf8hUo7A9THhyv7BuSHry71EbLVtBtkWpoiYtkm Lxbfk0Io63jsygs6vWbFU7x22HHFv8pIGgPMYkv1lyFXWbgPJqvCodSioUqa关键字格式:“*initial conditions, type=stress, input=bb.dat”上面的关键字,即绿色部分,全部插于*STEP语句之前(如下图),两语句之间不能有空格。
施加预应力场只是initial conditions关键字的一个应用,详见abaqus6.8帮助文档,《ABAQUS Analysis User’s Manual》的第28.2节“initial conditions”。
实例:点焊所产生的焊点中存在着残余应力,本文就是教大家如何完成焊点残余应力的模拟。
原理说明:先在模型上施加一个任意载荷(记为状态1),可得出此载荷作用下模型上的等效节点载荷,然后通过keywords让这个等效节点载荷作用于模型上,使它与之前施加在模型上的载荷相平衡,便得到了一个位移为0的初始状态(记为状态2),该状态下,模型中含有的应力场与状态1相同,只是模型不发生变形而已。
具体操作步骤:1、建立有限元模型,部件类型为轴对称2、设置材料常数(自己任意设)3、分析步,设置两个分析步4、设置任意一个自己需要的载荷,此载荷即为与初始应力对应的载荷。
让该载荷从分析步2开始作用。
分析步1空着。
原因不明。
(让载荷从step2开始,得到的分析结果图显示很光滑。
但若是让载荷从step1开始并延续到step2,或者从没有step2 的step1开始,得到的分析图都不是光滑,原因未知。
)设置边界条件,从状态Initial开始。
5、mesh:网格算法为网格类型为CAX4R6、Job模块下,创建工作名称为“Job-1-1”的名字,提交分析。
7、将分析得到的应力场保存为一个文本文件,ipn格式。
1.按ABAQUS入门指南操作。
用define cutting plane把部件一分为二,显示dependent part instances cannot be本人菜鸟1级。
按ABAQUS入门指南操作。
用define cutting plane把部件一分为二,显示dependent part instances cannot be edited or assigned mesh attributes.用别的par tition cell也是一样,请教大侠帮忙。
这是因为你在进入assembly后创建的是dependent part instances,要建成independent part ;还有另外一种方法就是,如图所示在模型树上找到这个part右键点然后有变成independent选项,你试试吧,应该可以的2.abaqus中出现,Dependent part instances connot be edited ,怎么办?如何设置independent选项呢?谢谢左边模型树里面,Assembly-instances里面的子选项右击,点Make independent 就行了选择make independent,这是在creat instance里面出现的问题,对以后的mesh 有影响你就用part的方式编辑就可以了3. ABAQUS中混凝土梁中的箍筋如何加入abaqus加钢筋的办法有几种:1。
采用rebar layer 的办法,在part里面画一个面,然后在property里面定义一个surface 为rebar layer,把这个surface的属性赋给前面的part里面的那个面。
然后在interation中embed中把钢筋层embed到混凝土实体中去。
2。
采用桁架单元的办法,在part里面建好纵筋和箍筋的钢筋骨架,在property中分别赋予截面和属性,在interation中的embed把钢筋骨架embed到混凝土的实体中去。
太 阳 能第08期 总第352期2023年08月No.08 Total No.352Aug., 2023SOLAR ENERGY0 引言随着中国明确提出2030年碳达峰与2060年碳中和的宏伟目标,光伏市场再次得到了广泛而深入的发展。
根据国家能源局发布的信息显示,2022年中国光伏发电新增装机规模再创新高,已达到87.41 GW 。
近些年,集中式地面光伏电站大规模建设,使其应用模式必须谋求多元化空间发展。
柔性光伏支架系利用高强柔索[1]做为光伏阵列的主载体,并结合光伏阵列荷载相对偏小的特征,实现光伏阵列发电单元的大跨度布置,继而实现光伏与农业、畜牧业等的有机结合,并充分利用土地资源。
目前,柔性光伏支架多以单独柔索受力,虽额外配以辅助限位索,但对光伏阵列的稳定性控制仍不尽人意——“小风有微震,大风有晃动”,这会对光伏阵列的发电效率造成一定影响;还有一种索桁架柔性光伏支架,其应用较好,但成本偏高。
本文对“柔+刚”新概念柔性光伏支架进行研究,即基于ABAQUS 有限元受力分析充分的情况,以柔索为主要载体,在柔索上固定刚性梁的“刚柔并济”的柔性光伏支架,吸收“柔性”和“刚性”两种支架的优势特征,并对两种光伏支架的劣势特征进行弱化或解决,从而使光伏与农业、畜牧业等更好的结合,得到更好的发展。
对3种不同状态的简支梁的内力大小和变形模态分别进行数值计算模拟分析,再对分析结果进行对比,并得出结论。
1 条件设定3种状态不同的简支梁(均布荷载q 均为1 kN/m ,跨度均为6 m)的截面示意图如图1所示。
简支梁的截面为H 型钢,截面高度为200 mm ,翼缘宽度为75 mm ,腹板及翼缘厚度均为2 mm 。
需要说明的是,此截面尺寸不常用,本文只用来做理论计算,研究计算方法。
根据GB 50017—2017《钢结构设计标准》[2]和GB 50010—2010《混凝土结构设计规范》[3],选取简支梁的材料物理参数,如表1所示。
abaqus问答自己总结1、Q:预拉钢筋怎样施加预应力,请各位指点~~~~Q:我在文档里看到要在inp文件定义一个rebar,但是rebar只能用于shell, membrane, and solid elements 。
我现在想做的是一个预应力拉索,不是镶嵌在shell, membrane, and solid 这些单元里的,而是独立的一根拉锁。
拉索单元打算用truss,但是怎样在truss上使用rebar啊?请高手指点还有个问题,我看到别人的inp文件,如下:*rebar,element=continuum,material=rebar2,name=ubartop1,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,rebarelset,rebar name,所施加预应力的值,另prestress hold 为保持所施加的预应力的值不变,我的理解是防止别的构件吃掉所施加的预应力,造成所施加预应力的损失。
使用了这个命令之后就避免了这种损失,保证所施加的预应力都施加到了钢筋上。
A:谢谢指点,你所说的应该是把预应力加在rebar上面,但我发觉truss单元不能定义成rebar,其实是我多想了,truss本来就可以当拉索,实际工程中加预应力只是为了使钢绞线拉紧,起到张拉作用,而在abaqus里,truss本身就是拉紧的,不用施加预应力A:我知道模拟加强筋的时候需要用rebar,但钢筋确实可以直接用truss来模拟,而lz所说的预应力其实其实只是施工时的张力而并不是真正意义上的预应力,比如螺栓预应力之类的。
如果是索的话可能是要施加预应力的,仅个人看法。
ABAQUS中预应力Truss单元的两种实现方法例题:100 m 长钢缆水平放置从 x = 0到 x = 100。
两端固定。
无初始拉力,计算下垂量。
截面; A = 0.01539 m2,Density: r =7800 kg/m3, g = 9.8 m/s2,E=2.1e+11 N/m2 Analytical solution of maximum displacement (u2) at x = 50m :U2_max = -((3*r*g*L^4)/(64*E))^(1/3) = -1.194944005 m方法一. 沿 truss element 加沿长度方向初始拉应力 (see job-1.inp)此文件中使用了initial condition, type = stress方法加初始拉应力。
因工程上无此初应力, 更好的方法是使用降温法。
算完后再升温。
用降温法。
算完后再升温。
NOTE: 降温法施加预应力(激活钢绞线)。
温度=-力/(膨胀系数*弹模*钢绞线面积)1、第一步,在truss单元中施加一个初始应力,让计算处于初始平衡状态;初始应力设置过小,可能不收敛,应多次试算,找到一个合理的应力值。
一般情况下,这个初始值对最终值的影响不大,可以忽略。
2、第二步,施加truss单元的自重荷载,打开非线性开关(nlgeom=YES )考虑几何非线性问题;3、本例中初始值采用0.1Mpa。
自重作用下缆索的拉应力约为80Mpa。
最大位移为 -1.195 m,与理论计算值吻合得很好。
*HeadingCable appling gravity load with initail stressThe maximum Analytical displacement without initail stress (at x = 50 m)U2 = -1.194944005 meter**** Method 1. Using * initial condition,type = stress method***Preprint, echo=NO, model=NO, history=NO, contact=NO*Node1, 0., 0.101, 100., 0.*NGEN, NSET = NALL1, 101, 1*Element, type=T2D21, 1, 2*ELGEN, elset = ELALL1, 100, 1*ELSET,ELSET=EL_OUT1, 51, 100*Solid Section, elset=ELALL, material=steel0.01539,***Nset, nset=Left1,*Nset, nset=right101,*Nset, nset=mid51,**** MATERIALS***Material, name=steel*Density7800.,*Elastic2.1e+11, 0.3*initial condition, type = stress** Note: the solution will not converge as the initial stress < 100,000 N/m^2 ELALL, 100000*BoundaryLeft, 1, 2Right, 1,2*Step, name=Step-0, inc=1000Initial stress equilibrium*Static1, 1., 1e-05, 1.*Output, field, variable=PRESELECT*Output, history, variable=PRESELECT*Node print, nset = mid, freq = 1000U,*EL PRINT, ELSET=EL_OUT, freq = 1000S*END STEP** ----------------------------------------------------------------**** STEP: Step-1***Step, name=Step-1, nlgeom=YES, inc=1000Apply gravity load*Static0.01, 1., 1e-05, 1.** Name: GRAVITY-1 Type: Gravity*DloadELALL, GRAV, 9.8, 0., -1.**** OUTPUT REQUESTS***Restart, write, number interval=1, time marks=NO*Output, field, variable=PRESELECT*Output, history, variable=PRESELECT*Node print, nset = mid, freq = 1000U,*EL PRINT, ELSET=EL_OUT, freq = 1000S*End Step方法二. 使用 STABILIZE parameter on the *STATIC.(see job-2.inp)“stabilization”在结构上附加artificial viscous damping(粘滞阻尼),使得计算结果to go beyond the instability point。
ABAQUS建模如何施加预应力本文参考了百度文库中的文章:/link?url=dt_VLOGCUf8hUo7A9THhyv7BuSHry71EbLVtBtkWpoiYtkm Lxbfk0Io63jsygs6vWbFU7x22HHFv8pIGgPMYkv1lyFXWbgPJqvCodSioUqa关键字格式:“*initial conditions, type=stress, input=bb.dat”上面的关键字,即绿色部分,全部插于*STEP语句之前(如下图),两语句之间不能有空格。
施加预应力场只是initial conditions关键字的一个应用,详见abaqus6.8帮助文档,《ABAQUS Analysis User’s Manual》的第28.2节“initial conditions”。
实例:点焊所产生的焊点中存在着残余应力,本文就是教大家如何完成焊点残余应力的模拟。
原理说明:先在模型上施加一个任意载荷(记为状态1),可得出此载荷作用下模型上的等效节点载荷,然后通过keywords让这个等效节点载荷作用于模型上,使它与之前施加在模型上的载荷相平衡,便得到了一个位移为0的初始状态(记为状态2),该状态下,模型中含有的应力场与状态1相同,只是模型不发生变形而已。
具体操作步骤:1、建立有限元模型,部件类型为轴对称2、设置材料常数(自己任意设)3、分析步,设置两个分析步4、设置任意一个自己需要的载荷,此载荷即为与初始应力对应的载荷。
让该载荷从分析步2开始作用。
分析步1空着。
原因不明。
(让载荷从step2开始,得到的分析结果图显示很光滑。
但若是让载荷从step1开始并延续到step2,或者从没有step2 的step1开始,得到的分析图都不是光滑,原因未知。
)设置边界条件,从状态Initial开始。
5、mesh:网格算法为网格类型为CAX4R6、Job模块下,创建工作名称为“Job-1-1”的名字,提交分析。
abaqus在混凝土中加钢筋的两个办法1。
采用rebar layer 的办法,在part里面画一个面,然后在property里面定义一个surface为rebar layer,把这个surface的属性赋给前面的part里面的那个面。
然后在interation中embed中把钢筋层embed到混凝土实体中去。
2。
采用桁架单元的办法,在part里面建好纵筋和箍筋的钢筋骨架,在property中分别赋予截面和属性,在interation中的embed把钢筋骨架embed到混凝土的实体中去。
3。
如果是作构件的话,第二种办法建的比较精确,而且后处理比较方便,查看钢筋单元的应力比较直观,如果是作结构的话,第一种钢筋层的办法比较好,但是个人觉得钢筋层的办法纵筋和箍筋的位置定义的不是很明确。
用truss模拟钢筋,要在mesh的时候指定他是truss单元,过程如下: __________________________________| part 模块:用wire的方法画线||__________________________________|||\/__________________________________| property 模块:创建钢筋的section || property(在category里面选beam-> || truss ||__________________________________|||\/__________________________________| assembly 模块:建立instance, ||__________________________________|||\/__________________________________| Interaction 模块: 在constaint里面|| 选embedded ||__________________________________|||\/________________________________| mesh 模块:指定单元属性|| 钢筋单元必须为truss单元(T3D2等) ||__________________________________|附件是一个用REBAR在混凝土中加入钢筋的例子钢筋混凝土建模:混凝土模型+钢筋模型+混凝土和钢筋的相互作用●ABAQUS中分别定义混凝土本构和钢筋的本构关系。
ABAQUS中预应力Truss单元的两种实现方法例题:100 m 长钢缆水平放置从 x = 0到 x = 100。
两端固定。
无初始拉力,计算下垂量。
截面; A = 0.01539 m2,Density: r =7800 kg/m3, g = 9.8 m/s2,E=2.1e+11 N/m2 Analytical solution of maximum displacement (u2) at x = 50m :U2_max = -((3*r*g*L^4)/(64*E))^(1/3) = -1.194944005 m方法一. 沿 truss element 加沿长度方向初始拉应力 (see job-1.inp)此文件中使用了initial condition, type = stress方法加初始拉应力。
因工程上无此初应力, 更好的方法是使用降温法。
算完后再升温。
用降温法。
算完后再升温。
NOTE: 降温法施加预应力(激活钢绞线)。
温度=-力/(膨胀系数*弹模*钢绞线面积)1、第一步,在truss单元中施加一个初始应力,让计算处于初始平衡状态;初始应力设置过小,可能不收敛,应多次试算,找到一个合理的应力值。
一般情况下,这个初始值对最终值的影响不大,可以忽略。
2、第二步,施加truss单元的自重荷载,打开非线性开关(nlgeom=YES )考虑几何非线性问题;3、本例中初始值采用0.1Mpa。
自重作用下缆索的拉应力约为80Mpa。
最大位移为 -1.195 m,与理论计算值吻合得很好。
*HeadingCable appling gravity load with initail stressThe maximum Analytical displacement without initail stress (at x = 50 m) U2 = -1.194944005 meter**** Method 1. Using * initial condition,type = stress method***Preprint, echo=NO, model=NO, history=NO, contact=NO*Node1, 0., 0.101, 100., 0.*NGEN, NSET = NALL1, 101, 1*Element, type=T2D21, 1, 2*ELGEN, elset = ELALL1, 100, 1*ELSET,ELSET=EL_OUT1, 51, 100*Solid Section, elset=ELALL, material=steel0.01539,***Nset, nset=Left1,*Nset, nset=right101,*Nset, nset=mid51,**** MATERIALS***Material, name=steel*Density7800.,*Elastic2.1e+11, 0.3*initial condition, type = stress** Note: the solution will not converge as the initial stress < 100,000 N/m^2 ELALL, 100000*BoundaryLeft, 1, 2Right, 1,2*Step, name=Step-0, inc=1000Initial stress equilibrium*Static1, 1., 1e-05, 1.*Output, field, variable=PRESELECT*Output, history, variable=PRESELECT*Node print, nset = mid, freq = 1000U,*EL PRINT, ELSET=EL_OUT, freq = 1000S*END STEP** ---------------------------------------------------------------- **** STEP: Step-1***Step, name=Step-1, nlgeom=YES, inc=1000Apply gravity load*Static0.01, 1., 1e-05, 1.** Name: GRAVITY-1 Type: Gravity*DloadELALL, GRAV, 9.8, 0., -1.**** OUTPUT REQUESTS***Restart, write, number interval=1, time marks=NO*Output, field, variable=PRESELECT*Output, history, variable=PRESELECT*Node print, nset = mid, freq = 1000U,*EL PRINT, ELSET=EL_OUT, freq = 1000S*End Step方法二. 使用 STABILIZE parameter on the *STATIC.(see job-2.inp)“stabilization”在结构上附加artificial viscous damping(粘滞阻尼),使得计算结果to go beyond the instability point。
但计算结果必须验证,并必须保证 ALLSD 比内能ALLIE 小很多。
NOTE:1.第一步用*Static, stabilize=2E-10。
笫二步不用stabilize (相当于*Static, stabilize=0)。
2.使用nlgeom=YES in the step to apply the gravity load.3.最终拉应力 = 8E7 N/m2与方法一相等。
4.The maximum displacement (at node 51) equals the analytical result.5.Check the ALLSD and ALLIE. The ALLIE is greater than ALLSD. (See figure 1)6.The deformation shape of the cable can be examined by CAE. It may need to set the deformation scale factor to a large number (10 – 1000).7.使用此法必须极端谨慎。
稍微不慎,结果会完全不对。
For example, 用*Static, stabilize=2E-4 (default value of the stabilize parameter)重算此题。
其结果如下;. Check the ALLSD and ALLIE. The ALLIE is less than ALLSD. (See figure 2)The deformation shape with deformation scale factor 1000 is shown in figure 3. In the figure, only the first and last elements are deformed. The maximum deformation value is not correct.使用 stabilize parameter 学问很多,一般是越小越好。
因为stabilizeparameter 是 artificial value, 无法确定理论上的最佳值。
我是用试错法。
从开始 default value (2.0e-4) 往下减 (2.E-6, 2.0E-8,..),直到不收敛 (2.0E-12). 经过验证结果 (see the note 4, 5, and 6),我决定在计算中使用 2.0E-10。
*HeadingCable apply gravity load using stabilize parameterThe maximum Analytical displacement without initail stress (at x = 50 m) U2 = 1.194944005 meter**** Method 2. Using *Static, stabilize method***Preprint, echo=NO, model=NO, history=NO, contact=NO*Node1, 0., 0.101, 100., 0.*NGEN, NSET = NALL1, 101, 1*Element, type=T2D21, 1, 2*ELGEN, elset = ELALL1, 100, 1*ELSET,ELSET=EL_OUT1, 51, 100*Solid Section, elset=ELALL, material=steel 0.01539,***Nset, nset=Left1,*Nset, nset=right101,*Nset, nset=mid51,**** MATERIALS***Material, name=Steel*Density7800.,*Elastic2.1e+11, 0.3*BoundaryLeft, 1, 2Right, 1,2**** STEP: Step-1***Step, name=Step-1, nlgeom=YES, inc=1000 Apply gravity load** the default value of stabilize value is 2.0E-4**Static, stabilize=2E-4*Static, stabilize=2E-100.01, 1., 1e-05, 1.** Name: GRAVITY-1 Type: Gravity*DloadELALL, GRAV, 9.8, 0., -1.**** OUTPUT REQUESTS***Restart, write, number interval=1, time marks=NO *Output, field, variable=PRESELECT*Output, history, variable=PRESELECTALLSD , ALLIE*Node print, nset = mid, freq = 1000U,*EL PRINT, ELSET=EL_OUT, freq = 1000S*End Step***Step, name=Step-2, nlgeom=YES, inc=1000 Recovery*Static0.01, 1., 1e-05, 1.*End Step。