×

基于NURBS的复杂域边界单元改进方法

消耗积分:3 | 格式:rar | 大小:333 | 2008-12-03

麻酱

分享资料个

提出了基于非均匀有理B样条(NURBS)理论的边界单元改进方法。该方法以原分网后的节点和单元数据为依据,按照平均取值的原则获取边界单元边界边上的控制点,根据这些控制点利用NURBS理论表示边界边,并可通过不断调节权因子,使分网后的边界更加接近真实边界。最后以VC++为平台,通过算例验证了该改进方法的有效性,并推导了改进后边界单元的形函数。
关 键 词 非均匀有理B样条; 复杂域; 边界单元; 改进方法; 分网

有限元法(Finite Element Method,FEM)作为最有效的数值方法之一,已成功用于解决许多工程问题。有限元模型的建立,是进行有限元分析(Finite Element Analysis,FEA)的前提和基础。分网又是有限元建模中最耗时的环节,是有限元分析的关键,它对计算精度、计算规模等都会产生较大的影响。
为满足功能或性能要求,工程中存在大量由曲面和曲线组成的区域——复杂域组成的产品,如汽车覆盖件、车架、桥架等。这些复杂域能很好地满足产品性能和外观要求,但同时也给工程分析带来一定难度,主要体现在分网的精度上。
目前,在分网中,三角形单元应用较为广泛,但有些工程问题需要采用四边形单元。与三角形单元相比,四边形单元的数据输入和对输出应力的解释更容易,但不能很好地表达复杂域的实际边界,精度不能保证[1]。为解决该问题,可采用减小单元尺寸或提高单元阶次的方法来提高精度。但这些方法会使节点和单元数成倍增加,对计算机性能要求较高。
针对以上问题,本文提出了基于非均匀有理B样条(Non-Uniform Rational B-Spline,NURBS)的复杂域边界单元的改进方法。

1 NURBS基本理论
NURBS方法,是目前最佳的曲线、曲面造型方法。由于它能够精确描述包括二次曲线在内的各种曲线,无论是解析形状还是自由格式形状均有统一的表示参数,便于工程数据库的存取和应用,并可通过控制点和权因子方便地对曲线、曲面的形状进行局部调整。
NUBRS的提出,是为了将所有空间曲线、曲面用统一数学形式描述。国际标准组织(ISO)继美国的PDES后,于1991年颁布了关于工业产品数据交换的STEP国际标准,把NURBS作为定义工业产品几何形状的唯一数学方法。
NURBS对曲线的数学定义为:已知空间n+1个点di,(i=0,1,􀀢,n),则一条k次NURBS曲线可表示为分段有理多项式矢函数:
,,00()()/()nniiikiikiipuwdNuwNu===ΣΣ 1[,]kkuuu+∈ (1)
式中 di为控制点,i= 0,1,􀀢,n;wi为相应控制点的权因子(Weight Factor),i= 0,1,􀀢,n;为k次B样条基函数;[u)(,uNki0,u1,􀀢,un+k]为曲线节点矢量。
2 复杂域边界单元的改进
2.1 复杂域边界单元改进思想
采用一般四边形分网方法得到边界单元如图1所示。图中i、j为边界单元与复杂域边界的交点;为分析域边界的NURBS表示,i点、j点坐标分别为、。采用平均取值的方法,在边界边上的和间取m个控制点,即d()pu)(iup)(jup
)(iup
(jup0=,dpm−1= ,d)(jupj=。用式(1)表示边界单元的边界边。 ((()/ium+−
为保证在公共节点处连续,取首末权因子为w0=wm−1=1,其余的权因子wi],0[+∞∈可取不同的值。当wi取值适当时,边界单元边界边将在控制点处与边界实现最佳逼近。控制点越多,逼近程度越好。同时,为保证NURBS曲线不出现奇异,wi取值不小于0。

声明:本文内容及配图由入驻作者撰写或者入驻合作网站授权转载。文章观点仅代表作者本人,不代表电子发烧友网立场。文章及其配图仅供工程师学习之用,如有内容侵权或者其他违规问题,请联系本站处理。 举报投诉

评论(0)
发评论

下载排行榜

全部0条评论

快来发表一下你的评论吧 !

'+ '

'+ '

'+ ''+ '
'+ ''+ ''+ '
'+ ''+ '' ); $.get('/article/vipdownload/aid/'+webid,function(data){ if(data.code ==5){ $(pop_this).attr('href',"/login/index.html"); return false } if(data.code == 2){ //跳转到VIP升级页面 window.location.href="//m.obk20.com/vip/index?aid=" + webid return false } //是会员 if (data.code > 0) { $('body').append(htmlSetNormalDownload); var getWidth=$("#poplayer").width(); $("#poplayer").css("margin-left","-"+getWidth/2+"px"); $('#tips').html(data.msg) $('.download_confirm').click(function(){ $('#dialog').remove(); }) } else { var down_url = $('#vipdownload').attr('data-url'); isBindAnalysisForm(pop_this, down_url, 1) } }); }); //是否开通VIP $.get('/article/vipdownload/aid/'+webid,function(data){ if(data.code == 2 || data.code ==5){ //跳转到VIP升级页面 $('#vipdownload>span').text("开通VIP 免费下载") return false }else{ // 待续费 if(data.code == 3) { vipExpiredInfo.ifVipExpired = true vipExpiredInfo.vipExpiredDate = data.data.endoftime } $('#vipdownload .icon-vip-tips').remove() $('#vipdownload>span').text("VIP免积分下载") } }); }).on("click",".download_cancel",function(){ $('#dialog').remove(); }) var setWeixinShare={};//定义默认的微信分享信息,页面如果要自定义分享,直接更改此变量即可 if(window.navigator.userAgent.toLowerCase().match(/MicroMessenger/i) == 'micromessenger'){ var d={ title:'基于NURBS的复杂域边界单元改进方法',//标题 desc:$('[name=description]').attr("content"), //描述 imgUrl:'https://'+location.host+'/static/images/ele-logo.png',// 分享图标,默认是logo link:'',//链接 type:'',// 分享类型,music、video或link,不填默认为link dataUrl:'',//如果type是music或video,则要提供数据链接,默认为空 success:'', // 用户确认分享后执行的回调函数 cancel:''// 用户取消分享后执行的回调函数 } setWeixinShare=$.extend(d,setWeixinShare); $.ajax({ url:"//www.obk20.com/app/wechat/index.php?s=Home/ShareConfig/index", data:"share_url="+encodeURIComponent(location.href)+"&format=jsonp&domain=m", type:'get', dataType:'jsonp', success:function(res){ if(res.status!="successed"){ return false; } $.getScript('https://res.wx.qq.com/open/js/jweixin-1.0.0.js',function(result,status){ if(status!="success"){ return false; } var getWxCfg=res.data; wx.config({ //debug: true, // 开启调试模式,调用的所有api的返回值会在客户端alert出来,若要查看传入的参数,可以在pc端打开,参数信息会通过log打出,仅在pc端时才会打印。 appId:getWxCfg.appId, // 必填,公众号的唯一标识 timestamp:getWxCfg.timestamp, // 必填,生成签名的时间戳 nonceStr:getWxCfg.nonceStr, // 必填,生成签名的随机串 signature:getWxCfg.signature,// 必填,签名,见附录1 jsApiList:['onMenuShareTimeline','onMenuShareAppMessage','onMenuShareQQ','onMenuShareWeibo','onMenuShareQZone'] // 必填,需要使用的JS接口列表,所有JS接口列表见附录2 }); wx.ready(function(){ //获取“分享到朋友圈”按钮点击状态及自定义分享内容接口 wx.onMenuShareTimeline({ title: setWeixinShare.title, // 分享标题 link: setWeixinShare.link, // 分享链接 imgUrl: setWeixinShare.imgUrl, // 分享图标 success: function () { setWeixinShare.success; // 用户确认分享后执行的回调函数 }, cancel: function () { setWeixinShare.cancel; // 用户取消分享后执行的回调函数 } }); //获取“分享给朋友”按钮点击状态及自定义分享内容接口 wx.onMenuShareAppMessage({ title: setWeixinShare.title, // 分享标题 desc: setWeixinShare.desc, // 分享描述 link: setWeixinShare.link, // 分享链接 imgUrl: setWeixinShare.imgUrl, // 分享图标 type: setWeixinShare.type, // 分享类型,music、video或link,不填默认为link dataUrl: setWeixinShare.dataUrl, // 如果type是music或video,则要提供数据链接,默认为空 success: function () { setWeixinShare.success; // 用户确认分享后执行的回调函数 }, cancel: function () { setWeixinShare.cancel; // 用户取消分享后执行的回调函数 } }); //获取“分享到QQ”按钮点击状态及自定义分享内容接口 wx.onMenuShareQQ({ title: setWeixinShare.title, // 分享标题 desc: setWeixinShare.desc, // 分享描述 link: setWeixinShare.link, // 分享链接 imgUrl: setWeixinShare.imgUrl, // 分享图标 success: function () { setWeixinShare.success; // 用户确认分享后执行的回调函数 }, cancel: function () { setWeixinShare.cancel; // 用户取消分享后执行的回调函数 } }); //获取“分享到腾讯微博”按钮点击状态及自定义分享内容接口 wx.onMenuShareWeibo({ title: setWeixinShare.title, // 分享标题 desc: setWeixinShare.desc, // 分享描述 link: setWeixinShare.link, // 分享链接 imgUrl: setWeixinShare.imgUrl, // 分享图标 success: function () { setWeixinShare.success; // 用户确认分享后执行的回调函数 }, cancel: function () { setWeixinShare.cancel; // 用户取消分享后执行的回调函数 } }); //获取“分享到QQ空间”按钮点击状态及自定义分享内容接口 wx.onMenuShareQZone({ title: setWeixinShare.title, // 分享标题 desc: setWeixinShare.desc, // 分享描述 link: setWeixinShare.link, // 分享链接 imgUrl: setWeixinShare.imgUrl, // 分享图标 success: function () { setWeixinShare.success; // 用户确认分享后执行的回调函数 }, cancel: function () { setWeixinShare.cancel; // 用户取消分享后执行的回调函数 } }); }); }); } }); } function openX_ad(posterid, htmlid, width, height) { if ($(htmlid).length > 0) { var randomnumber = Math.random(); var now_url = encodeURIComponent(window.location.href); var ga = document.createElement('iframe'); ga.src = 'https://www1.elecfans.com/www/delivery/myafr.php?target=_blank&cb=' + randomnumber + '&zoneid=' + posterid+'&prefer='+now_url; ga.width = width; ga.height = height; ga.frameBorder = 0; ga.scrolling = 'no'; var s = $(htmlid).append(ga); } } openX_ad(828, '#berry-300', 300, 250);