然后现在是建模的原始文件:
def in_def
xx = in(msg+'('+'default:'+string(default)+'):')
if type(xx) = 3
in_def = default
else
in_def = xx
endif
end
;
def beam_canshu
default = 6.0
msg=' yuan gong banjing (r) '
r = in_def
a13=r
b13=0
a1=-a13
b1=0
a12=r*cos(pi/12)
b12=r*sin(pi/12)
a2=-a12
b2=b12
a11=r*cos(pi/6)
b11=r*sin(pi/6)
a3=-a11
b3=b11
a10=r*cos(pi/4)
b10=r*sin(pi/4)
a4=-a10
b4=b10
a9=r*cos(pi/3)
b9=r*sin(pi/3)
a5=-a9
b5=b9
a8=r*cos((pi*5)/12)
b8=r*sin((pi*5)/12)
a6=-a8
b6=b8
a7=0
b7=r
;
default = -5e4
msg=' cross force thickness1 (p0) '
p0 = in_def
;
default = -1e5
msg=' cross force xielv (l0) '
l0 = in_def
;
default = 0
msg=' vertical force xielv (c0) '
c0 = in_def
;
default = 5e04
msg=' horizontal crustal stress xiaojianju (m0) '
m0 = in_def
;
default = 1e06
msg=' horizontal crustal stress yuanlice (n0) '
n0 = in_def
;
f1=(p0+l0*r)*r*sin(pi/12)
f2=((p0+l0*r)*r*sin(pi/6))/2+(l0*sin(pi/12)*2)/3+(l0*(sin(pi/6)-sin(pi/12)))/3
f3=((p0+l0*r)*r*(sin(pi/4)-sin(pi/12)))/2+(l0*(sin(pi/6)-sin(pi/12))*2)/3+(l0*(sin(pi/4)-sin(pi/6)))/3
f4=((p0+l0*r)*r*(sin(pi/3)-sin(pi/6)))/2+(l0*(sin(pi/4)-sin(pi/6))*2)/3+(l0*(sin(pi/3)-sin(pi/4)))/3
f5=((p0+l0*r)*r*(sin((pi*5)/12)-sin(pi/4)))/2+(l0*(sin(pi/3)-sin(pi/4))*2)/3+(l0*(sin((pi*5)/12)-sin(pi/3)))/3
f6=((p0+l0*r)*r*(1-sin(pi/3)))/2+(l0*(sin((pi*5)/12)-sin(pi/3))*2)/3+(l0*(1-sin((pi*5)/12)))/3
f7=((p0+l0*r)*r*sin(pi/6))/2-(l0*sin(pi/12)*2)/3-(l0*(sin(pi/6)-sin(pi/12)))/3
f8=((p0+l0*r)*r*(sin(pi/4)-sin(pi/12)))/2-(l0*(sin(pi/6)-sin(pi/12))*2)/3-(l0*(sin(pi/4)-sin(pi/6)))/3
f9=((p0+l0*r)*r*(sin(pi/3)-sin(pi/6)))/2-(l0*(sin(pi/4)-sin(pi/6))*2)/3-(l0*(sin(pi/3)-sin(pi/4)))/3
f10=((p0+l0*r)*r*(sin((pi*5)/12)-sin(pi/4)))/2-(l0*(sin(pi/3)-sin(pi/4))*2)/3-(l0*(sin((pi*5)/12)-sin(pi/3)))/3
f11=((p0+l0*r)*r*(1-sin(pi/3)))/2-(l0*(sin((pi*5)/12)-sin(pi/3))*2)/3-(l0*(1-sin((pi*5)/12)))/3
;
c1=(c0*(1-cos(pi/12))*2+c0*(cos(pi/6)-cos(pi/12)))/3+(n0)*r(1-cos(pi/6))
c2=(c0*(cos(pi/6)-cos(pi/12))*2+c0*(cos(pi/4)-cos(pi/6)))/3+(n0)*r*(cos(pi/4)-cos(pi/12))
c3=(c0*(cos(pi/4)-cos(pi/6))*2+c0*(cos(pi/3)-cos(pi/4)))/3+(n0)*r*(cos(pi/3)-cos(pi/6))
c4=(c0*(cos(pi/3)-cos(pi/4))*2+c0*(cos((pi*5)/12)-cos(pi/3)))/3+(n0)*r*(cos((pi*5)/12)-cos(pi/4))
c5=(c0*(cos((pi*5)/12)-cos(pi/3))*2+c0*(1-cos((pi*5)/12)))/3+(n0)*r*cos(pi/3)
c6=-((c0*(cos((pi*5)/12)-cos(pi/3))*2+c0*(1-cos((pi*5)/12)))/3+(m0)*r*cos(pi/3))
c7=-((c0*(cos(pi/3)-cos(pi/4))*2+c0*(cos((pi*5)/12)-cos(pi/3)))/3+(m0)*r*(cos((pi*5)/12)-cos(pi/4)))
c8=-((c0*(cos(pi/4)-cos(pi/6))*2+c0*(cos(pi/3)-cos(pi/4)))/3+(m0)*r*(cos(pi/3)-cos(pi/6)))
c9=-((c0*(cos(pi/6)-cos(pi/12))*2+c0*(cos(pi/4)-cos(pi/6)))/3+(n0)*r*(cos(pi/4)-cos(pi/12)))
c10=-((c0*(1-cos(pi/12))*2+c0*(cos(pi/6)-cos(pi/12)))/3+(n0)*r(1-cos(pi/6)))
cc=((n0-m0)*(1-cos(pi/12)))/2
end
beam_canshu
sel beam id=1 begin=( a1,b1,0 ) end=( a2,b2,0 ) nesg=1
sel beam id=1 begin=( a2,b2,0 ) end=( a3,b3,0 ) nesg=1
sel beam id=1 begin=( a3,b3,0 ) end=( a4,b4,0 ) nesg=1
sel beam id=1 begin=( a4,b4,0 ) end=( a5,b5,0 ) nesg=1
sel beam id=1 begin=( a5,b5,0 ) end=( a6,b6,0 ) nesg=1
sel beam id=1 begin=( a6,b6,0 ) end=( a7,b7,0 ) nesg=1
sel beam id=1 begin=( a7,b7,0 ) end=( a8,b8,0 ) nesg=1
sel beam id=1 begin=( a8,b8,0 ) end=( a9,b9,0 ) nesg=1
sel beam id=1 begin=( a9,b9,0 ) end=( a10,b10,0 ) nesg=1
sel beam id=1 begin=( a10,b10,0 ) end=( a11,b11,0 ) nesg=1
sel beam id=1 begin=( a11,b11,0 ) end=( a12,b12,0 ) nesg=1
sel beam id=1 begin=( a12,b12,0 ) end=( a13,b13,0 ) nesg=1
plot set rot 90 0 0
sel beam id=1 prop emod=52.4e9 nu=0.30 xcarea=3.55e-3 xcj=400e-6 xciy=333e-6 xciz=67e-6
sel node fix z xr yr
sel node fix x y zr range id=1
sel node fix x y zr range id=13
sel node apply force=(cc,f1,0.0) range id=7
sel node apply force=(c1,f2,0.0) range id=8
sel node apply force=(c2,f3,0.0) range id=9
sel node apply force=(c3,f4,0.0) range id=10
sel node apply force=(c4,f5,0.0) range id=11
sel node apply force=(c5,f6,0.0) range id=12
sel node apply force=(c6,f11,0.0) range id=2
sel node apply force=(c7,f10,0.0) range id=3
sel node apply force=(c8,f9,0.0) range id=4
sel node apply force=(c9,f8,0.0) range id=5
sel node apply force=(c10,f7,0.0) range id=6
his sel node ydisp id=7
his sel beam force fx end1 cid=2
his sel beam force fx end1 cid=4
his sel beam force fx end1 cid=6
his sel beam force fx end1 cid=7
his sel beam force fx end1 cid=9
his sel beam force fx end1 cid=11
plot create geom_and_disp
plot set title text
Geometry of the modeled system
plot set cap size 25
plot set rot 90 0 0
plot add sel geom black black
plot add axes black
plot add sel disp blue
plot create shear_force
plot set title text
Shear force distribution in the beam
plot set cap size 25
plot set rot 90 0 0
plot add sel geom black black node=off
plot add sel beam force fy scale=0.06
plot create axis_force
plot set title text
axis force in the beam
plot set cap size 25
plot set rot 90 0 0
plot add sel geom black black node=off
plot add sel beam force fx
plot create moment
plot set title text
Bending moment distribution in the beam
plot set cap size 25
plot set rot 90 0 0
plot add sel geom black black node=off
plot add sel beam moment mz scale=0.06
print sel beam id
print sel beam local
print sel beam position
solve ratio=1e-7
print sel beam force
利用FLAC及正交实验实现支护荷载反分析
来源:互联网 发布日期:2011-09-16 10:46:52 浏览:5880次
导读:土工作者的家园,最具人气的岩土工程社区,本站点是一个面向广大岩土工程工作者的岩土工程综合站点,主要内容有岩土工程若干方面的技术问题,包括常用岩土工程规范...
相关内容
- “5G+工业互联网”将迎来升级版政策,专家建议大规模推进工业企业数字化转型
- 中国工程院院士邬贺铨:5.5G呼之欲出 目标并不是取代5G
- 华为携手联通完成5G-A技术落地 机构建议积极关注射频等投资机会
- 华星创业:持续探索5G和5.5G应用领域的发展机会
- 高通公司取得5G空中用户设备模式切换专利,实现UE标高超过阈值时切换到定向发送模式
- 刚刚,中东宣布携手迈向 5G-A 商用元年
- 立昂技术:公司将切入和研发基于5G应用的的云网融合产品
- 华为本月底将有重磅发布,5.5G建设共识度提升,产业链迎来新机遇
- 5G搭台,创新唱戏:爱立信使能差异化服务推动商业成功
- 政策鼓励先行先试 5G消息规模应用或提质提速
AiLab云推荐
最新资讯
本月热点
热门排行
-
华为5G-A实现规模组网,下载速度超过4Gbps,比5G快10倍
阅读量:56297
-
最坏时期已过,美国5G市场即将迎来回暖?
阅读量:56079
-
回望2023年5G发展:释放和提升价值成最大课题
阅读量:13077
-
全球吹响5G-A商用号角
阅读量:12143
-
5.5G商用来了!多了半个G,和你我有啥关系?
阅读量:6200
-
从5G迈向5G-A,运营商的增长新机会在哪?| MWC 2024
阅读量:6038
推荐内容
- 2024冻立方中国(厦门)冷冻冷藏食品展
- 2024第二十一届中国畜牧业展览会(畜博会CAHE)
- 2024第十三届中国(西安)糖酒食品交易会(西安糖酒会)
- 2024西部乡村振兴博览会暨农业机械装备展
- 2024年第28届伊朗国际石油、天然气、炼油和石化展览会
- 中东石油展 2024年第28届伊朗国际石油展
- 2024重庆国际餐饮(火锅)食材展览会
- 2024深圳国际渔业博览会
- 2024华晨第十一届中国(山东)秋季国际宠物水族展
- 2024中亚一带一路工博会
- 2024第七届雄安城市建设及绿色建筑博览会 第七届雄安装配式建筑及绿色建材展览会 2024中国绿色建筑与智能建造科技大会 第五届雄安金桥梁项目对接会
- 2024第七届中国(新疆)国际智能交通产业博览会
- 2024越南(河内)国际五金及工具展览会
- 2024山东连锁加盟展览会 开年山东首展
- 2024武汉国际电子元器件、材料及生产设备展览会
- 2024山东(济南)连锁加盟博览会助力行业发展
- CRFE国际连锁加盟展览会 -山东招商加盟专业展
- 2024第五届深圳国际老龄博览会
- 2024年山西国际水展览会
- 2024哈尔滨东北亚国际口腔设备及器材展览会
- 2024第二十一届中国(漯河)食品博览会
- 2024第十届中国(北京)国际警用装备及反恐技术装备展览会暨学术研讨会 (CIPATE)
- 2024孟加拉国际暖通、空调制冷、净化、楼宇自控展