经典示例：HPMESH划分箱梁、降温法施加预应力与DLOAD移动载荷、生死单元应用、加强筋

Lxcad

成功效验移动载荷通过箱梁（比如运梁车）：
HPMESH划分箱梁、降温法施加预应力与DLOAD移动载荷、生死单元应用我一步步演示ABAQUS涉及功能，希望大家喜欢，作为中秋前的奉献。全部hm文件inp与dload文件均提供。
（补充：其中还可以用初应变施加预应力）

[ 本帖最后由 Lxcad 于 2007-9-18 18:30 编辑 ]

Lxcad

总结经典示例：HPMESH划分箱梁、降温法施加预应力与DLOAD移动载荷、生死单元应用

思路：先用HPMESH对箱梁网格划分，设置边界条件set组，元素表面组，附8.0hm文件
杀死预应力钢绞线
加载重力
降温法施加预应力（激活钢绞线）。温度=-力/（膨胀系数*弹模*钢绞线面积）
DLOAD移动载荷

陈火红老师二次开发书中MARC中有一移动载荷的例子，我综合其他ABA的技术。以上的示例文件。大家能不能鼓掌鼓励一下。

[ 本帖最后由 Lxcad 于 2007-9-18 15:05 编辑 ]

Lxcad

降温法施加预应力（激活钢绞线）。温度=-力/（膨胀系数*弹模*钢绞线面积）
施加预应力后发生起拱现象。
钢绞线通过植入法。

[ 本帖最后由 Lxcad 于 2007-9-18 15:13 编辑 ]

Lxcad

将DLOAD改一下，就可以验证2个移动载荷了。   


       SUBROUTINE DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,COORDS,
     1                 JLTYP,SNAME)
C
      INCLUDE 'ABA_PARAM.INC'
C
      DIMENSION TIME(2),COORDS(3)
      CHARACTER*80 SNAME
C
      dis=TIME(1)
    x1min=80.d0-35.d0
    x1max=80.d0+35.d0
    x2min=-80.d0-35.d0
    x2max=-80.d0+35.d0
    ymin=12.5d0-dis
    ymax=-12.5d0-dis
      F=0.0d0
    if(COORDS(1).le.x1max.and.COORDS(1).ge.x1min.and.
    1COORDS(2).le.ymin.and.COORDS(2).ge.ymax) then
    F=75.d4/(24.d0*25.d0)
    endif
    if(COORDS(1).le.x2max.and.COORDS(1).ge.x2min.and.
    1COORDS(2).le.ymin.and.COORDS(2).ge.ymax) then
    F=75.d4/(24.d0*25.d0)
    endif
c     调试专用
c    if(TIME(1).ge.60.0d0.and.COORDS(1).le.xmax.and.
c    1COORDS(1).ge.xmin) then
c    ii=1
c    endif
      RETURN
      END

Lxcad

网格又分了，可以在箱梁随心所欲设置车辆。运行该inp大约30分钟！附所有文件inp,hpmesh,dload。

Lxcad

DLOAD移动载荷
运行后可以看到箱梁顶的载荷一步步移动，很棒的约！
有谁告诉我AVI文件制作为GIF动画，我定将它做成gif。

1. 
   
2. SUBROUTINE DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,COORDS,
   
3.      1                 JLTYP,SNAME)
   
4. C
   
5.       INCLUDE 'ABA_PARAM.INC'
   
6. C
   
7.       DIMENSION TIME(2),COORDS(3)
   
8.       CHARACTER*80 SNAME
   
9. C
   
10.       dis=TIME(1)
    
11.     xmin=-35.d0
    
12.     xmax=35.d0
    
13.     ymin=35.d0-dis
    
14.     ymax=-35.d0-dis
    
15.       F=0.0d0
    
16.     if(COORDS(1).le.xmax.and.COORDS(1).ge.xmin.and.
    
17.     1COORDS(2).le.ymin.and.COORDS(2).ge.ymax) then
    
18.     F=75.d4/(24.d0*50.d0)
    
19.     endif
    
20.     if(TIME(1).ge.60.0d0.and.COORDS(1).le.xmax.and.
    
21.     1COORDS(1).ge.xmin) then
    
22.     ii=1
    
23.     endif
    
24.       RETURN
    
25.       END
    

复制代码

Lxcad

思路：先用HPMESH对箱梁网格划分，设置边界条件set组，元素表面组，附8.0hm文件

Lxcad

箱梁hpmesh网格。

Lxcad

预应力网格。

Lxcad

加载重力

Lxcad

inp文件有点心痛呀！

Lxcad

去掉网格的inp文件
*Surface, type=ELEMENT, name=top_sur
top, S1
*Solid Section, elset=press, material=STEEL
3.68155,
** Region: (Section-2-SOLIDMAP:SOLIDMAP)
** Section: Section-2-SOLIDMAP
*Solid Section, elset=solidmap, material=CON
1.,
** Constraint: Embed-1
*Embedded Element
press
**
** MATERIALS
**
*Material, name=CON
*Density
2.6e-05,
*Elastic
3.5e+06, 0.1667
*Material, name=STEEL
*Density
7.85e-05,
*Elastic
1.95e+07, 0.3
*Expansion
1e-05,
** ----------------------------------------------------------------
**
** STEP: Step-1
**
*Step, name=Step-1
*Static
1., 1., 1e-05, 1.
*model change,type=element,remove
press
**
** BOUNDARY CONDITIONS
**
** Name: Disp-BC-1 Type: Displacement/Rotation
*Boundary
bom_l, 1, 1
** Name: Disp-BC-2 Type: Displacement/Rotation
*Boundary
bom_l, 2, 2
** Name: Disp-BC-3 Type: Displacement/Rotation
*Boundary
bom_l, 3, 3
** Name: Disp-BC-4 Type: Displacement/Rotation
*Boundary
bom_r, 2, 2
** Name: Disp-BC-5 Type: Displacement/Rotation
*Boundary
bom_r, 3, 3
**
** LOADS
**
** Name: Load-1   Type: Gravity
*Dload
, GRAV, 980., 0., 0., -1.
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End Step
** ----------------------------------------------------------------
**
** STEP: Step-2
**
*Step, name=Step-2
*Static
1., 1., 1e-05, 1.
*model change,add=with strain
press
**
** FIELDS
**
** Name: Field-1   Type: Temperature
*Temperature
pre_node, -764.
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End Step
** ----------------------------------------------------------------
**
** STEP: Step-3
**
*Step, name=Step-3, inc=80
*Static, direct
10., 800.,
**
*Dsload
top_sur, PNU, 1.
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End Step

kiang

第一个支持，down下来看看

useful

好文，鼓掌致谢！

jinchengsimu

鼓掌表示感谢！

Lxcad

谢谢DAVA给分和楼上2位的鼓励！

lxm9977

小乐不一般那，好贴好贴。比上月饼强多了！ ^_^

Lxcad

原帖由 lxm9977 于 2007-9-18 15:42 发表
小乐不一般那，好贴好贴。比上月饼强多了！ ^_^

六小斑竹、Q斑竹捧场，我中秋节太高兴了！

cip05yh

ding ding ding

garethz

小乐版主出手就是不一样啊！！
好帖！

Allen-2000

支持斑竹好贴！

wujian84

多谢帮助了