用PYTHON实现弹性地基梁

nbiao

我是新手的，之前潜水比较多，在坛子里学了很多东西。
从开始学习python也要快一个月了，深感python的在abaqus中的例子很少，帮助文档也不够特别清楚，abaqus自己的python类使用起来也特别麻烦。所以想和大家交流一下，共同提高

因为回答了论坛里面几个朋友关于python的一些问题，被版主zsq-w加分了，感谢版主对我的支持哈！~~
应版主的要求，现在用python做一个弹性地基梁的例子。
这个例子中几个可以看看的地方有：
1、python用连接单元实际单向受压的地基弹簧
2、python中批量选择wire赋连接单元属性
3、python实现集合set的合并
4、python建模的一般思路

现在比较晚了哈，我也有点懒，没有做视图的切换，大家切换一下就可以正常显示模型了。还有就是自己加上荷载就可以计算
附图，弹性地基梁的变形和地基反力




代码：

1. # -*- coding: cp936 -*-
   
2. from abaqus import *
   
3. from abaqusConstants import *
   
4. 
   
5. import connectorBehavior
   
6. import math
   
7. import mesh
   
8. import numpy as np
   
9. import os
   
10. import regionToolset
    
11. import step
    
12. 
    
13. #准备数据
    
14. Hei=0.6  #梁厚度
    
15. pos=0.2  #泊松比
    
16. ela=3.45E7 #弹模
    
17. beam_len = 10 #梁长度
    
18. stif = 1000 #单个地基弹簧刚度
    
19. elastic = (ela , pos)
    
20. #创建模型
    
21. md_1 = mdb.Model(name='Beam on Plastic Foundation')
    
22. #定义材料
    
23. mat = md_1.Material(name='mat')
    
24. mat.Elastic ( table = ( elastic,))
    
25. #定义截面
    
26. mat= md_1.RectangularProfile( name='beam' , a=1,b=Hei)
    
27. 
    
28. #===============================================================================
    
29. #开始画SKETCH
    
30. ##画弹性地基梁
    
31. sk_beam = md_1.ConstrainedSketch(name='beam', sheetSize = 100);
    
32. sk_beam.Line( point1=(-beam_len/2,0), point2=(beam_len/2,0))
    
33. ##细分地基梁
    
34. kkk = sk_beam.geometry.keys()
    
35. while max([sk_beam.geometry.getSize() for i in kkk])>1:
    
36.     for j in kkk:
    
37.         edg=sk_beam.geometry[j]
    
38.         sz = edg.getSize()
    
39.         if sz>1:
    
40.             v1,v2=edg.getVertices()
    
41.             cord1=v1.coords
    
42.             cord2=v2.coords
    
43.             print cord1, cord2
    
44.             mid_pt =( (cord1[0]+cord2[0])/2 , (cord1[1]+cord2[1])/2)
    
45.             sk_beam.Line ( point1=mid_pt, point2=cord2)
    
46.             sk_beam.Line ( point2=mid_pt, point1=cord1)
    
47.             sk_beam.delete( (edg,))
    
48.         kkk = sk_beam.geometry.keys()
    
49. #===============================================================================
    
50. #画完 sketch，添加到part中
    
51. pt_beam = md_1.Part(name='beam',dimensionality=TWO_D_PLANAR,type=DEFORMABLE_BODY)
    
52. pt_beam.BaseWire(sketch=sk_beam)
    
53. #赋截面
    
54. md_1.BeamSection(name='sec_beam', profile='beam',integration=DURING_ANALYSIS,
    
55.                  material = 'mat')
    
56. region = regionToolset.Region(edges = pt_beam.edges)
    
57. pt_beam.SectionAssignment ( region=region, sectionName ='sec_beam',
    
58.                                 offsetType=MIDDLE_SURFACE)
    
59. pt_beam.assignBeamSectionOrientation(region=region, method=N1_COSINES, n1=(0.0, 0.0,-1.0))
    
60. #组装Assembly
    
61. as_1  = md_1.rootAssembly
    
62. ins_1 = as_1.Instance(name='beam',part = pt_beam,dependent=ON)
    
63. ###===============================================================================
    
64. #画wire
    
65. #画AXIAL连接单元,定义受压不受拉的弹簧单元
    
66. rfid1=[]
    
67. rfid2=[]
    
68. mer=1
    
69. a = md_1.rootAssembly
    
70. sec_conn_2 = md_1.ConnectorSection(name='Spring', translationalType=AXIAL)
    
71. elastic_1 = connectorBehavior.ConnectorElasticity(components=(1, ), behavior=NONLINEAR,
    
72.                                                   table=((0.0, 1.0), (0.0, 0.0), (-stif, -1.0)))
    
73. elastic_1.ConnectorOptions()
    
74. mer=1
    
75. v11 = a.instances['beam'].vertices
    
76. csys=a.DatumCsysByThreePoints(origin=(0,0,0), point1= (0,10,0),
    
77.                                  point2=(10,10,0),name='Datum csys', coordSysType=CARTESIAN)
    
78. for v in v11:
    
79.     md_1.sections['Spring'].setValues(  behaviorOptions = (elastic_1, ) )
    
80.     p1on =(v.pointOn[0][0],v.pointOn[0][1]-1,v.pointOn[0][2])
    
81.     prf = a.ReferencePoint(point= p1on)
    
82.     as_1.WirePolyLine(points=((as_1.referencePoints[prf.id],v), ), mergeWire=False, meshable=False)
    
83.     edge1 = as_1.edges.findAt( (p1on, ) )
    
84.     as_1.Set ( name='tmp' , edges=edge1 )
    
85.     region = regionToolset.Region(edges=as_1.sets['tmp'].edges)
    
86.     csa=as_1.SectionAssignment(sectionName='Spring',region=region)
    
87.     dt=as_1.datums[csys.id]
    
88.     as_1.ConnectorOrientation( region=csa.getSet(), localCsys1 = dt)
    
89.     #建地基弹簧的集合
    
90.     if mer:
    
91.         as_1.Set( name='groundConnnector',edges=edge1 )
    
92.         mer=0
    
93.     else:
    
94.         as_1.SetByMerge( name='groundConnnector',sets=(as_1.sets['tmp'],as_1.sets['groundConnnector']))
    
95. del as_1.sets['tmp']
    
96. #给参考点建集合，加约束
    
97. refs = tuple( [as_1.referencePoints for i in as_1.referencePoints.keys()])
    
98. as_1.Set( referencePoints = refs, name='prfs')
    
99. md_1.DisplacementBC(name='BC-prfs',
    
100. createStepName='Initial',region=as_1.sets['prfs'], u1=SET, u2=SET, ur3=SET)
     
101. verts = v11.findAt( ((-beam_len/2,0,0),),((beam_len/2,0,0),) )
     
102. as_1.Set(vertices = verts, name='bc')
     
103. md_1.DisplacementBC(name='fix-bc',
     
104.     createStepName='Initial',region=as_1.sets['bc'], u1=SET, u2=UNSET, ur3=UNSET,
     
105.     amplitude=UNSET, distributionType=UNIFORM, fieldName='', localCsys=None)
     
106. #建立分析步
     
107. as_1.Set(edges=ins_1.edges , name='all')
     
108. st_1 = md_1.StaticStep(name='Step-1', previous='Initial',initialInc=0.1)
     
109. #场变量
     
110. reg=as_1.sets['all']
     
111. md_1.FieldOutputRequest(name='F-Output-3',
     
112. createStepName='Step-1', variables=('U', 'RF', 'SF'),region=reg )
     
113. #画网格
     
114. pt_beam.seedEdgeBySize(edges=pt_beam.edges, size=1, deviationFactor=0.1,
     
115. constraint=FINER)
     
116. elemType1 = mesh.ElemType(elemCode=B23, elemLibrary=STANDARD)
     
117. pt_beam.setElementType ( regions=(pt_beam.edges,),elemTypes=(elemType1,))
     
118. pt_beam.generateMesh()
     
119. #===============================================================================
     
120. #创建JOB
     
121. mdb.Job(name='Job-1', model='Beam on Plastic Foundation')

复制代码

.py文件

前面又看了一下代码，代码里面有一个import numpy as np，我机器上能运行，大家用的话应该会提示错误，把这句删掉就行，因为后面的代码里面没有用到numpy模块
现在提一个问题~~
abaqus能不能用uel实现弹性地基梁，哪位前辈可以给点建议，或者给个例子？

ffubei

很好！该例子可以深度挖掘。

panpan200304

进来学习下啊

guiwh

abaqus是做什么用的？

nbiao

4# guiwh

好可爱啊...来这个论坛逛

水凸凸

mergeWire=False运行时提示关键字错误怎么回事啊