flac3d50结构单元教程
flac 结构单元用法
（1）cable单元模拟锚杆（全长连接、非全长连接，预应力锚杆），土钉等
（2）beam单元模拟支撑、围檩、冠梁等
（3）liner单元模拟面层，地下连续墙
Fix x range x -0.1 0.1
Fix x range x 5.9 6.1
Fix y range y -0.1 0.1
Fix y range y 7.9 8.1
Fix z range z -0.1 0.1
set mech force=50
Gen zone brick size 6 8 8
Mode mohr
Prop bulk 1e8 shear 0.3e8 fric 35
Prop cohesion 1e10 tens 1e10 ;注意在此这个值给的很大。
Init dens 1000
Set gravity 0 0 -10
1、 得到初始应力的方法：
方法1、可以先给一些材料参数很大的值，进行初始求解，在计算之前再将材料参数设为正常值，即可。如在手册中给的第一个示例中就是这样做的。下面是例子，These are only initial values that are used during the development of gravitational stresses within the body. In effect, we are forcing the body to behave elastically during the development of the initial in-situ stress state.* This prevents any plastic yield during the initial loading phase of the analysis.
（新）Flac3D5.0笔记
（新）Flac3D5.0笔记FLAC3D5.0导入dat、txt文件的方法FLAC3D5.0导入.sav文件的方法输入命令流：Restore 文件名.savFLAC3D5.0导入.f3prj文件的方法页脚内容1视图窗口关闭后，如何再次呈现？输入命令plot zone页脚内容2当命令流输错可从这里撤销页脚内容3点亮相应命令流后才可撤销透明度的调整保存项目文件页脚内容4保存文件的重命名做切片页脚内容5倾向和倾角法向和过一个点删除切片Ctrl+R为查看模型后还原显示接触面，双击I nterface显示接触面页脚内容6接触面属性项被激活后，可以更改属性隐藏网格显示输出网格页脚内容7如果点击了透明，将不显示Information 只要不点击透明，则显示Information 页脚内容8等同于点击显示节点坐标页脚内容9fish简介@之前要有（空格）之前加@是Flac3D3.0和5.0区别而在之间没有（空格）按下列步骤取消fish语言前要加@页脚内容10点击fish页脚内容11软件的刷新如何显示软件最终的安全系数显示软件节点.查看应力及最大剪应力的步骤页脚内容12查看位移云图的步骤页脚内容13输出云图页脚内容14设置出图边宽设置分辨率可设为1080 Reflect镜像restore调用之前的计算结果Origin原点normal法向量（该点与原点的坐标差）页脚内容15对于结构单元如桩，若点击Zone将不出现任何模型按如下步骤才可显示使用该命令后保存的文件与说建模的文件保存在同一文件夹页脚内容16查看约束反力的步骤查看分组的步骤位移云图的显示颜色相反如何调整页脚内容17连接两个不同尺寸的单元页脚内容18页脚内容19。
flac3d5.0结构单元教程
1. FLAC3D V5.0界面操作 2. FLAC3D基本操作方法vs应用流程； 3. FLAC3D内置Fish语言的应用； 4. FLAC3D结构单元vs接触单元；
5. FLAC3D渗流模块
6. 其他
StructuralElement
FLAC3D结构单元
1. 2. 3. 4.
4.1、简支梁（beam单元）承受两个相等集中载 ; 荷
========================================================== ; Setup histories for monitoring behavior. history add id=10 sel node ydisp id=7 history add id=30 sel beamsel moment mz end2 cid=1 ; moment, right of SEL-1 history add id=31 sel beamsel moment mz end1 cid=2 ; moment, left of SEL-2 ; ======================================================== ; Bring the problem to equilibrium 如何设置结构单元的跟踪变量！ solve ratio=1e-7 save equal-concent-loads ; ======================================================== ; Print out beam responses. list sel beam force list sel beam moment list sel node disp range id=7 return
flacd结构单元教程ppt文档
• 土工格栅(geogrid)
通过连接实现岩土体或结构与其它结构发生
• 衬砌单元(liner)
相互作用。
➢ 注意：结构节点并不是简单地与实体网格的
节点(gridpoint)建立联系，也不能建立node
与gridpoint之间的link
2、结构单元的建模方法
➢ 梁单元
sel beam id 1 beg 4 0 -1 end 5 0 -2 nseg 4
sel beam id=1 begin=( 6, 0, 0) end=( 9, 0, 0) nseg=3
; ======================================================
; Assign beam properties
sel beam id=1 prop emod=2e11 nu=0.30 &
2、结构单元的建模方法—线型结构单元
先建立节点再联接成单元的方法；
2、结构单元的建模方法—壳型结构单元
➢ 壳单元
2、结构单元的建模方法—壳型结构单元
def set_vals global ptA = 25.0 * sin( 40.0*degrad ) ; global ptB = 25.0 * cos( 40.0*degrad ) end @set_vals generate zone cylinder p0=( 0.0, 0.0, 0.0 ) &
; Print out beam responses.
list sel beam force
list sel beam moment
list sel node disp range id=7
return
FLAC3D50模型及输入参数说明
1.1模型参数代码可参考manual中各个章节的command命令及说明,注意单位。
用prop 赋值。
16 ttable 塑性拉应变-抗拉强度的表号下列参数可以显示、绘图与通过fish访问1 es_plastic 塑性切应变2 et_plastic 塑性拉应变3 ff_count 检测切应变反向的数4 ff_cvd 体应变,εvd经典粘弹性模型的材料参数(Classical Viscoelastic (Maxwell Substance) –MODEL mechanical viscous)1 bulk 弹性体积模量,K2 shear 弹性剪切模量,G3 viscosity 动力粘度,η粘弹性模型的材料参数(Burgers Model –MODEL mechanical burgers)1 bulk 弹性体积模量,K2 kshear Kelvin弹性剪切模量,G K3 kviscosity Kelvin动力粘度,ηK4 mkshear Maxwell切边模量,G M5 mviscosity Maxwell动力粘度,ηM二分幂律模型的材料参数(Power Law –MODEL mechanical power)1 a_1 常数,A12 a_2 常数,A23 bulk 弹性体积模量,K4 n_1 指数,n15 n_2 指数,n26 rs_1 参考应力,σ1ref7 rs_2 参考应力,σ2ref8 shear 弹性剪切模量,G蠕变模型材料参数(WIPP Model –MODEL mechanical wipp)1 act_energy 活化能,Q2 a_wipp 常数,A3 b_wipp 常数,B4 bulk 弹性体积模量,K5 d_wipp 常数,D6 e_dot_star临界稳定状态蠕变率,7 gas_c 气体常数,R8 n_wipp 指数,n9 shear 弹性剪切模量,G10 temp 温度,T下列参数可以显示、绘图与通过fish访问1 e_prime 累积主蠕变应变2 e_rate 累积主蠕变应变率Burger、蠕变组合材料模型的材料参数(Burgers-Creep Viscoplastic Model –MODEL mechanical cvisc)1 bulk 弹性体积模量,K2 cohesion 内聚力,c3 density 密度,ρ4 dilation 剪胀角,Ψ5 friction 内摩擦角,Φ6 kshear Kelvin弹性剪切模量,G K7 kviscosity Kelvin粘度,ηK8 shear 弹性剪切模量,G9 tension 抗拉强度,σt10 mviscosity Maxwell动力粘度,ηM下列计算参数可以显示、绘图与通过fish访问1 es_plastic 累积塑性切应变2 et_plastic 累积塑性拉应变幂律模型的材料参数(Power-Law Viscoplastic Model –MODEL mechanical cpower)1 a_1 常数,A12 a_2 常数,A23 bulk 弹性体积模量,K4 cohesion 内聚力,c5 dilation 剪胀角,Ψ6 friction 内摩擦角,Φ7 n_1 指数,n18 n_2 指数,n29 rs_1 参考应力,σ1ref10 rs_2 参考应力,σ2ref11 shear 弹性剪切模量,G12 tension 抗拉强度,σt粘塑形模型的材料参数(WIPP-Creep Viscoplastic Model –MODEL mechanical pwipp)1 act_energy 活化能,Q2 a_wipp 常数,A3 b_wipp 常数,B4 bulk 弹性体积模量,K5 d_wipp 常数,D6 e_dot_star临界稳定状态蠕变率,7 gas_c 气体常数,R8 kshear 材料参数,KΦ9 n_wipp 指数,n10 kdil 材料参数,q k11 kvol 材料参数,qΦ12 shear 弹性切变模量,G13 temp 温度,T14 tension 抗拉强度,σt以下计算参数可以显示、绘图与通过fish访问1 e_prime 累积主蠕变应变2 e_rate 累积主蠕变应变率3 es_plastic 累积塑性切应变4 et_plastic 累积塑性拉应变碎盐变形模型的材料参数(Crushed-Salt Model –MODEL mechanical cwipp)1 act_energy 活化能,Q2 a_wipp 常数,A3 b_f 最终体积模量,K f4 b_wipp 常数,B5 b0 蠕变压实系数,B06 b1 蠕变压实系数,B17 b2 蠕变压实系数,B28 bulk 弹性体积模量,K9 d_f 最终密度,ρf10 d_wipp 常数,D11 e_dot_star临界稳定状态蠕变率,12 gas_c 气体常数,R13 n_wipp 指数,n14 rho 密度,ρ15 s_f 最终切变模量,G f16 shear 弹性切变模量,G17 temp 温度,T以下计算参数可以显示、绘图与通过fish访问1 frac_d 当前碎片密度,ρd2 s_g1 蠕变压实参数,G3 s_k1 蠕变压实参数,K均质流体模型的材料参数1 permeability 等方向渗透性,k2 porosity 孔隙率,n(默认时,n=0、5)各向异性流体模型的材料参数1 fdd k1-k2的平面倾向2 fdip k1-k2的平面倾角1.2模型适用说明遍布节理模型适用于Mohr-Coulomb材料来明确显示力在各个方向上的差异性。
2024版FLAC3D5.0培训
06
总结与展望
本次培训总结
培训内容丰富
涵盖了FLAC3D5.0的基本原理、 建模方法、分析步骤、后处理等 多个方面，使学员能够全面了解
并掌握该软件的使用。
培训方式多样
采用了理论讲解、案例分析、实 践操作等多种培训方式，使学员 在理论学习的基础上，通过实践 操作加深了对软件的理解和掌握。
培训效果显著
程中的应力、变形和稳定性。
02
支护结构设计与优化
根据隧道开挖模拟结果，设计合理的支护结构，如锚杆、喷射混凝土等，
并利用FLAC3D5.0对支护结构进行优化。
03
隧道施工风险评估
基于FLAC3D5.0的模拟结果，对隧道施工过程中可能出现的风险进行评
估，提出相应的应对措施。
基坑开挖与支护设计
基坑开挖过程模拟 利用FLAC3D5.0建立基坑三维模型，模拟基坑的开挖过程， 分析开挖过程中的应力、变形和稳定性。
高效建模技巧
利用对称性简化模型
对于具有对称性的结构，可以只建立一半或四分之一的模型，通过设置对称边界条件来模拟 整个结构，从而大大提高建模效率。
使用模板快速创建复杂模型
FLAC3D5.0提供了丰富的模板库，用户可以直接调用模板来创建复杂的模型，避免了繁琐的 建模过程。
批量修改模型参数
通过编写脚本或使用内置工具，可以实现对模型参数的批量修改，提高建模效率。
边界条件设置方法
根据实际问题的要求，设置合理的边界条件。对于 固定边界，可将其节点位移约束为零；对于自由边 界，可不施加任何约束。同时，还需考虑边界条件 的对称性和周期性等因素。
网格密度控制
根据计算精度和计算效率的要求，合理控制网格的 密度。在关键区域和应力集中区域可采用较密的网 格，以提高计算精度。
(2024年)FLAC3D5.0培训教程
精度和计算效率的需求。
2024/3/26
13
接触面处理及摩擦模拟
2024/3/26
接触面定义
01
支持定义不同材料之间的接触面，包括摩擦系数、刚度等参数
设置。
接触面行为模拟
02
能够模拟接触面的滑动、张开和闭合等行为，以及接触面间的
传热和传质过程。
动画展示技巧 在制作动画时，可采用一些技巧来提高动画的展示效果， 如使用透明度渐变来突出关键区域的变化、使用色彩对比 来区分不同物理量的分布情况等。
结果数据对比 在动画制作中，可将不同方案或不同时间步的计算结果进 行对比展示，以便更直观地评估不同方案的效果或观察模 型的动态响应过程。
25
06
总结与展望
21
05
数据可视化与后处理
2024/3/26
22
结果数据输出格式
文本文件输出
可将模型计算结果以文本文件形式输出，方便用户进行自定义处理 和分析。
Excel文件输出
可将模型计算结果直接导出到Excel文件中，便于用户进行数据整理、 分析和可视化。
图像文件输出
可将模型计算结果以图像形式输出，如等值线图、云图等，方便用户 进行直观分析和展示。
施方法
学习在FLAC3D中施加边界条件和 荷载的方法，确保模拟过程的真实 性。
11
03
高级功能与技巧
2024/3/26
12
复杂模型处理技术
复杂地形建模
利用地形数据生成三维地形模型， 包括不规则地形、断层、节理等。
复杂结构建模
支持多种结构单元，如梁、板、 壳等，实现复杂结构的精细化建
FLAC3D50模型及输入参数说明
材料模型及关键字说明（model）MODEL关键字说明力学模型anisotropic横向同性弹性模型cam-clay修正剑桥模型doubleyield D-Y模型drucker德鲁克普拉格模型elactic各向同性弹性模型mohr摩尔-库伦模型null空模型orthotropic正交各向异性弹性模型ss应变硬化/软化模型subliquitous双线性应变硬化/软化多节理模型ubiquitous多节理模型finn纯动力学模型蠕变模型burger Burger材料粘弹性模型cpower幂律材料模型cvisc Burger、蠕变组合材料模型cwipp岩盐变形模型power二分幂律模型pwipp粘塑性模型viscous经典粘弹性模型wipp蠕变模型流体模型fl_anistropic各向异性流体模型fl_isotropic各向同性流体模型fl_null空流体模型热力模型th_anistropic均质热导模型th_null空热导模型Example:model mohr range x=1,10 y=1,10 z=1,10 （group name）prop bulk 25.0e9 shear 5.2e9 coh 0.23e6 dil 20.0 frict 38.0 tension 0.0e61.1模型参数代码可参考manual中各个章节的command命令及说明，注意单位。
用prop 赋值。
各向同性弹性材料参数（Isotropic Elastic – MODEL mechanical elastic）1bulk体积模量，或松散系数K2shear切边模量，G横向同行弹性模型材料参数（Transversely Isotropic Elastic – MODEL mechanical anisotropic）1dd同性平面的倾向2dip同性平面的倾角3e1同性平面的弹性模量4e3垂直同性平面的弹性模量5g12切变模量6nu12同性平面内施力时的泊松比7nu13垂直同性平面内施力时的泊松比正交各向异性弹性模型正交各向异性弹性模型材料参数（Orthotropic Elastic – MODEL mechanical orthotropic）1dd轴1-‘2’所定义平面的倾向2dip轴1-‘2’所定义平面的倾角3e11’方向弹性模量4e22’方向弹性模量5e33’方向弹性模量6g12平行于轴1-‘2’平面的切变模量7g13平行于轴1-‘3’平面的切变模量8g23平行于轴2-‘3’平面的切变模量9nu12沿2’方向施力，1’方向的泊松比10nu13沿3’方向施力，1’方向的泊松比11nu23沿3’方向施力，2’方向的泊松比12nx轴1-‘2’所定义平面单位法线x分量13ny轴1-‘2’所定义平面单位法线y分量14nz轴1-‘2’所定义平面单位法线z分量15rot旋转角德鲁克-普拉格模型德鲁克-普拉格模型的材料参数（Drucker-Prager – MODEL mechanical drucker）1bulk弹性体积模量，K2ksnear材料参数，kφ3qdil材料参数，qφ4qvol材料参数，qΨ5shear弹性切变模量，K6tension抗拉强度，σt摩尔-库伦模型摩尔-库伦模型的材料参数（Mohr-Coulomb – MODEL mechanical mohr）1bulk弹性体积模量，K2cohesion内聚力，c3dilation剪胀角，Ψ4friction内摩擦角，Φ5shear弹性切边模量，G6tension抗拉强度，σt多节理模型的材料参数（Ubiquitous-Joint – MODEL mechanical ubiquitous）1bulk弹性体积模量，K2cohesion内聚力，c3dilation剪胀角，Ψ4friction内摩擦角，Φ5jcohesion节理内聚力，c j6jddirection弱面dip方向（倾向）7jdilation节理剪胀角，Ψj8jdip弱面dip角度（倾角）9jfriction节理摩擦角，Φj10jnx弱面单位法线x分量11jny弱面单位法线y分量12jnz弱面单位法线z分量13jtension抗拉强度，σsj14shear弹性切边模量，G15tension抗拉强度，σt应变硬化/软化模型应变硬化/软化模型的材料参数（Strain-Hardening/Softening – MODEL mechanical ssoftening）1bulk弹性体积模量，K2cohesion内聚力，c3ctable塑性剪切应变-内聚力的表号4diation剪胀角，Ψ5dtable塑性剪切应变-剪胀角的表号6friction内摩擦角，Φ7ftable塑性剪切应变-摩擦角的表号8shear弹性切边模量，G9tension抗拉强度，σs10ttable塑性剪切应变-抗拉强度的表号双线性应变硬化/软化多节理模型的材料参数（Bilinear Strain-Hardening/Softening Ubiquitous-Joint –MODEL mechanical subiquitous）1bijoint=0，为线性节理，默认=1，为双线性节理2bimatrix=0，线性矩阵=1，双线性矩阵3bulk弹性体积模量，K4c2table塑性剪切应变-内聚力c2的表号5cjtable节理塑性剪切应变-节理内聚力c j1的表号6cj2table节理塑性剪切应变-节理内聚力c j2的表号7cohesion内聚力，c18co2内聚力，c29ctable塑性剪切应变-内聚力c1的表号10d2table塑性剪切应变-剪胀角Ψ2的表号11di2剪胀角，Ψ212dilation剪胀角，Ψ113djtable节理塑性剪切应变-节理剪胀角Ψj1的表号14dj2table节理塑性剪切应变-节理剪胀角Ψj2的表号15dtable塑性剪切应变-剪胀角Ψ1的表号16f2table塑性剪切应变-摩擦角Φ2的表号17fjtable节理塑性剪切应变-节理摩擦角Φj1的表号18fj2table节理塑性剪切应变-节理摩擦角Φj2的表号19fr2摩擦角Φ220friction摩擦角Φ121ftable塑性剪切应变-摩擦角Φ1表号22jc2节理内聚力c j123jcohesion节理内聚力c j224jddirection弱面dip方向（倾向）25jdilation节理剪胀角Ψj126jdip弱面dip角度（倾角）27jd2节理剪胀角，Ψj128jfriction节理摩擦角，Φj129jf2节理摩擦角，Φj230jnx弱面单位法线x分量31jny弱面单位法线y分量32jnz弱面单位法线z分量33jtension节理抗拉强度，σsj34shear弹性切边模量，G35tension抗拉强度，σs36tjtable节理塑性剪切应变-节理抗拉强度σtj的表号37ttable节理塑性剪切应变-节理抗拉强度σt的表号下列参数可以显示、绘图或者fish访问1es_plastic塑性切应变2et_plastic塑性拉应变3etj_plastic节理塑性拉应变4esj_plastic节理塑性切应变D-Y模型的材料参数（Double-Yield – MODEL mechanical doubleyield）1bulk弹性体积模量，K2cap_pressure cap压力，p c3cohesion内聚力，c4cptable塑性体应变-cap压力的表号5ctable塑性切应变-内聚力的表号6dilation剪胀角，Ψ7dtable塑性切应变-剪胀角的表号8ev_plastic塑性体应变总量9friction内摩擦角，Φ10ftable塑性切应变-摩擦角的表号11multiplier当前塑性-cap模量与弹性体积和切变模量的倍数，R12shear最大弹性切变模量，G13tension抗拉强度，σs14ttable塑性拉应变-抗拉强度的表号下列计算参数可以显示、绘图和通过fish访问1es_plastic累积塑性切应变2et_plastic累积塑性拉应变3ev_plastic累积塑性体应变修正剑桥模型修正剑桥模型的材料参数（Modified Cam-Clay – MODEL mechanical cam-clay）1bulk_bcund最大的弹性体积模量，K max2cv初始容积，v c3kappa弹性膨胀线斜率，k4lamda常态固结线斜率，λ5mm摩擦常数，M6mpc预固结压力，p c07mp1预固结压力，p18mv_1指定在参考压力常态固结线的容积v A 9poisson泊松比，ν10shear弹性剪切模量，G下列参数可以显示、绘图以及fish访问1bulk体积模量，K2cam_cp当前平均有效应力3cam_ev累积总容积应变4camev_cp累积塑性容积应变5cq当前平均差分应力纯动力学模型的材料参数（config dynamic, model mech finn）1bulk弹性体积模量，K2cohesion内聚力，c3ctable弹性切应变-内聚力的表号4dilation剪胀角，Ψ5dtable塑性切应变-剪胀角Ψ的表号6ff_c1常量，c17ff_c2常量，c28ff_c3常量，c39ff_c4常量，c410ff_latensy反向之间的最小时间步数11ff_switch=0：Martin（1995）公式=1：Byme（1991）公式12friction内摩擦角，Φ13ftable塑性切应变-摩擦角的表号14shear最大弹性切变模量，G15tension抗拉强度，σt16ttable塑性拉应变-抗拉强度的表号下列参数可以显示、绘图和通过fish访问1es_plastic塑性切应变2et_plastic塑性拉应变3ff_count检测切应变反向的数4ff_cvd体应变，εvd经典粘弹性模型经典粘弹性模型的材料参数（Classical Viscoelastic （Maxwell Substance）– MODEL mechanical viscous）1bulk弹性体积模量，K2shear弹性剪切模量，G3viscosity动力粘度，η粘弹性模型粘弹性模型的材料参数（Burgers Model – MODEL mechanical burgers）1bulk弹性体积模量，K2kshear Kelvin弹性剪切模量，G K3kviscosity Kelvin动力粘度，ηK4mkshear Maxwell切边模量，G M5mviscosity Maxwell动力粘度，ηM二分幂律模型二分幂律模型的材料参数（Power Law – MODEL mechanical power）1a_1常数，A12a_2常数，A23bulk弹性体积模量，K4n_1指数，n15n_2指数，n26rs_1参考应力，σ1ref7rs_2参考应力，σ2ref8shear弹性剪切模量，G蠕变模型材料参数（WIPP Model – MODEL mechanical wipp）1act_energy活化能，Q2a_wipp常数，A3b_wipp常数，B4bulk弹性体积模量，K5d_wipp常数，D6e_dot_star临界稳定状态蠕变率，ϵ̇ss∗7gas_c气体常数，R8n_wipp指数，n9shear弹性剪切模量，G10temp温度，T下列参数可以显示、绘图和通过fish访问1e_prime累积主蠕变应变2e_rate累积主蠕变应变率Burger、蠕变组合材料模型Burger、蠕变组合材料模型的材料参数（Burgers-Creep Viscoplastic Model – MODEL mechanical cvisc）1bulk弹性体积模量，K2cohesion内聚力，c3density密度，ρ4dilation剪胀角，Ψ5friction内摩擦角，Φ6kshear Kelvin弹性剪切模量，G K7kviscosity Kelvin粘度，ηK8shear弹性剪切模量，G9tension抗拉强度，σt10mviscosity Maxwell动力粘度，ηM下列计算参数可以显示、绘图和通过fish访问1es_plastic累积塑性切应变2et_plastic累积塑性拉应变幂律模型幂律模型的材料参数（Power-Law Viscoplastic Model – MODEL mechanical cpower）1a_1常数，A12a_2常数，A23bulk弹性体积模量，K4cohesion内聚力，c5dilation剪胀角，Ψ6friction内摩擦角，Φ7n_1指数，n18n_2指数，n29rs_1参考应力，σ1ref10rs_2参考应力，σ2ref11shear弹性剪切模量，G12tension抗拉强度，σt粘塑形模型的材料参数（WIPP-Creep Viscoplastic Model – MODEL mechanical pwipp）1act_energy活化能，Q2a_wipp常数，A3b_wipp常数，B4bulk弹性体积模量，K5d_wipp常数，D6e_dot_star临界稳定状态蠕变率，ϵ̇ss∗7gas_c气体常数，R8kshear材料参数，KΦ9n_wipp指数，n10kdil材料参数，q k11kvol材料参数，qΦ12shear弹性切变模量，G13temp温度，T14tension抗拉强度，σt以下计算参数可以显示、绘图和通过fish访问1e_prime累积主蠕变应变2e_rate累积主蠕变应变率3es_plastic累积塑性切应变4et_plastic累积塑性拉应变碎盐变形模型碎盐变形模型的材料参数（Crushed-Salt Model – MODEL mechanical cwipp）1act_energy活化能，Q2a_wipp常数，A3b_f最终体积模量，K f4b_wipp常数，B5b0蠕变压实系数，B06b1蠕变压实系数，B17b2蠕变压实系数，B28bulk弹性体积模量，K9d_f最终密度，ρf10d_wipp常数，D11e_dot_star临界稳定状态蠕变率，ϵ̇ss∗12gas_c气体常数，R13n_wipp指数，n14rho密度，ρ15s_f最终切变模量，G f16shear弹性切变模量，G17temp温度，T以下计算参数可以显示、绘图和通过fish访问1frac_d当前碎片密度，ρd2s_g1蠕变压实参数，G3s_k1蠕变压实参数，K均质流体模型均质流体模型的材料参数1permeability等方向渗透性，k2porosity孔隙率，n（默认时，n=0.5）各向异性流体模型的材料参数1fdd k1-k2的平面倾向2fdip k1-k2的平面倾角3frot k1轴和倾角矢量的转角4h1k1方向的渗透性5h2k2方向的渗透性6h3k3方向的渗透性均质热导模型的材料参数1conductivity等方向传热系数，K2spec_heat比热容，C v1.2模型适用说明遍布节理模型适用于Mohr-Coulomb材料来明确显示力在各个方向上的差异性。
FLAC,FLAC3D基础与应用(结构单元)教学提纲
连接 SELs
● 错误！
11
连接 SELs
● 正确！
12
Liner结构单元
● 三节点扁平有限单元 ● 每个节点有6个自由度
3个移动，3个旋转
● 能够抵抗膜及弯矩荷载 ● 能够承受主方向的拉压应力 ● 能够模拟管片与土体之间的分离及随后的重新接触 ● 能够模拟管片与土体之间的摩擦相切互向作用
法向
● 全连接
相邻的Liner单元在连接处共用一 个节点，连接处重叠单元不能 发生移动或旋转
● 结点连接
即结点间的连接在6个方向的自 由度上用弹簧来模拟，每个自 由度都可具有一定的特性
通缝拼接
错缝拼接
16
连接方式的影响
● 变形结果
zr属性更改 注：变形放大200倍
6属性更改
17
结构的动力响应
config dyn
• Link可以与任何位置的
grid进行联系，而不一定
要与grid的坐标一致。 群桩 = 插秧
7
默认的连接属性
建模SEL结构模型时，程序自动建立结构 node与zone的连接 (node-zone links)
8
Node-Node Links
● SEL nodes 之间不会自动生成联系. ● 必须手动设置node之间的联系 (e.g., beam and cable) 这
CS_nGap off
sel node fix x y z xr yr zr ran id=1
sel set damp combined def f1
pile
whilestepping
f0=10000*sin(10*dytime)
np = nd_head
loop while np # null
Flac3D5.0教程
Flac3D5.0教程第⼆章完成第⼀个简单分析计算样例by typing the commands from the keyboard, pressingat the end of each command line, and seeing the results directly.newgen zone brick size 6 8 8（This command will create an initial mesh that is six zones in the x-direction, eight zones in the y-direction and eight zones in the z-direction. For our model, the z-axis is oriented in the verticaldirection.）Plot（the Plot Base/0>prompt will appear. As long as this prompt appears, any subsequentcommands will be associated with the PLOT command. The plot view is identified as Base/0,which is the default view.）create Trenchadd surface yellowadd axes blackshowNow, give the zones a material model and properties. For this example, we use the Mohr-Coulombelastic-plastic model. Go back to the Flac3D>prompt in the command window andtype in thefollowing command:modelmohrThis will specify the Mohr-Coulomb model. Every zone in the grid could conceivably have adifferent material model and property. However, by not specifying a range of zonesafter the MODEL command, FLAC3D assumes that all zones are to be Mohr-Coulomb material.prop bulk = 1e8 shear = 0.3e8 fric = 35propcoh = 1e10 tens = 1e10You see that very high cohesion and tensile strength values are given. These areonlyinitial values that are used during the development of gravitational stresses within the body.In effect, we are forcing the body to behave elastically during the development of theinitial insitustress state.* This prevents any plastic yield during the initial loading phase of theanalysis.For this problem, loading is due to gravity. To apply gravity, use the commands setgrav 0, 0, -9.81ini dens = 1000In order to develop a gravitational body force, themass density must also be initialized. The INI command is used to initialize the mass density to1000 kg/m3 for all zones in the model.Next, the boundary conditions for the problem are set. At the Flac3D>prompt, type fix x range x -0.1 0.1fix x range x 5.9 6.1fix y range y -0.1 0.1fix y range y 7.9 8.1fix z range z -0.1 0.1With these commands, roller boundaries are placed on five sides of the model. The boundaries are “fixed” only in the specified direction (i.e., no displacement or velocity is allowed). The FIX commands perform the following functions.1. The gridpoints along the boundary planes at x = 0 and x = 6 are fixed in the x-direction. These two planes fall within the coordinate ranges specified by the range keywords for the first two FIX commands.2. The gridpoints along the boundary planes at y = 0 and y = 8 are fixed in the y-direction. These planes fall within the ranges specified for the third and fourth FIX commands.3. The gridpoints along the bottom boundary (z = 0) are fixed in the z-direction. This plane falls within the range for the fifth FIX command.We wish to monitor the change in the values of selected variables in the model during the calculational stepping. A HISTORY command can assist in helping us determinewhether a stable equilibrium solution or unstable collapse is occurring. We type the following commands:hist n = 5histunbalhistgpzdisp 4,4,8We choose to monitor the change in variables every five calculation steps. It is always a good idea to monitor the maximum unbalanced force in a model. If the unbalanced force approaches a very small value and displacement histories become constant, this indicates that an equilibrium state is reached.To allow gravitational stresses to develop within the body, we timestep the simulation to equilibrium.Here the SOLVE command is used to detect equilibrium automatically.setmech force=50solveWhen the unbalanced force falls below the limiting value (a limiting force of 50 N is specified with the SET command), the run will stop.* The plots are updated, since they are still visible on the screen. Shutting down the plots will cause the model to cycle faster.plothist 1hist 2The unbalanced force history approaches zero, andthe displacement history becomes constant; both are indicators that an equilibriumstate has beenreached.Note that each history is numbered sequentially from 1 as it is entered via the HIST command. Return to the Flac3D> prompt and typeprinthistfor a listing of the histories and their corresponding numbers.Plotcreat trenchadd contour dispadd axes blackshowclearaddbcontourszzadd axesplot create GravVplot set plane dip=90 dd=0 origin=3,4,0plot set rot 15 0 20plot set center 2.5 4.2 4.0plot add bound behindplot add bcontszz planeplot add axesplot showThis sequence will create a view, which we have called GravV, and make it the current view. We then set a plane for that view oriented at a dip angle of 90? (from the xy-plane, assuming that negative-z is “down”), a dip direction of 0?(measured clockwise from the positive y-axis in the xy-plane) and with one point on the plane at (x = 3, y = 4, z = 0). We add a wire-frame boundary plotted behind the plane and a block contour plot of the vertical stress component, σzz, on the plane. Finally, the model axes are added. The block contour plot, as opposed to an interpolated contour plot, displays the value of the stress calculated at each zone centroid. The color of each zone corresponds directly to the zone-based stress.It is wise to save the initial state so that you can restore it at any time to perform parameter studies.savetrench.savWe can change the current view from GravV to Trench with the commandplot current TrenchNow we can excavate a trench in the soil. First, typepropcoh=1e3 tens=1e3This will set the cohesion and tensile strength for all zones to 1000 Pa. These valuesfor strength arehigh enough to prevent failure in our initial state (i.e., unexcavated), but you shouldalways checkfor possible failure in the initial state by performing a few calculation steps. Toexcavate the trench,entermodel null range x=2,4 y=2,6 z=5,10With a low cohesion and vertical unsupported trench walls, collapse should occur. Because wewant to examine this process realistically, the large-strain logic is specified. This isdone by typingset largeFor plotting purposes, we wish to see only the change in displacements from the trench excavation, and not from the previous gravitational loading, so we zero out the x-, y- and z-displacement components:inixdis=0 ydis=0 zdis=0We purposely set the cohesion low enough to result in failure, so we do not want to use the SOLVE command with a limit for out-of-balance force (which checks for equilibrium). Oursimulation willnever converge to the equilibrium state. Instead, we can step through the simulationprocess onetimestep at a time, and plot and print the results of the collapse as it occurs. This is thereal powerof the explicit method. The model is not required to converge to equilibrium at eachcalculationcycle because we never have to solve a set of linear algebraic equationssimultaneously, as is thecase for implicit codes, with which many engineers are familiar.we use the STEP command: step 2000第三章FLAC 3D基础知识gridpointzonehorizontal boundary stressstructural cables(tiebacks)model boundaryinternal boundaries(excavation)roller bottom boundaryNamed ObjectsMacro Object—Typically, such an object contains a long, complex string that may be used repeatedly in the model. The pre-processor compares a string of command tokens to the list of defined macros and replaces any matching macro object with its fully expanded contents.macro Pt0 ’p0 0 0 0’macro Pt1 ’p1 add 10 0 0’macro Pt2 ’p2 add 0 10 0’macro Pt3 ’p3 add 0 0 10’macroModel_Size ’size 4 5 6’macroBig_Brick ’zone brick Pt0 Pt1 Pt2 Pt3 Model_Size’.genBig_Brick.macro ’Pt0’ ’p0 15 15 15’genBig_BrickThis pre-processing has two effects:(1) macro objects may be nested (but not recursively); and(2) the macro object name is removed from the command string.Model Object—Model objects, such as ranges in space, groups of zones in a model, or plot views, can be given userdefined names. Those objects can then be referred to by their names.gen zone brick size 6 6 6group Tunnel range x 1 5 y 0 6 z 1 5modelmohr..modelnul range group Tunnel.The named range and the named group are two very different model objects. The range pertains toa specified volume of space (or range of values), whereas the group identifies acollection of finitedifference zones in the model.Sign ConventionsDIRECT STRESS —Positive stresses indicate tension; negative stresses indicate compression.第四章初级实体建模技术Grid generation is performed via the GENERATE command and associated keywords that both define the number of zones in a model and shape the grid to fit a specified problem region.The number of zones is specified by the size keyword.It is best to start with a grid that has few zones (say, 1000 to 1500) to perform simple test runs andmake refinements to the grid. Then, increase the number of zones to improve theaccuracy.Using actual coordinatesgen zone brick size 6,8,8 p0 -10, -10, -20 &p1 10, -10, -20 &p2 -10, 10, -20 &p3 -10, -10, 0plot surfOnly four corners are required to define a parallelepiped-shaped mesh. More corners can be specifiedto define an irregular surface. Example 2.14 shows how to make a sloping surface atthe top of themesh.gen zone brick size 6,8,8 p0 -10, -10, -20 &p1 10, -10, -20 p2 -10, 10, -20 &p3 -10, -10, 0 p4 10, 10, -20 &p5 -10, 10, 10 p6 10, -10, 0 &p7 10, 10, 10plot surfIn the tutorial example, we noted that the boundaries of the model were influencing the results (see Figure 2.6). The boundary must be placed far enough away from the excavation to reduce these effects. A gradually graded mesh can be created in FLAC3D to move the model boundaries farther out without significantly increasingthe number of zones. For example, the command GEN zone radbrick creates a radially graded mesh around a brick-shaped mesh. The command in Example 2.15 creates a 3 ×5 ×5 zone brick-shaped mesh surrounded by a 7-zone radially graded mesh.gen zone radbrick&p0 (0,0,0) p1 (10,0,0) p2 (0,10,0) p3 (0,0,10) &size 3,5,5,7 &ratio 1,1,1,1.5 &dim 1 4 2 fillplot surfThe x,y,z coordinate system in FLAC3D is always right-handed and, as a default, the z-axis is drawn in the vertical direction on the screen. In Example 2.14 we assumed the z-axis was pointing in the vertical direction. However, we do not have to interpret the z-axis to mean the up direction. For Example 2.15, we will assume that the y-axis is the vertical direction. As long as we define the model in a right-handed system, we can create the grid in any direction we desire.The size keyword, as used in Example 2.15, specifies the number of zones in the brick and the number radially surrounding the brick (see Figure 2.11). The keyword ratio is given to manipulate the grid spacing. The four values that follow ratio are geometric ratios between successive zone sizes. The first three values are ratios for the zones in the brick, and the fourth value is the ratio for the zones surrounding the brick. In Example 2.15 above, the ratio is 1.0 for the brick zones and1.5 for the zones surrounding the brick. This causes each successive ring of zones surrounding the brick to be 1.5 times larger in the radial direction. The dim keyword defines the dimensions of the brick region (i.e., 1 m ×4 m ×2 m). The fill keyword fills the brick region with zones. If fill is omitted, no zones will be generated within the brick region.Notice that size, ratio and dim all refer to local axes as defined by p0, p1, p2, p3 in Figure 2.11, not to the global x,y,z-axes.Creating a model by reflecting elements on planes of symmetrygen zone radbrick&p0 (0,0,0) p1 (10,0,0) p2 (0,10,0) p3 (0,0,10) &size 3,5,5,7 &ratio 1,1,1,1.5 &dim 1 4 2 fillgen zone reflect dip 0 dd 90gen zone reflect dip 90 dd 90plot surfdip 为倾⾓的意思，其值⼤⼩为⾯与xy⾯所夹的倾⾓，dd为倾向，其值⼤⼩为⾯的法线在xy平⾯上的投影与y轴的夹⾓。
04_FLAC3D5.0_结构单元和接触面
3、结构单元的参数取值
某些结构单元参数的取值要视具体情况 而定，根据经验且必要时调整参数通过 试算来确定。
4、结构单元实例分析
• 4.1、简支梁（beam单元）承受两个相等集 中载荷
• 4.2、简支梁（shell单元）承受两个相等集 中载荷
4.1、简支梁（beam单元）承受两个相等集中 载荷
Simple Beam – Two Equal Concentrated Loads
2、结构单元的建模方法—线型结构单元
先建立节点再联接成单元的方法；
2、结构单元的建模方法—壳型结构单元
壳单元
2、结构单元的建模方法—壳型结构单元
def set_vals global ptA = 25.0 * sin( 40.0*degrad ) ; global ptB = 25.0 * cos( 40.0*degrad ) end @set_vals generate zone cylinder p0=( 0.0, 0.0, 0.0 ) & p1=( @ptA, 0.0, @ptB ) & p2=( 0.0, 25.0, 0.0 ) & p3=( 0.0, 0.0, 25.0 ) & p4=( @ptA, 25.0, @ptB ) & p5=( 0.0, 25.0, 25.0 ) & size=(1, 2, 2) sel shell id=5 range cylinder end1=(0.0, 0.0,0.0) & end2=(0.0,25.0,0.0) radius=24.5 not plot add zg plot ad sel geom delete zones ; delete all zones sel node init zpos add -25.0
FLAC3D基本原理和应用特点第三讲FLAC3D动力分析、自定义本构以及结构单元
30
Friction=30, hr=0.43, Go=440
Damping Ratio (%)
20
10
0
0.0001
0.001
0.01
0.1
1
Shear Strain Amplitude (%)
14
14
等效线性方法
• 等效线性方法是岩土地震工程中模拟波的传播的最常用的方法。 • 假定土体是粘弹性体，参照实验室得到的切线模量及阻尼比与剪应变幅值的关
From Seed & Idriss (1970)
Go = small-strain
27
shear modulus
27
滞后阻尼
特定的曲线可以使用滞后阻尼 来描述，并在本构模型计算中 根据应变的大小获得不同的切 线模量。
表观应变是偏应变，根据先前 反转的点进行累计得到的。使 用堆栈的数据结构可以描述应 变反转现象。
19
19
弹塑性模型
简单的理想弹塑性本构模型仅仅在发生屈服时才会出现滞回特性
stress strain
注意：
即使这样粗糙的模型也能够作出连续的 阻尼比和模量衰减曲线。
在屈服条件下会产生 体积改变，但通常都 是剪胀。
20
20
FLAC3D中的阻尼比 1. 弹塑性本构模型使用瑞利（粘性）阻尼 2. 弹塑性本构模型使用滞后（HD）阻尼. 3. 复杂本构模型拥有连续的屈服应力应变关系和对应的加卸载
2. 不能计算永久变形。等效线性方法模型在加荷与卸荷时模量相同，不能计 算土体在周期荷载作用下发生的剩余应变或位移。
3. 塑形屈服模拟不合理。在塑性流动阶段，普遍认为应变增量张量是应力张 量的函数，称之为“流动法则”。然而，等效线性方法使用的塑性理论认为应 变张量（而不是应变增量张量）是应力张量的函数。因此，塑性屈服的模 拟不合理。
flac3d教程
flac3d教程
FLAC3D是一种常用的三维有限差分软件，用于地质工程、岩土力学和地下空间开发等领域的数值模拟。
该软件具有强大的土体和岩体模拟能力，可以模拟地表沉降、岩石崩塌、地下水渗流等复杂地质现象。
使用FLAC3D进行模拟需要按照以下步骤进行操作：
1. 创建模型：首先要创建一个FLAC3D模型文件，可以通过几何建模软件或文本编辑器创建一个文本文件，并使用FLAC3D的特定语法定义模型的几何形状和参数。
2. 设定材料参数：在模型中定义岩土体的物理和力学参数，例如密度、弹性模量、摩擦角等。
这些参数将在模拟过程中用于计算岩土体的应力和变形。
3. 定义边界条件：为模型设置边界条件，如固支、自由表面、初始应力等。
这些边界条件将在模拟中约束模型的行为。
4. 施加荷载：根据实际情况为模型施加相应的荷载，例如施加地震力、垂直载荷等。
可以根据需要在模拟过程中改变或删除荷载。
5. 运行模拟：使用FLAC3D软件运行模拟，计算模型在荷载作用下的应力和变形响应。
模拟可以在软件界面中进行，也可以通过命令行方式进行。
6. 分析结果：模拟完成后，可以通过FLAC3D软件提供的各种功能和工具来分析模型的结果。
例如，绘制应力云图、位移云图、剪切云图等，以及输出模型的计算数据。
需要注意的是，在使用FLAC3D进行模拟时，应根据具体问题进行合理的模型设计和参数设定，并且进行准确的边界条件设置。
同时，还需要对模拟结果进行合理分析和解释，以得出有关工程或地质现象的结论。
flac结构单元用法
flac结构单元用法本人没有做过桩锚或桩＋支撑的围护形式，不知道pile单元在靠近开挖面一侧当土体被开挖时，pile于土体的接触会是什么样子？（1）cable单元模拟锚杆（全长连接、非全长连接，预应力锚杆），土钉等（2）beam单元模拟支撑、围檩、冠梁等（3）liner单元模拟面层，地下连续墙（4）pile单元模拟钻孔灌注桩，SMW工法，超前支护（注浆钢管、微型桩），抗滑桩等。
一般的搅拌桩可用zone＋interface解决。
FLAC学习总结一个在使用flac3d的兄弟(QQ:65006196)1、得到初始应力的方法：方法1、可以先给一些材料参数很大的值，进行初始求解，在计算之前再将材料参数设为正常值，即可。
如在手册中给的第一个示例中就是这样做的。
下面是例子，These are only initial values that are used during the development of gravitational stresses within the body. In effect, we are forcing the body to behave elastically during the development of the initial in-situ stress state.* This prevents any plastic yield during the initial loading phase of the analysis.Gen zone brick size 6 8 8Mode mohrProp bulk 1e8 shear 0.3e8 fric 35Prop cohesion 1e10 tens 1e10 ;注意在此这个值给的很大。
Init dens 1000Set gravity 0 0 -10Fix x range x -0.1 0.1Fix x range x 5.9 6.1Fix y range y -0.1 0.1Fix y range y 7.9 8.1Fix z range z -0.1 0.1set mech force=50solve;---------------------- mode null ---------------------Prop coh 1e3 tens 1e3 ;改为正常值(在此例中我们故意给小值) Mode null range x 2 ,4 y 2 , 6 z 5, 10Set largeIni xdis 0 ydis 0 zdis 0 ;清零，不影响结果，为画图方便。
FLAC3D50模型及输入参数说明
1.1 模型参数代码可参考manual 中各个章节的command 命令及说明,注意单位。
用prop 赋值。
1.1.1 各向同性弹性模型1.1.2 横向同行弹性模型1.1.3 正交各向异性弹性模型1.1.4 德鲁克-普拉格模型1.1.5 摩尔 - 库伦模型1.1.6 多节理模型1.1.7 应变硬化/软化模型1.1.8 双线性应变硬化/软化多节理模型1.1.9 D-Y1.1.10 修正剑桥模型1.1.11 纯动力学模型1.1.12 经典粘弹性模型1.1.131.1.14 二分幂律模型1.1.151.1.16 Burger 、蠕变组合材料模型1.1.17 幂律模型1.1.18 粘塑形模型1.1.19 碎盐变形模型1.1.20 均质流体模型1.1.21 各向异性流体模型1.1.22 均质热导模型1.2 模型适用说明遍布节理模型适用于Mohr-Coulomb 材料来明确显示力在各个方向上的差异性。
双线性软化应变遍布节理模型综合了软化应变Mohr-Coulomb 模型与遍布节理模型,这种模型包含面向矩阵与遍布节理的一个双线性断裂点集。
改进的Cam-clay 模型反映了形变度与抗破坏能力对体积变化的影响。
Mohr-Coulomb 模型最适用于一般工程研究,同时,Mohr-Coulomb 的内聚力与摩擦角参数相对于地质工程材料的其它属性,更容易获得。
软化应变与遍布节理塑性模型实际上就是Mohr-Coulomb 模型的变形,这些模型如果在附加材料参数的值较高时将得出与Mohr-Coulomb 模型同样的结果。
Druck-Prager 模型就是一个相对于Mohr-Coulomb 模型的破坏标准的简化体,但就是它一般不适于用来描述地质工程材料的破坏情况。
它主要就是用来把FLAC3D 与其它一些有Druck-Prager 模型但却没有Mohr-Coulomb 模型的数学软件作比较。
在摩擦力为零的时候请注意,此时Mohr-Coulomb 模型退化为Tresca 模型,而Druck-Prager 模型退化为Von Mises 模型。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
3、结构单元的参数取值
➢ 梁单元
➢ 锚索单元
• emod——弹性模量，E
• emod——弹性模量， E
• nu——泊松比，ν
• xcarea——横截面积，A
• xcarea——横截面积，A
FLAC3D 5.0培训日程安排
1. FLAC3D V5.0界面操作 2. FLAC3D基本操作方法vs应用流程； 3. FLAC3D内置Fish语言的应用； 4. FLAC3D结构单元vs接触单元； 5. FLAC3D渗流模块 6. 其他
StructuralElement
FLAC3D结构单元
1. 结构单元的类型 2. 结构单元的建模方法 3. 结构单元的参数取值 4. 结构单元实例分析 5. 关于link
两种建模方式各有 各的优点，第二种 方式适合建立复杂 曲线结构单元（但 是要注意它不会自 动建立link！！若不 手动link就无任何作 用）
sel beamsel cid=3 id=1 node 3 4
➢ 桩单元
sel pile id 1 beቤተ መጻሕፍቲ ባይዱ 0 0 0 end 0 0 10 nseg 4
2、结构单元的建模方法
➢ 锚索单元 sel cable id 1 beg 4 0 -1 end 5 0 -2 nseg 4
2、结构单元的建模方法—线型结构单元
起始点坐标并给定分段数目的方法； 建立梁单元，并显示 单元坐标系！
2、结构单元的建模方法—线型结构单元
ID号相同，共用Node，ID不同，各个ID对应的结构单元有各自独立的node。除非设 置联系，否则即使节点位于同一位置也不会传递力。
p1=( @ptA, 0.0, @ptB ) & p2=( 0.0, 25.0, 0.0 ) & p3=( 0.0, 0.0, 25.0 ) & p4=( @ptA, 25.0, @ptB ) & p5=( 0.0, 25.0, 25.0 ) & size=(1, 2, 2) sel shell id=5 range cylinder end1=(0.0, 0.0,0.0) &
节点(gridpoint)建立联系，也不能建立node
与gridpoint之间的link
2、结构单元的建模方法
➢ 梁单元
sel beam id 1 beg 4 0 -1 end 5 0 -2 nseg 4
sel node id=1 0 0 0 sel node id=2 2 0 0 sel node id=3 4 0 -1 sel node id=4 5 0 -2 sel beamsel cid=1 id=1 node 1 2 ; sel beamsel cid=2 id=1 node 2 3
1、结构单元的类型
FLAC3D中包含六种形式的结构单元，可以分成两类：
➢ 线型结构单元：
• 梁单元(beam) • 锚索单元(cable) • 桩单元(pile)
➢ 壳型结构单元：
• 壳单元(shell)
➢ FLAC3D中的结构单元是岩土工程中实际结 构的一种“抽象”，即采用简单的单元形式 来模拟复杂的结构体。
• gr_coh——单位长度上水泥浆粘结力cg
• xciy——梁结构y轴惯性矩, Iy • gr_fric——水泥浆的摩擦角φg
• xciz——梁结构z轴惯性矩，Ix • gr_k——单位长度上水泥浆刚度kg
• xcij——极惯性矩，J
• gr_per——水泥浆外圈周长Pg
• density——密度，ρ
结构单元的显示！GUI操作和命令操作（manual）！ 调整好显示效果后可以将显示的命令文件另存出来，以备下次使用。（最适用于 几何模型相同，参数不同的，不同工况分析的比较）
2、结构单元的建模方法—线型结构单元
先建立节点再联接成单元的方法；
2、结构单元的建模方法—壳型结构单元
➢ 壳单元
2、结构单元的建模方法—壳型结构单元
➢ 结构单元由结构节点(node)和结构构件 (SELs)构成。
➢ 结构单元中的节点(node)可以与周围的实体 网格(zone)或其它结构节点建立连接(link)，
• 土工格栅(geogrid)
通过连接实现岩土体或结构与其它结构发生
• 衬砌单元(liner)
相互作用。
➢ 注意：结构节点并不是简单地与实体网格的
def set_vals global ptA = 25.0 * sin( 40.0*degrad ) ; global ptB = 25.0 * cos( 40.0*degrad ) end @set_vals generate zone cylinder p0=( 0.0, 0.0, 0.0 ) &
end2=(0.0,25.0,0.0) radius=24.5 not plot add zg plot ad sel geom delete zones ; delete all zones sel node init zpos add -25.0
2、结构单元的建模方法—壳型结构单元
通过附着在实体网格表面来生成shell单元。
The shells can then be repositioned if ecessary by using the SEL node init command
2、结构单元的建模方法—注意事项
➢ FLAC3D是岩土工程的专业软件，因此一般很少用来做专门的结构 分析。在涉及到结构单元的问题中，往往都要考虑结构与周围的实 体单元的相互作用。在结构单元的建模时要特别注意一个基本原则： 一个zone至多包含一个structure node!
• slide——大变形滑动标志
• pmoment——塑性矩，Mp • thexp——热膨胀系数，αt • ydirection——矢量Y
• slide_tol——大变形滑动容差 • ycomp——抗压强度(力) • density——密度
• thexp——热膨胀系数
3、结构单元的参数取值
某些结构单元参数的取值要视具体情况 而定，根据经验且必要时调整参数通过 试算来确定。