请教三维框架有限元计算边界处理的罚函数问题
浏览：118167 回答：1
各位专家，学习《工程中的有限元方法》（曾攀 译）关于三维框架有限元计算边界处理时遇到一个问题，该书采用了罚函数来处理边界条件，具体为：
罚函数C = 10^4 × max { 总体刚度矩阵中有约束的相应行的第一个元素 } ——说明：“有约束的行”根据边界约束对应的自由度确定
总体刚度矩阵中支座约束自由度所在行的第1个元素 = 自己 + 罚函数C
荷载数组中对应约束自由度的项 = 自己 + 惩罚参数 × 相应约束位移
VB中的算法如下：
惩罚参数 = 0
For i = 1 To 结构自由度
If 惩罚参数 < 总刚(i, 1) Then 惩罚参数 = 总刚(i, 1)
Next i
惩罚参数 = 惩罚参数 * 10000
For i = 1 To 给定位移数量
位移自由度 = 位移自由度数组(i)
总刚(位移自由度, 1) = 总刚(位移自由度, 1) + 惩罚参数
荷载数组(位移自由度) = 荷载数组(位移自由度) + 惩罚参数 * 位移数组(i)
Next i
用罚函数处理边界位移以前没有学过，请专家指导计算的原理，或推荐参考书。不胜感激！！
:-|
罚函数C = 10^4 × max { 总体刚度矩阵中有约束的相应行的第一个元素 } ——说明：“有约束的行”根据边界约束对应的自由度确定
总体刚度矩阵中支座约束自由度所在行的第1个元素 = 自己 + 罚函数C
荷载数组中对应约束自由度的项 = 自己 + 惩罚参数 × 相应约束位移
VB中的算法如下：
惩罚参数 = 0
For i = 1 To 结构自由度
If 惩罚参数 < 总刚(i, 1) Then 惩罚参数 = 总刚(i, 1)
Next i
惩罚参数 = 惩罚参数 * 10000
For i = 1 To 给定位移数量
位移自由度 = 位移自由度数组(i)
总刚(位移自由度, 1) = 总刚(位移自由度, 1) + 惩罚参数
荷载数组(位移自由度) = 荷载数组(位移自由度) + 惩罚参数 * 位移数组(i)
Next i
用罚函数处理边界位移以前没有学过，请专家指导计算的原理，或推荐参考书。不胜感激！！
:-|
