ABAQUS怎么施加位移场

zhanglei830412

         各位大神：
                 小弟有一个问题又要请教大家了。本人最近在做一个地层沉降发生后对隧道结构影响的研究。在数值模拟中，要怎么实现地层沉降呢？本人的想法是通过施加发生的位移场来实现地层的不均匀沉降，那么abaqus怎么施加这么大范围的位移场？求大神答疑解惑···

kxh86

使用Python程序根据每个节点的坐标创建相应的位移载荷或位移约束。也可能需要预定义场。

ilwj25

你怎么知道位移会有多大呢？隧道是被动承受地层沉降产生的作用力的，你分析的就是这个啊。具体怎么做建议先可以看看地应力平衡的帖子。

ilwj25

还有孔隙水压力等问题，要看固结，不然没沉降，呵呵

zhanglei830412

ilwj25 发表于 2012-9-21 08:49
你怎么知道位移会有多大呢？隧道是被动承受地层沉降产生的作用力的，你分析的就是这个啊。具体怎么做建议先 ...

地应力平衡的问题我都了解的，关键是这个项目的这个子问题是模拟地层沉降发生对隧道结构的影响。那么必须有一个自主的沉降位移。这个怎么模拟比较好呢？

zhanglei830412

ilwj25 发表于 2012-9-21 08:50
还有孔隙水压力等问题，要看固结，不然没沉降，呵呵

我要研究的不是隧道开挖引起的地表或地层位移沉降，这个研究的很多了，数值模拟已经没有意义了····

zhanglei830412

kxh86 发表于 2012-9-20 22:23
使用Python程序根据每个节点的坐标创建相应的位移载荷或位移约束。也可能需要预定义场。 ...

不知道大神能否具体点说明···还有对不均匀沉降发生对结构影响的研究，尤其是隧道结构，有什么模拟想法没/

kxh86

zhanglei830412 发表于 2012-9-21 14:10
不知道大神能否具体点说明···还有对不均匀沉降发生对结构影响的研究，尤其是隧道结构，有什么模拟想法 ...

对不起，前面的发言有误。好像是不能用预定义场的方法定义位移场，但是可以用定义BC的方法在初始分析步定义位移。预定义场也是定义在初始分析步的。如果用定义初始分析步的BC的方法做的话，可以用Python程序读取模型的所有节点，根据节点的坐标，创建相应的BC，这样要创建大量的BC。下面是以前自己编的一些Python的片段。
        self.myMdb=Mdb()         #创建Mdb，因为自己定义了类、方法、函数，在一个类的不同方法函数中使用同一个变量，使用了self
        self.myMdl=self.myMdb.Model('Model-1') #创建Model-1
        myMdl=self.myMdl
        myAss=myMdl.rootAssembly        
# Create Assembly; BC nodes
        myCsys=msPart.DatumCsysByThreePoints(coordSysType=CARTESIAN,name='Datum csys-1',
                                             origin=(0, 0, 0), point1=(1, 0, 0), point2=(0, 1, 0)) #创建用户坐标
        CsId=myCsys.id  #获取用户坐标的ID
        self.myDatum=msPart.datums[CsId] #获取用户坐标
        myAss=myMdl.rootAssembly
        myIns=myAss.Instance(dependent=ON, name='SWH', part=msPart) #创建转配实体
        nds=myIns.nodes #获取装配实体有限元模型的所有节点！！！！
        ndx=[] # x=0, Left surface  
        ndr=[] # x=20, Right surface
        ndyz=[] # Nodes for BC-YZ    #预定义三个列表，用于存放节点编号
        for i in nds: # get Node Label
            cd=i.coordinates            
            if cd[0]==-90:          #获取X坐标为-90的节点的编号
                ndx.append(i.label)
            if cd[0]==90:       #获取X坐标为90的节点的编号
                ndr.append(i.label)
            if round(cd[0],3)==-90 and round(cd[1],3)==0 and round(cd[2],3)==0: #获取X=-90、Y=0、Z=0的节点编号
                ndyz.append(i.label)
        ndsx=nds.sequenceFromLabels(tuple(ndx)) # get Nodes #获取节点
        ndsr=nds.sequenceFromLabels(tuple(ndr))   #获取节点
        ndsyz=nds.sequenceFromLabels(tuple(ndyz)) #获取节点
        setNdx=myAss.Set(name='Set-BCX', nodes=ndsx) # Node Set #使用上面的节点编号创建set
        setNdr=myAss.Set(name='Set-NdCp', nodes=ndsr)
        setNdyz=myAss.Set(name='Set-BCYZ', nodes=ndsyz)
        myMdl.DisplacementBC(amplitude=UNSET, createStepName='Initial', name='BC-YZ',                                                     #使用节点Set创建BC
              region=myAss.sets['Set-BCYZ'], u1=UNSET, u2=SET, u3=SET, ur1=UNSET, ur2=UNSET, ur3=UNSET)
        myMdl.DisplacementBC(amplitude=UNSET, createStepName='Initial', name='BC-X',
              region=myAss.sets['Set-BCX'], u1=SET, u2=UNSET, u3=UNSET, ur1=UNSET, ur2=UNSET, ur3=UNSET)
        myMdl.DisplacementBC(createStepName='Step-1', name='BC-RP', region=myAss.sets['Set-RP'],
                             u1=-1, u2=0, u3=0, ur1=0, ur2=0, ur3=0)

zhanglei830412

kxh86 发表于 2012-9-21 14:39
对不起，前面的发言有误。好像是不能用预定义场的方法定义位移场，但是可以用定义BC的方法在初始分析步定 ...

不知道大神是土木那个方向的？你觉得我要模拟不均匀沉降对隧道结构影响的思路，采用施加位移场是可行的吗？还有其他好的思路吗？还有python语音我还不是很明白，这方面可能还有很多问题要问你··方便留个联系方法不？

kxh86

zhanglei830412 发表于 2012-9-21 14:46
不知道大神是土木那个方向的？你觉得我要模拟不均匀沉降对隧道结构影响的思路，采用施加位移场是可行的吗 ...

我硕士用Abaqus做的复合材料结构渐进损伤强度分析，C919相关的；现在做的是太阳能电池板，给天宫一号配套的。对土木不太了解，只是从有限元原理上能理解一点。用Abaqus/CAE建模的过程会用Python语言保存在JNL文件中，jnl文件保存了所有建模步骤及参数设置，是过程化的程序，可以修改成结构化的程序已完成一些循环、判断等功能。将jnl文件的扩展名改为py，就是Python文件，可以用Abaqus run script直接运行。Run Script在菜单栏第一个菜单下面，也可在启动abaqus时运行py文件。在Abaqus二次开发板块应该有很多贴子介绍Python在Abaqus。有本书《Python语言在Abaqus中的应用》。

zhanglei830412

kxh86 发表于 2012-9-21 15:28
我硕士用Abaqus做的复合材料结构渐进损伤强度分析，C919相关的；现在做的是太阳能电池板，给天宫一号配套 ...

大神能不能说一下，这个施加位移场的思路。是不是获取到有限元模型中所有节点的坐标，然后施加相应的位移边界···那么如果我通过其它方法得到了所有节点的坐标，可以再inp文件中施加位移边界，施加初始位移吗？

kxh86

zhanglei830412 发表于 2012-9-21 18:33
大神能不能说一下，这个施加位移场的思路。是不是获取到有限元模型中所有节点的坐标，然后施加相应的位移 ...

在inp中也可以实现，我没用过。

zhanglei830412

求救········· 小弟有一个问题又要请教大家了。本人最近在做一个地层沉降发生后对隧道结构影响的研究。在数值模拟中，要怎么实现地层沉降呢？本人的想法是通过施加发生的位移场来实现地层的不均匀沉降，那么abaqus怎么施加这么大范围的位移场？求大神答疑解惑···

zhanglei830412

大神能不能说一下，这个施加位移场的思路

zhanglei830412

现在这个课题这方面的研究内容初步定为：已发生不均匀变形条件下地铁结构的受力和变形的计算。请问数值模拟该怎么施加这个不均匀变形···

zhanglei830412

kxh86 发表于 2012-9-21 19:02
在inp中也可以实现，我没用过。

问个小问题。麻烦赐教，我在第一步中施加了一个重力荷载，在第2步施加一个位移荷载，那么这个位移荷载是相对于初始位移的大小，还是相对于重力荷载产生的位移之后的大小？

kxh86

zhanglei830412 发表于 2012-9-23 22:08
问个小问题。麻烦赐教，我在第一步中施加了一个重力荷载，在第2步施加一个位移荷载，那么这个位移荷载是 ...

要看Step-2中有没有重力载荷作用。我的理解是：在Load Manage中看Step-2的状态，如果是Propagated，就是Step-1的力载荷在Step-2中继续作用；如果选择Propagated并单击Deactive，使Step-2状态为Inactive，则Step-1中的力载荷不作用于Step-2。我不确定是否准确，楼主可以做个只有几个单元的小模型验证。

zhanglei830412

kxh86 发表于 2012-9-23 22:24
要看Step-2中有没有重力载荷作用。我的理解是：在Load Manage中看Step-2的状态，如果是Propagated，就是S ...

谢谢，太感谢你了。也问了一下其他人，有说跟OP有关的。我自己做个拉杆试试。还有个问题。在后处理中，ABAQUS图形可以不显示单元划分，只在实体上显示应力和变形不？
跟下图相似，而没有一个个小单元，怎么设置，能教下么？

kxh86

在Visualization模块工具栏单击第一个工具Common Option，在Common Plot Options对话框中第一个标签页Basic中最右边的Visible Edges选择No edges。

kxh86

zhanglei830412 发表于 2012-9-23 22:32
谢谢，太感谢你了。也问了一下其他人，有说跟OP有关的。我自己做个拉杆试试。还有个问题。在后处理中，AB ...