作者:二狗家大叔
物体受力在单一应力状态下,内部各点也是单一状态。在部分科研理论分析中,需要在物体内部挖孔,利用逆推法推导物理量。例如下图所示,受Y方向某拉力作用,各点应力状态为:
![9d6cd7b6bc9528ef1c25c1cf35e874e8.png](https://img-blog.csdnimg.cn/img_convert/9d6cd7b6bc9528ef1c25c1cf35e874e8.png)
![93b02a89c215b8bccd8d173566d2cd10.png](https://img-blog.csdnimg.cn/img_convert/93b02a89c215b8bccd8d173566d2cd10.png)
在圆孔中心位置建立圆柱坐标系,该应力状态在圆柱坐标系下的公式为:
![743382b9bcd6c98cc2e7c03ca5dee646.png](https://img-blog.csdnimg.cn/img_convert/743382b9bcd6c98cc2e7c03ca5dee646.png)
在这种情况下反推物理量,需要对曲面施加基于圆柱坐标系的面力。
案例如下:在圆弧面基于圆柱坐标系施加等效于单向应力状态的面力。
![114b719f57413586517becf3ed5b667d.png](https://img-blog.csdnimg.cn/img_convert/114b719f57413586517becf3ed5b667d.png)
![d0f87781d55fa615c0647384cabddebe.png](https://img-blog.csdnimg.cn/img_convert/d0f87781d55fa615c0647384cabddebe.png)
加载前先建立圆柱坐标系(注意R轴方向为0度位置,T轴方向为角度增大方向,示意图见文后的加载图)
具体设置方法为:Load>Create Load>Mechanical>surface traction
选中中间曲面后,先设置径向力,按以下参数设置:
![d0f87781d55fa615c0647384cabddebe.png](https://img-blog.csdnimg.cn/img_convert/d0f87781d55fa615c0647384cabddebe.png)
![73a289413f08846cdd4d88cabcb6b0a8.png](https://img-blog.csdnimg.cn/img_convert/73a289413f08846cdd4d88cabcb6b0a8.png)
Distribution:应力分配,点击后面的f(x)创建一个基于圆柱坐标系的表达式,Local system 要选择圆柱坐标系,Th为角度变量。
![a325552000faf27e846fe09b26032465.png](https://img-blog.csdnimg.cn/img_convert/a325552000faf27e846fe09b26032465.png)
Traction:选择General,为一般力。
Vector:点击选择图标后,依次选择(0,0,0) (-1,0,0) ,坐标选择建立的圆柱坐标系。
注:面力方向矢量是基于所选坐标系,(-1,0,0)就是沿圆柱坐标系下的R轴反向。
Magnitude:选择应力大小为1。
然后在创建一个Load,设置切向力,如下图所示,也是基于圆柱坐标系。
![ad5c13e7cacdf87951e12dab4ce03cea.png](https://img-blog.csdnimg.cn/img_convert/ad5c13e7cacdf87951e12dab4ce03cea.png)
![fc1204464cd1bf932b44876d64e949ac.png](https://img-blog.csdnimg.cn/img_convert/fc1204464cd1bf932b44876d64e949ac.png)
再创建一个Load,在整体坐标系下对两侧的平面施加Y方向的面力,大小为1,同时对后面的面施加全约束。
最后加载形式为下图所示:
![272c1ac8799ef75a8a1c452f55afb15c.png](https://img-blog.csdnimg.cn/img_convert/272c1ac8799ef75a8a1c452f55afb15c.png)
求解结果如下图:
![d0f87781d55fa615c0647384cabddebe.png](https://img-blog.csdnimg.cn/img_convert/d0f87781d55fa615c0647384cabddebe.png)
![d0f87781d55fa615c0647384cabddebe.png](https://img-blog.csdnimg.cn/img_convert/d0f87781d55fa615c0647384cabddebe.png)
![652e357119db5a6e358eefab2b3e5834.png](https://img-blog.csdnimg.cn/img_convert/652e357119db5a6e358eefab2b3e5834.png)
大部分位置应力在0.99~1.01之间,为单向应力状态,加载方式正确。
本问题的关键是面力的方向问题,在选择面力的方向矢量时,是基于所选坐标系。对于圆柱坐标系,切向力矢量为(0,-1,0)时,即力的方向只沿着theta的反方向。